- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Matlab 一维插值interp1 例子 及 可视拟合界面
一维插值:
已知离散点上的数据集?,即已知在点集X上对应的函数值Y,构造一个解析函数(其图形为一曲线)通过这些点,并能够求出这些点之间的值,这一过程称为一维插值。
MATLAB命令:yi=interp1(X, Y, xi, method)
一些个人经验说明:
①关于拟合参数的,X必须是向量,行向量或列向量均可,不可以是复数
②Y是向量或矩阵.但必须满足行数与length(X)相同即size(Y,1)==length (X)
③针对以上说明的例子
??? function tux=[5 1 2 20 14 21]y=rand(6,2)%按列计算的
xi=linspace(0,21,100);yi=interp1(x,y,xi,cubic)plot(x,y,o,xi,yi)size(x)size(y,1)length(x)
结果
size(x)?????? 6???? 1???? size(y,1)???????? 6??? length(x)????? 6
该命令用指定的算法找出一个一元函数?,然后以?给出?处的值。xi可以是一个标量,也可以是一个向量,是向量时,必须单调,method可以下列方法之一:
‘nearest’:最近邻点插值,直接完成计算;?????????
‘spline’:三次样条函数插值;
‘linear’:线性插值(缺省方式),直接完成计算;???
‘cubic’:三次函数插值;
对于[min{xi},max{xi}]外的值,MATLAB使用外推的方法计算数值。
%-- 09-4-1 下午8:38 --%%已知数据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(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 onplot(x,yi_cubic,g,LineWidth,1.5);grid on;plot(x,yi_v5cubic,y,LineWidth,1);title(Cubic VS V5cubic );%创建新图形窗口figureyi_nearest=interp1(t,p,x,nearest);%最邻近插值法plot(t,p,ko);hold onplot(x,yi_nearest,g,LineWidth,1.5);grid on;title(Nearest Method);
%以下是根据拟合估计
msg=?????????? year????? Cubic????? Linear??? Nearest??? Spline;for i=0:8??? n=10*i;??? year=1905+n;??? pop(i+1,1)=year;??? pop(i+1,2)=yi_cubic((year-1900)/0.01+1);???? pop(i+1,3)=yi_linear((year-1900)/0.01+1);????? pop(i+1,4)=yi_nearest((year-1900)/0.01+1);?????? pop(i+1,5)=yi_spline((year-1900)/0.01+1);endP=round(pop);
disp(msg)disp(P)
结果如图:
由此可见,各种插值的优劣,在速度上,Nearest最快,然后是Linear再到Cubic,最慢的是Spline.但是精度和曲线的平滑度恰好相反,Nearest甚至不连续~~
系统默认的是Linear
当然也可以用图形界面 更易理解
附录: Matlab 样条工具箱(Spline ToolBox)【信息来源教师博客】? Matlab样条工具箱中的函数提供了样条的建立,操作,绘制等功能;
一. 样条函数的建
您可能关注的文档
- 飞行航线上的五维可视化气象信息.pdf
- 食品感官资料.doc
- 高一上学期期末考试(英语)10.doc
- 高三第一学期期末试题江苏省苏北四市(徐、淮、连、宿)高三3月联考试题(物理)WORD版.doc
- 香港著名旅游景点及地图.pdf
- 高中历史地图.docx
- 高中历史相关知识点 or地图.docx
- 高功率准连续波腔内和频全固态黄光激光器.pdf
- 高性能仿真网络电子地图网站的设计与实现.pdf
- 高斯-克吕格与兰勃特投影,地理坐标与直角坐标网.pdf
- 化学反应速率与平衡.ppt
- 企业会计信息质量问题及对策【精品发布】.doc
- 免疫学的应用高二上学期生物人教版选择性必修一.pptx
- 高中物理自由落体运动课题研究力学部分+课件+-2023-2024学年高一上学期物理人教版必修第一册.pptx
- 北师大版生物七年级上册 《细胞的基本结构和功能》第2课时CEP课件1.pptx
- Unit 4 Grammar and usage(教学课件)-高中英语牛津译林版(2020)选择性必修第三册.pptx
- 特异性免疫高二上学期生物人教版选择性必修1.pptx
- 期末复习课件专题9+物体的浮沉条件及应用(2)++2023-2024学年人教版八年级物理下册.pptx
- 其他植物激素高二上学期生物人教版选择性必修1.pptx
- 三大费用核算节税技巧【会计实务操作教程】.pptx
文档评论(0)