function y=lagrange(x0,y0,x)
%拉格朗日插值函数
%n 个节点数据以数组 x0, y0 输入(注意 Matlat 的数组下标从1开始),
%m 个插值点以数组 x 输入,输出数组 y 为 m 个插值
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
先将上面的复制到Matlab函数命令下,然后运行下方实例
x0=1:1:20;
y0=x0.^2-20*x0-5;
x=1:0.1:20;
z=lagrange(x0,y0,x);
plot(x,z,':',x0,y0,'ko');
拉格朗日插值算法.zip
需积分: 0 11 浏览量
2022-12-24
10:24:14
上传
评论
收藏 739B ZIP 举报
![avatar](https://profile-avatar.csdnimg.cn/default.jpg!1)
Echo_w01
- 粉丝: 3
- 资源: 10
最新资源
- python爬虫代码详解
- 力道图数据,是一篇论文里的
- 无损压缩音频格式ERAC
- google-chrome-stable_current_x86_64 (1)(1).rpm
- Java源代码案例 - 使用正则表达式解析sql语法树.zip
- Dijkstra算法和Floyd算法 C++源代码案例.zip
- visualstudio安装教程
- 基于Selenium的Java爬虫实战(内含谷歌浏览器Chrom和Chromedriver版本123.0.6272.2)
- 小熊派BearPi-Pico H2821 SLE通信OLED显示.zip
- Anaconda3-2023.09-0-Linux-aarch64.sh
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback-tip](https://img-home.csdnimg.cn/images/20220527035111.png)