這篇文章給大家分享的是有關(guān)C語言實現(xiàn)牛頓迭代法解方程的示例的內(nèi)容。小編覺得挺實用的，因此分享給大家做個參考，一起跟隨小編過來看看吧。

創(chuàng)新互聯(lián)2013年開創(chuàng)至今，是專業(yè)互聯(lián)網(wǎng)技術(shù)服務(wù)公司，擁有項目做網(wǎng)站、成都網(wǎng)站制作網(wǎng)站策劃，項目實施與項目整合能力。我們以讓每一個夢想脫穎而出為使命，1280元鼓樓做網(wǎng)站,已為上家服務(wù),為鼓樓各地企業(yè)和個人服務(wù),聯(lián)系電話:18982081108

 C語言實現(xiàn)牛頓迭代法解方程詳解

利用迭代算法解決問題，需要做好以下三個方面的工作：

一、確定迭代變量

在可以用迭代算法解決的問題中，我們可以確定至少存在一個可直接或間接地不斷由舊值遞推出新值的變量，這個變量就是迭代變量。

二、建立迭代關(guān)系式

所謂迭代關(guān)系式，指如何從變量的前一個值推出其下一個值的公式（或關(guān)系）。迭代關(guān)系式的建立是解決迭代問題的關(guān)鍵，通常可以使用遞推或倒推的方法來完成。

三、對迭代過程進(jìn)行控制

在什么時候結(jié)束迭代過程？這是編寫迭代程序必須考慮的問題。不能讓迭代過程無休止地執(zhí)行下去。迭代過程的控制通?？煞譃閮煞N情況：一種是所需的迭代次數(shù)是個確定的值，可以計算出來；另一種是所需的迭代次數(shù)無法確定。對于前一種情況，可以構(gòu)建一個固定次數(shù)的循環(huán)來實現(xiàn)對迭代過程的控制；對于后一種情況，需要進(jìn)一步分析得出可用來結(jié)束迭代過程的條件。

接下來，我介紹一種迭代算法的典型案例----牛頓-拉夫遜（拉弗森）方法

牛頓-拉夫遜（拉弗森）方法，又稱牛頓迭代法，也稱牛頓切線法：先任意設(shè)定一個與真實的根接近的值x0作為第一次近似根，由x0求出f(x0),過(x0，f(x0))點做f(x)的切線，交x軸于x1，把它作為第二次近似根，再由x1求出f(x1),過(x1，f(x1))點做f(x)的切線，交x軸于x2，……如此繼續(xù)下去，直到足夠接近（比如|x- x0|<1e-6時）真正的根x*為止。

而f '(x0)=f(x0)/( x1- x0) 

所以 x1= x0- f(x0)/ f ' (x0)。

我們來看一副從網(wǎng)上找到的圖:

接下來，我們來看一個例子:

我們還是直接上代碼:

例子:用牛頓迭代法求下列方程在值等于2.0附近的根：2x3-4x2+3x-6=0。

#include  
#include   
int main(void) 
{ 
  float x,x0,f,f1;  
  x = 2.0; 
  do{  
      x0=x; 
      f=2*x0*x0*x0-4*x0*x0+3*x0-6; 
      f1=6*x0*x0-8*x0+3; 
      x=x0-f/f1;  
  //函數(shù)fabs：求浮點數(shù)x的絕對值 
  //說明：計算|x|, 當(dāng)x不為負(fù)時返回 x，否則返回 -x    
  }while(fabs(x-x0)>=1e-5); 
   printf ("%f\n",x); 
  return 0 ; 
 }

執(zhí)行結(jié)果:

當(dāng)x=1.5時，方程2x3-4x2+3x-6=0。附近的根為2.000000 。

感謝各位的閱讀！關(guān)于“C語言實現(xiàn)牛頓迭代法解方程的示例”這篇文章就分享到這里了，希望以上內(nèi)容可以對大家有一定的幫助，讓大家可以學(xué)到更多知識，如果覺得文章不錯，可以把它分享出去讓更多的人看到吧！

另外有需要云服務(wù)器可以了解下創(chuàng)新互聯(lián)建站www.cdcxhl.com，海內(nèi)外云服務(wù)器15元起步，三天無理由+7*72小時售后在線，公司持有idc許可證，提供“云服務(wù)器、裸金屬服務(wù)器、高防服務(wù)器、香港服務(wù)器、美國服務(wù)器、虛擬主機、免備案服務(wù)器”等云主機租用服務(wù)以及企業(yè)上云的綜合解決方案，具有“安全穩(wěn)定、簡單易用、服務(wù)可用性高、性價比高”等特點與優(yōu)勢，專為企業(yè)上云打造定制，能夠滿足用戶豐富、多元化的應(yīng)用場景需求。

網(wǎng)站欄目：C語言實現(xiàn)牛頓迭代法解方程的示例-創(chuàng)新互聯(lián)
文章出自：http://weahome.cn/article/cdsidi.html