数学实验——线性代数方程组的数值解.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.00