MATLAB编程与应用系列-第3章 矩阵运算(2)

本系列教程来源于出版设计《基于MATLAB编程基础与典型应用书籍》,如涉及版权问题,请联系:156204968@qq.com。 出版社:人民邮电出版社, 页数:525。

本系列教程目前基于MATLABR2006a,可能对于更高级版本的功能和函数有差异,教程中如有问题,请联系:156204968@qq.com

3.1.5 矩阵函数

MATLAB提供了丰富的函数来实现对矩阵的各种运算,下面将逐步介绍。因篇幅关系,将不作数学方面的解释。表3.1所示为常用的矩阵运算函数以及函数相对应的功能描述。
表3.1 常用矩阵运算函数

函数名 功能
det(X) 计算方阵行列式
rank(X) 求矩阵的秩,得出的行列式不为零的最大方阵边长。
trace(X) 矩阵A的迹b,即A的对角线元素之和
expm(A) 使用Pade近似算法计算eA,这是一个内部函数,A为方阵
expm1(A) 使用一个M文件和内部函数相同的算法计算e^A
expm2(A) 使用泰勒级数计算e^A
expm3(A) 使用特征值和特征向量计算e^A
logm(X) 计算矩阵X的对数,它是expm(X)的反函数
funm(X, fun) 指定的函数fun计算方阵X的函数矩阵
sqrtm(X) 计算矩阵A的平方根A1/2,相当于X*X=A,求X
polyvalm(P, X) 按照矩阵运算规则计算多项式的值。其中,P为多项式系数向量,方阵X为多项式变量,返回值为多项式值
inv(X) 求矩阵的逆阵,当方阵X的det(X)不等于零,逆阵X-1才存在。X 与X-1相乘为单位矩阵。
pinv(X) 求解矩阵X的伪逆B
norm(X , ref) 求解矩阵或者向量的不同范数, ref指定了求解范数的类型
cond(X, p) 返回矩阵X的p-范数的条件数。如果p=2对应为2范数
[v,d]=eig(X) 计算矩阵特征值和特征向量。如果方程Xv=vd存在非零解,则v为特征向量,d为特征值。
diag(X) 产生X矩阵的对角阵
[l,u]=lu(X) 方阵分解为一个准下三角方阵和一个上三角方阵的乘积。l为准下三角阵,必须交换两行才能成为真的下三角阵。
[q,r]=qr(X) m×n阶矩阵X分解为一个正交方阵q和一个与X同阶的上三角矩阵r的乘积。方阵q的边长为矩阵X的n和m中较小者,且其行列式的值为1。
[u,s,v]=svd(X) m×n阶矩阵X分解为三个矩阵的乘积,其中u,v为n×n阶和m×m阶正交方阵,s为m×n阶的对角阵,对角线上的元素就是矩阵X的奇异值,其长度为n和m中的较小者。

注意:1. 当矩阵为长方阵时,方程AX=I和XA=I至少有一个无解,这时A的伪逆能在某种程度上代表矩阵的逆,若A为非奇异矩阵,则pinv(A)=inv(A)。

  1. 在上表中det(a)=0或det(a)虽不等于零但数值很小接近于零,则计算inv(a)时,其解的精度比较低, 用条件数(求条件数的函数为cond)来表示,条件数越大,解的精度越低,MATLAB会提出警告:“条件数太大,结果可能不准确”。

3.1.6 矩阵转置

矩阵转置的运算符为“'”,“A'”表示矩阵A的转置,若矩阵A的元素为实数,则与线性代数中矩阵的转置相同。若A为复数矩阵,则A转置后的元素由A对应元素的共轭复数构成。

“A.'”表示数组A的转置,如果数组A为复数数组,则仅表示复数矩阵数组的转置,而不进行矩阵或者数组的共轭操作。

