北航數(shù)值分析B第二題報(bào)告完整版資料

1、數(shù)學(xué)分析B 計(jì)算實(shí)習(xí)題目第二題 1 算法設(shè)計(jì)方案1. 總體思路：由于需要四種方法進(jìn)行細(xì)化，所以大體需要四個(gè)模塊，實(shí)際上采用在main函數(shù)里面進(jìn)行聲明和調(diào)用幾個(gè)子函數(shù)，f1(),f2(),f3(),分別是線性插值，二次插值和三次插值，prepare()函數(shù)用于求關(guān)于Mi的線性方程的參數(shù)，chase()函數(shù)用于采用追趕法求出Mi，three_moment()函數(shù)用于求出三彎矩法的y表達(dá)式，求出各種方法下的y的表達(dá)式后，進(jìn)而x就容易求了，利用final_data()函數(shù)將數(shù)據(jù)輸出，在matlab下畫出細(xì)化后的輪廓線。2. 線性插值：插值法均采用牛頓法，d用于存放一階差商值，這個(gè)值不僅在線性插值中要用

2、，在二次插值與三次插值中都需要它進(jìn)行計(jì)算，同理，二次插值中的二階差商也要用作三次插值中的計(jì)算。線性插值算出一階差商后基本就可以代入公式得到y(tǒng)的表達(dá)式了。3. 二次插值：每三個(gè)點(diǎn)為一個(gè)插值區(qū)域，共六個(gè)插值區(qū)域，利用差商表計(jì)算得出二階差商，展開牛頓二次插值多項(xiàng)式，得出結(jié)果。4. 三次插值：每四個(gè)點(diǎn)為一個(gè)插值區(qū)域，共四個(gè)插值區(qū)域，這時(shí)計(jì)算量就已經(jīng)有點(diǎn)大了，為簡化運(yùn)算，采用不單獨(dú)計(jì)算，而是在二次插值的基礎(chǔ)上進(jìn)行計(jì)算，利用二次插值的已有數(shù)據(jù)，但是有兩個(gè)二次插值在二次插值函數(shù)中并沒有得到計(jì)算，即f(x3,x4,x5)和f(x7,x8,x9)，所以在二次插值函數(shù)中要添加相關(guān)代碼。5. 樣條插值：計(jì)算最為復(fù)雜

3、，但其實(shí)并不比三次插值容易錯(cuò)，因?yàn)橹皇枪胶芏?，所以化為三個(gè)函數(shù)來完成，每個(gè)函數(shù)負(fù)責(zé)一部分計(jì)算。6. 輸出數(shù)據(jù)：前面的計(jì)算只是算出公式，對(duì)于不同的插值法，y的分段表達(dá)式不一樣，采用if，else if語句以及for循環(huán)并用到了pow()函數(shù)算冪，完成對(duì)表中數(shù)據(jù)x從1到190的細(xì)化。2 全部源程序#include#includedouble x13=0.00,4.74,9.50,19.00,38.00,57.00,76.00,95.00,114.00,133.00,152.00,171.00,190.00;double y13=0.00,5.32,8.10,11.97,16.15,17.10,1

4、6.34,14.63,12.16,9.69,7.03,3.99,0.00;double d12=0;/存放一階差商double da12=0;/存放線性插值常數(shù)項(xiàng)double dt8=0;/存放二階差商double dta8=0;/存放二階插值x項(xiàng)系數(shù)double dtb8=0;/存放二階插值常數(shù)項(xiàng)double ddt4=0;/存放二階差商用作計(jì)算三階差商double dth4=0;/存放三階差商double dtha4=0;/存放三階插值附加x2項(xiàng)系數(shù)double dthb4=0;/存放三階插值附加x項(xiàng)系數(shù)double dthc4=0;/存放三階插值附加常數(shù)項(xiàng)double h13=0;/存放

