数学建模线性代数matlab代码,数学建模matlab简介第三部分线性代数常用算法

数学建模matlab简介第三部分线性代数常用算法

MATLAB 简介数学建模第三部分 线性代数常用算法一、矩阵的创建1、直接输入法采用赋值命令来完成，整个矩阵必须以方括号[ ] 为其首尾，行与行之间用 “；”或回车分隔，元素间用“， ”或空格分隔。矩阵元素可以是包含已定义变量的任何表达式。例 1、创建矩阵 Ca=2.7358; b=33/79;C=[1,2*a,b*sqrt(a); sin(pi/6),a+5*b,3.5]例 2、分行输入矩阵 DD=[1，2，34 5 67 8 9注：利用下标可以获得矩阵元素或修改矩阵。例 3、D(3,2) %获得矩阵元素D(3,2)=2 %修改矩阵元素2、几种常见矩阵的生成命令eye(n) 产生 n 阶单位矩阵eye(size(A)) 产生与矩阵 A 同维的单位矩阵ones(n) 产生 n 阶全 1 矩阵ones(size(A)) 产生与矩阵 A 同维的全 1 矩阵zeros(n) 产生 n 阶零矩阵zeros (size(A)) 产生与矩阵 A 同维的零矩阵zeros(m,n) 产生 m 行 n 列的零矩阵company(A) 产生 A 的伴随矩阵二、矩阵运算命令A’ 矩阵转置A+B 矩阵相加A-B 矩阵相减A*B 矩阵相乘k*A 数乘矩阵inv(A) 矩阵求逆A^n 矩阵求幂det(A) 矩阵的行列式rank(A) 矩阵的秩eig(A) 矩阵的特征值MATLAB 简介数学建模例 4、计算 AB’A=[11 12 13 14 15; 16 17 18 19 20]; B=ones(size(A));A*B’例 5、计算 13,AA=[1 2 3;4 5 6;7 8 9]; A^3inv(A)三、解线性方程组左除法 A\B 求解矩阵方程 AX=B右除法 B/A 求解矩阵方程 XA=Brref(A) 求矩阵的行最简形null(A) 求以 A 为系数矩阵的齐次线性方程组的基础解系例 6、解方程组 4231yxA=[1 2;3 -2]; B=[1;4]; x=A\B例 7、解方程组 4231zyxA=[1 2 1;3 -2 1]; B=[1;4]; x=A\B例 8、解方程组 241yxA=[1 2;-2 -4]; B=[1;-2]; x=A\BA=[1 2;-2 -4;0 0]; B=[1;-2;0]; x=A\B例 9、求线性方程组的通解 12431xx方法一：A=[1 -1 1 -1；-1 1 1 -1；2 -2 -1 1]；b=[1;1;-1];[rank(A),rank([A,b])]rref([A,b])方法二：x0=A\b, x=null(A)MATLAB 简介数学建模四、特征值与特征向量[V，D]=eig(A) 得到方阵 A 的特征值与特征向量(其中 D 为特征值构成的对角阵，V 的每一列为相应的已单位化正交化的特征向量) 。例 10、求 A 的特征值与特征向量A=[1 2 2 ; 2 4 4 ; 2 4 4]; [V,D]=eig(A)