一维插值[实验].pptVIP

下一页 上一页 主 页 Gongqu * 用MATLAB作分段线性插值计算 分段线性插值命令: y=interp1(x0,y0,x) xi处的插值结果 插值节点 被插值点 例：在1-12的11小时内，每隔1小时测量一次温度，测得的温度依次为：5，8，9，15，25，29，31，30，22，25，27，24。试估计每隔1/10小时的温度值。 用MATLAB作分段线性插值计算 hours=1:12; temps=[5 8 9 15 25 29 31 30 22 25 27 24]; h=1:.1:12; t=interp1(hours,temps,h) plot(hours,temps,+,h,t) title(线性插值下的温度曲线), xlabel(Hour),ylabel(Degrees Celsius) 用MATLAB作分段线性插值计算 Method可取： ‘nearest’ ：最邻近插值；‘spline’ ：三次样条插值； ‘cubic’ ： 立方插值;缺省时： 分段线性插值。 注意：所有的插值方法都要求x是单调的，并且xi不能够超过x的范围。 用MATLAB作插值计算 一维插值函数： yi=interp1(x，y，xi，method) 插值方法 被插值点 插值节点 xi处的插值结果 例:在 上取11个点，作三次样条插值, 观察三次样条插值曲线与g(x)的误差. 用MATLAB作三次样条插值计算 x0=linspace(-5,5,11); y0=1./(1+x0.^2); x=linspace(-5,5,100); y=interp1(x0,y0,x,spline); x1=linspace(-5,5,100); y1=1./(1+x1.^2); plot(x1,y1,k,x0,y0,+,x,y,r); M文件yangcha.m 用MATLAB作三次样条插值计算 为例，作插值的比较 例：以 插值方法的比较 n=11;m=101; x=-5:10/(m-1):5; y=1./(1+x.^2); plot(x,y,k:), pause,hold off, y2=interp1(x0,y0,x);hold on, plot(x,y2,b), gtext(Pies.Lin.), pause,hold off, y3=interp1(x0,y0,x,spline); hold on,plot(x,y3,m), gtext(Spline) 范例1：地图绘制问题 某国的地图如上页图,其国土面积为156.6500万平方公里。对某一地区或国家，如何根据测绘部门测量的数据绘制这样一张地区或国家的地图。 测量出该国边界上的一系列点的平面坐标如下页表。 范例1：地图绘制问题 表1 下边疆采样坐标 注：表中数据为相对数字，比例尺：30 (数据单位)：100公里(实际单位)。 范例1：地图绘制问题 表2 上边疆采样坐标 范例1：地图绘制问题 简单地图再现 相邻两点依次用直线段连接，构成一个封闭的图形：地图 范例1：地图绘制问题 x0=x1;y0=y1; %下边界处理，插值点 n=length(x0) x=linspace(x0(1),x0(n),1000); y=interp1(x0,y0,x,spline); plot(x,y); hold; clear y x0 y0 x0=x2; y0=y2; % 上边界处理，插值点 n=length(x0); x=linspace(x0(1),x0(n),1000); y=interp1(x0, y0, x, spline); plot(x, y); MATLAB程序