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

数学实验报告实验5 线性代数方程组的数值解法 实验5 线性代数方程组的数值解法 分1 黄浩 2011011743 实验目的 学会用MATLAB软件数值求解线性代数方程组，对迭代法的收敛性和解的稳定性作初步分析； 通过实例学习用线性代数方程组解决简化的实际问题。 实验内容 《数学实验》第二版（问题1） 问题叙述： 通过求解线性方程组A1x=b1,A2x=b2，理解条件数的意义和方程组性态对解的影响，其中A1是n阶范德蒙矩阵，即 A1=1x0x02…x0n-11x1x12…x1n-1……………1xn-1xn-12…xn-1n-1 , xk=1+0.1k , k=0,1,…,n-1 A2是n阶希尔伯特矩阵，b1，b2分别是A1,A2的行和。 编程构造A1（A2可直接用命令产生）和b1，b2；你能预先知道方程组A1x=b1和A2x=b2的解吗？令n=5，用左除命令求解（用预先知道的解可验证程序）。 令n=5,7,9，…，计算A1和A2的条件数。为观察他们是否病态，做以下试验：b1，b2不变，A1和A2的元素A1(n,n)，A2(n,n)分别加扰动ε后求解；A1和A2不变，b1，b2的分量b1(n),b2(n)分别加扰动ε后求解。分析A与b的微小扰动对解的影响。ε取10^-10,10^-8,10^-6。 经扰动得到的解记做x，计算误差x-xx，与用条件数估计的误差相比较。 模型转换及实验过程： （1）小题. 由b1,b2为A1,A2的行和，可知方程组A1x=b1和A2x=b2的精确解均为n行全1的列向量。在n=5的情况下，用matlab编程（程序见四.1），构造A1,A2和b1，b2，使用高斯消去法得到的解x1,x2及其相对误差e1,e2（使用excel计算而得）为： xx1e1x2e211.0000000000000800.0000000000000801.0000000000000200.00000000000002010.999999999999463-0.0000000000005370.999999999999558-0.00000000000044211.0000000000009700.0000000000009701.0000000000020100.00000000000201010.999999999999288-0.0000000000007120.999999999996863-0.00000000000313711.0000000000001800.0000000000001801.0000000000015600.000000000001560由上表可见，当n=5时，所得的解都接近真值，误差在10^-12的量级左右。 （2）小题 分别取n=5,7,9,11,13,15，计算A1和A2的条件数c1和c2，（程序见四.2），结果如下： n579111315c13.574E+058.739E+072.274E+106.518E+122.059E+158.291E+17c24.766E+054.754E+084.932E+115.227E+146.263E+173.675E+17由上表可见，二者的条件数都比较大，可能是病态的。为证实A1和A2是否为病态，先保持b不变，对A(n,n)做扰动，得到该情况下的高斯消元解，（程序见四.3），结果如下： (为使结果清晰简洁，在此仅列出n=5,9,13的情况，n=7,11,15略去) ε=10^-10时： nxx1xx250.9999999284994020.9999999370003031.0000002510851691.0000012599939300.9999996704145830.9999943300273341.0000001916677021.0000088199574730.9999999583331440.999995590021263nxx1xx290.9999975548050940.9999829278035681.0000149244303261.0012291970545280.9999603223369260.9784890657946091.0000600154317051.1577467696177110.9999435079279020.4084498576816771.0000338876411232.2304238808731590.999987348370823-0.43549412