线性代数方程组的值解法.docxVIP

线性代数方程组的数值解法 【实验目的】 1. 学会用MATLAB 软件数值求解线性代数方程组，对迭代法的收敛性和解的稳定性作初步分析； 2. 通过实例学习用线性代数方程组解决简化的实际问题。 【实验内容】 【题目1】 通过求解线性方程组和，理解条件数的意义和方程组的性态对解的影响。其中是n阶范德蒙矩阵，即 =，， 是n阶希尔伯特矩阵，，分别是，的行和。 （1）编程构造（可直接用命令产生）和，；你能预先知道方程组和的解吗？令n=5，用左除命令求解（用预先知道的解可检验程序）。 （2）令n=5,7,9,…，计算，的条件数。为观察它们是否病态，做以下试验：，不变，和的元素，分别加扰动后求解；和 不变，，的分量，分别加扰动求解。分析A和b的微小扰动对解的影响。取，，。 （3）经扰动得到的解记做，计算误差，与用条件数估计的误差相比较。 1.1构造，和， 首先令n=5，构造出，和，。首先运行以下程序，输出。 function A1=fdm(n) k=[]; n=5 o=ones(n,1); for i=1:n a(i)=1+0.1*(i-1); end for i=1:n k=[k,a(i)]; end k=k; for i=1:n-1 o=[o,k.*o(:,i)]; end A1=o; 直接由命令产生A2： n