數(shù)值分析實驗課

實驗報告求解方程迭代法實驗（數(shù)值計算方法，2學時）一實驗目的掌握方程迭代法思想，會用Aitken法、牛頓法、弦截法求解方程。二實驗內(nèi)容1用下列方法求方程X3-3X-1=0在％=2附近的根，要求準確到四位有效數(shù)字。（1）牛頓法。（2）單點弦截法（應數(shù)）。（3）雙點弦截法（信科）。2用Aitken法求方程x3-x2-1=0在x°=1.5附近的根，精度要求為￡=10-4。三實驗步驟（算法）與結果1（1）實驗步驟：輸出孔C語言代碼:#include"stdio.h”#include"math.h”

#defineNUM100intmain(){doublex[NUM];doublep,f;inti;x[0]=2;p=pow(10,-4);printf("thenumberis%5.4f\n",x[0]);for(i=1;i<20;i++){x[i]=x[i-1]-(pow(x[i-1],3)-3*x[i-1]-1)/(3*pow(x[i-1],2)-3);if(x[i]>x[i-1])f=x[i]-x[i-1];elsef=x[i-1]-x[i];printf("thenumberis%5.4f\n”,x[i]);if(f<p){break;}}getch();return0;}運算結果：*?QMJ5BFskA￡G「\DeMktQE逑隹方析作業(yè)'受途果匕上⑴Nlewtar^eKe'ilnumboF±?2.0300il我皿人觸卜Isilnumberis1=8795numbei*is1.87^4

(2)C語言代碼：#include"stdio.h”#include"math.h”#defineNUM100intmain(){doublex[NUM];doublep,q,f;inti;x[0]=2;x[1]=1.5;f=pow(x[0],3)-3*x[0]-1;p=pow(10,-4);printf("thenumberis%5.4f\n”,x[0]);for(i=1;i<10;i++){x[i+1]=(x[0]*(pow(x[i],3)-3*x[i]-1)-x[i]*f)/(pow(x[i],3)-3*x[i]-1-f);if(x[i]>x[i-1])

q=x[i]-x[i-1];elseq=x[i-1]-x[i];printf("thenumberis%5.4f\n”,x[i+1]);if(q<p){break;}}getch();return0;}運算步驟：C:\User&\Acer\De如叩\數(shù) 作業(yè)\實發(fā)誤11住J5ecantexetliethetlieLliethethetlicnunbermmbernLinberriLinbernunbernunbernLinbarLSISLSLSLSLSZ.MUUt)1.84001.S7601.87911.87V41.87941.8794,口回?將曲=響解■外=機心由兌出擔，城代兩次盼三個相牝的迭代廿沃u，H-zi-Kn)fl=Ol'￡f5h=f『+Il'￡f5輸出七C語言代碼:#include"stdio.h”#include"math.h”#defineNUM100intmain(){doublex[NUM];doubley[NUM];doublez[NUM];doublep,f;inti;x[0]=1.5;y[1]=pow(x[0],3)-pow(x[0],2)+x[0]-1;z[1]=pow(y[1],3)-pow(y[1],2)+y[1]-1;p=pow(10,-4);for(i=1;i<10;i++){x[i]=(x[i-1]*z[i]-pow(y[i],2))/(x[i-1]-2*y[i]+z[i]);y[i+1]=pow(x[i],3)-pow(x[i],2)+x[i]-1;z[i+1]=pow(y[i+1],3)-pow(y[i+1],2)+y[i+1]-1;printf("thenumberis%5.4f\n”,x[i]);if(x[i]>x[i-1])f=x[i]-x[i-1];elsef=x[i-1]-x[i];if(f<p){break;}}getch();return0;}運算結果：四實驗收獲與教師評語收獲：學會通過編程，利用Aitken,Newton,Secant三種迭代解法，計算方程的解，有實際應用