实验数学三：MATLAB在线性代数中应用.ppt

* MATLAB在线性代数中的应用 3.1实验目的 矩阵是人们用数学方法解决实际问题的重要工具，而MATLAB具有强大的矩阵运算功能，本实验的目的是学会用MATLAB软件进行线性代数中一些运算，包括矩阵的基本运算、计算矩阵行列式、线性方程组求解、矩阵的特征值和特征向量、矩阵的分解和化二次型为标准型等。 3.2实验内容 1、矩阵的基本运算 MATLAB中矩阵基本运算的指令和意义如下： A’:矩阵的转置； A+B:矩阵相加； A-B:矩阵相减； A*B:矩阵相乘； s*B:矩阵的数乘。 例1 解 输入命令 ： A=[1 2 3;0 -4 3;6 4 9]; B=[9 4 2;5 8 -4;0 3 1]; A ans = 1 0 6 2 -4 4 3 3 9 A+B ans = 10 6 5 5 4 -1 6 7 10 2*A-3*B ans = -25 -8 0 -15 -32 18 12 -1 15 例2 解 输入命令 ： A=[1 2 3;0 -4 3;6 4 9]; B=[1 2;6 -9; 2 0]; A*B ans = 19 -16 -18 36 48 -24 2、几个特殊矩阵的建立 zeros(m,n):产生一个m行、n 列的零矩阵； ones(m,n): 产生一个m行、n 列的元素全为1的矩阵； eye(m,n): 产生一个m行、n 列的单位矩阵； 例3 产生3行4列的零矩阵、2行3列的元素全为1的矩阵及4行5列的单位矩阵。 解：输入命令 ： zeros(3,4) ans = 0 0 0 0 0 0 0 0 0 0 0 0 ones(2,3) ans = 1 1 1 1 1 1 eye(4,5) ans = 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 3、矩阵的块操作 A(i,:):取出A的第i行的所有元素； A(：,j):取出A的第j列的所有元素； A([i,j],:):取出A的第i,j行的所有元素； A(i:j,k:s):取出A的第i 行到第j行与第k列到第s列的交叉元素； A(i,:)=a:将A的第i行的所有元素用a替代； A(find(A==a))=b:将A的等于a的所有元素用b替代； A(i,:)=[]:删除A的第i行. [A B]或[A;B]：将矩阵A和B拼成新的矩阵。 例4已知 对矩阵A进行一些块操作。 解：对矩阵A进行一些块操作命令如下 A=[1 2 1 0;2 3 6 9;3 4 5 8;4 5 6 7]; (1)写出A的第2行元素： A(2,:) ans = 2 3 6 9 (2)写出A的第3列元素： A(:,3) ans = 1 6 5 6 (3)写出A的第2和3行元素： A([2,3],:) ans = 2 3 6 9 3 4 5 8 （4）写出A的第2 行到第3行与第3列到第4列的交叉元素： A(2:3,3:4) ans = 6 9 5 8 （5）将A的第2行的所有元素用9替代； A(2,:)=9 A = 1 2 1 0 9 9 9 9 3 4 5 8 4 5 6 7 (6)A(find(A==3))=0:将A的等于3的所有元素用0替代； A(find(A==3))=0 A = 1 2 1 0 2 0 6 9 0 4 5 8 4 5 6 7 (7)删除A的