非線性方程組的數(shù)值算法研究 畢業(yè)論文_第1頁(yè)
已閱讀1頁(yè),還剩20頁(yè)未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、<p>  南 昌 工 程 學(xué) 院</p><p><b>  畢 業(yè) 論 文</b></p><p>  理學(xué) 系(院) 信息與計(jì)算科學(xué) 專業(yè)</p><p>  畢業(yè)論文題目 非線性方程組的數(shù)值算法研究 </p><p>  學(xué)生姓名 張浩浩 </

2、p><p>  班 級(jí) 09信息與計(jì)算科學(xué) </p><p>  學(xué) 號(hào) 2009101533 </p><p>  指導(dǎo)教師 禹海雄 </p><p>  完成日期2013年04月13</p><p><b

3、>  目錄</b></p><p><b>  摘要3</b></p><p>  Abstract4</p><p><b>  第一章 緒論5</b></p><p>  第二章、求解非線性方程組的幾種方法6</p><p><b> 

4、 2.1、牛頓法6</b></p><p>  2.1.1牛頓法的引入與介紹6</p><p>  2.1.2牛頓法的算法8</p><p>  2.1.3牛頓法代碼程序編程8</p><p>  2.2、擬牛頓法11</p><p>  2.2.1擬牛頓法的引入與介紹11</p>

5、<p>  2.2.2擬牛頓法的算法12</p><p>  2.2.3擬牛頓法的題例分析12</p><p>  2.3、割線法14</p><p>  2.3.1割線法的引入與介紹14</p><p>  2.3.2割線法的總結(jié)陳述15</p><p>  3.3.3割線法題例分析16<

6、;/p><p><b>  結(jié)束語(yǔ)18</b></p><p><b>  參考文獻(xiàn)19</b></p><p><b>  致謝詞20</b></p><p>  非線性方程組的數(shù)值算法研究</p><p>  Study on numerical

7、algorithms for nonlinear equations</p><p>  總計(jì) 畢 業(yè) 論 文 21 頁(yè)</p><p>  表 格 2 個(gè)</p><p><b>  摘要</b></p><p>  論文講解的是非線性

8、方程的數(shù)值的求解方法,課本中我們接觸到求解線性方程的方法比較多,相對(duì)于求解非線性方程組數(shù)值方法比較繁瑣和計(jì)算量大,同時(shí)課本上只是簡(jiǎn)單介紹非線性方程組數(shù)值的求解方法。我在圖書館查閱了資料和在老師的指導(dǎo)下,仔細(xì)研究把思路整理如下:我在這邊論文就如何求解非線性數(shù)值的求解方法闡述了牛頓法、擬牛頓法、割線法三種方法來求解非線性方程組數(shù)值,并且通過了列舉了題例可以比較更鮮明的可以看出3種方法的聯(lián)系和特點(diǎn)。</p><p> 

9、 關(guān)鍵詞:Newton法、迭代法、擬Newton法。</p><p><b>  .</b></p><p><b>  Abstract</b></p><p>  The paper explained that the method of solving nonlinear equations numerically

10、, text we come into contact with the method of solving linear equation more, compared with the numerical method for solving nonlinear equations is complicated and large amount of calculation, at the same time the textboo

11、k only briefly method for solving nonlinear equations. I looked up information and under the guidance of the teacher in the library, a careful study of the ideas as follows: I how to solve nonlinear num</p><p&

12、gt;  Keywords: Newton method, iterative method, quasi Newton method.</p><p><b>  第一章 緒論</b></p><p>  我們先了解非線性方程組的一般形式如下:</p><p>  我們可以看出是在空間的實(shí)值函數(shù)。</p><p>  

13、我們?cè)儆孟蛄哭D(zhuǎn)換下可以得到:</p><p><b>  ,x=,0=</b></p><p>  我們此時(shí)可以把方程換成:</p><p>  。 (1)</p><p>  我們可以把F可以看做在區(qū)域內(nèi)展開的非線性映像,表示為</p>&

