非線性方程求解算法的程序設(shè)計(jì)及比對(duì)課程設(shè)計(jì)畢業(yè)設(shè)計(jì)-資料下載頁
2025-06-04 22:47本頁面
【正文】 [1] 林成森 .數(shù)值計(jì)算方法 .北京 :科學(xué)出版社 .2021年 [2] 金聰 .熊盛武 .數(shù)值分析 .武漢 :武漢理工大學(xué)出版社 .2021年 [3] 顏暉 .C語言程序設(shè)計(jì) .北京 :高等教育出版社 .2021年 [4] 張立科 .MATLAB 應(yīng)用 .北京：人民郵電出版社 .2021 年 [5] 賈得彬 .數(shù)值計(jì)算方法 .北京：水利水電出版社 .2021 年 [6] 李慶揚(yáng) ,王能超 ,易大義 .數(shù)值分析基礎(chǔ) .北京：清華大學(xué)出版 ,2021年 . [7] 徐士良 .數(shù)值方法與計(jì)算機(jī)實(shí)現(xiàn) .北京：清華大學(xué)出版社 .2021年 . 17 附錄 A 本附錄包含不動(dòng)點(diǎn)迭代的流程圖 . 10)( xx ?? 開 始 輸 入 x0,? ,N 1?k k+1?k x1? x0 輸出近似 根 x1 |x1 x0|ε ? 輸出迭代 失敗標(biāo)志 結(jié) 束 n kN? y n y 18 附錄 B 本附錄 含二分法 ,不動(dòng)點(diǎn)迭代 ,牛頓法 ,弦截法的的 C語言計(jì)算程序 . 二分法 C語言程序 include define eps /* 容許誤差 */ double jisuan(double x) { return pow(x,4)+2*x*xx3 。 } main() { double a,b,y,x,h,j。 int k=0。 printf(a=)。 scanf(%lf,amp。a)。 printf(b=)。 scanf(%lf,amp。b)。 if(jisuan(a)*jisuan(b)=0) /* 判斷 是否符合二分法使用的條件 */ { printf(Not bisect of bisect)。 return。 } do { x=(a+b)/2。 k++。 j=jisuan(x)。/* f(x)的結(jié)果 */ if(jisuan(a)*jisuan(x)0) /* 如果 f(a)*f(x)0,則根在區(qū)間的左半部分 */ b=x。 else /* 否則根在區(qū)間的右半部分 */ a=x。 h=fabs(ba)。 printf(k=%2d , a=%lf , b=%lf , x=%lf , h=%lf , j=%lf\n,k,a,b,x,h,j)。 }while(fabs(ba)eps)。/*判斷是否達(dá)到精度要求 ,若沒有達(dá)到 ,繼續(xù)循環(huán) */ x=(a+b)/2。 /* 取最后的小區(qū)間中點(diǎn)作為根的近似值 */ printf(\n The zui hou root is x=%lf, k=%d\n,x,k)。 } 不動(dòng)點(diǎn)迭代法 C語言程序 include define eps /* 容許誤差 */ double jisuan(double x) { double h 。 h=sqrt(sqrt(x+4)1)。 19 return h 。 } main() { double x0,x1,h。 int k=0。 printf(input 迭代初值 x1: )。 scanf(%lf,amp。x1)。 do { x0=x1。 x1=jisuan(x0)。 k++。 h=fabs(x1x0)。 printf(\n k=%2d , x1=%lf ,h=%lf ,k ,x1,h)。 }while(fabs(x1x0)eps)。/*判斷是否達(dá)到精度要求 ,若沒有達(dá)到 ,繼續(xù)循環(huán) */ printf(\n最后結(jié)果 The root is \n)。 printf(x=%lf, k=%d\n,x1,k)。 } 牛頓迭代法 C語言程序 include include double js1(double x) { double j。 j=pow(x,4)+2*pow(x,2)x3。 return j。 } double js2(double c) { double j。 j=4*pow(c,3)+4*c1。 return j 。 } void main() { double x1,x,f,f1,h。 int k=0 。 printf(chu zhi x :)。 scanf(%lf,amp。x1)。 do { x=x1。 k++。 f=js1(x1)。 20 f1=js2(x1)。 x1=xf/f1。 h=fabs(x1x)。 printf(k=%2d , x1=%lf , f=% , f1=% , h=%lf\n,k,x1,f,f1,h)。 }while(fabs(x1x)1e4)。 printf(最后結(jié)果 \n)。 printf(k=%2d,x=%lf\n,k,x1)。 } 牛頓迭代法 C語言程序 include define eps /* 容許誤差 */ define N 100 /* 最大迭代次數(shù) N */ double f(double x) /* 定義函數(shù) f(x) */ { double y。 y=pow(x,4)+2*pow(x,2)x3。 return(y)。 } main() { double x0,x1,x2,h。 int i。 printf(input x0,x1=)。 scanf(%lf,%lf,amp。x0,amp。x1)。 for(i=1。i=N。i++) { x2=x1(f(x1)*(x1x0))/(f(x1)f(x0))。 /* 弦截法迭代公式 */ if(fabs(x2x1)eps || fabs(f(x2))eps) /*滿足精度要求輸出近似根并退出 */ { h=fabs(x2x1)。 printf(\nRoot of equation is\n)。 printf(k=%2d , x=% , h=%lf\n,i,x2,h)。 return。 } h=fabs(x2x1)。 printf(\nk=%2d , x=% , h=%lf\n,i,x2,h)。 x0=x1。 /* 準(zhǔn)備下一次迭代的初值 */ x1=x2。 } printf(\nAfter %d repeat, no solved.\n,N)。 /* 輸出無解信息 */ }
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1