大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全家桶1年46,售后保障稳定
lagrange插值法:求拉格朗日插值多项式matlab实现(内附代码及例题)
关于拉格朗日插值法相关理论知识,在这里小编不在赘述,请不明白的小伙伴自行百度。小编只负责给出matlab源码。
**例题:**看下面例题(如图):
matlab代码:
%%%% 求拉格朗日多项式及基函数 %%%%
%%%% Liu Deping %%%%
%%%% 2020.06.14 %%%%
%输入的量:n+1个节点(x_i,y_i)(i = 1,2, ... , n+1)横坐标向量X,纵坐标向量Y
%输出的量:n次拉格朗日插值多项式L和基函数l
X=input('请输入横坐标向量X:\nX='); %输入的数据为一维数组,例如:[1,3,4,5](下同);
Y=input('请输入纵坐标向量Y:\nY=');
m = length(X);
L = ones(m,m);
for k = 1 : m
V = 1;
for i = 1 : m
if k ~= i
V = conv(V,poly(X(i))) / (X(k) - X(i));
end
end
L1(k, :) = V;
l(k, :) = poly2sym(V);
end
fprintf('基函数为:\n');
for k=1:m
fprintf('q%d(x)=%s\n',k,l(k));
end
L = Y * l;
fprintf('拉格朗日多项式为:\nP(x)=%s\n',L);
操作步骤:
step1:复制代码至matlab编辑器,并点击运行;
step2:按照提示,输入坐标向量(输入完坐标向量后敲回车键);
上面例题的运行结果如下:
最后,喜欢的小伙伴点个赞,有兴趣的可以加个关注,小编将持续更新~~
发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/206828.html原文链接:https://javaforall.cn
【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛
【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...