5、樣條插值的一些參數(shù)double a13=0;double b13=0;double c13=0;double e13=0;double p13=0;double q13=0;double z13=0;double m13=0;double ta12=0;double tb12=0;double tc12=0;double td12=0;double ans1190=0;double ans2190=0;double ans3190=0;double ans4190=0;void main()/調(diào)用幾個(gè)子程序void f1(); /線性插值void f2();/二次插值void f3();/三次

6、插值void prepare();/樣條插值計(jì)算關(guān)于Mn的線性方程組的有關(guān)參數(shù)void chase();/追趕法解線性方程組void three_moment();/三彎矩法void final_data();f1();f2();f3();prepare();chase();three_moment(); final_data(); void f1()/線性插值函數(shù)int i=0;printf(分段線性插值法：nn);for(i=0;i12;i+)di=(yi+1-yi)/(xi+1-xi);dai=yi-di*xi;printf(插值區(qū)間：(%6.2f,%6.2f)t ,xi,xi+1);p

7、rintf(插值函數(shù)：P1(x)=%6.2fx+%6.2fn ,di,dai);for(i=0;i=190;i+)/計(jì)算具體數(shù)據(jù)if(i5)ans1i=d0*i+da0;else if(i10)ans1i=d1*i+da1;else if(i20)ans1i=d2*i+da2;else if(i39)ans1i=d3*i+da3;else if(i58)ans1i=d4*i+da4;else if(i77)ans1i=d5*i+da5;else if(i96)ans1i=d6*i+da6;else if(i115)ans1i=d7*i+da7;else if(i134)ans1i=d8*i+d

8、a8;else if(i153)ans1i=d9*i+da9;else if(i172)ans1i=d10*i+da10;else if(i191)ans1i=d11*i+da11;void f2()/二次插值int i=0;printf(nn);printf(二次插值法;nn);for(i=0;i6;i+)dti=(d2*i+1-d2*i)/(x2*i+2-x2*i);dtai=d2*i-dti*(x2*i+x2*i+1);dtbi=y2*i-d2*i*x2*i+dti*x2*i*x2*i+1; printf(插值區(qū)間：(%6.2f,%6.2f) ,x2*i,x2*i+2);printf(插

9、值函數(shù)：P1(x)=(%6.5fx2)+(%6.3fx)+(%.2f)n ,dti,dtai,dtbi); dt6=(d4-d3)/(x5-x3);dta6=d3-dt6*(x3+x4);dtb6=y3-d3*x3+dt6*x3*x4;/計(jì)算以x3x4x5插值的函數(shù)，為三次插值計(jì)算做準(zhǔn)備dt7=(d10-d9)/(x11-x9);dta7=d9-dt7*(x9+x10);dtb7=y9-d9*x9+dt7*x9*x10;/計(jì)算以x7x8x9插值的函數(shù)，為三次插值計(jì)算做準(zhǔn)備for(i=0;i=190;i+)/計(jì)算具體數(shù)據(jù)if(i10)ans2i=dt0*pow(i,2)+dta0*i+dtb0;

