线性代数基础性实验.docVIP

线性代数基础性实验 一、实验目的 1.掌握MATLAB软件中矩阵和行列式的如下输入命令和基本运算命令： 矩阵的输入命令和产生特殊矩阵的命令； 矩阵的加减、数乘、乘法、转置运算和MATLAB中矩阵的除法运算； 求方阵的行列式的命令和求可逆矩阵的逆矩阵的命令； 利用矩阵的运算求解矩阵方程的方法。 构造已知矩阵的行（列）向量组和子矩阵的方法； 2.掌握MATLAB软件中关于向量组和线性方程组的如下命令和方法： 将矩阵化为行最简形矩阵的命令和求矩阵的秩的命令； 求向量组的最大无关组和将不属于最大无关组的向量用最大无关组线性表示的方法； 求齐次线性方程组的基础解系和通解的方法； 求非齐次线性方程组的特解和通解的方法。 3.掌握MATLAB软件中关于矩阵对角化和二次型的如下命令和方法： 求方阵的特征值和特征向量的命令； 将可对角化矩阵对角化的方法； 利用正交矩阵将实对称矩阵对角化的方法； 利用正交变换化二次型为标准形的方法。 二、相关函数（命令）及简介 详见教育在线中的辅助讲义——《MATLAB线性代数实验指导书》（重庆交通大学理学院线性代数课程组编），以下简称辅助讲义。 三、实验内容 要求： 1）编写一个M文件完成以下所有的实验题目，并通过教育在线中的提交作业系统进行提交。 2）M文件名为“LA学号”，如学号为100001的同学，M文件名应为“LA100001”。 没有按要求提交作业的同学本次实验成绩记为零分。 1.计算并化简符号矩阵或行列式。（参见辅助讲义P12，例3-11，） 计算行列式的值，并利用pretty命令进行化简。 解： %1. syms x y A=[x y x+y;y x+y x;x+y x y]; d1=det(A) d2=simple(d1) pretty(d2) 则显示结果如下： d1 = - 2*x^3 - 2*y^3 d2 = - 2*x^3 - 2*y^3 pretty(d2) 3 3 - 2 x - 2 y 例：3-11计算行列式 的值。 在Matlab编辑器中建立M文件：LX0711.m syms a b c d A=[1 1 1 1;a b c d;a^2 b^2 c^2 d^2;a^4 b^4 c^4 d^4]; d1=det(A) d2=simple(d1) %化简表达式 d1 pretty(d2) %让表达式 d2 符合人们的书写习惯 则显示结果如下： d1 = b*c^2*d^4-b*d^2*c^4-b^2*c*d^4+b^2*d*c^4+b^4*c*d^2-b^4*d*c^2-a*c^2*d^4+a*d^^4+a*b^2*d^4-a*b^2*c^4-a*b^4*d^2+a*b^4*c^2+a^2*c*d^4-a^2*d*c^4-a^2*b*d^4+a^2*b*4+a^2*b^4*d-a^2*b^4*c-a^4*c*d^2+a^4*d*c^2+a^4*b*d^2-a^4*b*c^2-a^4*b^2*d+a^4*b^2 d2 = (-d+c)*(b-d)*(b-c)*(-d+a)*(a-c)*(a-b)*(a+c+d+b) (-d+c)(b-d)(b-c)(-d+a)(a-c)(a-b) 2. 在MATLAB环境中生成特殊的矩阵。（利用命令linspace，或符号“：”） 已知：命令A = -10:10或A = linspace(-10, 10, 21)可以形成一个的行向量： 。 试借助于上面的命令采用2种方法构造出如下矩阵M： 。 3.用2种方法求解矩阵方程。（参见辅助讲义P9，例3-5和P12，例3-12） 设，且满足，试分别用上述2种方法求矩阵。 法一：利用逆矩阵； (A-2*E)*X=A X=A*(A-2*E)-1 解： %3（1）. A=[1 2 3; 2 2 1; 3 4 3] E=eye(3) B=A-2*E inv(B) format rat X=A*inv(B) 例：3-5 求的逆矩阵 A=[1 2 3; 2 2 1; 3 4 3]; inv(A)或 A^(-1) 则结果显示为 ans = 1.0000 3.0000 -2.0000 -1.5000 -3.0000 2.5000 1.0000 1.0000 -1.0000 法二：利用初等变换； 解： %3.（2） A=[1 2 3; 2 2 1; 3 4 3] E=eye(3) B=A-2*E C=[-1 2 3 1 0 0 ;2 0 1 0 1 0 ;3 4 1 0 0 1] D=rref(C) X=A*D(:, 4:6) 方法二：由增广矩阵进行初等行变换 在 Matlab