数值计算实验牛顿前插和后插插值.pdfVIP

数值计算实验二 姓名:方小开 学号:20060810202 班级:计科0602 一. 实验目的： 1、 差分的matlab 实现； 2、 Newton 插值的matlab 实现； 二. 实验原理： MATLAB 在线性代数，矩阵分析，数值及优化，数理统计和随机信号分析，电路系统， 系统动力学，信号与图像处理，控制理论分析和系统设计，过程控制，建模和仿真，通信系 统，等有广泛的应用。它具有功能强大，界面友好，语言自然即开放性等特点。 三．试验环境 MATLAB7.0 四. 试验过程及现象： 1、牛顿插值公式： 把下面的matlab 程序在matlab 中建立M-file 文件并保存； function [d,v1]=newtonjz(x,y,v) %d 插商表 v 是要插入x v1 是插入的y 值 n=length(x); d=zeros(n,n); d(:,1)=y; for j=2:n for i=j:n d(i,j)=(d(i,j-1)-d(i-1,j-1))/(x(i)-x(i-j+1)); end end w=1; v1=d(1,1); for i=2:n w=w*(v-x(i-1)); v1=v1+d(i,i)*w; end 分别给x,y 赋初值，并调用Newton 插值函数得到结果如下： x=[0.40,0.55,0.65,0.80,0.90,1.05] ; y=[0.41075,0.57815,0.69675,0.88811,1.02652,1.25382] ; [z,xy]=newtonjz(x,y,0.596); z = 0.4108 0 0 0 0 0 0.5782 1.1160 0 0 0 0 0.6967 1.1860 0.2800 0 0 0 0.8881 1.2757 0.3589 0.1973 0 0 1.0265 1.3841 0.4335 0.2130 0.0312 0 1.2538 1.5153 0.5249 0.2287 0.0314 0.0003 xy = 0.63192 Newton 前插公式： 把Newton 前插公式的matlab 程序写在matlab 中建立M-file 文件并保存； function [d,v1]=newtonBefore(x,y,t) n=length(x); d=zeros(n,n); d(:,1)=y; for j=2:n for i=1:n-j+1 d(i,j)=(d(i+1,j-1)-d(i,j-1)); end end w=1; m=1; v1=d(1,1); for i=2:n w=w*(t-i+2); m=m*(i-1); v1=v1+d(1,i)*(w/m); end 分别给x,y 赋初值，并调用Newton 前插函数得到结果如下； x=[1 1.05 1.10 1.15 1.20 1.25 1.30]; y=[1 1.0247 1.04881 1.07238 1.09544 1.11803 1.14017]; [z,qc]=newtonBefore (x,y,0.2); qc =1.00499226380800 3、Newton 后插公式： 把Newton 后插公式的matlab 程序写在matlab 中建立M-file 文件并保存； function [d,v1]=newtonAfter(x,y,t) n=lengt