数学实验——性代数方程组的数值解.docxVIP

数学实验报告 实验5 线性代数方程组的数值解法 实验5 线性代数方程组的数值解法 分1 黄浩 2011011743 实验目的 学会用MATLAB软件数值求解线性代数方程组，对迭代法的收敛性和解的稳定性作初步分析； 通过实例学习用线性代数方程组解决简化的实际问题。 实验内容 《数学实验》第二版（问题1） 问题叙述： 通过求解线性方程组A1x=b1,A2x=b A A2是n阶希尔伯特矩阵，b1，b2分别是A1,A 编程构造A1（A2可直接用命令产生）和b1，b2；你能预先知道方程组A1x=b1和A2x=b2 令n=5,7,9，…，计算A1和A2的条件数。为观察他们是否病态，做以下试验：b1，b2不变，A1和A2的元素A1(n,n)，A2(n,n)分别加扰动ε后求解；A1和A2不变，b1，b2的分量b1(n),b2(n)分别 经扰动得到的解记做x，计算误差x-xx 模型转换及实验过程： （1）小题. 由b1,b2为A1,A2的行和，可知方程组A1x=b1和A2x=b2的精确解均为n行全1的列向量。在n=5的情况下，用matlab编程（程序见四.1），构造A1,A2和b1，b x x1 e1 x2 e2 1 1.000000000000080 0.000000000000080 1.000000000000020 0.000000000000020 1 0.999999999999463 -0.000000000000537 0.999999999999558 -0.000000000000442 1 1.000000000000970 0.000000000000970 1.000000000002010 0.000000000002010 1 0.999999999999288 -0.000000000000712 0.999999999996863 -0.000000000003137 1 1.000000000000180 0.000000000000180 1.000000000001560 0.000000000001560 由上表可见，当n=5时，所得的解都接近真值，误差在10^-12的量级左右。 （2）小题 分别取n=5,7,9,11,13,15，计算A1和A2的条件数c1和c2，（程序见四.2） n 5 7 9 11 13 15 c1 3.574E+05 8.739E+07 2.274E+10 6.518E+12 2.059E+15 8.291E+17 c2 4.766E+05 4.754E+08 4.932E+11 5.227E+14 6.263E+17 3.675E+17 由上表可见，二者的条件数都比较大，可能是病态的。为证实A1和A2是否为病态，先保持b不变，对A(n,n)做扰动，得到该情况下的高斯消元解，（ (为使结果清晰简洁，在此仅列出n=5,9,13的情况，n=7,11,15略去) ε=10^-10时： n xx1 xx2 5 0.999999928499402 0.999999937000303 1.000000251085169 1.000001259993930 0.999999670414583 0.999994330027334 1.000000191667702 1.000008819957473 0.999999958333144 0.999995590021263 n xx1 xx2 9 0.999997554805094 0.999982927803568 1.000014924430326 1.001229197054528 0.999960322336926 0.978489065794609 1.000060015431705 1.157746769617711 0.999943507927902 0.408449857681677 1.000033887641123 2.230423880873159 0.999987348370823 -0.435494123652275 1.000002687862696 1.878873743709766 0.999999751193406 0.780281609027228 n xx1 xx2 13 1.000778403386929 0.999999876758327 0.993774398727431 1.000026118511664 1.022687241045961 0.998779479002650 0.950184025603932 1.023304442731751 1.073410603535381 0.767841043550333 0.923508105988653 2.366938459461343 1.05