第十三次作业解答.doc

插值与拟合作业 1.已知一个地区的边界点数据如下表所列,试估算该地区的边界线长及近似面积. 边界点数据表 7.0 10.5 13.0 17.5 34.0 40.5 44.5 48.0 56.0 44 45 47 50 50 38 30 30 34 44 59 70 72 93 100 110 110 110 61.0 68.5 76.5 80.5 91.0 96.0 101.0 104.0 106.5 36 34 41 45 46 43 37 33 28 117 118 116 118 118 121 124 121 121 111.5 118.0 123.5 136.5 142.0 146.0 150.0 157.0 158.0 32 65 55 54 52 50 66 66 68 121 122 116 83 81 82 86 85 68 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,x,x,y2,*); hold on x1=7.0:0.5:158; y11=spline(x,y1,x1); plot(x1,y11,.); hold on y21=spline(x,y2,x1); plot(x1,y21,.); hold on p=polyfit(x1,y11,26); x11=7.0:0.1:158; yi=polyval(p,x11); plot(x11,yi,ro); hold on p1=polyfit(x1,y21,26); yii=polyval(p1,x11); plot(x11,yii,ro); hold on t=sym(t); y1i=p.*[t^26,t^25,t^24,t^23,t^22,t^21,t^20,t^19,t^18,t^17,t^16,t^15,t^14,t^13,t^12,t^11,t^10,t^9,t^8,t^7,t^6,t^5,t^4,t^3,t^2,t^1,1,]; f1=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ]*y1i.; s1=int(f1,7,158); y2i=p1.*[t^26,t^25,t^24,t^23,t^22,t^21,t^20,t^19,t^18,t^17,t^16,t^15,t^14,t^13,t^12,t^11,t^10,t^9,t^8,t^7,t^6,t^5,t^4,t^3,t^2,t^1,1,]; f2=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ]*y2i.; s2=int(f2,7,158); s=abs(s1-s2); A=eval(s) ff1=sqrt(1+(diff(f1,t))^2); ff2=sqrt(1+(diff(f2,t))^2); L=int(ff1,7,158)+int(ff2,7,158); t1=7.0:0.1:158; yt1=spline(x,y1,t1); yt2=spline(x,y2,t1); L1=trapz(t1,sqrt(1+gradient(yt1,t1).^2))+trapz(t1,sqrt(1+gradient(yt2,t1).^2)); double(L1) 2.许多供水厂由于没有测量流入或流出水塔流量的设备，而只能测量水塔中的水位，试通过测得的某时刻水塔中水位的数据，估计在任意时刻（包括泵灌水期间）流出水塔的流量.下表给出了水塔水位与时间的原始数据，其中时间用小时(h),高度用米(m)，水塔为圆柱体，其直径为1