10、else if(i39)ans2i=dt1*pow(i,2)+dta1*i+dtb1;else if(i77)ans2i=dt2*pow(i,2)+dta2*i+dtb2;else if(i115)ans2i=dt3*pow(i,2)+dta3*i+dtb3;else if(i153)ans2i=dt4*pow(i,2)+dta4*i+dtb4;else if(i191)ans2i=dt5*pow(i,2)+dta5*i+dtb5;void f3()/三次插值int i=0;printf(nn); printf(三次插值法;nn);ddt0=(d2-d1)/(x3-x1);ddt1=(d4-d

11、3)/(x5-x3);ddt2=(d8-d7)/(x9-x7);ddt3=(d10-d9)/(x11-x9);dth0=(ddt0-dt0)/(x3-x0);dth1=(dt2-ddt1)/(x6-x3);dth2=(ddt2-dt3)/(x9-x6);dth3=(dt5-ddt3)/(x12-x9);for(i=0;i4;i+)dthai=(x3*i+x3*i+1+x3*i+2)*dthi;dthbi=(x3*i*x3*i+1+x3*i*x3*i+2+x3*i+1*x3*i+2)*dthi;dthci=(x3*i*x3*i+1*x3*i+2)*dthi; printf(插值區(qū)間：(%6.2f

12、,%6.2f)n ,x0,x3);printf(插值函數(shù)：P1(x)=(%.6fx3)+(%.2fx2)+(%.2fx)+(%.2f)n ,dth0,dt0-dtha0,dta0+dthb0,dtb0-dthc0); printf(插值區(qū)間：(%6.2f,%6.2f)n ,x3,x6); printf(插值函數(shù)：P1(x)=(%.6fx3)+(%.2fx2)+(%.2fx)+(%.2f)n ,dth1,dt6-dtha1,dta6+dthb1,dtb6-dthc1); printf(插值區(qū)間：(%6.2f,%6.2f)n ,x6,x9); printf(插值函數(shù)：P1(x)=(%.6fx3)+

13、(%.2fx2)+(%.2fx)+(%.2f)n ,dth2,dt3-dtha2,dta3+dthb2,dtb3-dthc2);printf(插值區(qū)間：(%6.2f,%6.2f)n ,x9,x12); printf(插值函數(shù)：P1(x)=(%.6fx3)+(%.2fx2)+(%.2fx)+(%.2f)n ,dth3,dt7-dtha3,dta7+dthb3,dtb7-dthc3);for(i=0;i=190;i+)/計(jì)算具體數(shù)據(jù)if(i20)ans3i=dth0*pow(i,3)+(dt0-dtha0)*pow(i,2)+(dta0+dthb0)*i+dtb0-dthc0;else if(i7

14、7)ans3i=dth1*pow(i,3)+(dt6-dtha1)*pow(i,2)+(dta6+dthb1)*i+dtb6-dthc1;else if(i134)ans3i=dth2*pow(i,3)+(dt3-dtha2)*pow(i,2)+(dta3+dthb2)*i+dtb3-dthc2;else if(i191)ans3i=dth3*pow(i,3)+(dt7-dtha3)*pow(i,2)+(dta7+dthb3)*i+dtb7-dthc3;void prepare()/樣條插值計(jì)算關(guān)于Mn的線性方程組的有關(guān)參數(shù)int i=0;int j=0;for(i=1;i13;i+) hi=

15、xi-xi-1;for(j=1;j=11;j+)aj=hj+1/(hj+hj+1);cj=1-aj;bj=(yj+1-yj)/hj+1-(yj-yj-1)/hj)*6/(hj+hj+1);/printf(a=%.3f,b=%.3f,c=%.3fn,aj,bj,cj);c12=0;b0=0;b12=0;a0=0;void chase()/追趕法解線性方程組(注意此模塊中a對(duì)應(yīng)追趕法中的c；c對(duì)應(yīng)追趕法中的d)int i=0;int j=0;p0=2;z0=b0/p0;for(i=0;i=0;j-) mj=zj-qj*mj+1; /printf(m=%.5f,mj);void three_mome