14、lt;p>  下面我們來介紹簡(jiǎn)單的邊值問題:</p><p>  , 。 (2)</p><p>  我們此時(shí)定義f在D=上二階可微連續(xù),</p><p>  現(xiàn)在我們求解(2)上x的數(shù)值。</p><p>  我們用差分方法離散化得到:</p><p>  , j=0,1,2,3,、、n+1 ,<

15、/p><p><b>  在得到:</b></p><p>  j+1,2,、、、n,</p><p>  我們?cè)谵D(zhuǎn)化矩陣又可以得到:</p><p><b>  A=</b></p><p><b>  在從映像轉(zhuǎn)換成:</b></p>&l

16、t;p><b>  ,</b></p><p><b>  方程(2)轉(zhuǎn)化為:</b></p><p><b>  Ax+=0</b></p><p>  本論文將介紹求解非線性方程組的牛頓法,迭代法,牛頓法,這是本人對(duì)非線性方程數(shù)值求解的認(rèn)識(shí),我會(huì)使用這些方法并為為開展進(jìn)一步研究。</p

17、><p>  第二章、求解非線性方程組的幾種方法</p><p><b>  2.1、牛頓法</b></p><p>  2.1.1牛頓法的引入與介紹</p><p>  我們?cè)趯W(xué)習(xí)中關(guān)于方程f(x)=0的求解這種題型接觸的太多了,f(x)作為線性方程函數(shù),解法多樣也很容易求解值。我們來比較一下牛頓法,牛頓法簡(jiǎn)單的來說其實(shí)也

18、是一種線性化方法,他的理念就是把非線性方程f(x)轉(zhuǎn)化成某種類型的線性方程求解x的值。非線性方程不過是線性方程的擴(kuò)展,非線性方程組就是在此基礎(chǔ)上加以延伸。</p><p>  下面我們來介紹了解一下牛頓法的理論:</p><p><b>  我們看下例題:</b></p><p><b>  (1)</b></p&g

19、t;<p>  從上面非線性方程組我們可以看出,... 是的多元函數(shù),這是我們也可以用向量把它轉(zhuǎn)化為</p><p>  我們同時(shí)把他轉(zhuǎn)化為:</p><p><b>  (2)</b></p><p>  我們可以看出時(shí),至少有一個(gè)變量是在的非線性函數(shù),我們這時(shí)(1)就可以看作非線性方程組,非線性方程組的求解實(shí)際上就是n=1求根

20、的應(yīng)用。也就是把單一變量的函數(shù)轉(zhuǎn)化為向量函數(shù),這個(gè)時(shí)候就可以用求解單變量的方法來求解非線性方程組(2)。若果知道方程組=0的一個(gè)近似根,再用函數(shù)的分量在用多元函數(shù)泰勒的方法展開,提取線性方程就可以得到:</p><p><b>  ,</b></p><p><b>  我們令,</b></p><p><b>

21、  得到:</b></p><p><b>  (3)</b></p><p>  其中 (4)</p><p>  我們這時(shí)可以把(4)作為雅克比矩陣,(3)的線性方程組的解我們記作為,就可以得到:</p><p> ?。╧=0,1,2,.....)。

22、 (5)</p><p>  這就是我們所說的求解非線性方程組(2)的牛頓法。</p><p>  下面我們來簡(jiǎn)單介紹非線性方程組求解牛頓法的算法:</p><p>  從上面的實(shí)例我們可以看得出牛頓法求解非線性方程的主要理論是用在(k=0,1,2,...)的基礎(chǔ)上進(jìn)行迭代計(jì)算。我們這時(shí)所要做的就是計(jì)算出F(x)的雅克比矩陣,通過得到它的逆,直到達(dá)到所需要

23、的精度的范圍內(nèi)才停止迭代。</p><p>  2.1.2牛頓法的算法</p><p><b>  牛頓法算法如下:</b></p><p>  1.首先我們把所要求解的非線性方程組定義為,并為之確定精度。</p><p>  2.把轉(zhuǎn)化為雅克比矩陣,得到。求解方法如下:</p><p>  3.

