數(shù)值分析作業(yè)(C語(yǔ)言編程實(shí)現(xiàn))電子教案

1、數(shù)值分析作業(yè)（C語(yǔ)言編程實(shí)現(xiàn)）學(xué)習(xí)好資料L(10分)用Romberg方法計(jì)算枳分二去的近似值.誤差不超過(guò)101#include#includedoublef(doublex)doubleans;ans=exp(x);returnans;voidmain()doublea=1,b=3,error=0.0001,t2020,h,c;inti,j,k,m,n;h=b-a;t00=h*(f(a)+f(b)/2;k=1;while(1)t0k=0;m=1;for(j=0;jk-1;j+)m=m*2;for(i=1;i=m;i+)t0k=t0k+h*f(a+(i-0.5)*h);tk=(tk+t0k-1)

2、/2;for(j=1;j=k;j+)c=1;for(n=0;nj;n+)c=c*4;tjk-j=(c*tj-1k-j+1-tj-1k-j)/(c-1);if(fabs(tk0-tk-10)error)printf(n積分結(jié)果I%lfn,tk0);break;elseh=hk+;積分結(jié)果Ia17.3672SE2. CO分)用簡(jiǎn)單迭代法求解方程Fy加尸必要求心用#include#includedoublef(doublet)doubleans;ans=pow(cos(t),1.0/3);returnans;voidmain()doublex=0,eslong=0.000001,x0;intN=20

3、,i;printf(n近似初值x0=%lfn,x);for(i=0;iN;i+)x0=x;x=f(x);printf(x%d=%lfn,i+1,x);if(fabs(x-x0)eslong)break;if(fabs(x-x0)eslong)printf(得到近似結(jié)果為x%lfnn,x,i);elseprintf(迭代失敗n);近似初值xQ=0.R000QBxl=i.aaeoaex2=0.814477=B.SS2B51x4=O.8&978Bx5=0.8&7394x6-B.B64B22x7-a.Bt5695xS=B.8&5399/9=S.8&54?9=0.865465xll-0.865477x1

4、2=0.8S5473x13=0-Sg5474x14=3.S65474得到近似結(jié)果為x=0.8654743. （10分）用四階RK方法求初值問(wèn)題的數(shù)值解.取步長(zhǎng)0.1.汁算結(jié)果取8位小數(shù).y=l+/0x|xo）=o#include#includedoublea=0,b=1,x,y=0,h=0.1,k1,k2,k3,k4;inti,N;doublef(doublet,doubles)doubleans;ans=1+t*t;returnans;voidmain()N=(b-a)/h;x=a;printf(n初值為(x0,y0)=(%.8f,%.8f)n,x,y);for(i=0;i1次輸出結(jié)果為 選

5、符果力 出結(jié)臬為第6次輸出結(jié)果為第9次匕結(jié)果為 第8次狗出結(jié)臬為 第9次送7孤里人我)C 0.10000000 .C 0.20030000 ,0-30000000 ,C目回國(guó)團(tuán)M ,0.10033333 0.2026666?)0.30900M0Q 0.42133333 30.50B00000 , 0.41Gt67 ( 0.60000000 . 9.7200000 ? =1.0000000081,33333333;4、（10分）已知函數(shù)表姐F，用三次插旗求吊卜5）的近似值.修L22.94,65.8fr.7此KM433.715丸3679.249&79#includevoidmain()double

6、datax4=1.2,2.9,4.6,5.8,datay10=14.84,33.71,58.36,79.24,l3,x=1.5,y;inti,j;y=0;for(i=0;i=3;i+)li=1;for(j=0;ji;j+)li=(x-dataxj)/(dataxi-dataxj)*li;for(j=i+1;j=3;j+)li=(x-dataxj)/(dataxi-dataxj)*li;y=y+datayi*li;printf(nf(x)在x=%f處的近似值為：y=%fn,x,y);R2在X-1.500000處的近似值為；V-17.750000精品資料5.。力分已知數(shù)據(jù)表如L試求其最小二乘:次擬

7、合名頂式即L223A.44.55.66.77X910.0IL714.872L44313g44.7361.4681.57105.11BL-91#includevoidmain()doubledatay9=11.7,14.87,21.44,31.39,44.73,61.46,81.57,105.11,131.91;intm=2,i,j,k;doublep,data94,a34,datax9=1.2,2.3,3.4,4.5,5.6,6.7,7.8,8.9,10.0;for(i=0;i9;i+)for(j=1;j2*m+1;j+)dataij=1;for(k=0;kj;k+)dataij=dataij

8、*dataxi;for(i=0;im+1;i+)for(j=0;jm+1;j+)aij=0;for(k=0;k9;k+)aij=aij+dataki+j;)a00=9;a0m+1=0;for(i=0;i9;i+)a0m+1=a0m+1+datayi;for(i=1;im+1;i+)aim+1=0;for(j=0;j9;j+)p=datayj;for(k=0;ki;k+)p=p*dataxj;aim+1=aim+1+p;)/代成m+1行，m+2列增廣矩陣/for(i=0;im+1;i+)/顯示方程組/for(j=0;jm+2;j+)if(j!=m+1)printf(%f)a%d,aij,j);i

