实验要求

通过n+1个节点的次数不超过n的Lagrange插值多项式为:

其中,Lagrange插值基函数 ,

k=0,1,…,n。

x 0 10 20 30 40 50 60 70 80 90 100 110 120
y 5 1 7.5 3 4.5 8.8 15.5 6.5 -5 -10 -2 4.5 7

利用Lagrange插值多项式 求被插值函数f(x)在点x=65处的近似值。建议:画出Lagrange插值多项式 的曲线。

代码实现

X=[0,10,20,30,40,50,60,70,80,90,100,110,120];
Y=[5,1,7.5,3,4.5,8.8,15.5,6.5,-5,-10,-2,4.5,7];
n=length(X);%得到点的个数
disp(n);%n是13
syms t;
f=0;
for k=1:nl=Y(k);for j=1:nif j==kcontinue;elsel=l*((t-X(j))/(X(k)-X(j)));endendf=f+l;
end
disp('f是');
u=simplify(f);%把f化简
disp(u);
s=eval(subs(f,'t',65));
disp('x的值为65时,表达式的值为');
disp(s);
t=0:0.001:120;%产生从0到120长度为0.001的数组
d=eval(f);%eval函数可以把一个字符串变成可执行的表达式
figure
plot(t,d,'r');
grid on;%显示网格
xlabel('x');
ylabel('function');
for k=1:ntext(X(k),Y(k),'*','color','g');%图像上标记某一个点
end
text(65,s,'*','color','k');

运行图像

Matlab实现Lagrange插值多项式相关推荐

  1. 解决龙格现象matlab,matlab实现Lagrange多项式插值观察龙格现象

    Matlab进行Lagrange多项式插值 拉格朗日插值法对函数y=1./(1+25*x.^2)在区间[-1,1]进行5次.10次.15次插值观察龙格现象 主程序 1.拉格朗日 function [c ...

  2. Lagrange插值多项式算法

    标题:Lagrange插值多项式算法 一.公式介绍 Lagrange多项式公式: L[n](u) = 累加符号 l[i](u) * y[i] {i = 0,1,...,n}Lagrange因子公式 : ...

  3. matlab中如何用lagrange插值,MATLAB实现Lagrange插值函数

    MATLAB实现Lagrange插值函数 MATLAB实现Lagrange插值函数 MATLAB实现Lagrange插值函数 首先我们绘制Lagrange基函数 首先给出一个Lagrange基函数比较 ...

  4. 拉格朗日Lagrange插值多项式

    分享一下我老师大神的人工智能教程!零基础,通俗易懂!http://blog.csdn.net/jiangjunshow 也欢迎大家转载本篇文章.分享知识,造福人民,实现我们中华民族伟大复兴! 本文转载 ...

  5. matlab:Lagrange插值函数构造

    L a g r a n g e Lagrange Lagrange基函数及 L a g r a n g e Lagrange Lagrange多项式的构造 构 造 L a g r a n g e 基 ...

  6. matlab用牛顿差值计算三次差值多项式,计算方法用Newton插值多项式求函数的近似值.docx...

    计算方法用Newton插值多项式求函数的近似值 计算方法课程设计题 目: 用Newton插值多项式 处理磁化曲线学 院: 理学院班 级: 学 生 姓 名: 学 生 学 号: 指 导 教 师: 2017 ...

  7. 求拉格朗日多项式matlab,拉格朗日插值多项式积分求圆周率近似Matlab实现

    Lagrange 插值多项式积分求圆周率近似 摘要: 公式1:y1=4/(1+x^2) 公式2:y2=4*sqrt(1-x^2) 分别对公式1.公式2求其拉格朗日插值多项式,再对其求0-1上的定积分来 ...

  8. 数值分析上机题matlab线性方程组,数值分析上机实验报告 - 线性方程组部分实验题1...

    s=A(i,(i+1):n)*x((i+1):n,1); else s=0; end x(i,1)=(b(i)-s)/A(i,i);end %Cholosky分解方法***************** ...

  9. 实验Matlab数值运算,MATLAB数值实验一(数据的插值运算及其应用完整版

    <MATLAB数值实验一(数据的插值运算及其应用完整版>由会员分享,可在线阅读,更多相关<MATLAB数值实验一(数据的插值运算及其应用完整版(6页珍藏版)>请在人人文库网上搜 ...

最新文章

  1. Linux服务器oracle数据库重启服务、重启监听方法,oracle数据库sysdba管理员登录方法
  2. day17——其他内置函数
  3. CONVERT_DATE_INPUT
  4. 电离辐射防护与辐射源安全基本标准_辐射防护与安全机考难点考点解析辐射防护标准...
  5. Python:eval的妙用和滥用
  6. 查看linux系统核数
  7. 2020美国纽约大学计算机科学排名,2020美国纽约大学排名第几
  8. 程序员拒绝一个合理需求的15个方法!
  9. php写时间 自动,ThinkPHP6.0 自动写入时间戳
  10. 发送的消息无法订阅_微信服务号和订阅号的如何选择
  11. 网络干货,无论是运维还是开发都要知道的网络知识系列之(八)
  12. 堆排序(C++实现,非常简洁明了)
  13. H.266/VVC代码学习笔记15:VTM6.0中的xCheckRDCostMergeTriangle2Nx2N()函数
  14. java int 32位_JAVA基本数据类型所占字节数是多少?(32位系统)
  15. 解决“gettools.exe 已停止工作”问题
  16. .NetCore对接各大财务软件凭证API——用友系列(2)
  17. 拓展——RGB颜色查询对照表
  18. 关于Linux备份文件和应用的几个命令:tar和cp
  19. java公路车的气嘴,5分钟了解自行车常见的三种气嘴
  20. IAR下载程序只有提示音,没有任何反应。或报错The configuration does not have debuggable output.(A debug-only project shoul

热门文章

  1. ccsds ldpc matlab,CCSDS标准的LDPC编译码仿真
  2. php做手机网站,帝国CMS如何制作手机网站
  3. 银行信用评分卡建模原理
  4. 多边形套索及魔棒工具
  5. ai人工智能在手机的应用_常识在人工智能中仍然不常见
  6. 海康威视:单季盈利下滑与华为并无关系,2019 年将进入 AI 落地阶段...
  7. 【解决方案】助力电子商务平台建设TSINGSEE青犀视频实现AI全流程监管
  8. Linux嵌入式开发——vim编辑器
  9. 战国策评析(三)---齐
  10. 猫猫学习ios 之第三方登录友盟实现