→比较梯形积分法法和基于插值的方法: E >>x0=0: pi/10: pi; yO=sin(xo) I1=trapz(xo, yO) I2=quadspln(xO, yo, 0, pi) 给定5个不均匀分布的样点: 虚>xO=[0,0.4,12,pi];yo=sin(x0); plot(xo, yo, xo, yo,2o) I=quadspln(xo, yO, 0, pi) I1=trapz(xo, yo) 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 比较梯形积分法法和基于插值的方法： 给定5个不均匀分布的采样点：
样条插值的结果与理论之间的比较: >>x=[0:0.01:pi,pi];yoa=sin(x); y=interp1(xo, yo, x, ',) plot(xo, yo, x, y,':',x, yOa, x0, yo,'0,) 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 样条插值的结果与理论之间的比较：
例 8.5 →已知其中的150个数据点,用 quads1n( 函数计算出该定积分的值 3丌/2 cos(15.)da 0 求数值解: >>x=[0:3*p立/2/200:3*pi/2]; y=cs(15*x); I=quadspln(x, y, 0, 3*pi/2) 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 例 8.5 已知其中的150个数据点，用quadspln() 函数计算出该定积分的值 求数值解：
绘制曲线: e>x0=[0:3*pi/2/1000:3*pi/2]; yO=cos(15*x0); y1=interp1(x, y, xO, ',) plot(x,y, xo, y1,':') 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 绘制曲线：
81.3二维网格数据的插值问题 二维插值的函数: 21interp2(Co, y0, 20, 31, 91, method) 其中,x0,yo,x为已知二维数据,x1,y为 插值点构成的新的网格参数,a1矩阵为在 所选插值网格点处的函数近似值 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 8.1.3 二维网格数据的插值问题 二维插值的函数： 其中， 为已知二维数据， 为 插值点构成的新的网格参数， 矩阵为在 所选插值网格点处的函数近似值