?? lagrange.m
字號:
function s=Lagrange(x,y,x0)
%Lagrange插值,x和y為已知的插值點及其函數值
%x0為要求的插值點的x值
nx=length(x);
ny=length(y);
if nx~=ny
% warning('矢量x與y的長度應該相同')
return
end
m=length(x0);
%按照公式,對要求的插值點矢量x0的每個元素進行計算
for i=1:m
t=0.0;
for j=1:nx
u=1.0;
for k=1:nx
if k~=j
u=u*(x0(i)-x(k))/(x(j)-x(k));
end
end
t=t+u*y(j);
end
s(m)=t;
end
return
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -