大家好,又见面了,我是你们的朋友全栈君。
1、一维插值:
已知离散点上的数据集 ,即已知在点集X= 上的函数值Y=
,构造一个解析函数(其图形为一曲线)通过这些点,并能够求出这些点之间的值,这一过程称为一维插值。
MATLAB命令:yi=interp1(X, Y, xi, method)
该命令用指定的算法找出一个一元函数 ,然后以 给出
处的值。xi可以是一个标量,也可以是一个向量,是向量时,必须单调,method可以下列方法之一:
‘nearest’:最近邻点插值,直接完成计算;
‘spline’:三次样条函数插值;
‘linear’:线性插值(缺省方式),直接完成计算;
‘cubic’:三次函数插值;
对于[min{xi},max{xi}]外的值,MATLAB使用外推的方法计算数值。
例1:已知某产品从1900年到2010年每隔10年的产量为:75.995, 91.972, 105.711, 123.203,
131.699, 150.697, 179.323, 203.212, 226.505, 249.633, 256.344,
267.893,计算出1995年的产量,用三次样条插值的方法,画出每隔一年的插值曲线图形,同时将原始的数据画在同一图上。
解:程序如下
year=1900:10:2010;
product=[75.995, 91.972, 105.711, 123.203, 131.699, 150.697,
179.323, 203.212, 226.505, 249.633, 256.344, 267.893]
p1995=interp1(year,product,1995)
x=1900:2010;
y=interp1(year,product,x,’cubic’);
plot(year,product,’o’,x,y);
计算结果为:p1995=252.9885。
2、二维插值
已知离散点上的数据集 ,即已知在点集 上的函数值
,构造一个解析函数(其图形为一曲面)通过这些点,并能够求出这些已知点以外的点的函数值,这一过程称为二维插值。
MATLAB函数:Zi=interp2(X,Y,Z,Xi,Yi,method)
该命令用指定的算法找出一个二元函数 ,然后以 给出 处的值。返回数据矩阵 ,Xi,Yi是向量,且必须单调,
和meshgrid(Xi,Yi)是同类型的。method可以下列方法之一:
‘nearest’:最近邻点插值,直接完成计算;
‘spline’:三次样条函数插值;
‘linear’:线性插值(缺省方式),直接完成计算;
‘cubic’:三次函数插值;
例2:已知1950年到1990年间每隔10年,服务年限从10年到30年每隔10年的劳动报酬表如下:
表:某企业工作人员的月平均工资(元)
年份 1950 1960 1970 1980 1990
服务年限
10 150.697 179.323 203.212 226.505 249.633
20 169.592 195.072 239.092 273.706 370.281
30 187.652 250.287 322.767 426.730 598.243
发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/139513.html原文链接:https://javaforall.cn