解方程组基本思想.docVIP

. .. 四：基本方法基本思路将在解题的过程中得到体现。1．（求线性方程组的唯一解或特解），这类问题的求法分为两类：一类主要用于解低阶稠密矩阵 —— 直接法；一类是解大型稀疏矩阵 —— 迭代法。 1.1利用矩阵除法求线性方程组的特解（或一个解）方程：AX=b，解法：X=A\b，(注意此处’\’不是’/’)例1-1??? 求方程组 的解。解:??? A = ; = ；b=(1,0,0,0,1)’由于rank(A)=5,rank( )=5??? %求秩，此为R（A）=R（ ）=n的情形，有唯一解。??????X= A\b????? %求解 X =(2.2662, -1.7218, 1.0571,-0.5940, 0.3188)’ 或用函数rref求解，sv=rref(A:b);所得sv的最后一列即为所要求的解。1．2 利用矩阵的LU、QR和cholesky分解求方程组的解 这三种分解，在求解大型方程组时很有用。其优点是运算速度快、可以节省磁盘空间、节省内存。I) LU分解又称Gauss消去分解，可把任意方阵分解为下三角矩阵的基本变换形式（行交换）和上三角矩阵的乘积。即A=LU，L为下三角阵，U为上三角阵。则：A*X=b????????? 变成L*U*X=b所以X=U\(L\b)???? 这样可以大大提高运算速度。命令??? [L，U]=lu (A)