MATLAB 数据插值
零件加工问题:
x=[0,3,5,7,9,11,12,13,14,15];
y=[0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.0,1.6];
x1=0:0.1:15;
y1=interp1(x,y,x1,'spline');
plot(x1,y1)
数据插值的实现方法:
为什么这两种插值方法都用3次多项式而不用更高次的?
多项式次数并非越高越好。次数越高,越容易产生震荡而偏离原函数,这种现象称为龙格现象。
interp 2():二维插值函数
调用格式:Z1=interp2((X,Y,Z,X1,Y1,method)
粮储仓的通风控制问题:
分析:粮食水分吸收和蒸发量相等,这个湿度称为平衡点湿度。平衡点湿度根粮温,粮湿有关。
x=20:10:90;%粮湿
y=(0:5:20)';%粮温
z=[8.9,10.32,11.3,12.5,13.9,15.3,17.8,21.3;
8.7,10.8,11,12.1,13.2,14.8,16.55,20.8;
8.3,9.65,10.88,12,13.2,14.6,16.4,20.5;
8.1,9.4,10.7,11.9,13.1,14.5,16.2,20.3;
8.1,9.2,10.8,12,13.2,14.8,16.9,20.9];
xi=20:90;
yi=(0:20)';
zi=interp2(x,y,z,xi,yi);
surf(xi,yi,zi)