线性方程组直接实验.docVIP

3.2 线性方程组直接实验 实验目的 1．掌握解线性方程组的几种基本常用的直接法，并能比较它们各自的优缺点； 2．能熟练地对方阵A进行三角分解； 3．利用LU分解法求解线性方程组，能够体会LU分解法的优点； 4．会用追赶法求解三对角方程组，明白用追赶法解三对角方程组的优点。 实验题目 1、用列主元法求解线性方程组： 程序如下： A [0.001 2.000 3.000;-2.000 1.072 5.643;-1.000 3.172 4.623]; b [1.000 3.000 2.000] b 1 3 2 x magauss2 A,b ;x 结果为： ans -0.4646 -0.0702 0.3803 2、分别用列主元法和顺序高斯消去法求解下面的线性方程组，分析对结果的影响： 3、用追赶法分别对n 10，n 100，n 1000解方程组Ax b，其中 , . 再用LU分解法解此方程组，并对二者进行比较； 4、线性方程组直接实验-----LU分解的优点。 实验原理 1、列主元gauss 1.输入系数矩阵A，右端B，置k: 1; 2.对k 1, ……，n-1进行如下操作： （1）选列主元，确定k（r）,使a r k k 的K次导 max | a ik 的K次导|，k i n。若a r k k 0,则停止计算，否则，进行下一步； （2）若r k k,交换 A k ,b k 的第k，b k 两行； （3）消元：对i,j k+1, ……，n,计算 m ik [a il /a kk ] 的K次导,a ik 的K+1次导为0， a ik 的K+1次导 a ij 的K次导-m ik *[a kj 的K次导]， b i 的K+1次导 b i 的K次导-m ik *[b k 的K次导]； c.回代：X n [b n 的n次导]/[a nn 的n次导], 对K n-1, ……,1计算 n X k b k 的K次导 -∑[a kj 的K次导]*x j /[a kk 的K次导]J k+1 2、追赶法 1.输入对称正定矩阵A和右端向量b; 2.cholesky分解： d（1）t 11 a 11 ,l i1 a i1 /d 1 ,i 2, …,n, 对K 2, …,n计算： k-1 d k a kk - ∑t kj *l kj , j 1 k-1 t ik a ik - ∑t ij *l kj ,l ik t ik /d k ,i k+1, …,n; j-1 3.用向前消去法解下三角方程组Ly b： y 1 b 1 , k-1 对k 2, …,n计算y k b k - ∑l kj *y j ; j 1 4.解对角形方程组Dz y: 对K 1, …,n,计算z k y k /d k ; T 5.用回代法解上三角方程组L x z: X n z n , 对K n-1, …,1 计算x k z k -l jk *x j . 3、LU分解法 4、顺序高斯消去法 设 A, x, b 考察线性方程组 Ax b . 1 顺序高斯消去法的消元过程: 对方程组的增广矩阵 作初等行变换如下: 第一步, 记, , , . 设 , 对计算 , 第二步, 设, 对计算 , , 如此继续, 第步结束, 得矩阵 n 然后进行回代, 即可解出: , , . 实验内容 1、用列主元法求解线性方程组 1.1实验步骤： 新建一个M文件，编写程序，运行程序，记录结果 1.2实验程序： clc A [0.001 2.000 3.000;-2.000 1.072 5.643;-1.000 3.712 4.623] b [1.000 3.000 2.000] x magauss2 A,b 1.3实验设备： matlab软件。 2、分别用列主元法和顺序高斯消去法求解下面的线性方程组 2.1实验步骤： 新建一个M文件，编写程序，运行程序，记录结果 2.2实验程序： clear clc A [0.3*10^ -16 59.14 3 1;1 2 1 1;11.2 9 5 2;5.291 -6.13 -1 2] b [51.97 2 1 46.78]’ x1 magauss2 A,b x2 magauss A,b 2.3实验设备： matlab软件。 3、用追赶法分别对n 10，n 100，n 1000解方程组Ax b 3.1实验步骤： （1）新建一个M文件，编写程序（如下程序1），运行程序，记录结果。 （2）把程序中的n 10改为n 100, 运行程序，记录结果。 （3）把程序中的n 10改为n 1000, 运行程序，记录结果。 （4）再用LU分解法解此题，新建一个M文件