matlab差商chaS,牛顿插值MATLAB算法.doc
MATLAB程序设计期中作业
——编程实现牛顿插值
成员:刘 川(P091712797) 签名_____
汤 意(P091712817) 签名_____
王功贺(P091712799) 签名_____
班级:2009信息与计算科学
学院:数学与计算机科学学院
日期:2012年05月02日
牛顿插值的算法描述及程序实现
一:问题说明
在我们的实际应用中,通常需要解决这样的问题,通过一些已知的点及其对应的值,去估算另外一些点的值,这些数据之间近似服从一定的规律,于是,这就引入了插值法的思想。
插值法是利用函数f (x)在某区间中若干点的函数值,作出适当的特定函数,在这些点上取已知值,在区间的其他点上用这特定函数的值作为函数f (x)的近似值。如果这特定函数是多项式,就称它为插值多项式。利用插值基函数很容易得到拉格朗日插值多项式,公式结构紧凑,在理论分析中甚为方便,但当插值节点增减时全部插值基函数均要随之变化,整个公式也将发生变化,这在实际计算中是很不方便的,为了克服这一缺点,提出了牛顿插值。
二:算法分析
newton插值多项式的表达式如下:
其中每一项的系数ci的表达式如下:
即为f (x)在点处的i阶差商,(,),由差商的性质可知:
牛顿插值的程序实现方法:
第一步:计算。
第二步:计算牛顿插值多项式中,,得到n个多项式。
第三步:将第二步得到的n个多项式相加,得到牛顿插值多项式。
第四步:利用所得到的插值多项式,估算取其它值时的值。
第五步:作出所求多项式在插值结点周围的函数图像。
三:编程实现
function [p2,z]=newTon(x,y,t)
%输入参数中x,y为元素个数相等的向量,t为待估计的点,可以为数字或向量。
%输出参数中p2为所求得的牛顿插值多项式,z为利用多项式所得的t的函数值。
n=length(x);
chaS(1)=y(1);
for i=2:n
x1=x;y1=y;
x1(i+1:n)=[];
y1(i+1:n)=[];
n1=length(x1);
s1=0;
for j=1:n1
t1=1;
for k=1:n1
if k==j
continue;
else
t1=t1*(x1(j)-x1(k));
end
end
s1=s1+y1(j)/t1;
end
chaS(i)=s1;
end
b(1,:)=[zeros(1,n-1) chaS(1)];
cl=cell(1,n-1);
for i=2:n
u1=1;
for j=1:i-1
u1=conv(u1,[1 -x(j)]);
cl{i-1}=u1;
end
cl{i-1}=chaS(i)*cl{i-1};
b(i,:)=[zeros(1,n-i),cl{i-1}];
end
p2=b(1,:);
for j=2:n
p2=p2+b(j,:);
end
if length(t)==1
rm=0;
for i=1:n
rm=rm+p2(i)*t^(n-i);
end
z=rm;
else
k1=length(t);
rm=zeros(1,k1);
for j=1:k1
for i=1:n
rm(j)=rm(j)+p2(i)*t(j)^(n-i);
end
z=rm;
end
end
plot(t,z,y,x,y,*r)
四:实例验证
clc
clear
x=[0.4 0.55 0.65 0.80 0.90 1.05];
y=[0.41075 0.57815 0.69675 0.88811 1.02652 1.25386];
t=0.4:0.1:1.05;
[u,v]=newTon(x,y,t)
执行结果:
u =
0.0085 0.0032 0.1587 0.0073 0.9971 0.0004
v =
0.4108 0.5211 0.6367 0.7586 0.8881 1.0265 1.1752
1.3356 1.5095
则所求得的牛顿多项式为:
牛顿多项式的函数图像及已知节点在坐标中的显示如下:
五:结果分析
本程序给出了计算牛顿插值多项式的函数,通过调用函数可以求得牛顿多项式与待估算点的值,作出了节点及待求多项式的函数图像,能够比较清晰的通过图像显示出来,总体来说,计算结果是比较理想的,达到了我们的目的。
然而程序在实现过程中,依旧存在着一些不足之处,总体反应在灵活性方面,参数的输入必须为三个,否则程序会出错,有时候我们仅需要得到牛顿多项式,而不需要去估算某个具体的x对应的函数值。另外,本程序的函数在实现过程中会给出多项式函数的函数图像,没有设置参数对其判断是否需要。因此,这些都是有待改进的地方。
展开阅读全文
matlab差商chaS,牛顿插值MATLAB算法.doc相关推荐
- 重读微积分(六):差商与牛顿插值
一 理解极限 二 微分学 本系列所有代码皆用R语言完成. 4 差商与牛顿插值 如果采取间隔不等的采样,差商会变得稍显复杂,对于x0,x1,-,xnx_0,x_1,\ldots,x_nx0,x1,- ...
- matlab全域基函数,多项式函数插值:全域多项式插值(一)单项式基插值、拉格朗日插值、牛顿插值 [MATLAB]...
全域多项式插值指的是在整个插值区域内形成一个多项式函数作为插值函数.关于多项式插值的基本知识,见"计算基本理论". 在单项式基插值和牛顿插值形成的表达式中,求该表达式在某一点处的值 ...
- 牛顿插值法 matlab m文件,牛顿插值法matlab程序
<牛顿插值法matlab程序>由会员分享,可在线阅读,更多相关<牛顿插值法matlab程序(3页珍藏版)>请在人人文库网上搜索. 1.计算方法数值实验报告班级090712学号0 ...
- matlab自带拉格朗日插值,MATLAB实现拉格朗日插值
<MATLAB实现拉格朗日插值>由会员分享,可在线阅读,更多相关<MATLAB实现拉格朗日插值(14页珍藏版)>请在人人文库网上搜索. 1.实用标准文档文案大全数值分析上机报告 ...
- 计算方法 差商与牛顿插值
题目 牛顿插值 求差商代码 #include <iostream> #include <algorithm> #include <vector> #include ...
- 牛顿插值实验报告C语言,牛顿插值C语言.doc
<数值分析>实验报告 学号: 20120921101 姓名: 孙腾 班级: 计算1201 日期:2012/10/20 题目:Newton插值法的函数表如下 x 0.0 0.1 0.2 0. ...
- 艾特肯法方程解matlab程序,牛顿迭代法matlab代码
牛顿法 迭代公式: x(k1) xk [2 f (x(k) )]1f (x(k) ) Matlab 代码: function [x1,k] =newton(x1,eps) hs=inline('(x ...
- 基于MATLAB的三维数据插值拟合与三次样条拟合算法(附完整代码)
目录 一. 三维插值 例题1 二. 高维度插值拟合 格式一 格式二 格式三 格式四 格式五 例题2 三. 单变量三次样条插值 例题3 例题4 四. 多变量三次样条插值 例题6 一. 三维插值 首先三维 ...
- 牛顿插值与正交多项式曲线拟合的Matlab实现
一.相关知识提要 1.Newton插值法 2. 正交多项式的方法作曲线拟合 二.代码实现 1.Newton插值法 function [] = Newton_Chazhi(x ...
最新文章
- 《我想进大厂》之Spring夺命连环10问
- python平方数迭代器_对python中的高效迭代器函数详解
- 特斯拉Model3交付危机:周产5000辆仍需排队一年,造空者花式补刀
- 设置CodeRush Xpress的线条颜色
- oracle开放查询表权限_oracle 查询当前用户的表和其他用户的表
- 0-1背包问题暴力递归
- 24小时临时邮箱_实用网站 | 临时邮箱,悄悄来悄悄去~
- 一块电路板从“出生”到“成年”
- Android之Bitmap的内存优化方案总结
- redux-4-ways
- 算法十——深度优先搜索和广度优先搜索
- androidstudio build tools安装_Android Studio4.0 安装及配置
- drupal ajax 分页,在drupal中使用Ajax时要保持一个404
- pass by value pass by reference
- deepin系统安装成功了之后重启电脑没有deepin启动选项的简单解决办法
- Django的数据迁移
- iso27001认证怎么申请?(iso27001认证基本条件和详细流程)
- centos查看yum上jdk的版本
- PHPwamp开发环境搭建详解
- 基于张量的多元多阶马尔科夫多模态预测方法