MATLAB实现Lagrange插值函数

MATLAB实现Lagrange插值函数

MATLAB实现Lagrange插值函数

首先我们绘制Lagrange基函数

首先给出一个Lagrange基函数比较复杂的设法:

function y=lagrange(x0,y0,x)

n=length(x0);m=length(x);

for i=1:m

z=x(i);

s=0.0;

for k=1:n

p=1.0;

for j=1:n

if j~=k

p = p*(z-x0(j))/(x0(k)-x0(j));

end

end

s=p*y0(k)+s;

end

y(i)=s;

end

罗里吧嗦我也看不太懂自己在干嘛了/擦汗

然后是简单版,但是这里使用了几个可能大家不太熟悉的函数(我自己不太熟悉)

其中有

conv函数:多项式乘积(向量形式)

polyder函数:求解多项式微分表达式(求导)

polyval函数:多项式求值

deconv函数:多项式作除

w=1;l=zeros(9);t=1:0.01:9;

for i=1:9

x(i)=i;

w=conv(w,[1,-x(i)]);

end

dw=polyder(w);

for i=1:9

dwx(i)=polyval(dw,x(i));

l(i,:)=deconv(w,conv([1,-x(i)],dwx(i)));

end

for i=1:9

yi=polyval(l(i,:),t');

plot(t,yi,'Color',[rand(),rand(),rand()]);

hold on;

end

%基函数的表达式为[2,0,1,7,1,1,0,0,0]*[l1(x),l2(x),l3(x),l4(x),l5(x),l6(x),l7(x),l8(x),l9(x)]'

%我们可以通过直接调用基函数进行绘制langange函数,也可以直接编写langange函数,在此为充分熟悉langange函数,采用第二种实现方式。

x0 = [1,2,3,4,5,6,7,8,9];

y0 = [2,0,1,7,1,1,0,0,0];

x = linspace(1,9,1000);

%我们调用lagrange函数

y1 = lagrange(x0,y0,x);

plot(x,y1,'g')

hold on

plot(x0,y0,'*')

图像如图所示:

在MATLAB中有现成的可以调用求lagrange函数的函数:

p=polyfit(x,y,n)

根据拟合数据x,y求n阶多项式p(x)的系数,p是系数

我们通过以上实际上是加深自己对lagrange函数的理解,实际使用时本人还是觉得直接用现成的最方便。

MATLAB实现Lagrange插值函数相关教程

matlab中如何用lagrange插值,MATLAB实现Lagrange插值函数相关推荐

  1. matlab样条插值如何用,三次样条插值matlab实现

    %三次样条差值-matlab通用程序 - zhangxiaolu2015的专栏 - CSDN博客 https://blog.csdn.net/zha %[图文]三次样条插值算法详解_百度文库 http ...

  2. matlab中如何用数组绘图,matlab 基础学习(1)——数组,矩阵操作和绘图基础学习...

    一维数组(创建,赋值,取值,分割取值,四则运算等操作) 创建一维数组 >> A=[] %创建空的数组 >> B=[1 2 3 4 5] %两种创建行向量,数组的方式 >& ...

  3. 将z的第二列除以根号三MATLAB,matlab中如何用牛顿法求根号2的近似值

    如何用matlab求出图中各条直线的斜率 1.瞬时斜率:求导.方法:d(k)=(u(k)-u(k-1))/T(k);d(0)=0;其中d就是瞬时斜率2.平均斜率:求拟合.方法:p=polyfit(x, ...

  4. matlab中的timer模块,[转载]Matlab中Timer的使用

    Matlab中Timer的使用 鉴于Matlab中缺乏多线程机制,使用Timer无疑是一个很重要的工具,Matlab中Timer是一个Java对象. (1)Timer 的定义 t=timer(); 设 ...

  5. m 文件 dll matlab 中调用_如何在matlab中调用python程序

    现在python很火,很多代码都是python写的,如果你和我一样,习惯了使用matlab,还想在matlab中调用Python的代码,应该怎么办呢?其中一条思路:首先在matlab中调用系统脚本命令 ...

  6. matlab中m文件是什么,MATLAB中M文件的使用

    运行方法: 使用matlab的时候,可以在"Command Window"内直接书写matlab代码,也可以将代码保存到M文件中,然后运行该文件.使用matlab主界面菜单&quo ...

  7. matlab 如何axis,在matlab中axis是什么意思,matlab中axis的用法

    在matlab中axis是什么意思轴在matlab中是什么意思,matlab中的Axis是指在绘图中通常使用axis函数来设置坐标值的范围,在matlab的命令窗口中输入doc轴或help轴即可获得该 ...

  8. 在MATLAB中使用数学符号,在matlab中怎么输入特殊符号 function在MATLAB中怎么用

    导航:网站首页 > 在matlab中怎么输入特殊符号 function在MATLAB中怎么用 在matlab中怎么输入特殊符号 function在MATLAB中怎么用 相关问题: 匿名网友: 一 ...

  9. matlab中a2=poly(p2),插值与拟合matlab实现

    插值与拟合的Matlab实现 王正盛编写 在科技工程中,除了要进行一定的理论分析外,通过实验.观测数据,做分析.处理也是必不可少的一种途径.由于实验测定实际系统的数据具有一定的代表性,因此在处理时必须 ...

最新文章

  1. 电源稳定性测试软件,电源稳定性测试
  2. Python面试必知(四)
  3. JAXB做错了; 尝试Xembly
  4. python 运算及注释
  5. 在Kubernetes上部署Hyperledger Fabric
  6. 【windows10】CMD命令行隐藏窗口运行
  7. Python——输入若干行,每行若干行整数,求所有整数的最大值
  8. 如何进行数据安全管理体系建设?
  9. 华为rh2288v2服务器系统,扩展能力强 华为RH2288H V2服务器解析
  10. GIt上传到远程库Gitee出现的问题,需要身份表明和私人邮箱的问题
  11. java Date与 double 互转
  12. 使电动机反转的matlab仿真图,关于电机正反转电路、电机顺序启动、逆序停止的电路仿真图分析...
  13. 从一个资深面试官角度看Java面试:面试了20家大厂与面试了上百人之后,发现这样介绍项目经验,显得项目很牛逼!
  14. 波罗的海国家同意支持分布式账本技术开发
  15. Python基础五十四题19~36题
  16. 2021 谷歌浏览器等无法下载android studio
  17. linux rmvb转mp4,linux系统下实现播放rmvb文件
  18. 详说MOS管开关时的米勒效应
  19. JointJS动态流程图升级
  20. PCL 非线性最小二乘法拟合圆柱

热门文章

  1. Katalon自动化测试基础教程
  2. 【转】WIN32编程中的HRESULT类型和LRESULT类型的区别
  3. Django的初步学习及线上部署
  4. Element ‘web-app’ cannot have character [children], because the type’s content type is element-onl
  5. 【干货】ArcGIS自定义符号库
  6. Vue初始化创建项目
  7. 自然语言处理入门学习系列一
  8. 面向对象的主要优点有哪些?
  9. win10 uwp 入门
  10. 数据结构——顺序表的逆置