#incl">
數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)報(bào)告new

數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)報(bào)告new

ID:17478944

大小:75.50 KB

頁(yè)數(shù):7頁(yè)

時(shí)間:2018-09-02

數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)報(bào)告new_第1頁(yè)
數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)報(bào)告new_第2頁(yè)
數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)報(bào)告new_第3頁(yè)
數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)報(bào)告new_第4頁(yè)
數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)報(bào)告new_第5頁(yè)
資源描述:

《數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)報(bào)告new》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。

1、實(shí)習(xí)報(bào)告課程名稱:數(shù)值計(jì)算方法上機(jī)學(xué)院:電氣信息工程學(xué)院姓名:李勇學(xué)號(hào):31005010632011年12月30日數(shù)值計(jì)算方法上機(jī)習(xí)題1.用Newton迭代法求解方程在2.0附近的實(shí)根,。源程序:#include#includemain(){doublex,f,f1;scanf("%lf",&x);for(;;x=x-f/f1){f=xxx-3x-1;f1=3xx-3;if(fabs(f)<0.5e-5)break;}printf("r=%lf",x);getch();r

2、eturn0;}輸入:2.0運(yùn)行結(jié)果:r=1.8793852.設(shè)(1)(2)試用高斯-約旦消去法求方程組的解。(1)源程序:#include#includemain(){floata[5][5]={0,0,0,0,0,0,5,7,9,10,0,6,8,10,9,0,7,10,8,7,0,5,7,6,5},b[5]={0,1,1,1,1},s,t,am;inti,j,l,k;clrscr();for(k=1;k<=4;k++){am=fabs(a[k][k]);l=k;for(

3、i=k;i<=4;i++)if(fabs(a[i][k])>am){am=fabs(a[i][k]);l=i;}if(a[l][k]==0){printf("qiyijuzhen.");break;}if(l!=k){for(j=k;j<=4;j++){t=a[k][j];a[k][j]=a[l][j];a[l][j]=t;}s=b[l];b[l]=b[k];b[k]=s;}for(j=k+1;j<=4;j++)a[k][j]=a[k][j]/a[k][k];b[k]=b[k]/a[k][k];for(i=1

4、;i<=4;i++){if(i==k)continue;b[i]=b[i]-a[i][k]b[k];for(j=k+1;j<=4;j++)a[i][j]=a[i][j]-a[i][k]a[k][j];}}printf("x1=%f",b[1]);printf("x2=%f",b[2]);printf("x3=%f",b[3]);printf("x4=%f",b[4]);}運(yùn)行結(jié)果:x1=19.999809x2=-11.999882x3=-4.999953x4=2.999972(1)源程序:#i

5、nclude#includemain(){floata[4][4]={0,0,0,0,0,2,1,-3,0,3,6,2,0,4,5,-3},b[4]={0,-5,9,5},s,t,am;inti,j,l,k;clrscr();for(k=1;k<=3;k++){am=fabs(a[k][k]);l=k;for(i=k;i<=3;i++)if(fabs(a[i][k])>am){am=fabs(a[i][k]);l=i;}if(a[l][k]==0){printf("qiyij

6、uzhen.");break;}if(l!=k){for(j=k;j<=3;j++){t=a[k][j];a[k][j]=a[l][j];a[l][j]=t;}s=b[l];b[l]=b[k];b[k]=s;}for(j=k+1;j<=3;j++)a[k][j]=a[k][j]/a[k][k];b[k]=b[k]/a[k][k];for(i=1;i<=3;i++){if(i==k)continue;b[i]=b[i]-a[i][k]b[k];for(j=k+1;j<=3;j++)a[i][j]=a[i][j]

7、-a[i][k]a[k][j];}}printf("x1=%f",b[1]);printf("x2=%f",b[2]);printf("x3=%f",b[3]);}運(yùn)行結(jié)果:x1=-11.000001x2=8.000000x3=-3.0000003.已知函數(shù)的數(shù)表xi0.320.340.360.38Sinxi0.31456660.33348710.35227420.3709205試用拉格朗日插值法求時(shí)的函數(shù)值。源程序:#includevoidmain(){floatx[4]={0

8、.32,0.34,0.36,0.38},y[4]={0.3145666,0.3334871,0.3522742,0.3709205};floatu,v,l;inti,j;clrscr();scanf("%f",&u);v=0;for(j=0;j<=3;j++){l=1;for(i=0;i<=3;i++){if(i==j)continue;l=l(u-x[i])/(x[j]-x[i]);}v=v+ly

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動(dòng)畫的文件,查看預(yù)覽時(shí)可能會(huì)顯示錯(cuò)亂或異常,文件下載后無此問題,請(qǐng)放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫(kù)負(fù)責(zé)整理代發(fā)布。如果您對(duì)本文檔版權(quán)有爭(zhēng)議請(qǐng)及時(shí)聯(lián)系客服。
3. 下載前請(qǐng)仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時(shí)可能由于網(wǎng)絡(luò)波動(dòng)等原因無法下載或下載錯(cuò)誤,付費(fèi)完成后未能成功下載的用戶請(qǐng)聯(lián)系客服處理。