矩陣運(yùn)算的基礎(chǔ)

關(guān)于矩陣運(yùn)算的基礎(chǔ)第一頁，共四十二頁，2022年，8月28日MATLAB提供了一種計(jì)算機(jī)高級(jí)編程語言——M語言MATLAB提供了不同類型的數(shù)據(jù)MATLAB專門以矩陣作為基本的運(yùn)算單位MATLAB提供了關(guān)于數(shù)組和矩陣不同的運(yùn)算方法第二頁，共四十二頁，2022年，8月28日一.概述在M語言中最常用的數(shù)據(jù)類型表現(xiàn)手段和形式就是變量和常量M語言的基本處理單位是數(shù)值矩陣或者數(shù)值向量回顧有關(guān)概念變量和常量數(shù)組向量矩陣第三頁，共四十二頁，2022年，8月28日變量和常量變量：程序運(yùn)行過程中需要改變數(shù)值的量每一個(gè)變量都具有一個(gè)名字變量在內(nèi)存中占據(jù)一定的空間變量必須以字母開頭，后面可以是字母、數(shù)字或者下劃線的組合MATLAB僅識(shí)別前面N個(gè)字符，在不同的操作系統(tǒng)下可以識(shí)別的字符個(gè)數(shù)不同常量：在程序運(yùn)行的過程中不需要改變數(shù)值的量常量具有名字在M語言中不存在常量的定義，只在MATLAB中提供一些常用的常數(shù)作為常量第四頁，共四十二頁，2022年，8月28日數(shù)組是有序數(shù)據(jù)的集合數(shù)組的每一個(gè)成員（元素）都屬于同一種數(shù)據(jù)類型，它們使用同一個(gè)數(shù)組名稱和不同的下標(biāo)來唯一確定數(shù)組中的成員（元素）。在MATLAB中元胞數(shù)組比較特殊，數(shù)組中的元素可以是不同的數(shù)據(jù)類型。第五頁，共四十二頁，2022年，8月28日向量從編程語言的角度上看，向量其實(shí)就是一維數(shù)組從數(shù)學(xué)的角度上看，向量就是1×N或者N×1的矩陣，即行向量或列向量

b1，1b2，1B=b3，1和B=[b1，1b1，2b1，3······b1，n]∶∶bn，1第六頁，共四十二頁，2022年，8月28日矩陣是用一對(duì)圓括號(hào)或方括號(hào)括起來，符合一定規(guī)則的數(shù)學(xué)對(duì)象

b11b12b13B=b21b22b23b31b32b33對(duì)于編程語言，矩陣就是二維的數(shù)組第七頁，共四十二頁，2022年，8月28日二.矩陣和向量創(chuàng)建1.向量的創(chuàng)建（1）在命令窗口逐個(gè)輸入元素

>>X=[13pi3+5i]（2）利用冒號(hào)運(yùn)算符創(chuàng)建向量

X=J:INC:KJ為向量的第一個(gè)元素，K為向量的最后一個(gè)元素，INC為向量元素遞增的步長J、INC、K之間必須用“:”間隔若忽略INC，則默認(rèn)的遞增步長為1INC可以為正數(shù)，也可以為負(fù)數(shù)

>>X=1:10>>X=1:0.01:1.1

第八頁，共四十二頁，2022年，8月28日(3)定數(shù)線性采樣法：在設(shè)定的“總點(diǎn)數(shù)”下，均勻采樣生成向量（一維“行”數(shù)組）使用函數(shù)linspace和logspacelinspace是用來創(chuàng)建線性間隔向量的函數(shù)linspace

的基本語法

X=linespace(X1，X2，n)X1為向量的第一個(gè)元素，X2為向量的最后一個(gè)元素，n為向量具有的元素個(gè)數(shù)，函數(shù)將根據(jù)n的數(shù)值平均計(jì)算元素之間的間隔，間隔計(jì)算公式為若在表達(dá)式中忽略參數(shù)n，則系統(tǒng)默認(rèn)地將向量設(shè)置為100個(gè)元素第九頁，共四十二頁，2022年，8月28日使用linspace函數(shù)創(chuàng)建向量>>X=linspace(1，2，5)X=1.00001.25001.50001.75002.0000第十頁，共四十二頁，2022年，8月28日logspace是用來創(chuàng)建對(duì)數(shù)空間的向量logspace

的基本語法

