线性代数实验二行列式及方程组及求解.docVIP

实验二　　行列式与方程组的求解 一、实验目的 1掌握matlab软件求行列式的命令； 2 掌握matlab软件对矩阵进行初等行变换的命令； 3掌握matlab软件求满秩线性方程组的各种方法； 4掌握matlab软件的符号变量的应用. 二、实验指导 １.与本实验相关的Matlab 命令 clear-----------清除工作空间的各种变量 clc-------------清除工作窗口 [m,n]=size(A) ------ m,n分别为Ａ的行数与列数 det(A) ---------- 计算A的行列式 rank(A) -----------计算A的秩 U=rref(A) ----对矩阵Ａ进行初等行变换，Ｕ为Ａ的行最简型矩阵 A(:,i) =b------------将b赋值给A的第i列 A(:,1：4) ----------取 矩阵A的第1列到第4列 A(i,j)--------------引用矩阵A的第i行第j列的元素 syms x-----------定义x为符号变量 ＝＝---------关系运算符号：等于 factor(D) ------------对符号变量多项式D进行因式分解 solve(D) -----------求符号变量多项式D＝0的解 三、实验内容 例1 已知非齐次线性方程组： 解法一：初等行变换法 A=[6 2 3 4 5;2 -3 7 10 13;3 5 11 -16 21;2 -7 7 7 2;7 3 -5 3 10] %A为系数矩阵 b=[80;59;90;22;85] u=rref([A,b]) 解法二：求逆矩阵 A=[6 2 3 4 5;2 -3 7 10 13;3 5 11 -16 21;2 -7 7 7 2;7 3 -5 3 10] %A为系数矩阵 b=[80;59;90;22;85] x=inv(A)*b 解法三：矩阵左除法 A=[6 2 3 4 5;2 -3 7 10 13;3 5 11 -16 21;2 -7 7 7 2;7 3 -5 3 10] %A为系数矩阵 b=[80;59;90;22;85] x=A\b 例2 已知，用不用方法求矩阵A的逆。 分析：1）命令inv(A) 2）幂运算 A^-1 3）因为（右除） 4）因为（左除） 5) 利用伴随矩阵求逆（这个方法后面再讲） 程序如下： ％ 逆矩阵的各类求法 clear A=[-7 -2 -6 4 6;1 3 -6 3 11;3 -11 9 5 -2;-3 0 -2 9 -3;7 30 -18 11 4] An1=inv(A) % 1. 命令inv(A) An2=A^ -1 %2. 幂运算 A^-1 An3=eye(5)/A % 3. 右除 An4=A\eye(5) % 4. 左除 例三 解方程 解：％ 求解符号行列式方程 clear all syms x A=[3 2 1 1;3 2 2-x^2 1;5 1 3 2;7-x^2 1 3 2] D=det(A) %计算A的行列式 F=factor(D) %将A的行列式D因式分解 X=solve(D) %求方程D=0的解 四、实验习题 1.求解非齐次线性方程组 2. 设，求其逆. 3. 已知，求符号变量的行列式，并将结果因式分解.