matlab实现newton插值

男娘i 2022-01-28 15:01 398阅读 0赞

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2h1YW5namlhYWFhYQ_size_16_color_FFFFFF_t_70]  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2h1YW5namlhYWFhYQ_size_16_color_FFFFFF_t_70 1]  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2h1YW5namlhYWFhYQ_size_16_color_FFFFFF_t_70 2]  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2h1YW5namlhYWFhYQ_size_16_color_FFFFFF_t_70 3]

newton.m

function s=newton(x,y,t)
    %计算输入的x y是否长度相等
    if(length(x)==length(y))
        n=length(x);
    else
        disp('x和y的维数不相等！');
        return;  
    end
     
    syms p;
    s=y(1);
    dd=zeros(1,100);
    dxs=1;
     
     
    for i=1:n-1
        for j=i+1:n
            dd(j)=(y(j)-y(i))/(x(j)-x(i));  %用来构造均差
        end
        temp(i)=dd(i+1);
        dxs=dxs*(p-x(i));
        s=s+temp(i)*dxs;
        y=dd;     %第一次将一阶均差赋给y，然后通过上面的dd(j)来构造二阶均差
        s=simplify(s);
        s=subs(s,'p',t);
        s=vpa(s,6);
    end

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2h1YW5namlhYWFhYQ_size_16_color_FFFFFF_t_70]: /images/20220128/30631c59f6044652b9f0f3764d81074a.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2h1YW5namlhYWFhYQ_size_16_color_FFFFFF_t_70 1]: /images/20220128/cf46278fb3b9404a8d5b7f9241804a74.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2h1YW5namlhYWFhYQ_size_16_color_FFFFFF_t_70 2]: /images/20220128/11766fc0f7454b3fa7cd4cf9f9063fe1.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2h1YW5namlhYWFhYQ_size_16_color_FFFFFF_t_70 3]: /images/20220128/944edc1141ca4b44b67306de85c5c027.png