2019-吉林-《计算方法》大作业.docVIP

《计算方法》大作业 一、构造次数不超过三次的多项式P3（X），使满足： P3（0）= 1;P3（1）=0； P3′（0）=P3′（1）=0。 （10分） 解：求解这种带有导数的多项式逼近，我们考虑埃尔米特（Hermite）插值法。 采用构造差值基函数的方法，设 上式中，均为插值基函数， 其中，代入P3（0）= 1;P3（1）=0；P3′（0）=P3′（1）=0。求得，a=1，b=2； 所以综上 二、设f(xi)=i(i=0,1,2),构造二次式p2(x)，使满足： p2(xi)=f(xi)(i=0,1,2) （10分） 解：本题没有明确给出xi的取值，为便于计算，我们取xi=i（i=0,1,2）。 依据公式 化简，利用MATLAB进行快速计算，我们得到多项式为 MATLAB代码如下： a=[0,1,2];b=[0,1,2]; syms x y y1 y2 y=0; for i=1:length(a) y1=1;y2=1; for j=1:length(a) if i~=j y1=y1*(x-a(j)); y2=y2*(a(i)-a(j)); end end y=y1*y2*b(i)+y; end collect(y) 三、设节点xi=i(i=0,1,2,3),f(0)=1,f(1)=0,f(2)=-7,f(3)=26,构造次数不超过3次的多项式p3(x),满足 p3(xi)=f(xi),i=0,1,2,3 （10分） 解：依据公式 化简，利用MATLAB进行快速计算，我们得到多项式为 MATLAB代码如下： a=[0,1,2,3];b=[1,0,-7,26]; syms x y y1 y2 y=0; for i=1:length(a) y1=1;y2=1; for j=1:length(a) if i~=j y1=y1*(x-a(j)); y2=y2*(a(i)-a(j)); end end y=y1*y2*b(i)+y; end collect(y) 对于上题的问题，构造Newton插值多项式。 （10分） 解：采用newton插值法，如下。 由均差计算公式，计算结果如下表所示 表1 第四题newton插值法 i 一阶均差 二阶均差 三阶均差 0 0 1 1 1 0 -1 2 2 -7 -7 -3 3 3 26 33 20 23/3 易得P3(x)=1-x-3x(x-1)+23/3x(x-1)(x-2) =(23*x^3)/3 - 26*x^2 + (52*x)/3 + 1; 五、构造三次多项式P3（X）满足：P3（0）= P3（1）=0， P3′（0）=P3′（1）=1。 （10分） 解：本题与第一问基本类似，考虑埃尔米特（Hermite）插值法。 采用构造差值基函数的方法，设 上式中，均为插值基函数， 其中，代入P3（0）= 1;P3（1）=0；P3′（0）=P3′（1）=0。求得，c=1； 所以综上 在19题的插值条件上，另加上p4(2)=1,试用构造满足插值条件的四次插值多项式。 （10分） 解：作业中并未给出19题的具体内容，故而此题不可作答。 七、利用Doolittle分解法解方程组Ax=b即解方程组 （10分） 解：Doolittle分解法是对矩阵进行LU分解，其中L为单位下三角矩阵，U为上三角矩阵。 首先令该矩阵为A，利用MATLAB分解后的结果如下。 如上所示，再令 故x的求解公式为 求得方程解为 -4.4000 2.2732 1.6098 -0.0244 八、基于迭代原理证明 （10分） 解：构造数列{xn}如下： 易知待证的等式左边=lim xn； 于是题目变为求解该数列的极限。 易知该数列单调递增。 该数列有界。 x12,易得x22，同理，由归纳法易得，xn2; 综上，单调有界数列必有极限。故取该极限为x，则有 解得x=2.即