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

实验5 线性代数方程组的数值解法分1 黄浩 2011011743实验目的学会用MATLAB软件数值求解线性代数方程组，对迭代法的收敛性和解的稳定性作初步分析； 通过实例学习用线性代数方程组解决简化的实际问题。实验内容《数学实验》第二版（问题1）问题叙述：通过求解线性方程组，理解条件数的意义和方程组性态对解的影响，其中是n阶范德蒙矩阵，即是n阶希尔伯特矩阵，b1，b2分别是的行和。编程构造（可直接用命令产生）和b1，b2；你能预先知道方程组和的解吗？令n=5，用左除命令求解（用预先知道的解可验证程序）。令n=5,7,9，…，计算和的条件数。为观察他们是否病态，做以下试验：b1，b2不变，和的元素，分别加扰动后求解；和不变，b1，b2的分量b1(n),b2(n)分别加扰动后求解。分析A与b的微小扰动对解的影响。取10^-10,10^-8,10^-6。经扰动得到的解记做，计算误差，与用条件数估计的误差相比较。模型转换及实验过程：（1）小题.由b1,b2为,的行和，可知方程组和的精确解均为n行全1的列向量。在n=5的情况下，用matlab编程（程序见四.1），构造,和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，计算和的条件数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由上表可见，二者的条件数都比较大，可能是病态的。为证实和是否为病态，先保持b不变，对做扰动，得到该情况下的高斯消元解，（程序见四.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.4354941236522751.0000026878626961.8788737437097660.9999997511934060.780281609027228nxx1xx2131.0007784033869290.9999998767583270.9937743987274311.0000261185116641.0226872410459610.9987794790026500.9501840256039321.0233044427317511.0734106035353810.7678410435503330.9235081059886532.3669384594613431.057790678168565-4.0983444092691000.96809962810653313.49836909493934