青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

Tauruser

Enjoy Every Day
posts - 34, comments - 95, trackbacks - 0, articles - 5
  C++博客 :: 首頁 :: 新隨筆 :: 聯系 :: 聚合  :: 管理

Gauss消去法直接求解方程組(附例程)

Posted on 2006-06-04 20:22 Tauruser 閱讀(3965) 評論(7)  編輯 收藏 引用 所屬分類: 數值計算

3.2 Gauss消去法

3.2.1 順序消去法

  Gauss消去法就是將方程組(3.1.1)通過(n-1)步消元,將(3.1.1)轉化為上三角方程組
       (3.2.1)
再回代求此方程組的解.
  下面記增廣矩陣,即
   
  第1步 設,計算l,記為,若用第一行加到第i行,可消去,用Gauss變換矩陣表示
    
令   
其中
  一般地,假定已完成了(k-1)步消元,即已將轉化為以下形式:
   
  第k步,假定,計算
       (3.2.2)
,則
   
其中
     (3.2.3).
當k=1,2,…,n-1則可得到,即方程組(3.2.1).

  直接回代解(3.2.1)得,
     (3.2.4)
并且有,由以上順序消去過程可得如下定理.
  定理2.1 設非奇異,則通過兩行互換總可使,k=1,2,…,n-1.可將方程組(3.1.1)轉化為(3.2.1)并求得方程組(3.1.1)的解為(3.2.4),且有.
  如果不做行交換,則使的條件如下.

  定理2.2 非奇異,且各階順序主子式, 則,k=1,2,…,n-1.
  證明 用歸納法,當,故.現假設(k-1)成立,即,對i=1,2,…,k-1已推出,故Gauss消去法能進行(k-1)步消元,A已約化為,即
   

     

對k=1,2,…,n均成立,證畢.
  在整個消去法消元過程中,k從1到(n-1)共需乘除法運算次數為
     
加減法次數為
     
回代過程中由公式(3.2.4)可知乘除法次數為,加減法次數為,于是Gauss消去法的乘除法總次數為,加減法次數為
  例3.4 用Gauss消去法解方程組
      
并求detA.
  解 消元得
   
再由(3.2.4)回代,得解

講解:

  Gauss 消去法是將方程組AX=b,通過消元轉化為上三角方程組(3,2,1)求解,消元第一步做完后有
         
  用矩陣表示
  第K-1步完成后得到
  當,可做K步,得到
  得到,對應的方程組就是(3.2.1),利用公式(3.2.4)就可求得解。
  定理2.2給出了進行順序消去法的條件,即A的所有順序生子式,而方程(3.1.1)解存在唯一的條件是

好了,原理講完了,貼我的例程。

#include? < iostream >
#include?
< vector >
#include?
< cmath >
using ? namespace ?std;
class ?CGAUSSSOLVEEQU
{
private :
????vector
< vector < double >> ?m_equset;
????vector
< double > ?m_answer;
????unsigned?
int ?m_n;
public :
????
void ?inputEquSet();
????
void ?solveEquSet();
????
void ?outputAnswer();
}
;
void ?CGAUSSSOLVEEQU::inputEquSet()
{
????
double ?dtemp;
????vector
< double > ?vtemp;
????cout
<< " 請輸入你的方程個數: " ;
????cin
>> m_n;
????cout
<< " 請按照向量的形式輸入各變量的系數。最后一位為b。每個方程一行: " << endl;
????
for (unsigned? int ?i( 0 );i < m_n; ++ i)
????
{
????????m_equset.push_back(vtemp);
????????
for (unsigned? int ?j( 0 );j <= m_n; ++ j)
????????
{????
????????????cin
>> dtemp;
????????????m_equset[i].push_back(dtemp);
????????}

????????
if (m_equset[i].size() != m_n + 1 )
????????
{
????????????cout
<< " 輸入有誤,請重新輸入上一個方程。 " << endl;
????????????
-- i;
????????}

????}

????
}


