三个关于matlab机算的基于Matlab的线性代数机算2009.5章节.pptVIP

基于MATLAB的线性代数机算 内容概要 Matlab的初识 基本要求 线性代数机算示例 线性代数机算小结 Matlab的初识 基本要求 会矩阵的基本运算，包括：赋值，加、减、数乘，乘法，转置，幂等运算 会用命令 inv 计算矩阵的逆 会用命令 det 计算行列式 会用命令 rank 计算矩阵的秩 会用命令 rref 把矩阵变为行最简型 会用命令 rref 计算矩阵的逆 会用命令 rref 解方程组的解 会用命令 rref 找出向量组的最大无关组 会用命令 null 计算齐次线性方程组的基础解系 会用 左除‘\’ 运算计算非齐次方程组的特解 会用命令 orth 把向量组正交规范化 会用命令 eig 计算矩阵的特征值和特征向量 会用命令 eig 把二次型标准化 会用命令 eig 判断二次型的正定性 线性代数机算示例1----矩阵基本运算(I) 线性代数机算示例1----矩阵基本运算(II) 已知 线性代数机算示例2----逆矩阵的多种求法 线性代数机算示例3----求解方程组的唯一解 线性代数机算示例4----求解方程组的通解 线性代数机算示例5---求解向量组的极大无关组 已知向量组 求它的最大无关组，并用该最大无关组来线性表示其它向量 解：在MATLAB命令窗口输入： a1=[1;1;0;2;2]; % 输入5个列向量 a2=[3;4;0;8;3]; a3=[2;3;0;6;1]; a4=[9;3;2;1;2]; a5=[6;-2;2;-9;2]; A=[a1,a2,a3,a4,a5]; % 构造矩阵A [R,s]=rref(A) % R 即为行最简梯形矩阵 % s 为最大无关组的列号 线性代数机算示例6----向量组规范正交化 设向量组： 求由这三个向量生成的子空间V的一个标准正交基。 线性代数机算示例7----方阵的特征值和特征向量 lamad=eig(A) 线性代数机算示例8----二次型标准化 用正交变换法将二次型 化为标准形。 Institute of Intelligent Control and Image Engineering Institute of Intelligent Control and Image Engineering 朱娟娟 西安电子科技大学 智能控制与图像工程研究所ICIE MATLAB程序执行方式： MATLAB的启动 MATLAB的窗口 MATLAB的命令窗口 MATLAB的退出 解释执行(立即执行) rix oratory 矩阵实验室 Mat Lab 矩阵的输入 1、基本输入法： 在MATLAB命令窗口输入： A=[1,2,3;2,3,4] 2、命令输入法： A=rand(2,3) % rand 产生0~1之间按均匀概率密度分布的数 B=randn(2,3) % randn 产生以0为中心的正态分布的数 C=round(10*rand(3)) % *10表示最终的数在0~10之间 % 且 round 表示四舍五入取整 C=round(rand(3,5)*10) 计算： A+B 2. A.*B (A*B?) 3. A*C 4. det(D) 5. A’ 6. inv(D) 7. D^2 (D^50) 8. rank(D) 9. eig(D) 10. D的正定性 矩阵的各类运算 1.命令法： An1=inv(A)2.幂运算法： An2=A^-13.左除法： An3=A\eye(5) % eye(5)为5阶单位矩阵4.右除法： An4=eye(5)/A 5.初等行变换法：B=rref([A,eye(5)]); % 对矩阵[A,I] 进行初等行变换 % B为矩阵A的最简行阶梯矩阵if(rank(B(:,1:5))==5) % 判断最简行阶梯矩阵B的前5列是否为单位阵 An5=B(:,6:10) % 取出矩阵的后5列，并显示else disp(A不可逆);end A=[-7,-2,-6,4,6;1,3,-6,3,11;3,-11,