《实验3 插值与数值积分new》由会员分享，可在线阅读，更多相关《实验3 插值与数值积分new（8页珍藏版）》请在金锄头文库上搜索。
1、数学实验作业郭明钊 2012011880 化21一、机翼剖面积问题1、 问题分析：求剖面积，可以根据上下两边界线分别与坐标轴围成的面积之差来求。先插值，然后再用梯形公式求积分。2、matlab编写如下x1=0 3 5 7 9 11 12 13 14 15;y1=0 1.8 2.2 2.7 3.0 3.1 2.9 2.5 2.0 1.6;x2=0 3 5 7 9 11 12 13 14 15;y2=0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6;xa=0:0.1:15;%按0.1的间隔在x方向产生插值点ya=spline(x1,y1,xa);%在y方向计算三次样条插值
2、xb=0:0.1:15;yb=spline(x2,y2,xb);xa ya yb %输出插值后的结果subplot(1,2,1),plot(x1,y1,x2,y2);subplot(1,2,2),plot(xa,ya,xb,yb);s1=trapz(xa,ya) %上线与坐标轴围成的面积s2=trapz(xb,yb)%下线与坐标轴围成的面积s=s1-s2 %机翼的剖面积3、 运行结果如下插值前后图像对比插值后数组结果输出（加工所需数据）ans = 0 0 0 0.1000 0.1089 0.0499 0.2000 0.2134 0.0990 0.3000 0.3137 0.1474 0.400
3、0 0.4097 0.1951 0.5000 0.5018 0.2421 0.6000 0.5898 0.2884 0.7000 0.6740 0.3340 0.8000 0.7545 0.3788 0.9000 0.8314 0.4230 1.0000 0.9047 0.4665 1.1000 0.9747 0.5094 1.2000 1.0413 0.5515 1.3000 1.1047 0.5930 1.4000 1.1651 0.6338 1.5000 1.2225 0.6739 1.6000 1.2770 0.7134 1.7000 1.3287 0.7523 1.8000 1.37
4、78 0.7904 1.9000 1.4244 0.8280 2.0000 1.4685 0.8649 2.1000 1.5104 0.9012 2.2000 1.5499 0.9368 2.3000 1.5874 0.9719 2.4000 1.6229 1.0063 2.5000 1.6565 1.0401 2.6000 1.6884 1.0732 2.7000 1.7185 1.1058 2.8000 1.7471 1.1378 2.9000 1.7742 1.1692 3.0000 1.8000 1.2000 3.1000 1.8245 1.2302 3.2000 1.8480 1.2
5、599 3.3000 1.8704 1.2889 3.4000 1.8918 1.3174 3.5000 1.9125 1.3454 3.6000 1.9325 1.3727 3.7000 1.9519 1.3995 3.8000 1.9708 1.4258 3.9000 1.9894 1.4515 4.0000 2.0076 1.4767 4.1000 2.0258 1.5014 4.2000 2.0439 1.5255 4.3000 2.0620 1.5491 4.4000 2.0803 1.5722 4.5000 2.0989 1.5947 4.6000 2.1179 1.6168 4.
6、7000 2.1374 1.6383 4.8000 2.1575 1.6594 4.9000 2.1784 1.6799 5.0000 2.2000 1.7000 5.1000 2.2225 1.7196 5.2000 2.2459 1.7387 5.3000 2.2700 1.7573 5.4000 2.2948 1.7754 5.5000 2.3201 1.7930 5.6000 2.3459 1.8102 5.7000 2.3720 1.8269 5.8000 2.3984 1.8430 5.9000 2.4249 1.8588 6.0000 2.4515 1.8740 6.1000 2
7、.4781 1.8887 6.2000 2.5045 1.9030 6.3000 2.5307 1.9168 6.4000 2.5566 1.9301 6.5000 2.5821 1.9430 6.6000 2.6071 1.9553 6.7000 2.6315 1.9672 6.8000 2.6552 1.9786 6.9000 2.6780 1.9895 7.0000 2.7000 2.0000 7.1000 2.7210 2.0100 7.2000 2.7411 2.0195 7.3000 2.7602 2.0285 7.4000 2.7786 2.0370 7.5000 2.7961
8、2.0450 7.6000 2.8130 2.0525 7.7000 2.8291 2.0595 7.8000 2.8446 2.0660 7.9000 2.8595 2.0719 8.0000 2.8739 2.0773 8.1000 2.8878 2.0822 8.2000 2.9013 2.0865 8.3000 2.9144 2.0902 8.4000 2.9272 2.0933 8.5000 2.9397 2.0959 8.6000 2.9520 2.0979 8.7000 2.9641 2.0994 8.8000 2.9761 2.1002 8.9000 2.9881 2.1004
9、 9.0000 3.0000 2.1000 9.1000 3.0119 2.0990 9.2000 3.0238 2.0974 9.3000 3.0355 2.0952 9.4000 3.0469 2.0925 9.5000 3.0578 2.0893 9.6000 3.0683 2.0857 9.7000 3.0782 2.0815 9.8000 3.0873 2.0770 9.9000 3.0956 2.0721 10.0000 3.1029 2.0668 10.1000 3.1092 2.0611 10.2000 3.1143 2.0552 10.3000 3.1181 2.0490 1
10、0.4000 3.1206 2.0425 10.5000 3.1215 2.0358 10.6000 3.1209 2.0289 10.7000 3.1185 2.0219 10.8000 3.1143 2.0147 10.9000 3.1082 2.0074 11.0000 3.1000 2.0000 11.1000 3.0897 1.9924 11.2000 3.0772 1.9841 11.3000 3.0626 1.9742 11.4000 3.0459 1.9621 11.5000 3.0269 1.9469 11.6000 3.0059 1.9280 11.7000 2.9826
11、1.9046 11.8000 2.9573 1.8759 11.9000 2.9297 1.8413 12.0000 2.9000 1.8000 12.1000 2.8682 1.7516 12.2000 2.8342 1.6970 12.3000 2.7984 1.6377 12.4000 2.7606 1.5749 12.5000 2.7211 1.5099 12.6000 2.6798 1.4442 12.7000 2.6370 1.3790 12.8000 2.5927 1.3157 12.9000 2.5470 1.2556 13.0000 2.5000 1.2000 13.1000 2.4518 1.1501 13.2000 2.4026 1.1063 13.3000 2.3527 1.0687 13.4000 2.3021 1.0377 13.5000 2.2513 1.0134 13.6000 2.2004 0.9960 13.7000 2.1496
