- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Matlab线性插值
Matlab线性插值??
已知离散点上的数据集,即已知在点集X上对应的函数值Y,构造一个解析函数(其图形为一曲线)通过这些点,并能够求出这些点之间的值,这一过程称为一维插值。
MATLAB命令:yi=interp1(X, Y, xi, method)
该命令用指定的算法找出一个一元函数,然后以给出xi处的值。xi可以是一个标量,也可以是一个向量,是向量时,必须单调,method可以下列方法之一:
nearest:最近邻点插值,直接完成计算;
spline:三次样条函数插值;
linear:线性插值(缺省方式),直接完成计算;
cubic:三次函数插值;
对于[min{xi},max{xi}]外的值,MATLAB使用外推的方法计算数值。下面是例子:
时间 6:00 8:00 10:00 12:00 14:00 16:00 18:00 室内温度()各点的近似温度,并对四种插值方法下的结果进行比较。
t=6:2:18;
T=[18.0 20.0 22.0 25.0 30.0 28.0 24.0];
tt=6.5:2:17.5;%使用不同的方法行一插
yi_spline=interp1(t,T,tt,spline) %三次条插
yi_cubic=interp1(t,T,tt,cubic) %三次多式插yi_linear=interp1(t,T,tt) %线性插yi_nearest=interp1(t,T,tt,nearest) %最近邻点插%绘制像比
plot(t,,ko);
grid on;hold on;
plot(tt,yi_spline,b,tt,yi_cubic,g,tt,yi_linear,r-,tt,yi_nearest,k);
legend(原始数据,三次条插三次xlabel(时间/t);ylabel(室内温度(T)/摄氏度);
title(四种插值方法的比较 )
结果:
yi_spline =
18.5020 20.4986 22.5193 26.3775 30.2051 26.8178
yi_cubic =
18.5000 20.4812 22.6305 26.3086 29.8125 27.2344
yi_linear =
18.5000 20.5000 22.7500 26.2500 29.5000 27.0000
yi_nearest =
18 20 22 25 30 28
图片:
下面是例子:
t=1900:10:1990;
p=[75.995,91.972,105.711,123.203,131.669,150.697,179.323,203.212,226.505,249.633];
x=1900:0.01:1990;%使用不同的方法进行一维插值
yi_linear=interp1(t,p,x); %线性插值
yi_spline=interp1(t,p,x,spline);%三次样条插值
yi_cubic=interp1(t,p,x,cubic);%三次多项式插值
yi_v5cubic=interp1(t,p,x,v5cubic);%matlab5中使用的三次多项式插值
%绘制图像对比
%?subplot是将多个图画到一个平面上的工具。其中,m表示是图排成m行,n表示图排成n列,也就是整个figure中有n个图是排成一行的,一共m行,如果第一个数字是2就是表示2行图。p是指你现在要把曲线画到figure中哪个图上,最后一个如果是1表示是从左到右第一个位置。
subplot(2,1,1);
plot(t,p,ko);
hold on;
plot(x,yi_linear,g,LineWidth,1.5);grid on;
plot(x,yi_spline,y,LineWidth,1.5);
title(Linear VS Spline )
?
subplot(2,1,2);
plot(t,p,ko);
hold on
plot(x,yi_cubic,g,LineWidth,1.5);grid on;
plot(x,yi_v5cubic,y,LineWidth,1);
title(Cubic VS V5cubic );
%创建新图形窗口
figure
yi_nearest=interp1(t,p,x,nearest);%最邻近插值法
plot(t,p,ko);
hold on
plot(x,yi_nearest,g,LineWidth,1.5);grid on;
title(Nearest Method);
%以下是根据拟合估计
文档评论(0)