9、f(j!=m)printf(+);elseprintf(=%fn,aij);for(i=0;im;i+)/高斯消去法/if(aii!=0)for(j=i+1;jm+1;j+)aji=aji/aii;for(k=i+1;km+2;k+)ajk=ajk-aik*aji;elsebreak;if(amm!=0&i=m)amm+1=amm+1/amm;for(i=2;i=m+1;i+)for(j=1;ji;j+)am+1-im+1=am+1-im+1-am+1-im+1-j*am+1-jm+1;am+1-im+1=am+1-im+1/am+1-im+1-i;printf(方程組的解為：n);for(j

10、=0;ja.0+224B00a1*C33548方程組的解為:居-12.076541kl=-2.H00515a2=1,3987H7擬合多項(xiàng)式為二P2(x)=+C-2-000515x*x26、(20分)分別用列主元法和Ja8bi迭代法求解方程組.再-2不一工=3.2x1410_r.3%-15一上i-2x2+5工3=10列主元素法#include#includevoidmain()doublea34=1,-2,-13-2,10,-3,15,-1,-2,5,10,mov,comp;inti,j,k,nrow;for(i=0;i2;i+)comp=fabs(aii);for(k=i;k=comp)nro

11、w=k;comp=fabs(aki);for(j=0;j=3;j+)mov=a皿;aij=anrowj；anrowj=mov;printf(方程第d行互換位置后如下n,i+1);for(j=0;j3;j+)printf(%f)x1+(%f)x2+(%f)x3=%fn,aj0,aj1,aj2,aj3);if(aii!=0)for(j=i+1;j3;j+)aji=aji/aii；for(k=i+1;k=3;k+)a皿k=ajk-aik*a皿i；aji=0；printf(方程經(jīng)d次消元如下n,i+1);for(j=0;j3;j+)printf(%f)x1+(%f)x2+(%f)x3=%fn,aj0,

12、aj1,aj2,aj3);elsebreak;if(a22!=0&i=2)printf（方程化簡(jiǎn)得n）;for(i=0;i3;i+)printf(%f)x1+(%f)x2+(%f)x3=%fn,ai0,ai1,ai2,ai3);a23=a23/a22;for(i=2;i=3;i+)for(j=1;ji;j+)a3-i3=a3-i3-a3-i3-j*a3-j3;a3-i3=a3-i3/a3-i3-i;printf(方程組的解為：n);for(j=0;jx2+=3,0C0OS0(-1.000000)x14x2*-10.060000片程經(jīng)該消元如下xl+x2*=10.5G0OOBxl+x2+C6,G

13、00000x3-2.500600片程第2行互換位置后如下.2000009x1+10.090000x2*3-060000x3=151600000x2*k3=2.500066(0.000000xl+(3.009900x2*x3-10.500000方程經(jīng)2次消元如F-2-080009)x1+M2+C3-000000?x3=15.0OGOSS(0,008008)x1+(-?.0O0006x2+C6.500600xJ=2.500008(0.00000a)xl+(0.O009Oax2+0.2857147x3-11-57142質(zhì)程化簡(jiǎn)得C-2.060000)x1+-C10.0030S0x2+C-3-回回回團(tuán)

14、)箕3=15.008866xl+C-7.0O0006x2+三2.50S6880.00阿目ONI+(血酶的%M2+0-285714x3=11.571429方程組的解為:xl=118.9B0909x2=37.250000x3-40.509303Jacobi迭代法#include#includevoidmain()doublea37=1,-2,-1,3,-2,10,-3,15,-1,-2,5,10,error=0.000001,norm;intN=423,i,j,k;a04=0,a14=0,a24=0;for(i=0;i3;i+)把a(bǔ)矩陣轉(zhuǎn)化為b矩陣ai6=aii；for(j=0;j3;j+)aij

15、=-aij/ai6;ai3=ai3/ai6;aii=0;printf(化為b矩陣如下n);for(i=0;i3;i+)printf(%f%f%f%fn,ai0,ai1,ai2,ai3);for(i=1;iN;i+)for(j=0;j3;j+)a皿5=0;for(k=0;k3;k+)aj5=ak4*a皿k+aj5;a皿5=a皿5+aj3;norm=0;for(k=0;k3;k+)norm=norm+fabs(ak4-ak5);if(normerror)break;elsefor(k=0;k3;k+)ak4=ak5;if(normerror)printf（計(jì)算結(jié)果為n）;for(i=0;i3;i+

16、)printf(x%d=%.3fn”,i+1,ai5);elseprintf(迭代失敗n);K為隨陣如下0.6000392.0S000S1-0000803.Q8BS0Se.2000000.0000080.3600801.500000計(jì)算結(jié)果為xl=118.009x2=37.253x3=4930現(xiàn)場(chǎng)考試題目1#includestdio.h#includemath.hdoublef(doublex)doubleans;ans=exp(x);return(ans);voidmain()doublea=-1,b=1,error=0.0001,m=1,h,T0,T,F;intk;h=(b-a)/2;T0=h*(f(a)+f(b);while(1)F=0;for(k=1;k=pow(2.0,m-1);k+)F=F+f(a+(2*k-1)*h);T=T0/2+h*F;if(fabs(T-T0)error)break;m+;h=h/2;T0=T;printf(積分結(jié)果為I=%fn,T);題目2#includestdio.hdoublef(doubl