数值分析实验报告模板(1).docxVIP

数学与应用数学专业 学生实验报告 学 号： 20054053017 姓名： 王泽豆 班 级： B20数学2班 实验名称： 多项式插值计算及其收敛实验 数值分析实验课程实验报告 实验名称 多项式插值计算及其收敛实验 编号 　02 姓名 王泽豆 班级 B20数学2班 学号 20054053017 同组人姓名 无 同组人学号 无 任课教师 王小霞 指导教师 王小霞 实验地点 第三教学楼J3201 实验时间 2022年4月2日 一、实验环境： 1、实验硬件：windows10 2、实验软件：Matlab 二、实验目的 1.使学生深入理解Langrage插值法和Newton插值法及两者之间的异同，能用Matlab语言编写实现Langrage插值法和Newton插值法的程序，并用所编写的程序进行插值计算、验证Runge现象、分析插值多项式的收敛性; 2.?使学生深入理解教材介绍的两种分段低次插值法，熟悉掌握函数interp1的使用;3.?调用函数interp1,采用不同方法进行插值计算,对教材介绍的几种分段低次插值法进行分析比较。 实验内容（包括所建数学模型、算法、编写程序、运行结果（图和表形式）） 1.根据Matlab语言特点，描述Langrage插值法和Newton插值法。2.用Matlab语言编写按Langrage插值法和Newton插值法计算插值的程序3.对y=11+x2,x∈[-5.5]，?分别取3个，5个、9个、11个等距节点，用编写的程序进行插值计算并画图，以验证Runge现象、分析插值多项式的收敛性。4.?用函数interp1,?对y=11+x2，x∈[-5.5]，用n=11个节点(等分)作分段线性插值、分段 Langrage算法描述： 若n-1次多项式lj(x)(j=1,……,n)在n个节点x1x2 li(xk)=fx=1, k=j 就称这n个n-1次多项式l1(x),...ln(x)为节点x1, （2）插值基函数为 lk=x- 显然它满足条件（1）。于是，满足条件（1）的插值多项式Ln(x) Ln(x)=k=1n Ln(x)就称为拉格朗日Lagrange Newton插值多项式算法描述： Newton插值多项式的表达式如下: Nnx=c ci=fx0 (1)计算(4)计算f(0.,-.?)2+fx.-.)(5)计算f(x,5,-.mx.) Lagrange的算法程序： function y=lagrange(x0,yo,x) n=length(x0); m=length(x); for i=1:m z=x(i); s=0.0; for k=l:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end Newton插值法如下： function f = Newton(x,y, x0) syms t; if (length(x) == length(y)) n= length(x); c(1:n)=0; c(1:n) = 0; else display(???ùê?è?μ?xó?yμ???êy2??àμè! ); return end f = y(1); y(1)=0; 1=1; for(i=1:n-1) for (j=i+1:n) y1(j) = (y(j)-y(i)/(x(j)-x(i)); end C(i) = y1(i+1); 1 = 1*(t-x(i)); f=f+ c(i)*1; simplify(f); if(i==n-1) if(nargin == 3) f = subs(f,t,x0); else f = collect(f); f = vpa(f, 6); end end end 插值计算 Langrange算法 1.取3个等距点进行插值计算： 编写m函数如下： x0=linspace(-5,5,3); y0=1./(1+x0.^2); x=-5:0.1:5; y=1./(1+x.^2); y1=lagrange(x0,y0,x); plot(x,y,r,x,y1,b) 2.取5个等距点进行插值计算： x0=linspace(-5,5,5); y0=1./(1+x0.^2); x=-5:0.1:5; y=1./(1+x.^2); y1=lagrange(x0,y0,x); plot(x,y,r,x,y1,b) 3.取9个等距点进行插值计算 x0=linspace(-5,5,9); y0=1./(1+x0