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

计算方法实验报告 年级班级： 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,试由测量数据统计该国国土的近似面积，与它的精确值41288km^2作比较。 图（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 on plot(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 on yy2=interp1(x,y2,xx,spline); plot(xx,yy2,r) hold on syms x; s=int(-0.0107*x^2+1.8693*x+41.1473-0.0035*x^2+0.4648*x-51.3154,7,158); s=vpa(s); s=s/18^2*40^2 s =43819.988724279838764947279042147 3、下表给出的x，y数据位于机翼剖面的轮廓线上，y1和y2分别对应轮廓的上下线。假设需要得到x坐标每改变0.1时的y坐标。试完成加工所需数据，画出曲线，求机翼剖面的面积。 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 表 （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=spline(x,y1,b); m2=spline(x,y2,b); plot(b,m1,b,m2); m3=m1-m2; trapz(b,m3); ans= 11.3444