24、重復(fù)第二步方法,求解雅克比矩陣的逆。另外把乘以單位矩陣,我們可以用單位矩陣轉(zhuǎn)換求解的逆用來保存。</p><p><b>  4.與的相乘</b></p><p>  5.再用 (k=0,1,2,.....)來迭代。</p><p>  6.最后我們注意的時(shí)精度,其精度時(shí),我們需要重復(fù)2—5次,一直使精度達(dá)到最?。ň龋r(shí)停止迭代,最后的迭代結(jié)

25、果為。</p><p>  2.1.3牛頓法代碼程序編程</p><p>  最后我們介紹代碼的編程:</p><p>  #include <iostream.h></p><p>  #include <stdlib.h></p><p>  #include <math.h>&

26、lt;/p><p>  #include <conio.h></p><p>  #define f0(x1,x2) (x1+2*x2-3)</p><p>  #define f1(x1,x2) (2*x1*x1+x2*x2-5)</p><p>  #define x_ 0.000001</p><p>

27、  #define matrixNum 2</p><p>  double *matrixF2(double *x);</p><p><b>  int y=0;</b></p><p>  void main()</p><p><b>  {</b></p><p>

28、  int i,j,n;</p><p>  double p,*x;</p><p>  double *b;</p><p>  double *matrixF; //矩陣F</p><p>  double *matrixF_; //矩陣F的雅可比矩陣的逆 </p><p>  b=(double *

29、)malloc(matrixNum);</p><p>  matrixF=(double *)malloc(matrixNum);</p><p>  matrixF_=(double *)malloc(matrixNum*matrixNum);</p><p>  cout<<"請(qǐng)輸入初值:";</p><p&

30、gt;  for(i=0;i<matrixNum;i++)</p><p>  cin>>*(x+i);</p><p><b>  do</b></p><p><b>  {</b></p><p><b>  p=0.0;</b></p>

31、<p>  for(i=0;i<matrixNum;i++)</p><p><b>  *(b+i)=0;</b></p><p>  *matrixF=f0(*x,*(x+1));</p><p>  *(matrixF+1)=f1(*x,*(x+1));</p><p>  matrixF_=mat

32、rixF2(x);</p><p>  for(i=0;i<matrixNum;i++)</p><p><b>  {</b></p><p>  for(j=0;j<matrixNum;j++)</p><p>  *(b+i)+=*(matrixF_+i*matrixNum+j)*(*(matrixF+

33、j));</p><p>  *(x+i)=*(x+i)-*(b+i);</p><p>  cout<<*(x+i)<<" ";</p><p><b>  }</b></p><p>  cout<<endl;</p><p>  for

34、(i=0;i<matrixNum;i++)</p><p>  p+=pow(*(b+i),2);</p><p><b>  y++;</b></p><p>  }while(sqrt(p)>x_);</p><p>  cout<<"停止迭代,最終迭代結(jié)果為"<&l

35、t;*x<<','<<*(x+1)<<""<<endl;</p><p>  delete [] matrixF;</p><p>  delete [] matrixF_;</p><p><b>  getch();</b></p><p

36、><b>  }</b></p><p>  double *matrixF2(double *x)</p><p><b>  {</b></p><p><b>  int i,j;</b></p><p>  double t; </p><p&

37、gt;  double *matrixF1; //矩陣F的雅可比矩陣</p><p>  double *matrixF2; //矩陣F的雅可比矩陣的逆 </p><p>  matrixF1=(double *)malloc(matrixNum*matrixNum);</p><p>  matrixF2=(double *)malloc(matrixNu

38、m*matrixNum);</p><p>  for(i=0;i<matrixNum;i++)</p><p>  for(j=0;j<matrixNum;j++)</p><p><b>  if(i==j)</b></p><p>  *(matrixF2+i*matrixNum+j)=1;</p

39、><p>  else *(matrixF2+i*matrixNum+j)=0;</p><p>  *matrixF1=(f0((*x+x_),*(x+1))-f0(*x,*(x+1)))/x_;</p><p>  *(matrixF1+1)=(f0(*x,(*(x+1)+x_))-f0(*x,*(x+1)))/x_;</p><p>  *

