ch6 实验2MATLAB在线性代数中的应用.ppt

MATLAB在线性代数的应用 matlab在线性代数中的运用 （1）??? matlab的矩阵初等运算 （2）??? 向量组的线性相关性分析 （3）??? 行列式与线性方程组的解 （4）??? 特征值与特征向量 （5）??? 矩阵的相似对角化 （6）??? 化二次型为标准型 目标要求 会给矩阵赋值 会进行矩阵的基本运算，包括：加、减、数乘、乘法、转置、幂等运算 会用命令inv计算矩阵的逆 会用命令rank计算矩阵的秩 会用命令rref把矩阵变为行最简型 会用命令rref计算矩阵的逆 会用命令rref接方程组的解 会用命令rref找出向量组的最大无关组 会用命令null计算齐次线性方程组的基础解系 会用左除运算计算非齐次方程组的特解 会用命令orth把向量组正交规范化 会用命令eig计算矩阵的 特征值和特征向量 会用命令eig把二次型标准化 会用命令eig判断二次型的正定性 矩阵的赋值 赋值语句一般形式 变量=表达式（或数） 如： 》A=[1 2 3;4 5 6 ;7 8 9] 显示A=1 2 3 4 5 6 7 8 9 》X=[-1.2 sqrt(3) (1+2+3)/5*4] 显示X=-1.2000 1.7321 4.8000 规则：矩阵元素放在括号中，元素之间可以空格或逗号分隔，不同行以分号分隔，语句结尾用回车将显示结果 基本赋值矩阵 ones(m,n),zeros(m,n),eye(n),rand(m,n) 如：》A=ones(3,2) 》B=zeros(2,3) 显示A=1 1 显示B=0 0 0 1 1 0 0 0 1 1 C=eye(3) rand(m,n)是产生随机矩阵， 显示C=1 0 0 0 1 0 0 0 1 矩阵初等运算 A.^3 %“.*”是矩阵对应元素的乘积% ans = 1728 39304 -64 39304 343 658503 27 274625 343 A/B ans = 16.4000 -13.6000 7.6000 35.8000 -76.2000 50.2000 67.0000 -134.0000 68.0000 B\A ans = 109.4000 -131.2000 322.8000 -53.0000 85.0000 -171.0000 -61.6000 89.8000 -186.2000 大家思考一下左除和右除有区别吗??? 第二节 向量组的线性相关性分析 向量组的最大无关组 解：在mtlab中新建一个m文件程序为： A=[2 -1 3 5;4 -3 1 3;3 -2 3 4; 4 -1 15 17;7 -6 -7 0 ]; A=A; format rat; rref(A) 运行结果为ans = 1 0 0 2 1 0 1 0 -3 5 0 0 1 4 -5 0 0 0 0 0 相关命令：rref(A):对矩阵A进行初等行变换 det(A):计算矩阵A的行列式 inv(A):求矩阵A的逆矩阵 思考：如何用求逆矩阵或初等变换法求方程的解？ 第三节行列式与线性方程组的解 解：在matlab中输入 A=[1 –5 3 –3;2 0 1 –1;3 1 –1 2;4 1 3 -1]; D=det(A) ans= -55 B = 9 -4 1 X = 9*k -4*k k 用Matlab的eig命令可以求出矩阵的特征值与特征向量。 命令eig的使用方法有两种：只求矩阵A的特征值用命令eig(A)，求特征值和特征向量用命令[v,d]=eig(A)。 第四节 特征值和特征向量 解：在matlab中的程序为： clear A=[-29 6 18;20 5 12;-8 8 5] [V,D]=eig(A) V = 698/979 395/1409 989/3619 -859/1412 -2430/3089