16、nt()/三彎矩法int i=0;int j=0;printf(nn);printf(三次樣條插值結(jié)果：);for(j=1;j=12;j+)taj-1=yj-1/hj-mj-1*hj/6;tbj=yj/hj-mj*hj/6;tcj-1=mj-1/(6*hj);tdj=mj/(6*hj);for(i=1;i=12;i+)printf(插值區(qū)間：(%6.2f,%6.2f)n ,xi-1,xi);printf(插值函數(shù)s(x)=%.3f(%.2f-x)3+%.3f(x-%.2f)3+%.3f(%.2f-x)+%.3f(x-%.2f)n,tcj-1,xi,tdj,xi-1,tai-1,xi,tbi,x

17、i-1);for(i=0;i=190;i+)/計(jì)算具體數(shù)據(jù)if(i5) ans4i=tc0*pow(x1-i,3)+td1*pow(i-x1,3)+ta0*(x1-i)+tb1*(i-x0); else if(i10)ans4i=tc1*pow(x2-i,3)+td2*pow(i-x1,3)+ta1*(x2-i)+tb2*(i-x1);else if(i20)ans4i=tc2*pow(x3-i,3)+td3*pow(i-x2,3)+ta2*(x3-i)+tb3*(i-x2);else if(i39)ans4i=tc3*pow(x4-i,3)+td4*pow(i-x3,3)+ta3*(x4-i

18、)+tb4*(i-x3);else if(i58)ans4i=tc4*pow(x5-i,3)+td5*pow(i-x4,3)+ta4*(x5-i)+tb5*(i-x4);else if(i77) ans4i=tc5*pow(x6-i,3)+td6*pow(i-x5,3)+ta5*(x6-i)+tb6*(i-x5);else if(i96)ans4i=tc6*pow(x7-i,3)+td7*pow(i-x6,3)+ta6*(x7-i)+tb7*(i-x6);else if(i115)ans4i=tc7*pow(x8-i,3)+td8*pow(i-x7,3)+ta7*(x8-i)+tb8*(i-x

19、7);else if(i134)ans4i=tc8*pow(x9-i,3)+td9*pow(i-x8,3)+ta8*(x9-i)+tb9*(i-x8);else if(i153)ans4i=tc9*pow(x10-i,3)+td10*pow(i-x9,3)+ta9*(x10-i)+tb10*(i-x9);else if(i172)ans4i=tc10*pow(x11-i,3)+td11*pow(i-x10,3)+ta10*(x11-i)+tb11*(i-x10);else if(i191)ans4i=tc11*pow(x12-i,3)+td12*pow(i-x11,3)+ta11*(x12-i

20、)+tb12*(i-x11);void final_data()/輸出最終數(shù)據(jù)int j=0;printf(y值線性插值二次插值三次插值樣條插值n);printf(x值n);for(j=0;j=190;j+)printf(%d%.2f %.2f %.2f %.2f n,j,ans1j,ans2j,ans3j,ans4j);三全部數(shù)據(jù)分段線性插值法：插值區(qū)間：( 0.00, 4.74) 插值函數(shù)：P1(x)= 1.12x+ 0.00 插值區(qū)間：( 4.74, 9.50) 插值函數(shù)：P1(x)= 0.58x+ 2.55 插值區(qū)間：( 9.50, 19.00) 插值函數(shù)：P1(x)= 0.41x+

21、4.23 插值區(qū)間：( 19.00, 38.00) 插值函數(shù)：P1(x)= 0.22x+ 7.79 插值區(qū)間：( 38.00, 57.00) 插值函數(shù)：P1(x)= 0.05x+ 14.25 插值區(qū)間：( 57.00, 76.00) 插值函數(shù)：P1(x)= -0.04x+ 19.38 插值區(qū)間：( 76.00, 95.00) 插值函數(shù)：P1(x)= -0.09x+ 23.18 插值區(qū)間：( 95.00,114.00) 插值函數(shù)：P1(x)= -0.13x+ 26.98 插值區(qū)間：(114.00,133.00) 插值函數(shù)：P1(x)= -0.13x+ 26.98 插值區(qū)間：(133.00,152

22、.00) 插值函數(shù)：P1(x)= -0.14x+ 28.31 插值區(qū)間：(152.00,171.00) 插值函數(shù)：P1(x)= -0.16x+ 31.35 插值區(qū)間：(171.00,190.00) 插值函數(shù)：P1(x)= -0.21x+ 39.90二次插值法;插值區(qū)間：( 0.00, 9.50) 插值函數(shù)：P1(x)=(-0.05667x2)+( 1.391x)+(0.00) 插值區(qū)間：( 9.50, 38.00) 插值函數(shù)：P1(x)=(-0.00657x2)+( 0.595x)+(3.04) 插值區(qū)間：( 38.00, 76.00) 插值函數(shù)：P1(x)=(-0.00237x2)+( 0.

23、275x)+(9.12) 插值區(qū)間：( 76.00,114.00) 插值函數(shù)：P1(x)=(-0.00105x2)+( 0.090x)+(15.58) 插值區(qū)間：(114.00,152.00) 插值函數(shù)：P1(x)=(-0.00026x2)+(-0.065x)+(22.99) 插值區(qū)間：(152.00,190.00) 插值函數(shù)：P1(x)=(-0.00132x2)+( 0.265x)+(-2.85)三次插值法;插值區(qū)間：( 0.00, 19.00) 插值函數(shù)：P1(x)=(0.002330x3)+(-0.09x2)+(1.50x)+(0.00) 插值區(qū)間：( 19.00, 76.00) 插值函

