实验三
第三节  矩阵基本函数运算与矩阵元素的提取(第二章)
一、 矩阵基本函数运算
此运算是矩阵运算中最实用的部分,其基本命令如下:
命令集 矩阵的大小、行列式、逆、特征值、秩、迹、范数
size(A)  给出包含A的维数的一个行向量.在这个返回向量中的第一个元素是行数,随后是列数 .
[ mn ]=size(A)  给出A的维数、m行数和n列数,即两个标量.
length(x) 给出一个向量的长度,即x分量个数.
sum(A)  A是矩阵,给出一个行向量,其每个分量表示A相应的列和;若A是向量,给出此向量的分量和.
det(A)  求矩阵A的行列式.   
eig(A)  求包含矩阵A的特征值的向量.
[X,D]=eig(A)  求包含矩阵A的特征值对应的对角阵D和以相应特征向量为列的矩阵.
inv(A)A ^ (-1) 求矩阵A的逆矩阵.
rank(A)  求矩阵A的秩.
trace(A)  求矩阵A的迹(对角线元素之和).
norm(A1)  矩阵A1—范数或列和范数,定义如下.
norm(A2)  矩阵A2—范数.
norm(Ainf)  矩阵A的∞范数.
norm(x1)  向量x1—范数或列和范数,定义如下.
norm(x2)  向量linspace函数调用的格式为x2—范数.
norm(xinf)  向量x的∞范数.
范数定义如下:
,则相应范数定义如下
二、矩阵元素的提取
MATLAB中还有利用已存在的矩阵建立新矩阵的命令.以下假设矩阵 Am×n的矩阵,x是个有n个元素的向量.
1. 对角阵与三角阵的生成
命令集10 
diag(A)  生成一个由矩阵A主对角线元素组成的列向量.主对角线总是从矩阵左上角开始.对于方阵来说它结束于矩阵的右下角.
diag(x)  生成一个n维的方阵,它的主对角线元素值取自向量 x,其余元素的值都为0.
diag(A , k) 生成一个由矩阵Ak条对角线的元素组成的列向量. k= 0为主对角线;k< 0为下第k对角线;k> 0为上第k对角线.
diag(x , k)  生成一个(n+ a b s (k) )×(n+ a b s (k) )维的矩阵,该矩阵的第k条对角线元素取自向量x,其余元素都为零.关于参数k可参考上个命令.
triu(A)    生成一个和A大小相同的上三角矩阵.该矩阵的主对角线及以上元素取自A中相应元素,其余元素都为零.
triu(A , k)  生成一个和A大小相同的上三角矩阵.该矩阵的第k条对角线及以上元素取自A
相应元素,其余元素都为零.
命令t r i u ( A , 0 )等同于命令t r i u ( A ).
tril(A)  生成一个和A大小相同的下三角矩阵.该矩阵的主对角线及以下元素取自A中相应元素,其余元素都为零.
tril(A , k) 生成一个和A大小相同的下三角矩阵.该矩阵的第 k条对角线及以下元素取自A中相应元素,负数k表示主对角线下的对角线.其余元素都为零.命令t r i l ( A , 0 )等同于命令t r i l ( A ).
2. 向量和子矩阵的生成
MATLAB中可以使用冒号‘:’来代表一系列数值.有时也使用它来定义一个子矩阵.
命令集11
i : k  创建从i开始、步长为1、到k结束的数字序列,即i ,i+1, i+2, . . . , k .如果i>kMATLAB则返回一个空矩阵,也就是 [ ].数字ik不必是整数,该序列的最后一个数是小于或等于k.
i : j : k  创建从i开始、步长为j、到k结束的数字序列,即i, i+j, i+ 2j, . . ., k .对于j= 0,则返回一个空矩阵.数字ijk不必是整数,该序列的最后一个数是小于或等于k.
linspace(a , b)  在区间[a, b]上创建一个有1 0 0个元素的向量,这1 0 0个数把整个区间线性分隔.其中a是第一个元素,b是最后一个.
linspace(a, b, n)  在区间[a, b]上创建一个有n个元素的向量.这个命令和冒号表示形式相近,但是它直接定义了数据的个数,其步长为(b-a/(n-1) .
命令集12  定义子阵
A ( i , j )    返回矩阵A中第ij元素的值.
A ( : , j )    返回矩阵A中第j列列向量.
A ( i , : )      返回矩阵A中第i行行向量.
A ( : , j : k )  返回由矩阵A中的第j列,第j+ 1列,直到第k列列向量组成的子阵.
A ( i : k , : )  返回由矩阵A中的第i行,第i+ 1行,直到第k行行向量组成的子阵.
A ( i : k , j : l )  返回由二维矩阵A中的第i行到第k行行向量和第j列到第l 列列向量组成的子阵.
A ( : )        将矩阵A中的每列合并成一个长的列向量.
A(j:k)        返回一个行向量,其中的元素为A ( : )中的从第j个元素到第k个元素.
A([j1 j2 . . . ] )  返回一个行向量,其中的元素为A ( : )中的第j1j2 元素.
A(:,[j1 j2 . . .] )  返回矩阵A的第j1列、第j2列等的列向量.
A([i1 i2. . .  ] : , )  返回矩阵A的第i1行、第i2行等的行向量.
A([i1 i2. . . ] , [j1 j2. . .  ] )  返回矩阵第i1行、第i2行等和第j1列、第j2列等的元素.
二、矩阵元素的增减
MATLAB中可以通过增加元素、行和列将一个矩阵或者向量进行扩展.由于MATLAB可以自动地改变矩阵的大小,所以使用已存在的矩阵的一部分来创建一个新矩阵是很容易的这
在许多应用中都很有用.从已存在的矩阵中建立一个矩阵就和定义一个新矩阵一样.元素用空格或逗号分隔,行用分号或回车分隔.
1 假设下列矩阵已经定义为:
(a) 有几种方式可以将向量x扩展成1×4.假设想要的新向量是:
xnew=(9 10 0 5)
下列的三种方法都可以给出想要的结果:
  xnew=x; xnew(3)=0xnew(4)=5;      xnew=[x 0 5];

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。