数值分析MATLAB上机实验文档格式.docx

1、包含大量计算算法的集合，拥有600多个工程中要用到的数学运算函数。4出色的图像处理功能，可以方便地输出二维图像，便于我们绘制函数图像。目录1 第一题4 1.1 实验目的4 1.2 实验原理和方法4 1.3 实验结果51.3.1 最佳平方逼近法51.3.2 拉格朗日插值法71.3.3 对比82 第二题92.1 实验目的92.2 实验原理和方法102.3 实验结果102.3.1 第一问102.3.2 第二问112.3.3 第三问113 第三题123.1 实验目的123.2 实验原理和方法123.3 实验结果124 MATLAB程序141第一题某过程涉及两变量x和y，拟分别用插值多项式和多项式拟合给

2、出其对应规律的近似多项式，已知xi与yi之间的对应数据如下：12345678910xy34.658840.371914.6448-14.2721-13.357024.823475.2795103.574397.484778.2392请用次数分别为3，4，5，6的多项式拟合并给出最好近似结果f（x）。请用插值多项式给出最好近似结果。1.1实验目的：学习逼近和插值的原理和编程方法，由给出的已知点构造多项式，在某个范围内近似代替已知点所代表的函数，以便于简化对未知函数的各种计算。1.2试验原理和方法：实验原理：拉格朗日插值法中先构造插值基础函数：lkx=j=0jknx-xixk-xi k=0,1,2

3、,n，然后构造出拉格朗日多项式：pnx=k=0nj=0jknx-xixk-xifxk。最佳平方逼近中，设逼近函数Pnx=a0+a1x+a2x2+anxn，逼近函数和真实函数之差r=Pnx-y，r1r2rn=11x1x2x1nx2n1xnxnna0a1an-y1y2yn，即：r=Xa-Y，根据最小二乘准则令i=0nri2=min，可以得到a=XTX-1XTY。实验方法：逼近法采用最佳平方逼近，依据最小二乘原则：i=0nri2=min，由已知条件采用离散型。插值法采用拉格朗日插值法。在逼近法中，由于是离散型的，所以法方程系数阵设计成求和。分别求出3、4、5、6次的多项式，逼近结果和真实值有一定差距

4、，最小二乘正是让这些差距达到最小，理论上多项式次数越高结果和真实值差距越小。拉格朗日插值法中“la=la*（p-x（j）/（x（k）-x（j）”语句实现的是我们通常书写的连乘形式拉格朗日插值多项式，但是表示不方便，而如果用“s=collect（s）”函数将其展开成降幂排列多项式以后，由于余项问题结果会和原本的多项式有偏差，这种偏差随着x的增大而增大。求出多项式后和题目中给出的参考点进行比较。最后，选择六次最佳平方逼近多项式和拉格朗日插值多项式（九次）进行比较，选取xi=a+ih=1+0.2*i（i=0,1,45），分别绘制两者的图像进行比较。1.3试验结果1.3.1最佳平方逼近法三次多项式：-

5、 1.033*x3 + 19.33*x2 - 94.48*x + 131.8拟合结果：55.617011.8960-5.5610-2.952013.525037.672063.291084.184094.153087.0000四次多项式：- 0.3818*x4 + 7.368*x3 - 42.14*x2 + 73.53*x+ 0.74539.121232.080210.0852-5.5638-2.730021.560261.1172100.5882115.457272.0450五次多项式：0.09807*t5 - 3.079*t4 + 34.5*t3 - 163.5*t2 + 304.7*t

6、- 139.533.219145.77429.0320-16.5003-8.906326.908370.9835100.073899.416474.5000六次多项式：0.01936*t6 - 0.5408*t5 + 5.114*t4 - 16.9*t3 - 0.867*t2 + 66.38*t - 18.734.505641.149413.2700-13.9486-12.225023.711473.9500105.169496.345678.4000对比可知，六次多项式拟合结果最好。1.3.2拉格朗日插值法插值多项式5.353*10（-5）*x9 - 0.003088*x8 + 0.0722

7、9*x7 - 0.8792*x6 + 5.932*x5 - 22.41*x4 + 50.11*x3 - 86.47*x2 + 113.5*x - 25.2注：此多项式为拉格朗日多项式的近似式，当x=10的时候偏差可以达到23以上。对比数据：1.50001.90002.30002.70003.10003.50003.90004.30004.700042.149841.462035.118224.385211.2732-1.7813-12.3006-18.1566-17.9069111213141516175.10005.50005.90006.30006.70007.10007.50007.90

8、00-11.02262.028419.854940.362661.084079.568893.7700102.3677插值结果：42.384041.494735.074224.360111.2792-1.7683-12.2977-18.1626-17.9118-11.02102.033319.856540.358461.079479.570993.7788102.3713其中红点表示参考点。1.3.3比较选取xi=a+ih=1+0.2*i（i=0,1,45），分别绘制六次多项式拟合和拉格朗日插值结果图：其中绿线表示拉格朗日插值多项式图像，蓝线表示六次多项式拟合图像。两者效果近似但后者比前者低三

9、次。2第二题用雅格比法与高斯赛德尔迭代法解下列方程组Ax=b1或Ax=b2，研究其收敛性。上机验证理论分析是否正确，比较它们的收敛速度，观察右端项对迭代收敛有无影响。（1）A行分别为A1=6,2,-1,A2=1,4,-2,A3=-3,1,4； b1=-3,2,4T；b2=100,-200,345T。（2） A行分别为A1=1,0,8,0.8,A2=0.8,1,0.8,A3=0.8,0.8,1；b1=3,2,1 T； b2=5,0,-10T。（3）A行分别为A1=1,3,A2=-7,1；b1=4,6T。2.1试验目的学习jacobi迭代法和GuassSeidel迭代法的原理和编程方法，研究方程组

10、系数阵和右边项对方程的解及其收敛性的影响，判断迭代法的收敛条件。2.2实验原理和方法将方程组系数阵A分解为A=D+L+U,其中D为对角阵，L为减去D的下三角阵，U为减去D的上三角阵。Jacobi迭代法中构造如下迭代公式：xk+1=-D-1L+Uxk+D-1b而Gauss-Seidel迭代法的迭代公式为：xk+1=-D+L-1Uxk+D+L-1b初始值直接选取为0。在判断其收敛性时，分别求解其迭代矩阵的谱半径G，G=max1inli， li为迭代矩阵的特征值。分别编写jacobi迭代及其收敛判别函数和Seidel迭代及其收敛判别函数。如果在初试迭代步数之内还未收敛就进行收敛判别，收敛判别的依据是迭代矩阵的谱半径是否小于1。比较同一方程组的jacobi迭代法和Seidel迭代法的结果是否相同，在达到精度要求后比较两种方法的迭代次数，比较哪一个的效率更高。比较方程组系数阵和等号右边的变化会对方程的解和收敛速度造成什么影响。如果迭代不收敛，那么考虑为什么不收敛，如果把方程组系数阵进行强对角占优处理，是否会收敛。2.3实验结果规定误差界：1e-42.