line_int_wise.m

来自「为配套王能超版本数值分析简明教程编写,方便学生教师演示!」· M 代码 · 共 24 行

function yi=line_int_wise(x,y,xi)
%分段线性插值公式
%x为向量，全部插值节点
%y为向量，插值节点处函数值
%xi为标量，被估计函数自变量
n=length(x);m=length(y);
if n~=m
error('The length of X and Y must be equal!');
return;
end

for k=1:n-1
    
            %输入的插值节点必须互异
            if abs(x(k)-x(k+1))<eps
                error('the DATA is error')
                return;
            end
            if x(k)<=xi&xi<=x(k+1)
                yi=(xi-x(k+1))/(x(k)-x(k+1))*y(k)...
                    +(xi-x(k))/(x(k+1)-x(k))*y(k+1);
                return;
            end
        end