【例3.17】矩阵和数组转置运算。
在命令窗口输入以下内容,生成矩阵AA和矩阵A1:
>> AA=magic(3); %生成3阶的魔方阵
>> A1=eye(3); %生成3阶的单位对角阵
由A和A1组成复数矩阵,在命令窗口输入:
>> x3=AA+A1*i % 由AA和A1组成复数矩阵
生成以下复数矩阵x3:
x3 = 8.0000 + 1.0000i 1.0000 6.0000 3.0000 5.0000 + 1.0000i 7.0000 4.0000 9.0000 2.0000 + 1.0000i
①若矩阵AA进行转置运算,可在命令窗口中输入:
>> AA' %对实数矩阵进行转置运算
矩阵AA转置后如下:
ans = 8 3 4 1 5 9 6 7 2
②若对复数矩阵x3进行转置运算,可在命令窗口中输入:
>> x3' % 对复数矩阵进行转置运算
转置后的复数矩阵如下:
ans = 8.0000 - 1.0000i 3.0000 4.0000 1.0000 5.0000 - 1.0000i 9.0000 6.0000 7.0000 2.0000 - 1.0000i
③若对复数矩阵x3进行数组转置运算,可在命令窗口中输入:
>> x3.' % 对复数矩阵进行数组转置运算
进行转置运算后如下:
ans = 8.0000 + 1.0000i 3.0000 4.0000 1.0000 5.0000 + 1.0000i 9.0000 6.0000 7.0000 2.0000 + 1.0000i

3.1.7 矩阵和数组的数学函数

MATLAB基本数学函数库中还有很多基本的数学函数运算可以对矩阵或数组进行计算,各函数的调用格式和功能描述详见表3.2。
表3.2 基本函数

函数名 含义 函数名 含义
abs 绝对值或者复数模 Rat 有理数近似
sqrt 平方根 Mod 模除求余
real 实部 Round 4舍5入到整数
imag 虚部 Fix 向最接近0取整
conj 复数共轭 Floor 向最接近-∞取整
sin 正弦 Ceil 向最接近+∞取整
cos 余弦 Sign 符号函数
tan 正切 Rem 求余数留数
asin 反正弦 Exp 自然指数
acos 反余弦 Log 自然对数
atan 反正切 log10 以10为底的对数
atan2 第四象限反正切 pow2 2的幂
sinh 双曲正弦 Bessel 贝赛尔函数
cosh 双曲余弦 Gamma 伽吗函数
tanh 双曲正切

【例3.18】数组算术运算函数的使用。
首先在MATLAB工作空间创建以下的数组B:
>> B=[5 3 4 8 5;6 8 7 6 9;2 6 8 2 1;2 7 9 3 9;4 5 6 7 8]
读者可以在MATLAB命令行窗口中输入以下命令,不同数组算术运算函数的运算结果:
>> y1=sin(B) %求取数组的正弦
>> y2=sqrt(y1) %求取数组的方根
>> y3=imag(y2) %求取数组的虚部

【例3.19】接上例,将数组y1进行圆整操作。
在MATLAB命令行窗口中输入以下命令,对数组B的正弦函数计算结果进行不同的圆整操作:
>> B1=floor(y1) % 向最接近-∞取整
>> B2=ceil(y1) % 向最接近+∞取整
>> B3=round(y1) % 4舍5入到整数
>> B4=fix(y1) % 向最接近0取整

作者:德特数据
联系方式:156204968@qq.com

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,732评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,496评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,264评论 0 338
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,807评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,806评论 5 368
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,675评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,029评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,683评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 41,704评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,666评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,773评论 1 332
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,413评论 4 321
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,016评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,978评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,204评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,083评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,503评论 2 343

推荐阅读更多精彩内容

  • 今天群里一个朋友问 怎么用indexOf判断一维数组是否存在二维数组中 因为数组是引用类型的,所以不能直接判断--...
    反者道之动001阅读 4,054评论 0 0
  • 一张机,鸳鸯戏水映春趣,亭台蝶舞百花聚,闲愁几许,相思漫忆,月坠显芳绿。 二张机,细雨伶仃乱风絮,伊人独醉思卿意,...
    夏悠羽季阅读 330评论 0 3
  • 一位干信息化的老板,与几个朋友喝酒,醉倒后不省人事。被抬回家后,老婆试着用各种办法给他醒酒,都无济于事,于是打电话...
    zhixuehuisen阅读 192评论 0 0
  • 去了这么多地方,这是第一次正儿八经的写攻略。一个原因是因为舟山的景色实在没什么好写,另一个也是舟山委实过于混乱,没...
    上官愫茗阅读 961评论 2 3
  • 第一阶段的计划(到4月份之前): 一有时间就背单词 中午回寝室看高数教程 晚上21.30回寝室看高数教程 英语每天...
    一个有点理想的文艺青年阅读 177评论 0 0