Matlab中顺序查找、二分搜索及一二三次样条的一些总结

无论是哪种方法查表，本质上都是对数值的大小进行比较。首先看顺序查找

function []=TableLook(x,y,xx)
%对数组中的数据进行顺序查表，进行指定量与向量中的元素进行依次的比较
n=length(x);
if xx<x(1)||xx>x(n)
    error('这个数在给定数组范围外')
end
i=1;
while(1)
    if xx<=x(i+1)
        break
    else
        i=i+1;
    end
end
yi=y(i)+((y(i+1)-y(i))/(x(i+1)-x(i)))*(xx-x(i)) %一次样条
end

在这里主要是将输入数值xx与输入数组x、y中的数字进行逐一比较，数值较多的时候比较繁琐。

二分查找：

function []=TableLookBin(x,y,xx)
%二分搜索实现查表操作
n=length(x);
if xx<x(1)||xx>x(n)
    error('估值数据不在给定范围中')
end
iL=1;iU=n;
while(1)
if iU-iL<=1
    break
end
iM=fix((iL+iU)/2);
if x(iM)<xx
    iL=iM;
else
    iU=iM;
end
end
yi=y(iL)+((y(iL+1)-y(iL))/(x(iL+1)-x(iL)))*(xx-x(iL))  %一次样条
end

对于二三次样条来说，基本原理一致，式子给的是三次样条，但推到过程写的是二次样条，值得注意的是，样条的基本条件是建立在你已经找到了那个数据点在哪个区间的基础上，也就是说不管你用什么查表方式，都要先搜索到数据。

自然三次样条基本式子

 其中hi=X（i+1）-Xi，f【xi，xj】是一个有限差分，即f【xi，xj】=（fi-fj）/（xi-xj），推导过程就是四个条件。

我们以二次样条的推导过程为例子。

（1）函数必须经过所有结点（连续性条件），在这个条件下fi=ai+bi（x-xi）+ci（x-xi）^2，可以代入下限xi的值将ai算出来，即fi=ai。也就是说，高次样条下的常数项必须等于每个区间的左端点数值。

此时需要的条件数减少到2（n-1）

（2）可以从样条图形中看出，相邻多项式