X=logspace(X1，X2，n)該函數(shù)創(chuàng)建的向量第一個(gè)元素值為10X1，而最后一個(gè)元素的數(shù)值為10X2，n為向量的元素個(gè)數(shù)，元素彼此之間的間隔按照對(duì)數(shù)空間的間隔設(shè)置若在表達(dá)式中忽略參數(shù)n，則系統(tǒng)默認(rèn)地將向量設(shè)置為50個(gè)元素第十一頁，共四十二頁，2022年，8月28日例使用logspace函數(shù)創(chuàng)建向量>>X=logspace(1，3，3)X=101001000第十二頁，共四十二頁，2022年，8月28日創(chuàng)建列向量使用分號(hào)作為元素與元素之間的間隔使用轉(zhuǎn)置運(yùn)算符“'”例2-6：>>A=[1；2；3；4；5；6]或>>A=(1:6)'第十三頁，共四十二頁，2022年，8月28日2.創(chuàng)建矩陣矩陣的元素可以為任意MATLAB數(shù)據(jù)類型的數(shù)值或?qū)ο髣?chuàng)建矩陣的方法直接輸入法使用數(shù)組編輯器

第十四頁，共四十二頁，2022年，8月28日(1)直接輸入法規(guī)則：整個(gè)矩陣的元素必須用[]括住同一行的矩陣元素之間必須用逗號(hào)或空格分隔在[]內(nèi)矩陣的行與行之間必須用分號(hào)分隔，也可以在需要分行的地方用回車鍵間隔矩陣元素可以是任何MATLAB表達(dá)式，可以是實(shí)數(shù)，也可以是復(fù)數(shù)，復(fù)數(shù)用i，j輸入

例2-7：

>>A=[1，2，3；4，5，6；7，8，9]>>X=[2pi/2；sqrt(3)3+5i]

第十五頁，共四十二頁，2022年，8月28日逗號(hào)和分號(hào)的作用逗號(hào)和分號(hào)可作為指令間的分隔符，MATLAB允許多條語句在同一行出現(xiàn)。分號(hào)如果出現(xiàn)在指令后，屏幕上將不顯示結(jié)果。只要是賦過值的變量，不管是否在屏幕上顯示過，都會(huì)存儲(chǔ)在工作空間中，以后可隨時(shí)顯示或調(diào)用。變量名盡可能不要重復(fù)，否則會(huì)覆蓋。例2-8：

>>A=[1,2,3;4,5,6;7,8,9],X=[2pi/2;sqrt(3)3+5i]

？

>>A=[1,2,3;4,5,6;7,8,9];X=[2pi/2;sqrt(3)3+5i]

？第十六頁，共四十二頁，2022年，8月28日例：

>>A=[1,2,3;4,5,6;7,8,9],X=[2pi/2;sqrt(3)3+5i]A=123456789X=2.00001.57081.73213.0000+5.0000i>>A=[1,2,3;4,5,6;7,8,9];X=[2pi/2;sqrt(3)3+5i]X=2.00001.57081.73213.0000+5.0000i第十七頁，共四十二頁，2022年，8月28日冒號(hào)的作用用于生成等間隔的向量，默認(rèn)間隔為1。例2-9：>>y=[1:3；4:6；7:9]第十八頁，共四十二頁，2022年，8月28日矩陣的基本運(yùn)算函數(shù)基本數(shù)學(xué)運(yùn)算規(guī)則數(shù)組的運(yùn)算運(yùn)算函數(shù)運(yùn)算指令三.基本運(yùn)算第十九頁，共四十二頁，2022年，8月28日函數(shù)說明zeros產(chǎn)生元素全為0的矩陣ones產(chǎn)生元素全為1的矩陣eye產(chǎn)生單位矩陣rand產(chǎn)生均勻分布的隨機(jī)數(shù)矩陣，數(shù)值范圍（0，1）randn產(chǎn)生均值為0，方差為1的正態(tài)分布隨機(jī)數(shù)矩陣diag獲取矩陣的對(duì)角線元素，也可生成對(duì)角矩陣tril產(chǎn)生下三角矩陣triu產(chǎn)生上三角矩陣pascal產(chǎn)生帕斯卡矩陣magic產(chǎn)生幻方陣第二十頁，共四十二頁，2022年，8月28日矩陣生成函數(shù)示例>>A=zeros(3)A=000000000>>A=ones(3)A=111111111>>A=eye(3)A=100010001>>A=rand(3)A=0.95010.48600.45650.23110.89130.01850.60680.76210.8214>>A=randn(3)A=-0.43260.28771.1892-1.6656-1.1465-0.03760.12531.19090.3273第二十一頁，共四十二頁，2022年，8月28日矩陣生成函數(shù)示例>>A=magic(3)A=816357(15)492>>A=magic(4)A=162313511108(34)97612414151第二十二頁，共四十二頁，2022年，8月28日

基本矩陣運(yùn)算運(yùn)算命令說明A’矩陣轉(zhuǎn)置A^n矩陣求冪，n可以為任意實(shí)數(shù)A*B矩陣相乘A/B矩陣右除（一般的除法，A/B=A÷B）A\B矩陣左除（一種倒置的除法，A\B=B÷A）A+B矩陣相加A-B矩陣相減inv矩陣求逆，注意不是所有的矩陣都有逆矩陣det求方陣的行列式rank求矩陣的秩eig求矩陣的特征向量和特征值svd對(duì)矩陣進(jìn)行奇異值分解norm求矩陣的范數(shù)第二十三頁，共四十二頁，2022年，8月28日1.矩陣加、減運(yùn)算（A＋B、A－B）規(guī)則：相加、減的兩矩陣必須有相同的行和列，兩矩陣對(duì)應(yīng)元素相加減。MATLAB允許參與運(yùn)算的兩矩陣之一是標(biāo)量，標(biāo)量與矩陣的所有元素分別進(jìn)行加減操作。例：A=[123；456]B=[345；789]C=3A+B=[468；111315]A+C=[456；789]B+C=[678；101112]第二十四頁，共四十二頁，2022年，8月28日2.矩陣乘運(yùn)算

A*B：A矩陣的列數(shù)必須等于B矩陣的行數(shù)。

s*A或A*s：標(biāo)量可與任何矩陣相乘，標(biāo)量s分別與矩陣A每個(gè)元素相乘。例：>>A=[123;456;780];B=[1;2;3];>>C=A*BC=143223>>D=[-1;0;2];>>F=pi*DF=-3.141606.2832

基本矩陣運(yùn)算（續(xù)）第二十五頁，共四十二頁，2022年，8月28日

3.矩陣除運(yùn)算及線性方程組的解在線性代數(shù)中沒有矩陣的除運(yùn)算，只有矩陣逆的運(yùn)算，在MATLAB中有兩種矩陣除運(yùn)算。

A/B—矩陣右除，相當(dāng)于Ainv(B) A\B—矩陣左除，相當(dāng)于inv(A)B

因此，x=A\B是線性方程組Ax=B的解。例：求解方程組3x1+x2-x3=3.6x1+2x2+4x3=2.1-x1+4x2+5x3=-1.4>>A=[31-1;124;-145];>>B=[3.6;2.1;-1.4];

>>x=A\Bx= 1.4818 -0.4606 0.3848基本矩陣運(yùn)算（續(xù)）第二十六頁，共四十二頁，2022年，8月28日4.矩陣乘方

A^n——A自乘n次冪

例>>a=[1,2,3;4,5,6;7,8,9];>>a^2ans=303642668196102126150

方陣>1的整數(shù)基本矩陣運(yùn)算（續(xù)）第二十七頁，共四十二頁，2022年，8月28日數(shù)組運(yùn)算指元素對(duì)元素的算術(shù)運(yùn)算，與通常意義上的由符號(hào)表示的線性代數(shù)矩陣運(yùn)算不同。1. 數(shù)組加減（+，-）運(yùn)算規(guī)則：相加、減的兩數(shù)組必須有相同的行和列，兩數(shù)組對(duì)應(yīng)元素相加減。MATLAB允許參與運(yùn)算的兩數(shù)組之一是標(biāo)量，標(biāo)量與數(shù)組的所有元素分別進(jìn)行加減操作

A+B A-B

基本數(shù)組（元素群）運(yùn)算

與矩陣加減運(yùn)算等效，數(shù)組之一也可為標(biāo)量。第二十八頁，共四十二頁，2022年，8月28日2.數(shù)組乘()

運(yùn)算AB

A，B兩數(shù)組必須有相同的行和列，兩數(shù)組相應(yīng)元素相乘。

sA或As

標(biāo)量與數(shù)組相乘，標(biāo)量s分別與數(shù)組A每個(gè)元素相乘，與sA或As相同。例16：>>A=[123;456;789];>>B=[246;135;7910];>>A.*Bans= 2818 41530 497290>>A=[123;456;789];>>B=[246;135;7910];>>A*Bans= 253746 5585109 85133172第二十九頁，共四十二頁，2022年，8月28日3.數(shù)組除(/，\)運(yùn)算C=A./B——數(shù)組右除

C(i,j)=A(i,j)/B(i,j)C=A.\B——數(shù)組左除

C(i,j)=B(i,j)/A(i,j)A./B=B.\AA./s=s.\A—A的元素分別被標(biāo)量s除s./A=A.\s—標(biāo)量s分別被A的元素除例：>>A=[123];>>B=[456];>>C1=A./BC1=0.25000.40000.5000>>C2=B.\AC2=0.25000.40000.5000>>C3=A.\BC3=4.00002.50002.0000>>A=[123];B=[456];>>A/Bans=0.4156>>A\Bans=0000001.33331.66672.0000第三十頁，共四十二頁，2022年，8月28日4.數(shù)組乘方（.^）A.^n——A的每個(gè)元素自乘n次A.^p——對(duì)A各元素分別求非整數(shù)冪p.^A——以p為底,分別以A的元素為指數(shù)求冪值C=A.^B——元素對(duì)元素的冪C(i,j)=A(i,j).^B(i,j)例：>>A=[123];B=[456];>>X=A.^2X=1.004.009.00>>Y=A.^0.5Y=1.00001.41421.7321

>>C=3.^BY=81.00243.00729.00

3^43^53^6>>Z=A.^BZ=1.0032.00729.001^42^53^6第三十一頁，共四十二頁，2022年，8月28日5.數(shù)組轉(zhuǎn)置（.’）例：>>A=[135;246]A=135246>>A'ans=123456>>A.'ans=123456結(jié)論：對(duì)于實(shí)數(shù)矩陣，矩陣轉(zhuǎn)置和數(shù)組轉(zhuǎn)置的計(jì)算結(jié)果是一致的。第三十二頁，共四十二頁，2022年，8月28日例：>>A=A*iA=0+1.0000i0+3.0000i0+5.0000i0+2.0000i0+4.0000i0+6.0000i>>A'ans=0-1.0000i0-2.0000i0-3.0000i0-4.0000i0-5.0000i0-6.0000i>>A.'ans=0+1.0000i0+2.0000i0+3.0000i0+4.0000i0+5.0000i0+6.0000i結(jié)論：對(duì)于復(fù)數(shù)矩陣，矩陣轉(zhuǎn)置和數(shù)組轉(zhuǎn)置的計(jì)算結(jié)果不一致。矩陣轉(zhuǎn)置運(yùn)算——共軛轉(zhuǎn)置數(shù)組轉(zhuǎn)置運(yùn)算——非共軛轉(zhuǎn)置第三十三頁，共四十二頁，2022年，8月28日函數(shù)的主要類別三角函數(shù)指數(shù)運(yùn)算函數(shù)復(fù)數(shù)運(yùn)算函數(shù)圓整和求余函數(shù)函數(shù)在處理參數(shù)時(shí)，是按照數(shù)組運(yùn)算的規(guī)則進(jìn)行的

基本數(shù)學(xué)函數(shù)

第三十四頁，共四十二頁，2022年，8月28日三角函數(shù)函數(shù)說明函數(shù)說明函數(shù)說明sin正弦函數(shù)tanh雙曲正切函數(shù)csch雙曲余割函數(shù)sinh雙曲正弦函數(shù)atan反正切函數(shù)acsc反余割函數(shù)asin反正弦函數(shù)atan2四象限反正切函數(shù)acsch反雙曲余割函數(shù)asinh反雙曲正弦函數(shù)atanh反雙曲正切函數(shù)cot余切函數(shù)cos余弦函數(shù)sec正割函數(shù)coth雙曲余切函數(shù)cosh雙曲余弦函數(shù)sech雙曲正割函數(shù)acot反余切函數(shù)acos反余弦函數(shù)asec反正割函數(shù)acoth反雙曲余切函數(shù)acosh反雙曲余弦函數(shù)asech雙曲反正割函數(shù)tan正切函數(shù)csc余割函數(shù)第三十五頁，共四十二頁，2022年，8月28日指數(shù)運(yùn)算函數(shù)函數(shù)說明函數(shù)說明exp指數(shù)函數(shù)realpow實(shí)數(shù)冪運(yùn)算函數(shù)log自然對(duì)數(shù)函數(shù)reallog實(shí)數(shù)自然對(duì)數(shù)函數(shù)log10常用對(duì)數(shù)函數(shù)realsqrt實(shí)數(shù)平方根函數(shù)log