MATLAB仿真说明.docVIP

MATLAB仿真指导说明 一、矩阵的建立和基本运算实验指导 1. 矩阵的运算 运算符：+（加）、-（减）、*（乘）、/（右除）、\（左除）、^（乘方）、’（转置）等； 常用函数：det（行列式）、inv（逆矩阵）、rank（秩）、eig（特征值、特征向量）、rref（化矩阵为行最简形） 例5： A=[2 0 –1;1 3 2]; B=[1 7 –1;4 2 3;2 0 1]; M = A*B % 矩阵A与B按矩阵运算相乘 det_B = det(B) % 矩阵A的行列式 rank_A = rank（A） % 矩阵A的秩 inv_B = inv（B） % 矩阵B的逆矩阵 [V,D] = eig(B) %矩阵B的特征值矩阵V与特征向量构成的矩阵D X = A/B % A/B = A*B-1，即XB=A，求X Y = B\A % B\A = B-1*A，即BY=A，求Y 2. 线性方程的组的求解 我们将线性方程的求解分为两类：一类是方程组求唯一解或求特解，另一类是方程组求无穷解即通解。可以通过系数矩阵的秩来判断： 若系数矩阵的秩r=n（n为方程组中未知变量的个数），则有唯一解 若系数矩阵的秩rn，则可能有无穷解。 线性方程组的无穷解=对应齐次方程组的通解+非齐次方程组的一个特解；其特解的求法属于解的第一类问题，通解部分属第二类问题。 利用矩阵除法求线性方程组的特解（或一个解） 方程：AX=b 解法：X=A\b 例:求方程组的解 解：在Matlab编辑器中建立M文件：LX0716.m A=[5 6 0 0 0 1 5 6 0 0 0 1 5 6 0 0 0 1 5 6 0 0 0 1 5]; B=[1 0 0 0 1]; R_A=rank(A) %求秩 X=A\B %求解 运行后结果如下 R_A = 5 X = 2.2662 -1.7218 1.0571 -0.5940 0.3188 这就是方程组的解。 例:求方程组 的一个特解 解：在Matlab编辑器中建立M文件：LX0717.m A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; B=[1 4 0]; X=A\B X = 0 0 -0.5333 0.6000 3. 求线性齐次方程组的通解 在Matlab中，函数null用来求解零空间，即满足A·X=0的解空间，实际上是求出解空间的一组基（基础解系）。 格式：z = null % z的列向量为方程组的正交规范基，满足 % z的列向量是方程AX=0的有理基 例:求解方程组的通解： 解：在Matlab编辑器中建立M文件：LX0719.m A=[1 2 2 1;2 1 -2 -2;1 -1 -4 -3]; format rat %指定有理式格式输出 B=null(A,r) %求解空间的有理基 运行后显示结果如下： B = 2 5/3 -2 -4/3 1 0 0 1 写出通解： syms k1 k2 X=k1*B(:,1)+k2*B(:,2) %写出方程组的通解 pretty(X) %让通解表达式更加精美 运行后结果如下： X = [ 2*k1+5/3*k2] [ -2*k1-4/3*k2] [ k1] [ k2] %下面是其简化形式 [2k1 + 5/3k2 ] [ ] [-2k1 - 4/3k2] [ ] [ k1 ] [ ] [ k2 ] 4. 求非齐次线性方程组的通解 非齐次线性方程组需要先判断方程组是否有解，若有解，再去求通解。因此， 步骤为：第一步：判断AX=b是否有解，若有解则进行第二步 第二步：求AX=b的一个特解 第三步：求AX=0的通解 第四步：AX=b的通解= AX=0的通解+AX=b的一个特解。 例:求解方程组 解：在Matlab编辑器中建立M文件：LX0720.m A=[1 -2 3 -1;3 -1 5 -3;2 1 2 -2];b=[1 2 3];B=[A b];n=4; R_A=rank(A),R_B=rank(B),format rat ifR_A==R_BR_A==n %判断有唯一解 X=A\