40、(matrixF1+2)=(f1((*x+x_),*(x+1))-f1(*x,*(x+1)))/x_;</p><p>  *(matrixF1+3)=(f1(*x,(*(x+1)+x_))-f1(*x,*(x+1)))/x_;</p><p>  //for(i=0;i<matrixNum;i++)</p><p>  // cout<<*(

41、x+i)<<endl; </p><p>  cout<<"矩陣F在["<<*x<<','<<*(x+1)<<"]的雅可比矩陣"<<endl;</p><p>  for(i=0;i<matrixNum;i++) </p>

42、<p><b>  {</b></p><p>  for(j=0;j<matrixNum;j++)</p><p>  cout<<*(matrixF1+i*matrixNum+j)<<" ";</p><p>  cout<<endl;</p><p

43、><b>  }</b></p><p>  //求矩陣F的雅可比矩陣的逆 </p><p>  t=*matrixF1;</p><p>  for(i=0,j=0;j<matrixNum;j++) </p><p><b>  {</b></p>&

44、lt;p>  *(matrixF1+i*matrixNum+j)/=t;</p><p>  *(matrixF2+i*matrixNum+j)/=t;</p><p><b>  }</b></p><p>  t=*(matrixF1+1*matrixNum);</p><p>  for(i=1,j=0;j&

45、lt;matrixNum;j++)</p><p><b>  {</b></p><p>  *(matrixF1+i*matrixNum+j)-=*(matrixF1+j)*t;</p><p>  *(matrixF2+i*matrixNum+j)-=*(matrixF2+j)*t;</p><p><b&g

46、t;  }</b></p><p>  t=*(matrixF1+1*matrixNum+1);</p><p>  for(i=1,j=0;j<matrixNum;j++)</p><p><b>  {</b></p><p>  *(matrixF1+i*matrixNum+j)/=t;</

47、p><p>  *(matrixF2+i*matrixNum+j)/=t;</p><p><b>  } </b></p><p>  t=*(matrixF1+1);</p><p>  for(i=i,j=0;j<matrixNum;j++)</p><p><b>  {

48、</b></p><p>  *(matrixF1+j)-=*(matrixF1+i*matrixNum+j)*t;</p><p>  *(matrixF2+j)-=*(matrixF2+i*matrixNum+j)*t;</p><p><b>  }</b></p><p>  for(i=0;i<

49、;matrixNum;i++) </p><p><b>  {</b></p><p>  for(j=0;j<matrixNum;j++)</p><p>  cout<<*(matrixF1+i*matrixNum+j)<<" ";</p><p&g

50、t;  cout<<endl;</p><p><b>  }</b></p><p>  for(i=0;i<matrixNum;i++) </p><p><b>  {</b></p><p>  for(j=0;j<matrixNum;j++)&l

51、t;/p><p>  cout<<*(matrixF2+i*matrixNum+j)<<" ";</p><p>  cout<<endl;</p><p><b>  }</b></p><p>  cout<<"第"<<y

52、<<"次迭代結(jié)果為"<<*x<<','<<*(x+1)<<""<<endl;</p><p><b>  getch();</b></p><p>  return matrixF2;</p><p>  delete

53、 [] matrixF1;</p><p>  delete [] matrixF2;</p><p><b>  }</b></p><p>  最后總結(jié):我們可以從上面的實(shí)例可以得到,牛頓法是求解非線性方程組最簡(jiǎn)單的一種線性方法,它的構(gòu)想是通過非線性方程組以線性方程組轉(zhuǎn)化,從而來形成一種迭代形式然后迭代達(dá)到迭代次數(shù)來逼近,最終來求解。牛頓法

54、的迭代方式通常都是最少2次迭代以上,并且收斂速度快。因此可以說是最常用的求解非線性方程組的方法</p><p><b>  2.2、擬牛頓法</b></p><p>  2.2.1擬牛頓法的引入與介紹</p><p>  上面我們?cè)敿?xì)介紹了牛頓法求解非線性方程組數(shù)值,我們仔細(xì)留一下,我們會(huì)發(fā)現(xiàn)牛頓法雖然有很好的收斂性,你有沒有發(fā)現(xiàn)牛頓法對(duì)它的初

