7.1.1 一维内插
线性内插是假设在二个已知数据中的变化为线性关系,因此可由已知二点的座标(a, b)去计算通过这二点的 斜线,公式如下:
其中 a<b<c 在上式的 b 点即是代表要内插的点,f(b) 则是要计算的内插函数值。下图即是一个以二种内插 法的比较
\pcxfile[12cm,5cm]{fig9_1.pcx}
\caption{线性式与 spline 函数的曲线契合}
线性内插是最简单的内插方法,但其适用范围很小;如果原来数据的函数f有极大的变化,假设其数据点之
间为线性变化并不合理。所以我们可以用二次、三次方程式或是另一种称为spline函数来近似原来数据的函
数。MATLAB的一维内插函数是interp1,其语法为interp1(x,y,xi),interp1(x,y,xi,'method');其中的x,y是原已知的
数据的x,y值,而xi则是要内插的数据点,另外method可以设定内插方法有
linear,cubic,spline,分别是一次、三
次方程式和spline函数,其中预设方法是linear。如果数据的变化较大,以
spline函数内插所形成的曲线最平滑
,所以效果最好。而三次方程式所得到的内插曲线平滑度,则介于线性与spline函数之间。
我们以下面的例子说明。假设有一个汽车引擎在定转速下,温度与时间(单位为sec)的三次量测值如下
其中温度的数据从 20oC变化到 503oC,如果要估计在t=2.6, 4.9 sec 的温度,可以下列指令计算
>> x=[0 1 2 3 4 5]'; % 键入时间
>> y=[0 20 60 68 77 110]'; % 键入第一组时间
>> y1=interp1(x,y,2.6) % 要内插的数据点为 2.6
y1 = % 对应 2.6 的函数值为 64.8
64.8
>> y1=interp1(x,y,[2.6 4.9]) % 内插数据点为 2.6, 4.9,注意用[ ]将多个内插点放在其中
y1 =
64.8
106.7
>> y1=interp1(x,y,2.6,'cubic') % 以三次方程式对数据点 2.6 作内插
y1 = % 对应 2.6 的函数值为 66.264
66.264
>> y1=interp1(x,y,2.6,'spline') % 以spline函数对数据点 2.6 作内插
y1 = % 对应 2.6 的函数值为 66.368
66.368
以下的例子还配合绘图功能,用以比较不同内插方法的差异。
>> h=1:12;
>> temp=[5 8 9 15 25 29 31 30 22 25 27 24]; % 这组温度数据变化较大
>> plot(h,temp,'--',h,temp,'+') % 将线性内插结果绘图
>> h_3=1:0.1:12 % 要每0.1小时估计一次温度值
>> t_3=interp1(h,temp,h_3,'cubic') % 以三次方程式做内插
>> t_s=interp1(h,temp,h_3,'spline') % 以spline函数做内插
>> hold on
>> subplot(1,2,1)
>> plot(h,temp,'--',h,temp,'+',h_3,t_3) % 将线性及三次方程式内插绘图
>> subplot(1,2,2)
>> plot(h,temp,'--',h,temp,'+',h_3,t_s) % 将线性方程式及spline内插绘图
>> hold off
上一页 下一页 讲义大纲