数值分析试验一.doc

数值分析第一次实验报告 姓名： 学号： 实验1： 实验项目的性质和任务 通过上机实验，使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步理解。 2．教学内容和要求 1）对高阶多多项式 编程求下面方程的解 并绘图演示方程的解与扰动量的关系。（实验2.6） 2）对，生成对应的Hilbert矩阵，计算矩阵的条件数；通过先确定解获得常向量b的方法,确定方程组 最后，用矩阵分解方法求解方程组，并分析计算结果。（第三章，实验题4） 3）对函数 的Chebyshev点 编程进行Lagrange插值，并分析插值结果。（第四章 实验1） 项目涉及核心知识点 病态方程求解、矩阵分解和方程组求解、Lagrange插值。 重点与难点 算法设计和matlab编程。 1）a．实验方案： 先创建一个20*50的零矩阵X，然后利用Matlab中的roots（）和poly（）函数将50个不同的ess扰动值所产生的50个解向量分别存入X矩阵中。然后再将ess向量分别和X的20个行向量绘图。即可直观的看出充分小的扰动值会产生非常大的偏差。即证明了这个问题的病态性。 b．编写程序： X=zeros(20,50); ve=zeros(1,21); ess=linspace(0,0.00001,50);k=1; while k=50 ve(2)=ess(k); X(1:20,k)=roots(poly(1:20)+ve); k=k+1; end m=1; while m=20 figure(m),plot(ess,X(m,:)); m=m+1; end C．实验结果分析和拓展 由上面的实验结果可以看出一个充分小的扰动值可以让方程的解产生非常大的偏差，而且这个偏差随着ess的变大偏差也随即变大。但可以看出在相对小的根处根比较稳定，也就是说这些根关于ess并不敏感，而在较大根处时，根很不稳定，即这些解关于ess的变化是敏感的。这就说明了这个问题本身就是一个病态问题，与算法好坏无关。 若扰动在x^18处，只要把程序中的ve（2）改为ve（3）即可，其图形和此类似。 d．实验结论： 高次多项式扰动求方程解问题是一个病态问题。 2）a．实验方案： 先创建一个20*20的零矩阵A，再通过给定解x和Hilbert矩阵求出列向量b，然后通过LU分解法求出方程HX=b的解X，然后将x-X’这一行向量存入A矩阵中，形成一循环，最后，如果Hilbert矩阵非病态的话，则可输出一个20*20的对角矩阵。 b.编写程序： n=2; A=zeros(20,20); while n=20 x=1:n; H=hilb(n); b=H*x; [L U]=lu(H); y=L\b;X=U\y; A(n,1:n)=x-X; n=n+1; end Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 4.455948e-017. Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 7.948463e-017. Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.798429e-016. Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 7.626119e-018. Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 6.040620e-017. Warning: Matrix is close to singular or b