计算方法实验报告二答案许晶.doc

计算方法实验报告年级班级： 2013级一班学 号： 201335010110姓 名： 许 晶日 期： 2015年11月24日计算方法实验（二）实验目的：1、利用MATLAB插值命令进行计算2、利用MATLAB数值积分命令进行计算实验内容：1、从1点12点的11小时内，每隔1小时测量一次温度，测得的温度的数值依次为：5，8，9，15，25，29，31，30，22，25，27，24试估计每隔1/20小时的温度值并画出图形。答:hours=1:12;temps=5 8 9 15 25 29 31 30 22 25 27 24;h=1:0.05:12;t=interp1(hours,temps,h,spline);plot(hours,temps,+,h,t,hours,temps,r:)2、图（1）是欧洲一个国家的地图，为计算出它的国土面积，首先对该国地图做出如下测量：以由西向东方向为x轴，由南向北方向为y轴，选择方便的原点，并从最西边界点到最东边界点在x轴上的区间适当的划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2, 得到表格1（单位mm）根据地图的比例尺我们知道18mm相当于40km,试由测量数据统计该国国土的近似面积，与它的精确值41288km2作比较。 图（1）答： 表（1）x=7.0 10.5 13.0 17.5 34.0 40.5 44.5 48.0 56.0 61.0 68.5 76.5 80.5 91.0 96.0 101.0 104.0 106.5 111.5 118.0 123.5 136.5 142.0 146.0 150.0 157.0 158.0;y1=44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68;y2=44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68;plot(x,y1,o)hold onplot(x,y2,rp)hold on xi=10:10:100;yi=interp1(x,y1,xi,spline);yi=interp1(x,y2,xi,spline); xx=min(x):0.1:max(x);yy1=interp1(x,y1,xx,spline);plot(xx,yy1) hold onyy2=interp1(x,y2,xx,spline);plot(xx,yy2,r)hold onsyms x;s=int(-0.0107*x2+1.8693*x+41.1473-0.0035*x2+0.4648*x-51.3154,7,158); s=vpa(s);s=s/182*402s =43819.9887242798387649472790421473、下表给出的x，y数据位于机翼剖面的轮廓线上，y1和y2分别对应轮廓的上下线。假设需要得到x坐标每改变0.1时的y坐标。试完成加工所需数据，画出曲线，求机翼剖面的面积。X035791112131415Y101.82.22.73.03.12.92.52.01.6Y201.21.72.02.12.01.81.21.01.6表 （2）答：x=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 ;y2=0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6;b=0:0.1:15;m1=