【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时的值

注:此情况下,内插精度较高

代码:

主代码:

matlab 复制代码
%线性插值只有内插时精度是最好的
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*')

插值函数:

matlab 复制代码
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);
相关推荐
studyer_domi9 小时前
matlab质子磁力仪传感器线圈参数绘图
人工智能·matlab
青橘MATLAB学习15 小时前
模糊综合评价法:原理、步骤与MATLAB实现
开发语言·算法·数学建模·matlab·分类
studyer_domi15 小时前
matlab 三维时频图绘制
开发语言·matlab
studyer_domi19 小时前
matlab飞行姿态pid控制
matlab
Matlab仿真实验室19 小时前
基于Matlab实现信道估计仿真(源码)
开发语言·matlab·信道估计仿真
studyer_domi20 小时前
matlab 汽车abs的pid控制仿真
开发语言·matlab·汽车
studyer_domi20 小时前
matlab 汽车abs的模糊pid和pid控制仿真
开发语言·matlab·汽车
freexyn1 天前
Matlab自学笔记四十七:如何把日期时间型数据作为横坐标进行绘图
开发语言·笔记·matlab
见你背影1 天前
matlab数据处理:创建网络数据
matlab
chinakq2 天前
C++学习计划(三周) 有matlab和python基础,如何快速入门C++(利用《C++ primer plus》)
c++·学习·matlab