现在时间是:
当前位置:首 页 >> 信息论>> 教学区>> 文章列表

10211081 作业3

作者:10211081   发布时间:2012-09-29 21:01:35   浏览次数:314

 

先算出6x值对应的y值,程序如下:
X=[1 2 3 4 5 6]
Y=zeros(1,6);
for i=1:6
    Y(1,i)=X(i)^3-4*X(i);
end
Y
算出的y值如下:
-3     0    15    48   105   192
接着利用以下程序:
function [C,D]=newpoly(X,Y)
n=length(X);
D=zeros(n,n);
D(:,1)=Y';
for j=2:n
    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
    C=conv(C,poly(X(k)));
    m=length(C);
    C(m)=C(m)+D(k,k);
End
command window中输入
X=[1 2 3 4 5 6]
X =
 
     1     2     3     4     5     6
>> Y=[-3 0 15 48 105 192]
Y =
 
    -3     0    15    48   105   192
 输入:
[C,D]=newpoly(X,Y)
答案为:
C =
     0     0     1     0    -4     0
 
D =
    -3     0     0     0     0     0
     0     3     0     0     0     0
    15    15     6     0     0     0
    48    33     9     1     0     0
  105    57    12     1     0     0
   192    87    15     1     0     0
牛顿插值公式为:N6x=-3+3x-1+6x-1)(x-2+1x-1)(x-2)(x-3
 






上一篇:没有了    下一篇:没有了

Copyright ©2019    计算数学达人 All Right Reserved.

技术支持:自助建站 | 领地网站建设 |短信接口 |燕窝 版权所有 © 2005-2019 lingw.net.粤ICP备16125321号 -5