【Matlab】线性插值

例：

有一组函数，x=0:10:100，y=f(x)，想线性插值的方式求函数值，例如：

x=2.5，根据x=0、f(0)、x=10、f(10)插值求x=2.5时的值

x=-10，根据x=0、f(0)、x=10、f(10)插值求x=-10时的值

x=55，根据x=50、f(50)、x=60、f(60)插值求x=55时的值

x=101，根据x=90、f(90)、x=100、f(100)插值求x=101时的值

注：此情况下，内插精度较高

代码：

主代码：

%线性插值只有内插时精度是最好的
clc
clear
close all
%%
x=0:1:100;
y=sin(x);
x0=9.05;
y0=ChaZhi(x,y,x0)
%%
figure(1)
plot(x,y,'b-')
hold on
plot(x0,y0,'r*')

插值函数：

function y0=ChaZhi(x,y,x0)%x是行矩阵
%%
xx=sort([x,x0]);
label=find(xx==x0);%找到x0的标签
label=label(1);
%% 找最近的两个数，超范围了就取最小/大的两个值
if label==1
    x1=xx(2);y1=y(1);
    x2=xx(3);y2=y(2);
elseif label==length(xx)
    x1=xx(length(xx)-1);y1=y(length(y));
    x2=xx(length(xx)-2);y2=y(length(y)-1);
else
    x1=xx(label-1);y1=y(label-1);
    x2=xx(label+1);y2=y(label);
end
%% 线性插值
y0=y2-(y2-y1)/(x2-x1)*(x2-x0);