55、值要求的什么嚴(yán)格,每步迭代都要計(jì)算,是個(gè)偏導(dǎo)數(shù)值建立的矩陣,我們不可能遇到的都是簡(jiǎn)單的數(shù)據(jù),假如我們遇到每個(gè)數(shù)值都很復(fù)雜,這個(gè)時(shí)候我們將無法進(jìn)行計(jì)算。例如n數(shù)值很大時(shí),我們不僅要浪費(fèi)時(shí)間,同時(shí)每步迭代都要求解線性方程組,計(jì)算工作量太大。同時(shí)還有其它問題,假如迭代過程中有一步處有奇異,那么這個(gè)時(shí)候牛頓法將無法計(jì)算。為了克服以上缺點(diǎn)我們下面來介紹擬牛頓法。</p><p>  我們聽到擬牛頓法就知道是對(duì)牛頓法的改進(jìn),

56、例如我們用矩陣來近似的轉(zhuǎn)換代替,這時(shí)我們就可以看到這樣形式的迭代法:</p><p>  k=0,1,2,3,.....................。</p><p>  我們把定義為非奇異的。</p><p>  我們?yōu)榱撕?jiǎn)單化,追求計(jì)算簡(jiǎn)單,我們就不多次計(jì)算逆矩陣,直接定義無限接近的逆矩陣,那么迭代就可以轉(zhuǎn)化為:,。</p><p>

57、  2.2.2擬牛頓法的算法</p><p>  現(xiàn)在我們來說下擬牛頓法的算法過程:</p><p>  首先我們要確定的初值,數(shù)值的精確度:,,并定義初始矩陣為:;</p><p>  其次求解的數(shù)值,假如,那么就令,停止;</p><p><b>  把進(jìn)行迭代計(jì)算;</b></p><p>

58、<b>  在求解;</b></p><p>  ,代入上面的(2),來循環(huán)計(jì)算。</p><p>  下面我為大家介紹一個(gè)實(shí)例來說明我的觀點(diǎn):</p><p>  2.2.3擬牛頓法的題例分析</p><p>  我們首先看下面的非線性方程組:</p><p><b>  定義:,精度

59、。</b></p><p>  對(duì)于上面的非線性方程組我用了擬牛頓法算法的源程序進(jìn)行迭代計(jì)算得到了以下數(shù)據(jù),我用了圖表表(1)表示:</p><p><b>  表1的迭代數(shù)值</b></p><p>  我們通過看上面的迭代結(jié)果可以得出:,,,||||的數(shù)值變化不是很大,精度取到取0.00000001,這是我們看到理論值于迭代值幾

60、乎相同。擬牛頓法可以很簡(jiǎn)單將不管是奇異還是接近于奇異的非線性方程組求解。這就說明了擬牛頓法的計(jì)算量小,不復(fù)雜,同時(shí)收斂性好。</p><p><b>  2.3、割線法</b></p><p>  2.3.1割線法的引入與介紹</p><p>  上面我們已經(jīng)介紹了求解非線性方程組的兩種方法:牛頓法和擬牛頓法?,F(xiàn)在我們來介紹另一種方法割線法,它

61、是通過計(jì)算函數(shù)值的一種迭代方法,這種方法簡(jiǎn)單而且計(jì)算有效方便。</p><p>  下面我們討論一個(gè)非線性方程組:</p><p><b>  (3 .1)</b></p><p>  我們知道建立非線性方程組迭代法,有一個(gè)特別簡(jiǎn)單的方法就是轉(zhuǎn)化成線性方程</p><p><b> ?。? .2)</b&

62、gt;</p><p>  我是通過無限接近(3 .1)求解的??梢钥吹贸鲆粋€(gè)非奇異矩陣,,我們已知F(x)在某點(diǎn)上的值,就可以利用插值得的方法來進(jìn)一步確定方程(3 .2)中的和,從而得到這一類不用到算導(dǎo)數(shù)就可以求數(shù)值的方法,這就是我準(zhǔn)備所要介紹的割線法:</p><p>  我們假設(shè)在中n+1個(gè)上面有互不相同的點(diǎn)(j=0.,1,...n)上所對(duì)應(yīng)的函數(shù)值F(),這是我們可以得到:<