24、數(shù)：P1(x)=(0.000037x3)+(-0.01x2)+(0.62x)+(3.04) 插值區(qū)間：( 76.00,133.00) 插值函數(shù)：P1(x)=(0.000018x3)+(-0.01x2)+(0.58x)+(0.38) 插值區(qū)間：(133.00,190.00) 插值函數(shù)：P1(x)=(-0.000014x3)+(0.01x2)+(-0.94x)+(65.55)三次樣條插值結(jié)果：插值區(qū)間：( 0.00, 4.74) 插值函數(shù)s(x)=0.000(4.74-x)3+1.122(x-0.00)3+0.000(4.74-x)+1.256(x-0.00)插值區(qū)間：( 4.74, 9.50)

25、插值函數(shù)s(x)=0.000(9.50-x)3+1.122(x-4.74)3+1.252(9.50-x)+1.704(x-4.74)插值區(qū)間：( 9.50, 19.00) 插值函數(shù)s(x)=0.000(19.00-x)3+1.122(x-9.50)3+0.858(19.00-x)+1.286(x-9.50)插值區(qū)間：( 19.00, 38.00) 插值函數(shù)s(x)=0.000(38.00-x)3+1.122(x-19.00)3+0.682(38.00-x)+0.876(x-19.00)插值區(qū)間：( 38.00, 57.00) 插值函數(shù)s(x)=0.000(57.00-x)3+1.122(x-3

26、8.00)3+0.876(57.00-x)+0.914(x-38.00)插值區(qū)間：( 57.00, 76.00) 插值函數(shù)s(x)=0.000(76.00-x)3+1.122(x-57.00)3+0.914(76.00-x)+0.867(x-57.00)插值區(qū)間：( 76.00, 95.00) 插值函數(shù)s(x)=0.000(95.00-x)3+1.122(x-76.00)3+0.867(95.00-x)+0.779(x-76.00)插值區(qū)間：( 95.00,114.00) 插值函數(shù)s(x)=0.000(114.00-x)3+1.122(x-95.00)3+0.779(114.00-x)+0.6

27、37(x-95.00)插值區(qū)間：(114.00,133.00) 插值函數(shù)s(x)=0.000(133.00-x)3+1.122(x-114.00)3+0.637(133.00-x)+0.513(x-114.00)插值區(qū)間：(133.00,152.00) 插值函數(shù)s(x)=0.000(152.00-x)3+1.122(x-133.00)3+0.513(152.00-x)+0.371(x-133.00)插值區(qū)間：(152.00,171.00) 插值函數(shù)s(x)=0.000(171.00-x)3+1.122(x-152.00)3+0.371(171.00-x)+0.222(x-152.00)插值區(qū)間

28、：(171.00,190.00) 插值函數(shù)s(x)=0.000(190.00-x)3+1.122(x-171.00)3+0.222(190.00-x)+0.000(x-171.00)y值 線性插值 二次插值 三次插值 樣條插值x值0 0.00 0.00 0.00 0.001 1.12 1.33 1.41 1.572 2.24 2.56 2.65 2.633 3.37 3.66 3.74 3.804 4.49 4.66 4.70 5.035 5.47 5.54 5.52 5.546 6.06 6.31 6.24 6.277 6.64 6.96 6.87 6.898 7.22 7.50 7.41

29、7.419 7.81 7.93 7.88 7.8810 8.30 8.33 8.30 8.3211 8.71 8.79 8.68 8.7612 9.12 9.23 9.04 9.2013 9.53 9.66 9.38 9.6214 9.93 10.08 9.73 10.0415 10.34 10.49 10.09 10.4516 10.75 10.88 10.48 10.8517 11.16 11.25 10.91 11.2418 11.56 11.62 11.41 11.6119 11.97 11.97 11.97 11.9720 12.19 12.31 12.30 12.3121 12.4

