matlab关于牛顿插值法的简单应用

2021-11-02 12:58:50 阅读：247 来源： 互联网

标签：end 插值法牛顿 0.5 0.4 0.7 0.6 zeros matlab

差商表

差值多项式


X=[1,2,3,4,5,6]; %X矩阵
Y=X.^3-4.*X;     %Y矩阵
N=6      %总共六个结点
f=zeros(N,N); %建立一个N维零方阵

for k = 1 : N
    f(k,1)=Y(1,k); %先给第一列赋值Y

end
 
for i = 2:N   %列
    for k = i:N %行
        
        f(k,i)=(f(k,i-1)-f(k-1,i-1))/(X(k)-X(k-i+1));  
        %这个X（k-i+1）要仔细琢磨一下其余的直接抄公式
        

    end
end
syms x;
y=Y(1)+f(2,2)*(x-1)+f(3,3)*(x-1)*(x-2)+f(4,4)*(x-1)*(x-2)*(x-3);

再举一个例子
用牛顿差值公式计算ln(0.54)

X=[0.4,0.5,0.6,0.7,0.8];
Y=[-0.916291,-0.693147,-0.510826,-0.357765,-0.223144];
N=5
f=zeros(N,N);

for k = 1 : N
f(k,1)=Y(1,k);

end

for i = 2:N %列
for k = i:N %行

    f(k,i)=(f(k,i-1)-f(k-1,i-1))/(X(k)-X(k-i+1));
    

end

end
x=0.54;
y=Y(1)+f(2,2)(x-0.4)+f(3,3)(x-0.4)(x-0.5)+f(4,4)(x-0.4)(x-0.5)(x-0.6)+…
f(5,5)(x-0.4)(x-0.5)(x-0.6)(x-0.7);

标签：end,插值法,牛顿,0.5,0.4,0.7,0.6,zeros,matlab
来源： https://blog.csdn.net/JIAHAORA/article/details/121097966

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

matlab关于牛顿插值法的简单应用