63、/p><p>  , i=0 ,.......n。 (3. 3)</p><p>  那么我們就可以用插值的方法確定和的數(shù)值。</p><p>  此時(shí)再將(3.1)的K次近似解,記作為,同時(shí)取n個(gè)輔助點(diǎn),... ,這時(shí)我們可以從(3.3)求解得</p><p><b>  和</b></p&g

64、t;<p>  ,i=1,.......n,</p><p>  為個(gè)達(dá)到分割,我們令兩式相減得到:</p><p>  ,i=1,.......n, (3.4)</p><p><b>  (3.5)</b></p><p>  可以從(3.4)求解得到:</p>&l

65、t;p><b>  ,</b></p><p>  我們?cè)儆汕蠼猓?.5)求的,</p><p>  將結(jié)果聯(lián)合(3.2)得到</p><p><b>  (3.6)</b></p><p>  通過(3.4)我們可以知道,這時(shí)給出的n+1個(gè)點(diǎn),,只要滿足一定的條件就可以求出,</p&g

66、t;<p>  2.3.2割線法的總結(jié)陳述</p><p>  下面看看我的總結(jié)陳述如下:</p><p>  我們定義為在中如果任何的一個(gè)n+1組成的n個(gè)向量,, i=1,2,...,n ,他們是線性相關(guān),這時(shí)候我們可以把這組點(diǎn){}在一般位置上。</p><p>  1,非奇異矩陣()可以滿足條件;.</p><p>  2

67、,n個(gè)向量,(i)可以組成的一組基</p><p>  3,在 y區(qū)間,都不存在不全為0的常數(shù),或或。</p><p>  現(xiàn)在來闡述一下割線法的幾何意義:</p><p>  我們知道雙點(diǎn)割線法是將過點(diǎn)和的割線與軸交點(diǎn)的橫坐標(biāo)作為方程的根的近似值,可以重復(fù)此過程進(jìn)行運(yùn)算,將和這兩點(diǎn)的割線與軸交點(diǎn)的橫坐標(biāo)理解為方程解的近似值。簡(jiǎn)單的說單點(diǎn)割線法是用過點(diǎn)和的割線與軸交

68、點(diǎn)的橫坐標(biāo)作為方程的根的近似值。</p><p><b>  割線法提例:</b></p><p>  3.3.3割線法題例分析</p><p><b>  用雙點(diǎn)割線法求方程</b></p><p><b>  在區(qū)間上的根,</b></p><p>

69、  在MATLAB命令窗口執(zhí)行。</p><p>  進(jìn)行迭代6次后得到計(jì)算結(jié)果,下圖表2</p><p><b>  表2 </b></p><p>  (3)下面我們來介紹割線法的收斂速度</p><p>  我們以方程為例,假設(shè)它的根為,那么假設(shè)在附近有二階連續(xù)導(dǎo)數(shù),,那么初值就會(huì)無限接近,此時(shí)雙點(diǎn)割線法的迭代

70、過程收斂,收斂速度的計(jì)算方法為:</p><p><b>  這說明,</b></p><p>  由此我們從上面得出結(jié)論,單點(diǎn)割線法時(shí)線性收斂的而雙點(diǎn)割線法是超線性收斂的。</p><p>  我們把輔助點(diǎn)定位,i=1,...,n, 其此時(shí)作為第i個(gè)分量為1的坐標(biāo)向量,是的已知向量,得到</p><p><b&g

71、t;  =</b></p><p><b>  =</b></p><p>  我們?nèi)绻阉涀鳛椋海?lt;/p><p><b>  其中</b></p><p>  此時(shí)這種割線法可以轉(zhuǎn)化為:</p><p><b> ?、?lt;/b></

