?? newpoly.m
字號(hào):
function [C,D]=newpoly(x,y) %牛頓插值
n=length(x);
D=zeros(n);
D(:,1)=y';
for j=2:n %計(jì)算差商矩陣
for k=j:n
D(k,j)=(D(k,j-1)-D(k-1,j-1))/(x(k)-x(k-j+1));
end
end
C=D(n,n);
for k=(n-1):-1:1 %構(gòu)造插值多項(xiàng)式
C=conv(C,poly(x(k)));
m=length(C);
C(m)=C(m)+D(k,k);
end
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -