非線性方程求解算法的程序設計及比對課程設計畢業(yè)設計(參考版)
2025-01-21 17:37本頁面
【正文】 } printf(\nAfter %d repeat, no solved.\n,N)。 x0=x1。 } h=fabs(x2x1)。 printf(k=%2d , x=% , h=%lf\n,i,x2,h)。 /* 弦截法迭代公式 */ if(fabs(x2x1)eps || fabs(f(x2))eps) /*滿足精度要求輸出近似根并退出*/ { h=fabs(x2x1)。i=N。x1)。 scanf(%lf,%lf,amp。 int i。 return(y)。 }牛頓迭代法C語言程序include define eps /* 容許誤差 */define N 100 /* 最大迭代次數(shù)N */double f(double x) /* 定義函數(shù)f(x) */{ double y。 printf(最后結果\n)。 printf(k=%2d , x1=%lf , f=% , f1=% , h=%lf\n,k,x1,f,f1,h)。 x1=xf/f1。 f=js1(x1)。 do { x=x1。 scanf(%lf,amp。 int k=0 。 return j 。 } double js2(double c) { double j。 j=pow(x,4)+2*pow(x,2)x3。 printf(x=%lf, k=%d\n,x1,k)。 }while(fabs(x1x0)eps)。 h=fabs(x1x0)。 x1=jisuan(x0)。x1)。 printf(input 迭代初值 x1: )。 }main(){ double x0,x1,h。 h=sqrt(sqrt(x+4)1)。 /* 取最后的小區(qū)間中點作為根的近似值 */ printf(\n The zui hou root is x=%lf, k=%d\n,x,k)。 }while(fabs(ba)eps)。 h=fabs(ba)。/* f(x)的結果 */ if(jisuan(a)*jisuan(x)0) /* 如果f(a)*f(x)0,則根在區(qū)間的左半部分 */ b=x。 k++。 return。b)。 printf(b=)。 scanf(%lf,amp。 int k=0。 x0y|x1 x0|ε?ynkN?n輸出近似根x1輸出迭代失敗標志結 束 附錄 B本附錄 含二分法,不動點迭代,牛頓法,弦截法的的C語言計算程序.二分法C語言程序include define eps /* 容許誤差 */double jisuan(double x){ return pow(x,4)+2*x*xx3 。k k+1222。迭代法是一種逐次逼近的方法,起著把根的精確值一步一步算出來的作用。 計算迭代公式 下面用二分法的程序來（程序詳見附錄B的二分法程序）計算（21）式的在的根,誤差度在,初始值計算把結果寫在表 21表 21 二分法的計算結果備注11. 2. 234567891011121314最后的結果 ：= ,由表 21二分法的最后結果與比較,我們可以看出與Mathematica 5求得結果在誤差范圍之內(nèi)幾乎一致,所以此算法是正確的,是可以應用到下一小節(jié),進行算法比較.下面用不動點迭代法的C語言的程序來（程序詳見附錄B不動點迭代法C語言程序）計算（22）式的在的根,誤差度在,初始迭代值為,計算把結果寫在表 ,來得出算法的正確與否.表 22 不動點迭代的計算結果備注11. 2. 3. 初始值2345最后結果: = , k=5由表 22不動點迭代法的最后結果與比較,我們可以看出與Mathematica 5求得
環(huán)評公示相關推薦
鄂ICP備17016276號-1