![MATLAB Simulink 2020系统仿真从入门到精通](https://wfqqreader-1252317822.image.myqcloud.com/cover/221/44175221/b_44175221.jpg)
2.5 MATLAB下矩阵的运算
2.5.1 矩阵的代数运算
矩阵的代数运算包括加、减、乘、数乘、点乘、乘方、左除、右除等。其中加、减、乘与大家所学的线性代数中的定义是一样的,相应的运算符为“+”“-”“*”。
对于上述的运算,需要注意的是,矩阵的加、减、乘运算对维数要求与线性代数中的要求一致。
1.矩阵的加减运算
设A=(aij),B=(bij)都是m× n矩阵,矩阵A与B的和记成A+B,规定为:
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/39_02.jpg?sign=1739105884-IOA3kqfisuNwcaC6YTLDJe1NTdB76lKe-0-7a798fceeaf059d2f27fd173839778ba)
1)交换律:A+B=B+A。
2)结合律:(A+B)+C=A+(B+C)。
例2-16:验证加法法则。本实例验证矩阵加法的交换律与结合律。
解:MATLAB程序如下。
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/39_03.jpg?sign=1739105884-8YMd7V6PxtSEUsZZleOjDLCnjFpRrnWR-0-3e84c41b65e1bbf5aa6c07803479b645)
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/39_04.jpg?sign=1739105884-XfJnhGdCQhlyzEaXyYOUeXSzEVyTg84U-0-2f462c95e1e365ad2c16ab47d643174e)
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/40_01.jpg?sign=1739105884-dg8RG1PpmUCcqbiKLTxwWx2eAxiqi6f9-0-c3b3b2990d2a42d21fbe8f4bcff06eb1)
减法运算法则为:A-B=A+(-B)。
例2-17:矩阵求差。本实例求矩阵的减法运算。
解:MATLAB程序如下。
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/40_02.jpg?sign=1739105884-gLHHKhpWbSpUZ5Tx8vMGYscyaht3fhQj-0-876e30389999d13573376c207191bcfa)
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/40_03.jpg?sign=1739105884-rrYZAUxmivqU4SKUuxZmuIC6duImQ0Bk-0-8b02e13d08c0008cd23454c347a79ca9)
2.矩阵的乘法运算
(1)数乘运算
数λ与矩阵A=(aii)m×n的乘积记成λA或Aλ,规定为:
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/40_04.jpg?sign=1739105884-Pfa8JNr5LEhdKr0xDEzwH5wfttVAodTh-0-478b611686171df0c1c6a76af82e6bba)
同时,矩阵还满足下面的规律:
λ(μA)=(λμ)A
(λ+μ)A=λA+μA
λ(A+B)=λA+λB
其中,λ,μ为数,A,B为矩阵。
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/41_01.jpg?sign=1739105884-UpLO3a2rcvVQ9l3RNpnfInvSMUa359U9-0-62799eafa41c04162c7cfc0cd8441093)
(2)乘运算
若3个矩阵有相乘关系,设A=(aij)是一个m× s矩阵,B=(bij)是一个s× n矩阵,规定A与B的积为一个m× n矩阵C=(cij):
cij=ai1b1j+ai2b2j+…+aisbsji=1,2,…,m; j=1,2,…,n
即C=A*B,需要满足以下3种条件。
◆ 矩阵A的列数与矩阵B的行数相同。
◆ 矩阵C的行数等于矩阵A的行数,矩阵C的列数等于矩阵B的列数。
◆ 矩阵C的第m行n列元素值等于矩阵 A的 m行元素与矩阵 B的 n列元素对应值积的和。
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/41_02.jpg?sign=1739105884-BhBWAXgZOyW9KjmIZZ3u1qWeBG8t2Qix-0-4d7f2ef92340bbc5d678fbc7a5e1eda1)
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/41_03.jpg?sign=1739105884-wQXrzEZ7fwwXOKrRnQObwuEDffNlvF26-0-ca8d17fa385acdfbb346c3c266643c61)
注意:
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/41_05.jpg?sign=1739105884-Z8bGrHcVRmnD8ivutbqvGCmacgPy7PSU-0-9fd83fb6b6e1fa6fa285c64436512e42)
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/42_01.jpg?sign=1739105884-fGM0BUPSMS5RMUuL3O9wIQLO4s9rG7ZP-0-9744c6d151cc2de4a7fb9aa1da1b3366)
若矩阵A、B满足AB=0,未必有A=0或B=0的结论。
3.点乘运算
点乘运算是指将两矩阵中相同位置的元素进行相乘运算,将积保存在原位置组成新矩阵。
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/42_02.jpg?sign=1739105884-VsI8ylJuGIpp8EkrJ0z20BijByhxu894-0-6d923d1b0d8073b2460c2a7e0e84fc39)
例2-18:矩阵点乘运算。
解:MATLAB程序如下。
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/42_03.jpg?sign=1739105884-sJuUIUvSPrqr2VzVNi9DrMeKaRwV0Wz0-0-9dbe4b3e4d14c5a00bdd8567e9fe4826)
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/42_04.jpg?sign=1739105884-EG3u6asZPxGJrG8zuOSCNe1YiJjPvEoE-0-ef007b2446ea288f6ff0f8f86bc3492f)
4.矩阵的除法运算
由于矩阵的特殊性计算左除A\ B时,A的行数要与B的行数一致,计算右除A/B时,A的列数要与B的列数一致。
(1)左除运算
A*B通常不等于B*A,除法也一样。因此除法要区分左除和右除。
线性方程组D*X=B,如果 D非奇异,即它的逆矩阵inv(D)存在,则其解用MTLAB表为:
X=inv(D)*B=D\B
符号“\”称为左除,即分母放在左边。
左除的条件:B的行数等于D的阶数(D的行数和列数相同,简称阶数)。
例2-19:求解矩阵左除。
解:MATLAB程序如下。
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/43_01.jpg?sign=1739105884-sWsOjxPtx1sSNxtQLkV7MLNiSvmLGKZl-0-318a3b0232e57b6709d08e9f08bfc18b)
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/43_02.jpg?sign=1739105884-CopGTBSjvnHeEHmuFpnWSvdEARnYfgZr-0-09158cfb46c5b2142aa18a2ef75b4839)
(2)右除运算
若方程组表示为X*D1=B1,D1非奇异,即它的逆阵inv(D1)存在,则其解为:
X=B1*inv(D1)=B1/D1
符号“/”称为右除。
右除的条件:B1的列数等于D的阶数(D的行数和列数相同,简称阶数)。
例2-20:验证矩阵的右除。
解:MATLAB程序如下。
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/43_03.jpg?sign=1739105884-M79bNCF1pMkW3fWQISTSI5FiU1Gk5eG0-0-520001cbfa30fa1d398b83bb63f6436f)
![](https://epubservercos.yuewen.com/9B02D1/23446275609890206/epubprivate/OEBPS/Images/43_04.jpg?sign=1739105884-Dz1h94RrgZwRuqBEiiQaaVpFwTsPSDH8-0-ebd0d9dec7dcc160544ef57ad79266f3)