30、1 12.63 12.61 12.6422 12.63 12.95 12.91 12.9523 12.85 13.24 13.19 13.2424 13.07 13.53 13.47 13.5225 13.29 13.80 13.73 13.7926 13.51 14.06 13.98 14.0427 13.73 14.31 14.22 14.2828 13.95 14.54 14.45 14.5129 14.17 14.76 14.67 14.7230 14.39 14.97 14.87 14.9231 14.61 15.16 15.07 15.1132 14.83 15.34 15.25

31、15.2933 15.05 15.51 15.43 15.4634 15.27 15.66 15.59 15.6235 15.49 15.81 15.74 15.7636 15.71 15.93 15.89 15.9037 15.93 16.05 16.02 16.0338 16.15 16.15 16.15 16.1539 16.20 16.24 16.27 16.2640 16.25 16.33 16.38 16.3741 16.30 16.41 16.48 16.4642 16.35 16.49 16.57 16.5543 16.40 16.57 16.65 16.6344 16.45

32、16.63 16.73 16.7145 16.50 16.70 16.80 16.7746 16.55 16.76 16.86 16.8347 16.60 16.81 16.91 16.8948 16.65 16.86 16.96 16.9349 16.70 16.91 17.00 16.9850 16.75 16.95 17.03 17.0151 16.80 16.98 17.06 17.0452 16.85 17.02 17.08 17.0653 16.90 17.04 17.09 17.0854 16.95 17.06 17.10 17.0955 17.00 17.08 17.11 17

33、.1056 17.05 17.09 17.11 17.1057 17.10 17.10 17.10 17.1058 17.06 17.10 17.09 17.0959 17.02 17.10 17.07 17.0860 16.98 17.09 17.05 17.0761 16.94 17.08 17.03 17.0562 16.90 17.07 17.00 17.0263 16.86 17.04 16.97 17.0064 16.82 17.02 16.94 16.9665 16.78 16.99 16.90 16.9366 16.74 16.95 16.86 16.8967 16.70 16

34、.91 16.82 16.8568 16.66 16.87 16.77 16.8069 16.62 16.82 16.72 16.7570 16.58 16.76 16.67 16.7071 16.54 16.71 16.62 16.6572 16.50 16.64 16.57 16.5973 16.46 16.57 16.51 16.5374 16.42 16.50 16.46 16.4775 16.38 16.42 16.40 16.4176 16.34 16.34 16.34 16.3477 16.25 16.27 16.28 16.2778 16.16 16.20 16.22 16.2

35、079 16.07 16.12 16.15 16.1380 15.98 16.04 16.08 16.0581 15.89 15.96 16.01 15.9782 15.80 15.88 15.93 15.9083 15.71 15.80 15.85 15.8184 15.62 15.71 15.76 15.7385 15.53 15.62 15.67 15.6486 15.44 15.53 15.58 15.5587 15.35 15.44 15.49 15.4688 15.26 15.35 15.39 15.3789 15.17 15.25 15.29 15.2790 15.08 15.1

36、5 15.18 15.1791 14.99 15.05 15.08 15.0792 14.90 14.95 14.97 14.9693 14.81 14.85 14.86 14.8594 14.72 14.74 14.75 14.7495 14.63 14.63 14.63 14.6396 14.50 14.52 14.51 14.5197 14.37 14.41 14.39 14.3998 14.24 14.29 14.27 14.2799 14.11 14.17 14.15 14.15100 13.98 14.05 14.02 14.02101 13.85 13.93 13.90 13.90102 13.72 13.81 13.77 13.77103 13.59 13.68 13.64 13.64104 13.46 13.55 13.51 13.50105 13.33 13.42 13.38 13.37106 13.20 13.29 13.24 13.24107 13.07 13.16 13.11 13.10108 12.94 13.02 12.98 12.9710