72、p><p>  我們通過割線法轉(zhuǎn)化得到的也可以說是離牛頓法。</p><p><b>  結(jié)束語(yǔ)</b></p><p>  現(xiàn)在的科學(xué)研究中,面對(duì)很多實(shí)際問題都無法用線性表達(dá)式有規(guī)律的計(jì)算出結(jié)果,而很多問題實(shí)際上都是非線性問題,非線性問題相比較線性問題要麻煩的多,我們常常需要構(gòu)造一個(gè)非線性方程通過對(duì)數(shù)值的研究計(jì)算與探考,求出結(jié)果。然而隨著科學(xué)的發(fā)

73、展,現(xiàn)在非線性的問題已經(jīng)應(yīng)用到在科學(xué)計(jì)算以及工程領(lǐng)域等多個(gè)方面,因此,研究非線性方程對(duì)科學(xué)計(jì)算和工程應(yīng)用等領(lǐng)域有很高的價(jià)值和意義,這也是這篇論文探考求解非線性方程數(shù)值的方法。</p><p>  本文主要研究了非線性方程迭代法的相關(guān)運(yùn)算以及Newton法,主要介紹了求解非線性方程目前比較常用的幾種迭代方法,牛頓法、割線法、。通過對(duì)幾種方法的計(jì)算精度和收斂性等作出比較,而得出結(jié)論,在這一章還介紹了一種求解非線性方程

74、的新的迭代法,相比其他傳統(tǒng)的迭代法,新的迭代法有其迭代收斂速度更快、精度更高等特點(diǎn)</p><p>  最后我再次對(duì)自己隨寫的論文坐下總結(jié),對(duì)于非線性方程組的求解數(shù)值的方法,論文介紹了牛頓法,擬牛頓法,割線法這三種最常見的方法,并把其算法步棸進(jìn)行了闡述,形象簡(jiǎn)單,容易讓人接受。</p><p><b>  參考文獻(xiàn)</b></p><p>  

75、[1] 徐萃薇、孫繩武、計(jì)算方法引論 。北京:高教出版社 2001.</p><p>  [2] 曾金平、數(shù)值計(jì)算方法。長(zhǎng)沙:湖南大學(xué)出版社 2004. </p><p>  [3]曲建明、求解非線性方程的拋物線迭代。</p><p>  [4] 林成森、數(shù)值計(jì)算方法。北京:科學(xué)出版社,1988.</p><p>  學(xué)出版社,2005、

76、</p><p>  [5]李慶樣、王能超、易大義。北京:清華大學(xué)出版社,2001.8</p><p>  [6]關(guān)冶,陸金蒲,數(shù)值分析基礎(chǔ)。高等教育出版社 1988.5</p><p>  [7]鄧建中,葛仁杰,程正興,計(jì)算方法,西安交通大學(xué)</p><p>  [8]王則柯,計(jì)算的負(fù)責(zé)性,湖南教育出版社</p><p

77、><b>  致謝詞</b></p><p>  本人2個(gè)月的時(shí)間終于將這篇論文寫完,在論文的寫作過程中遇到了無數(shù)的困難,都在同學(xué)和老師的幫助下度過了。尤其要特別感謝我的論文指導(dǎo)老師—禹海雄老師,她對(duì)我進(jìn)行了無私的指導(dǎo)和幫助,不厭其煩的幫助進(jìn)行論文的修改和改進(jìn)。另外,在校圖書館查找資料的時(shí)候,里面有我需要的各種資料和體例。在此向幫助和指導(dǎo)過我的各位老師表示最衷心的感謝!</p&g

78、t;<p>  感謝這篇論文所涉及到的各位學(xué)者。本文引用了數(shù)位學(xué)者的研究文獻(xiàn),如果沒有各位學(xué)者的研究成果的幫助和啟發(fā),我將很難完成本篇論文的寫作。感謝我的同學(xué)和朋友,在我寫論文的過程中給予我了很多參考素材,還在論文的撰寫和排版等過程中提供熱情的幫助。</p><p>  由于我的學(xué)術(shù)水平有限,所寫論文難免有不足之處,懇請(qǐng)各位老師和學(xué)友批評(píng)和指正同時(shí)衷心地感謝在百忙之中評(píng)閱論文和參加答辯的各位專家、教

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 眾賞文庫(kù)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論