数值分析实验报告二插值.docVIP

姓名 郭飞 评分 实验报告 实验报告 课程名称： 数值分析 课题名称： 比较不同插值方法的稳定性及运算量 专 业： 勘查技术与工程 姓 名： 郭飞 班 级： 061132班 完成日期： 2015 年 11 月 10 日 实验报告 一、实验名称 比较牛顿插值拉格朗日插值分段埃尔米特以及三次样条插值四种插值方法的稳定性以及运算量。 二、 实验目的 （1）培养编程与上机调试能力 （2）比较牛顿插值拉格朗日插值分段埃尔米特以及三次样条插值四种插值方法的稳定性以及运算量。 三、 实验要求 利用matlab编写算法，比较四种插值方法的优劣性。 四、 实验原理 用多项式近似逼近被差值函数，选用一些特殊点近似求得 五、 实验题目 比较四种种算法得出结论 六、 实验步骤 （）利用Matlab编程牛顿插值拉格朗日插值分段埃尔米特差值算法 （3）利用MATLAB程序自带的三次样条算法 七、 实验整体流程图或算法 八、 程序及其运行结果 (一)拉格朗日差值 function [ result,Ln] = lagelangri( X,Y,x0 ) if(length(X)~=length(Y)) disp(2??μ?úμ???êyó??úμ??μ??êy2?í?￡?????D?ê?è?￡?); return; end syms x l ; n=length(X); Ln=0; for i=1:n l=1; for j=1:n if i~=j l=l*(x-X(j))/(X(i)-X(j)); end end Ln=Ln+l*Y(i); end Ln=simplify(Ln) result=subs(Ln,x,5) 2,3,4 ans = 2 ans = 3 ans = 4 x=[1 2 3] x = 1 2 3 Unt159led4 Ln = x + 3 result = 8 ans = 8 Unt159led4 Ln = x + 3 result = 8 ans = 8 （二）牛顿差值 function [result, Nn] = niudun(X,Y,a) if(length(X)~=length(Y)) disp(2??μ?úμ???êyó??úμ??μ??êy2?í?￡?????D?ê?è?￡?); return; end syms x Nn; w=1; n=length(X); Nn=0; juncha(:,1)=Y; for j=2:n m=diff(juncha(:,j-1)); for k=1:(j-1) juncha(:,j)=diff(juncha(:,j-1)); for i=1:(n-j+1) juncha(i,j)=juncha(i,j)/(X(i+j-1)-X(i)); end end for j=1:n Nn=Nn+juncha(1,j)*w; w=w*(x-X(j)); end Nn=simplify(Nn) result=subs(Nn,x,a) end （三）艾尔米特差值 function f=Hermite(x,y,y_1,x0) syms t; f=0.0; if(length(x)==length(y)) if(length(y)==length(y_1)) n=len