基3_Matlab背景基础之三 数值分析与其MATLAB实现.pptVIP

0 1.0000 1.0000 1.0000 0.5000 0.8000 0.7500 0.8205 1.0000 0.5000 0.5000 0.5000 1.5000 0.3077 0.3500 0.2973 2.0000 0.2000 0.2000 0.2000 2.5000 0.1379 0.1500 0.1401 3.0000 0.1000 0.1000 0.1000 3.5000 0.0755 0.0794 0.0745 4.0000 0.0588 0.0588 0.0588 4.5000 0.0471 0.0486 0.0484 5.0000 0.0385 0.0385 0.0385 例 对 在[-5, 5]上, 用n=11个等距分点作分段线性插值和三次样条插值, 用m=21个插值点作图,比较结果. x y y1 y2 linear Spline 解 在命令窗口输入: 例 在一天24h内, 从零点开始每间隔2h测得的环境温度为 12, 9, 9, 10, 18, 24, 28, 27, 25, 20, 18, 15, 13 (单位: ) 推测在每1s时的温度. 并描绘温度曲线. t=0:2:24 T=[12 9 9 10 18 24 28 27 25 20 18 15 13] plot(t,T,*) ti=0:1/3600:24 T1i=interp1(t,T,ti) plot(t,T,*,ti,T1i,r-) T2i=interp1(t,T,ti,spline) plot(t,T,*,ti,T1i,r-,ti,T2i,g-) 例 在飞机的机翼加工时, 由于机翼尺寸很大, 通常在图纸上只能标出部分关键点的数据. 某型号飞机的机翼上缘轮廓线的部分数据如下: x 0 4.74 9.05 19 38 57 76 95 114 133 y 0 5.23 8.1 11.97 16.15 17.1 16.34 14.63 12.16 6.69 x 152 171 190 y 7.03 3.99 0 例 在飞机的机翼加工时, 由于机翼尺寸很大, 通常在图纸上只能标出部分关键点的数据. 某型号飞机的机翼上缘轮廓线的部分数据如下: x=[0 4.74 9.05 19 38 57 76 95 114 133 152 171 190] y=[0 5.23 8.1 11.97 16.15 17.1 16.34 14.63 12.16 9.69 7.03 3.99 0] xi=[0:0.001:190] yi=interp1(x,y,xi,spline) plot(xi,yi) 例 天文学家在1914年8月份的7次观测中, 测得地球与金星之间距离(单位: m), 并取其常用对数值与日期的一组历史数据如下所示, 试推断何时金星与地球的距离(单位: m)的对数值为 9.9352. 日期 18 20 22 24 26 28 30 距离对数 9.9618 9.9544 9.9468 9.9391 9.9312 9.9232 9.9150 解 由于对数值 9.9352 位于 24 和 26 两天所对应的对数值之间, 所以对上述数据用三次样条插值加细为步长为1的数据: 解 先绘图看变化规律 x=[18:2:30] y=[9.9618 9.9544 9.9468 9.9391 9.9312 9.9232 9.9150] plot(x,y,’*’); 对数值 9.9352 位于 24 和 26 两天所对应的对数值之间 解 由于对数值 9.9352 位于 24 和 26 两天所对应的对数值之间, 所以对上述数据用三次样条插值加细为步长为1的数据: x=[18:2:30] y=[9.9618 9.9544 9.9468 9.9391 9.9312 9.9232 9.9150] plot(x,y,’*’); xi=[18:1:30] yi=interp1(x,y,xi,spline) A=[xi;yi] A = 18.0000 19.0000 20.0000 21.0000 22.000