MATLAB实现牛顿插值的源程序

function yi=Newton(x,y,xi)
%Newton插值方法，给定一系列插值的点(x,y)，得到在x=xi处的，牛顿插值多项的值yi
n=length(x);
m=length(y);
if m~=n
    error('x,y的长度不一样,请重新输入!');
    return
end
A=zeros(n);    %定义差商表
A(:,1)=y;      %差商表第一列为y
for j=2:n               %j为列标
    for i=1:(n-j+1)     %i为行标
        A(i,j)=(A(i+1,j-1)-A(i,j-1))/(x(i+j-1)-x(i));   %计算差商表
    end
end
%根据差商表,求对应的牛顿插值多项式在x=xi处的值yi
N(1)=A(1,1);
for j=2:n
    T=1;
    for i=1:j-1
        T=T*(xi-x(i));
    end
    N(j)=A(1,j)*T;
end
yi=sum(N);   %将x=xi带入牛顿插值多项式，得到的yi的值
end

MATLAB实现牛顿插值的源程序相关推荐

matlab全域基函数,多项式函数插值：全域多项式插值（一）单项式基插值、拉格朗日插值、牛顿插值 [MATLAB]...
全域多项式插值指的是在整个插值区域内形成一个多项式函数作为插值函数.关于多项式插值的基本知识,见"计算基本理论". 在单项式基插值和牛顿插值形成的表达式中,求该表达式在某一点处的值 ...
matlab差商chaS,牛顿插值MATLAB算法.doc
MATLAB程序设计期中作业 --编程实现牛顿插值成员:刘川(P091712797) 签名＿＿＿＿＿汤意(P091712817) 签名＿＿＿＿＿王功贺(P091712799) 签名＿＿＿＿＿ ...
Matlab实现线性插值、抛物插值、牛顿插值、拉格朗日插值、分段抛物插值、分段线性插值
目录线性插值原理流程图代码抛物插值原理流程图代码拉格朗日插值代码牛顿插值原理代码分段线性插值代码线性插值原理流程图单个点的线性插值代码 X=[0.2 0.4]; ...
牛顿插值与正交多项式曲线拟合的Matlab实现
一.相关知识提要 1.Newton插值法 2. 正交多项式的方法作曲线拟合二.代码实现 1.Newton插值法 function [] = Newton_Chazhi(x ...
艾特肯法方程解matlab程序,牛顿迭代法matlab代码
牛顿法迭代公式: x(k1) xk [2 f (x(k) )]1f (x(k) ) Matlab 代码: function [x1,k] =newton(x1,eps) hs=inline('(x ...
hermit插值 matlab,埃尔米特(Hermite)插值及其MATLAB程序
%hermite.m %求埃尔米特多项式和误差估计的MATLAB主程序 %输入的量:X是n+1个节点(x_i,y_i)(i = 1,2, ... , n+1)横坐标向量,Y是纵坐标向量, %以f'(x ...
matlab用牛顿差值计算三次差值多项式,计算方法用Newton插值多项式求函数的近似值.docx...
计算方法用Newton插值多项式求函数的近似值计算方法课程设计题目: 用Newton插值多项式处理磁化曲线学院: 理学院班级: 学生姓名: 学生学号: 指导教师: 2017 ...
MATLAB之多项式插值
MATLAB之多项式插值一.算法原理 1.插值问题定义当精确函数 y = f(x) 非常复杂或未知时,在区间[a,b]上一系列节点x0,x1,x2,......xn处测得函数值f(x0).f(x1 ...
matlab自带拉格朗日插值,MATLAB实现拉格朗日插值
<MATLAB实现拉格朗日插值>由会员分享,可在线阅读,更多相关<MATLAB实现拉格朗日插值(14页珍藏版)>请在人人文库网上搜索. 1.实用标准文档文案大全数值分析上机报告 ...

MATLAB实现牛顿插值的源程序

MATLAB实现牛顿插值的源程序相关推荐

最新文章

热门文章