📄 拉格朗日插值多项式拟合.m
字号:
%拉格朗日插值多项式拟合
%___________________________________________________________
%调用合验证:
%调用函数lagran(X,Y)即得到拟合多项式C
%然后用polyval(C,z0)去验证,即求的是当x=zo的时候,多项式的值
%___________________________________________________________
% X是自变量,Y是对应的函数值
%拟合问题即是给定两组向量,分别是x量合y量
%然后用相应的多项式去拟合
%___________________________________________________________
function [C,L]=lagran(X,Y)
N=length(X)
n=N-1;
L=zeros(N,N);
for k=1:N
V=1;
for j=1:N
if k~=j
%其中poly(z0)得到的是一个多根为z0的多项式,比如poly(2),得到的是【1 -2】即x-2
V=conv(V,poly(X(j)))/(X(k)-X(j));
end
end
L(k,:)=V;
end
C=Y*L;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -