第3講矩陣的運算1matlab運算2矩陣分析3字符串

1、第3講 矩陣的運算3.1 MATLAB運算3.2 矩陣分析3.3 字符串,3.1 MATLAB運算,3.1.1算術運算,1. 基本算術運算 MATLAB的基本算術運算有：＋(加)、－(減)、*(乘)、/(右除)、\(左除)、^(乘方)。 注意，運算是在矩陣意義下進行的，單個數(shù)據(jù)的算術運算只是一種特例。,(1) 矩陣加減運算 假定有兩個矩陣A和B，則可以由A+B和A-B實現(xiàn)

2、矩陣的加減運算。 運算規(guī)則是：若A和B矩陣的維數(shù)相同，則可以執(zhí)行矩陣的加減運算，A和B矩陣的相應元素相加減。如果A與B的維數(shù)不相同，則MATLAB將給出錯誤信息，提示用戶兩個矩陣的維數(shù)不匹配。,(2) 矩陣乘法 假定有兩個矩陣A和B，若A為m×n矩陣，B為n×p矩陣，則C=A*B為m×p矩陣。,(3) 矩陣除法 在MATLAB中，有兩種矩陣除法運算：\和/，分

3、別表示左除和右除。 如果A矩陣是非奇異方陣，則A\B和B/A運算可以實現(xiàn)。 A\B等效于A的逆左乘B矩陣，也就是inv(A)*B，而B/A等效于A矩陣的逆右乘B矩陣，也就是B*inv(A)。,對于含有標量的運算，兩種除法運算的結果相同，如 3/4和 4\3有相同的值，都等于0.75。又如，設a=[10.5 B,25]，則a/5=5\a=[2.1000 5.0000]。 對于矩陣來說，左除和

4、右除表示兩種不同的除數(shù)矩陣和被除數(shù)矩陣的關系。對于矩陣運算，一般A\B≠B/A。,(4) 矩陣的乘方 一個矩陣的乘方運算可以表示成A^x，要求A為方陣，x為標量。,2. 點運算 在MATLAB中，有一種特殊的運算，因為其運算符是在有關算術運算符前面加點，所以叫點運算。 點運算符有 .*、./、.\和 .^。兩矩陣進行點運算是指它們的對應元素進行相關運算，要求兩矩陣的維數(shù)相同。,3.1.2

5、 關系運算 MATLAB提供了6種關系運算符：(大于)、>=(大于或等于)、==(等于)、～=(不等于)。 它們的含義不難理解，但要注意其書寫方法與數(shù)學中的不等式符號不盡相同。,關系運算符的運算法則為： (1) 當兩個比較量是標量時，直接比較兩數(shù)的大小。若關系成立，關系表達式結果為1，否則為0。 (2) 當參與比較的量是兩個維數(shù)相同的矩陣時，比較是對兩矩陣相同位置的元

6、素按標量關系運算規(guī)則逐個進行，并給出元素比較結果。最終的關系運算的結果是一個維數(shù)與原矩陣相同的矩陣，它的元素由0或1組成。,(3) 當參與比較的一個是標量， 而另一個是矩陣時，則把標量與矩陣的每一個元素按標量關系運算規(guī)則逐個比較，并給出元素比較結果。最終的關系運算的結果是一個維數(shù)與原矩陣相同的矩陣，它的元素由0或1組成。,例3-1 產(chǎn)生5階隨機方陣 A，其元素為 [10,90]區(qū)間的隨機整數(shù)，然后判斷A的元素是否能被3整除。

7、 (1) 生成5階隨機方陣A： A=fix((90-10+1)*rand(5)+10) (2) 判斷A的元素是否可以被3整除： P=rem(A,3)==0 其中，rem(A,3)是矩陣A的每個元素除以3的余數(shù)矩陣。此時，0被擴展為與A同維數(shù)的零矩陣，P是進行等于(==)比較的結果矩陣。,3.1.3 邏輯運算 MATLAB提供

8、了3種邏輯運算符：&(與)、|(或)和～(非)。 邏輯運算的運算法則為： (1) 在邏輯運算中，確認非零元素為真，用1表示，零元素為假，用0表示。,邏輯運算的運算法則為： (1) 在邏輯運算中，確認非零元素為真，用1表示，零元素為假，用0表示。 (2) 設參與邏輯運算的是兩個標量a和b，那么： a&b a,b全為非零時，運算結果為1，否則

9、為0。 a|b a,b中只要有一個非零，運算結果為1。 ～a 當a是零時，運算結果為1；當a非零時，運算結果為0。,(3) 若參與邏輯運算的是兩個同維矩陣，那么運算將對矩陣相同位置上的元素按標量規(guī)則逐個進行。最終運算結果是一個與原矩陣同維的矩陣，其元素由 1或 0組成。 (4) 若參與邏輯運算的一個是標量，一個是矩陣，那么運算將在標量與矩陣中的每個元素之間按標量規(guī)則逐個進行。最

10、終運算結果是一個與矩陣同維的矩陣，其元素由 1或 0組成。,(5) 邏輯非是單目運算符，也服從矩陣運算規(guī)則。 (6) 在算術、關系、邏輯運算中，算術運算優(yōu)先級最高，邏輯運算優(yōu)先級最低。,例3-2 建立矩陣A，然后找出大于4的元素的位置。 (1) 建立矩陣A： A=[4,-65,-54,0,6;56,0,67,-45,0] (2) 找出大于4的元素的位置：

11、 find(A>4),3.2 矩陣分析,3.2.1 對角陣與三角陣,1. 對角陣 只有對角線上有非 0 元素的矩陣稱為對角矩陣，對角線上的元素相等的對角矩陣稱為數(shù)量矩陣，對角線上的元素都為 1 的對角矩陣稱為單位矩陣。,(1) 提取矩陣的對角線元素 設A為m×n矩陣，diag(A)函數(shù)用于提取矩陣A主對角線元素，產(chǎn)生一個具有min(m,n)個元素的列向量。

12、 diag(A)函數(shù)還有一種形式diag(A,k)，其功能是提取第k條對角線的元素。 (2) 構造對角矩陣 設V為具有 m 個元素的向量，diag(V) 將產(chǎn)生一個m×m對角矩陣，其主對角線元素即為向量V的元素。 diag(V) 函數(shù)也有另一種形式 diag(V,k)，其功能是產(chǎn)生一個n×n(n=m+)對角陣，其第k條對角線的元素即為向量V的元素。,例3-3

13、先建立 5×5矩陣A，然后將A的第一行元素乘以1，第二行乘以2，…，第五行乘以5。 A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;11,18,25,2,19]; D=diag(1:5); D*A %用D左乘A，對A的每行乘以一個指定常

14、數(shù),2. 三角陣 三角陣又進一步分為上三角陣和下三角陣，所謂上三角陣，即矩陣的對角線以下的元素全為 0 的一種矩陣，而下三角陣則是對角線以上的元素全為0的一種矩陣。,(1) 上三角矩陣 求矩陣A的上三角陣的MATLAB函數(shù)是triu(A)。 triu(A) 函數(shù)也有另一種形式 triu(A,k)，其功能是求矩陣A的第k條對角線以上的元素。例如，提取矩陣A的第2條對角線以上的元素，形成新的

15、矩陣B。 (2) 下三角矩陣 在MATLAB中，提取矩陣A的下三角矩陣的函數(shù)是tril(A)和tril(A,k)，其用法與提取上三角矩陣的函數(shù) triu(A)和 triu(A,k)完全相同。,3.2.2 矩陣的轉置與旋轉1. 矩陣的轉置 轉置運算符是單撇號(’)。2. 矩陣的旋轉 利用函數(shù)rot90(A,k)將矩陣A旋轉90º的k倍，當k為1時可省略。,

16、3. 矩陣的左右翻轉 對矩陣實施左右翻轉是將原矩陣的第一列和最后一列調換，第二列和倒數(shù)第二列調換，…，依次類推。 MATLAB對矩陣A實施左右翻轉的函數(shù)是fliplr(A)。4. 矩陣的上下翻轉 MATLAB對矩陣A實施上下翻轉的函數(shù)是flipud(A)。,3.2.3 矩陣的逆與偽逆1. 矩陣的逆 對于一個方陣A，如果存在一個與其同階的方陣B，使得：

17、 A·B=B·A=I (I為單位矩陣) 則稱B為A的逆矩陣，當然，A也是B的逆矩陣。 求一個矩陣的逆是一件非常煩瑣的工作，容易出錯，但在MATLAB中，求一個矩陣的逆非常容易。 求方陣A的逆矩陣可調用函數(shù)inv(A)。,例3-4 用求逆矩陣的方法解線性方程組。 Ax=b 其解為： x=A-1b 命令為：

18、 a=[2,-3,1;8,3,2;45,1,-9]; b=[4;2;17]; x=inv(a)*b,2. 矩陣的偽逆 如果矩陣A不是一個方陣，或者 A是一個非滿秩的方陣時，矩陣A沒有逆矩陣，但可以找到一個與A的轉置矩陣 A’同型的矩陣B，使得： A·B·A=A B·A·B=B 此時稱矩陣B

19、為矩陣A的偽逆，也稱為廣義逆矩陣。 在MATLAB中，求一個矩陣偽逆的函數(shù)是pinv(A)。,3.2.4 方陣的行列式 把一個方陣看作一個行列式，并對其按行列式的規(guī)則求值，這個值就稱為矩陣所對應的行列式的值。 在MATLAB中，求方陣A所對應的行列式的值的函數(shù)是det(A)。,3.2.5 矩陣的秩與跡1. 矩陣的秩 矩陣線性無關的行數(shù)與列數(shù)稱為矩陣的秩。在MATLAB

20、中，求矩陣秩的函數(shù)是rank(A)。2. 矩陣的跡 矩陣的跡等于矩陣的對角線元素之和，也等于矩陣的特征值之和。在 MATLAB 中，求矩陣的跡的函數(shù)是trace(A)。,3.2.6 向量和矩陣的范數(shù) 矩陣或向量的范數(shù)用來度量矩陣或向量在某種意義下的長度。范數(shù)有多種方法定義，其定義不同，范數(shù)值也就不同。,1. 向量的3種常用范數(shù)及其計算函數(shù) 在MATLAB中，求向量范數(shù)的函數(shù)為：

21、 (1) norm(V)或norm(V,2)：計算向量V的2—范數(shù)。 (2) norm(V,1)：計算向量V的1—范數(shù)。 (3) norm(V,inf)：計算向量V的∞—范數(shù)。,2. 矩陣的范數(shù)及其計算函數(shù) MATLAB提供了求3種矩陣范數(shù)的函數(shù)，其函數(shù)調用格式與求向量的范數(shù)的函數(shù)完全相同。,3.2.7 矩陣的條件數(shù) 在MATLAB中，計算矩陣A的3種條件數(shù)的

22、函數(shù)是： (1) cond(A,1) 計算A的1—范數(shù)下的條件數(shù)。 (2) cond(A)或 cond(A,2) 計算A的 2—范數(shù)數(shù)下的條件數(shù)。 (3) cond(A,inf) 計算 A的 ∞—范數(shù)下的條件數(shù)。,3.2.8 矩陣的特征值與特征向量 在MATLAB中，計算矩陣 A 的特征值和特征向量的函數(shù)是eig(A)，常用的調用格式有3種：

23、 (1) E=eig(A)：求矩陣 A 的全部特征值，構成向量E。 (2) [V,D]=eig(A)：求矩陣A的全部特征值，構成對角陣 D，并求 A的特征向量構成V的列向量。,(3) [V,D]=eig(A,‘nobalance’)：與第2種格式類似，但第 2種格式中先對A作相似變換后求矩陣A的特征值和特征向量，而格式 3直接求矩陣A的特征值和特征向量。,例3-5 用求特征值的方法解方程： 3x5-7x

24、4+5x2+2x-18=0 p=[3,-7,0,5,2,-18]; A=compan(p); %A的伴隨矩陣 x1=eig(A) %求A的特征值 x2=roots(p) %直接求多項式p的零點,3.3 字符串,在MATLAB中，字符串是用單撇號(‘)括起來的字符序列。 MATLAB 將字符串當作一個行向量，每個

25、元素對應一個字符，其標識方法和數(shù)值向量相同。也可以建立多行字符串矩陣。,字符串是以ASCII碼形式存儲的。abs和double函數(shù)都可以用來獲取字符串矩陣所對應的ASCII碼數(shù)值矩陣。 相反，char函數(shù)可以把ASCII碼矩陣轉換為字符串矩陣。,例3-6 建立一個字符串向量，然后對該向量做如下處理： (1) 取第1～5個字符組成的子字符串。 (2) 將字符串倒過來重新排列。

26、 (3) 將字符串中的小寫字母變成相應的大寫字母，其余字符不變。 (4) 統(tǒng)計字符串中小寫字母的個數(shù)。,命令如下：ch='ABc123d4e56Fg9';subch=ch(1:5) %取子字符串revch=ch(end:-1:1) %將字符串倒排k=find(ch>='a'&ch<='z'); %找