void ?CGAUSSSOLVEEQU::solveEquSet()
{
????vector
< vector < double >> ::iterator?iter;
????iter
= m_equset.begin();
????
for (unsigned? int ?m( 0 );m < m_n - 1 ; ++ m)
????
{
????????
// 將絕對值最大的主元素移上去。此舉是為了減少誤差
???????? for (vector < vector < double >> ::iterator?iter2 = iter + 1 ;iter2 != m_equset.end(); ++ iter2)
????????
{
????????????
if (fabsl(iter -> front()) < fabsl(iter2 -> front()))
????????????
{
????????????????swap(
* iter, * iter2);
????????????}

????????}

????????
// 進行消元
???????? for (unsigned? int ?i = m + 1 ;i < m_n; ++ i)
????????
{
????????????
double ?dm;
????????????dm
= m_equset[i][m] / m_equset[m][m];
????????????
for (unsigned? int ?j = m;j < m_n + 1 ; ++ j)
????????????
{
????????????????m_equset[i][j]
-= dm * m_equset[m][j];
????????????}

????????}

????????
++ iter;
????}

????
// 初始化m_answer向量
???? for (unsigned? int ?i( 0 );i < m_n; ++ i)?m_answer.push_back( 0 );
????
// 求解答案
????m_answer[m_n - 1 ] = m_equset[m_n - 1 ][m_n] / m_equset[m_n - 1 ][m_n - 1 ];

????
for ( int ?i = m_n - 2 ;i >= 0 ; -- i)
????
{
????????m_answer[i]
= m_equset[i][m_n];
????????
for ( int ?j = m_n - 1 ;j > i; -- j)
????????????m_answer[i]
-= m_answer[j] * m_equset[i][j];
????????m_answer[i]
/= m_equset[i][i];
????}


????
}


void ?CGAUSSSOLVEEQU::outputAnswer()
{
????
for (unsigned? int ?i( 1 );i <= m_n; ++ i)
????
{
????????cout
<< " x( " << i << " )= " << m_answer[i - 1 ] << endl;
????}

}


int ?main()
{
????CGAUSSSOLVEEQU?myEqu;
????myEqu.inputEquSet();
????myEqu.solveEquSet();
????myEqu.outputAnswer();
????
return ? 0 ;
}


// Power?By?Tauruser?2006.6.4

Feedback

# re: Gauss消去法直接求解方程組(附例程)  回復  更多評論   

2006-11-17 13:38 by tim
ding

# re: Gauss消去法直接求解方程組(附例程)  回復  更多評論   

2007-06-01 11:12 by Snowfox
請問LZ:
能發下具體的結果嗎?我不明白你的輸入格式,還有能解幾元幾次的方程組。謝謝了!勞駕!
期待中。。。

# re: Gauss消去法直接求解方程組(附例程)  回復  更多評論   

2007-06-10 12:02 by 星夢情緣
我的BLOG集齊了各重解線性方程的方法,任意方程都可以解的

# re: Gauss消去法直接求解方程組(附例程)  回復  更多評論   

2008-06-04 20:11 by
有一個小錯誤

# re: Gauss消去法直接求解方程組(附例程)  回復  更多評論   

2008-06-04 20:11 by
for ( int i = m_n - 2 ;i >= 0 ; -- i)

i 重復定義了,不過謝謝你呀!

# re: Gauss消去法直接求解方程組(附例程)  回復  更多評論   

2008-06-04 20:12 by
好程序

# re: Gauss消去法直接求解方程組(附例程)  回復  更多評論   

