傻大方摘要：【计算|计算国土面积数值分析课程设计( 二 )|国土面积|数值|分析|课程设计】计算 s1 的关键代码： s10=trapz(x10,y1i)+trapz(x30,y3i)+trapz(x40,y4i) +trapz(x60,y6i); %用复化梯形公式 。 10、求除特殊段以外的各段与 X 轴围成的面积并求和 计算结果...

计算 s1 的关键代码： s10=trapz(x10,y1i)+trapz(x30,y3i)+trapz(x40,y4i) +trapz(x60,y6i);
%用复化梯形公式 。

10、求除特殊段以外的各段与 X 轴围成的面积并求和 计算结果为： 数值分析课程设计论文 9 s1 =7.4597e+004 s2 =2.1257e+005 4计算 S 从以上的过程已经得出 s=s2-s1 ， 还要将 s 按比例转化所得就是国土面积 S 。
s=s2-s1； S=s*100/9； S 计算结果为： S =1.5331e+006（平方公里） 5.计算相对误差 t 相对误差的求解公式为： 0065665 . 1 0065665 . 1 t e Se 计算结果为： t =0.0213 6误差分析 误差的产生主要来源于数据点的个数有点少以及数据点之间不够均匀 ， 这样三次 样条插值后的数据作图就会与 。

11、实际地图和地形相比有较大的误差 。
7分析和总结 由于 t 的值较小 ， 所以以上的计算结果 S=1.5331e+006 平方公里可以作为国土面 积的近似值 。
在以后的不规则图形面积的计算中 ， 此方法可以考虑选择使用 。
由于梯 形公式只有一次代数精确度代数精确度（ （3 3） ， 所以会产生计算的结果不够精确 。
但要很精确求出不 规则图形面积 ， 就必须测量出更多的数据点以及选用代数精确度更高的算法 ， 才能够 更加减少插值和计算上的误差 。
至此 ， 国土面积计算完成 。
数值分析课程设计论文 10 注注释：释： （1）三次样条插值：参照 参考文献2 4650 页 （2）复化梯形公式：参照 参考文献2 9091 页 （3）代 。

12、数精确度：参照 参考文献2 88 页 参考文献：参考文献： 1李玉莉等 ， MATLAB 函数速查手册 ，北京：化学工业出版社 ，2010 2袁东锦 ， 计算方法数值分析 ，南京：南京师范大学出版社 ，2007 3蒲俊 吉家锋 伊良忠 ， MATLAB6.0 数学手册 ， 上海浦东：浦东电子出版社 ，2002 附录附录： 1对附件数据的分段 下边疆: A03:B10;
A11:B13;
%(特殊段 s11) A13:B31;
A32:B39;
A39:B40;
%(特殊段 s12) A40:B56;
上边疆: A03:B11;
A11:B14;
%(特殊段 s21) A14:B40;
A40:B41;
%(特殊段 s 。

13、22) A41:B49;
A50:B52;
A53:B64;
2求国土面积及相对误差的完整代码 A1=xlsread(第一题数据,下边疆,A03:B10);
A2=xlsread(第一题数据,下边疆,A11:B13);
A3=xlsread(第一题数据,下边疆,A13:B31);
A4=xlsread(第一题数据,下边疆,A32:B39);
A5=xlsread(第一题数据,下边疆,A39:B40);
A6=xlsread(第一题数据,下边疆,A40:B56);
数值分析课程设计论文 11 %导入下边疆的实验数据并分好计算的数据段 x1=A1(:,1);
x2=A2(:,1);
x3=A3(:,1 。

14、);
x4=A4(:,1);
x5=A5(:,1);
x6=A6(:,1);
y1=A1(:,2);
y2=A2(:,2);
y3=A3(:,2);
y4=A4(:,2);
y5=A5(:,2);
y6=A6(:,2);
%给相应的变量赋值 x10=min(x1):0.001:max(x1);
%对每一段数据点按连结顺序进行点横坐标的加密处理 x20=max(x2):-0.001:min(x2);
x30=min(x3):0.001:max(x3);
x40=min(x4):0.001:max(x4);
x50=max(x5):-0.001:min(x5);
x60=min(x6):0.001:max(x 。

15、6);
y1i=interp1(x1,y1,x10,spline);
%对数据进行三次样条插值 y2i=interp1(x2,y2,x20,spline);
y3i=interp1(x3,y3,x30,spline);
y4i=interp1(x4,y4,x40,spline);
y5i=interp1(x5,y5,x50,spline);
y6i=interp1(x6,y6,x60,spline);
x=x10 x20 x30 x40 x50 x60;
y=y1i y2i y3i y4i y5i y6i;
%对三次样条插值后的数据按连结顺序合并 s10=trapz(x10,y1i)+trapz 。

【计算|计算国土面积数值分析课程设计】16、(x30,y3i)+trapz(x40,y4i) +trapz(x60,y6i);
%用梯形公式求除特殊段以外的各段与 X 轴围成的面积并求和 s11=-trapz(x20,y2i);
%对特殊段面积 s11 的计算 s12=-trapz(x50,y5i);
%对特殊段面积 s12 的计算 s1=s10-s11-s12;
%计算下疆界与 X 轴围成面积的精确值 plot(x,y) hold on;
A1=xlsread(第一题数据,上边疆,A03:B11);
A2=xlsread(第一题数据,上边疆,A11:B14);
A3=xlsread(第一题数据,上边疆,A14:B40);
A4=xlsrea 。

17、d(第一题数据,上边疆,A40:B41);
A5=xlsread(第一题数据,上边疆,A41:B49);
A6=xlsread(第一题数据,上边疆,A50:B52);
A7=xlsread(第一题数据,上边疆,A53:B64);
%导入上边疆的实验数据并分好计算的数据段 x1=A1(:,1);
x2=A2(:,1);
x3=A3(:,1);
x4=A4(:,1);
x5=A5(:,1);
x6=A6(:,1);
x7=A7(:,1);
y1=A1(:,2);
y2=A2(:,2);
y3=A3(:,2);
y4=A4(:,2);
y5=A5(:,2);
y6=A6(:,2);
y7=A7(:,2);
%给相应的变 。

18、量赋值 x10=min(x1):0.001:max(x1);
%对每一段数据点按连结顺序进行点横坐标的加密处理 x20=max(x2):-0.001:min(x2);
x30=min(x3):0.001:max(x3);
x40=max(x4):-0.001:min(x4);

稿源：(未知)

【傻大方】网址：/a/2021/0707/0022751417.html

标题：计算|计算国土面积数值分析课程设计( 二 )