创建于2019年6月3日,北京
现有一组数据点:
x=0.25 0.4 0.6 0.8 1
y=3.36 3.30 3.33 2.75 2.43
若直接采用matlab绘制曲线得到不是很完美的一组曲线:
为了更好的绘制该曲线,采用插值的方式来绘制。
matlab中常用的插值方式有以下几种:
'nearest'最近项插值; 'linear'线性插值 ;'spline' 立方样条插值;'cubic' 立方插值
在脚本文件中,插值的语法为
y=interp1(x0,y0,x,'method')
其中,x0,y0是已知的数据向量;x应以升序或者降序排列(所有的插值方法要求x0是单调的);method指定插值的方法,默认为线性插值,即linear,另可取另外三种插值方式
通过下面代码绘制四种方式结果对比:
clear all
x0=[0.25 0.4 0.6 0.8 1];
y0=[3.36 3.30 3.33 2.75 2.43];
x=linspace(x0(1),x0(end),100);
y1=interp1(x0,y0,x,'linear');
y2=interp1(x0,y0,x,'spline');
y3=interp1(x0,y0,x,'nearst');
y4=interp1(x0,y0,x,'cubic');
plot(x0,y0,'o',x,y1,'r-',x,y2,'b-',x,y3,'k-',x,y4,'g-','LineWidth',1.2,'markersize',4)
set(gca,'FontName','Times New Roman','FontSize',9);
legend('fontname{宋体}数据点','linear','spline','nearst','scubic','location','West')
legend boxoff;
结果可以看出最近项插值效果很差,样条插值在最后的点处出现一些误差,立方插值效果较好
个人观点:平时用spline样条插值应该是效果比较好的,本次算例出现问题应该是数据点特殊个边界条件的不合理
对一组数进行指数的拟合,采用ployfit和polyval函数,如下所示
clear all
a=xlsread('Case258(1).xlsx');%读取数据
x1=a(:,1);y1=a(:,2);y2=a(:,3);y3=a(:,4)%读取数据
p1=polyfit(x1(1:24),log(y1(1:24)),5);%生成多项式的系数,注意y值取了对数
y11=exp(polyval(p1,x1));%根据x1的值,对y重新赋值
plot(x1,y11,'-o')%绘图
set(gca,'XLim',[0.05 1.8]);%X轴的数据显示范围