《数值分析试验一》由会员分享，可在线阅读，更多相关《数值分析试验一（23页珍藏版）》请在金锄头文库上搜索。
1、 数值分析第一次实验报告 姓名： 学号：实验1： 1. 实验项目的性质和任务通过上机实验，使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步理解。2教学内容和要求 1）对高阶多多项式 编程求下面方程的解 并绘图演示方程的解与扰动量的关系。（实验2.6） 2）对，生成对应的Hilbert矩阵，计算矩阵的条件数；通过先确定解获得常向量b的方法,确定方程组 最后，用矩阵分解方法求解方程组，并分析计算结果。（第三章，实验题4）3）对函数 的Chebyshev点 编程进行Lagrange插值，并分析插值结果。（第四章 实验1）项目涉及核心知识点 病态方程求解、矩阵分解和方程组求解、Lagra
2、nge插值。重点与难点 算法设计和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 m=1; while m n=2; A=zeros(20,20); while n AA = 1.0
3、e+003 * Columns 1 through 10 0 0 0 0 0 0 0 0 0 0 -0.0000 0.0000 0 0 0 0 0 0 0 0 -0.0000 0.0000 -0.0000 0 0 0 0 0 0 0 -0.0000 0.0000 -0.0000 0.0000 0 0 0 0 0 0 0.0000 -0.0000 0.0000 -0.0000 0.0000 0 0 0 0 0 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0 0 0 0 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0
4、000 0.0000 0 0 0 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 0 0 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000
5、 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0001 -0.0003 0.0006 -0.0007 0.0005 0.0000 -0.0000 0.0000 -0.0001 0.0005 -0.0027 0.0096 -0.0223 0.0348 -0.0361 0.0000 -0.0000 0.0000 -0.0004 0.0030 -0.0098 0.0080 0.0593 -0.2570 0.5154 0.0000 -0.0000 0.0000 -0.0001 0.0005 -0.0029 0.0095 -0.0171 0.0086 0.0347 0
6、.0000 -0.0000 0.0000 -0.0000 0.0003 -0.0016 0.0059 -0.0133 0.0145 0.0094 0.0000 -0.0000 0.0000 -0.0001 0.0009 -0.0042 0.0118 -0.0182 0.0082 0.0185 0.0000 0.0000 -0.0000 0.0002 -0.0027 0.0187 -0.0762 0.1806 -0.2249 0.0813 0.0000 0.0000 -0.0000 0.0001 -0.0017 0.0120 -0.0497 0.1224 -0.1699 0.1064 0.0000 -0.0000 0.0000 -0.0003 0.0028 -0.0137 0.0371 -0.0464 -0.0164 0.1243 Columns 11 through 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