2009-03-25 22:06 by nan
能把改后的程序發上去嗎,謝謝
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            亚洲欧美欧美一区二区三区| 国内精品视频一区| 亚洲一区免费看| 亚洲一区二区三区精品在线| 欧美午夜精品久久久久久孕妇 | 国产精品久久久亚洲一区 | 亚洲伦理自拍| 亚洲视屏在线播放| 国产日韩免费| 欧美激情一区二区三区在线视频| 欧美www在线| 亚洲综合国产激情另类一区| 欧美在线看片| 亚洲最新视频在线| 亚洲影音先锋| 亚洲黄色av| 亚洲一区在线视频| 亚洲国产欧美一区二区三区久久 | 欧美大片一区二区| 性伦欧美刺激片在线观看| 久久国产精品久久久| 亚洲另类在线一区| 亚洲欧美日韩在线观看a三区| 伊人男人综合视频网| 99视频在线精品国自产拍免费观看 | 欧美一级日韩一级| 99re8这里有精品热视频免费 | 欧美ed2k| 久久久99久久精品女同性| 欧美激情第9页| 欧美一区激情视频在线观看| 欧美成人精品三级在线观看| 性色一区二区| 欧美日韩不卡视频| 免费观看久久久4p| 国产区精品在线观看| 最新中文字幕亚洲| 国产一区二区久久久| 亚洲精选91| 91久久精品www人人做人人爽 | 久久综合99re88久久爱| 免费欧美网站| 久久精品一区蜜桃臀影院| 欧美日韩一级片在线观看| 欧美va亚洲va国产综合| 国产日韩精品在线| 亚洲无限av看| 亚洲一区影音先锋| 欧美日韩中文在线| 日韩亚洲一区二区| 日韩视频不卡中文| 欧美顶级大胆免费视频| 蜜桃av噜噜一区| 精品福利免费观看| 久久成人精品无人区| 香蕉乱码成人久久天堂爱免费 | 国产区精品在线观看| 亚洲性视频h| 亚洲视频精选在线| 欧美视频中文字幕在线| 亚洲精品一区在线| 99精品欧美一区二区三区综合在线| 另类成人小视频在线| 免费成人黄色片| 在线日本成人| 蜜桃久久精品乱码一区二区| 欧美成人激情视频免费观看| 亚洲第一页在线| 欧美a级一区| 亚洲伦伦在线| 亚洲尤物视频在线| 国产日韩欧美不卡| 久久精品国产一区二区三| 老司机午夜精品视频| 在线看不卡av| 欧美激情国产日韩精品一区18| 最近看过的日韩成人| 在线视频亚洲| 国产欧美精品日韩精品| 久久久久久久网站| 亚洲激情二区| 亚洲欧美国产高清| 精品二区视频| 欧美精品国产精品日韩精品| 中文无字幕一区二区三区| 午夜电影亚洲| 一区二区三区自拍| 欧美激情综合五月色丁香| 亚洲小说欧美另类社区| 久久久91精品| 夜夜爽www精品| 国产精品一区二区三区免费观看| 久久久久久久久久久一区| 亚洲欧洲综合| 久久久噜噜噜久久| 99视频精品免费观看| 国产精品最新自拍| 欧美sm视频| 午夜精品999| 91久久视频| 久久频这里精品99香蕉| 中文av字幕一区| 在线电影国产精品| 国产精品九九| 欧美激情精品久久久久久久变态 | 美女日韩在线中文字幕| 日韩视频一区二区在线观看| 久久字幕精品一区| 午夜视频一区二区| 亚洲理伦电影| 伊人男人综合视频网| 国产精品进线69影院| 欧美承认网站| 久久精品国产精品亚洲精品| 国产精品99久久久久久久女警 | 欧美四级电影网站| 久久久免费精品视频| 亚洲欧美日韩精品久久| 99精品视频免费全部在线| 欧美成人激情视频免费观看| 久久国产88| 午夜免费日韩视频| 亚洲婷婷在线| 中文久久精品| 99re66热这里只有精品4| 亚洲第一天堂av| 伊人成人在线视频| 国产一区二区三区高清播放| 欧美午夜精品久久久久久孕妇| 欧美成人亚洲成人日韩成人| 麻豆成人在线播放| 久久阴道视频| 老司机成人网| 蜜臀99久久精品久久久久久软件 | 欧美在线精品免播放器视频| 亚洲影院在线| 亚洲欧美一区二区激情| 亚洲影院在线观看| 午夜精品久久久久久| 亚洲欧美日本国产有色| 亚洲免费在线视频| 亚洲欧美一区二区三区极速播放| 宅男精品导航| 亚洲欧美一区二区三区久久| 午夜影院日韩| 久久精品一区二区三区四区 | 香蕉尹人综合在线观看| 午夜在线不卡| 久久免费午夜影院| 欧美a级大片| 欧美日韩国产精品| 国产精品久久久久毛片大屁完整版 | 亚洲日本无吗高清不卡| 亚洲精品偷拍| 夜夜嗨av色一区二区不卡| 亚洲午夜精品国产| 欧美一二三区精品| 久久女同互慰一区二区三区| 老色鬼精品视频在线观看播放| 欧美国产一区在线| 欧美视频在线看| 国产真实久久| 日韩亚洲欧美一区| 欧美一区二区三区久久精品茉莉花 | 欧美大片在线看免费观看| 亚洲福利国产精品| 日韩视频在线观看一区二区| 欧美在线高清视频| 亚洲少妇中出一区| 欧美一级理论片| 免费欧美在线| 亚洲一区二区三区四区在线观看 | 亚洲黄色av一区| 亚洲女人天堂av| 另类成人小视频在线| 亚洲人成网站精品片在线观看| 亚洲影院免费观看| 女仆av观看一区| 国产精品一区久久久久| 亚洲美女尤物影院| 久久精品一区二区三区不卡| 亚洲三级国产| 久久午夜精品一区二区| 欧美视频在线观看一区二区| 亚洲福利国产| 久久av一区二区| 亚洲精品色图| 久久精品国产精品亚洲精品| 欧美a级一区| 美国十次成人| 国产精品欧美日韩| 亚洲精品久久久久久久久| 欧美在线视频一区| a4yy欧美一区二区三区| 久久亚洲私人国产精品va媚药| 国产精品美女久久久久av超清| 亚洲日本va午夜在线电影| 久久婷婷丁香| 欧美一区二区三区的| 国产精品久久久久久久7电影|