ch2 Matlab矩陣的生成與運(yùn)算.ppt

1、Matlab矩陣的生成與運(yùn)算,安徽工業(yè)大學(xué)數(shù)理學(xué)院 侯為根 ,【例1-1】表示矩陣,一、MATLAB矩陣的生成,1、直接輸入法,將矩陣的元素用方括號括起來，按矩陣行的順序輸入各元素，同一行的各元素之間用空格或逗號分隔，不同行的元素之間用分號分隔。,A=1,2,3;4,5,6;7,8,0,A=1,2,3;4 5,6;7,8 0,【例1-2】 試輸入復(fù)數(shù)矩陣,需要避免的語句,B=1+9*i,2+8*i,3+7*j; 4+6*j 5+5*i,6+4*i;7+3*i,8+2*i,i,B=1 +9*i,2+8*i,3+7*j; 4+6*j 5+5*i,6+4*i;7+3*i,8+2*i,i,2、利用M文

2、件建立矩陣,對于比較大且比較復(fù)雜的矩陣，可以為它專門建立一個M文件。,【例1-3】利用M文件建立mymat矩陣。,(1) 啟動有關(guān)編輯程序或MATLAB文本編輯器，并輸入待建矩陣.,(2) 把輸入的內(nèi)容以純文本方式存盤(設(shè)文件名為mymat.m)。,(3) 運(yùn)行該M文件，就會自動建立一個名為mymat的矩陣，可供以后使用。,A=4 10 1 6 2;8 2 9 4 7; 7 5 7 1 5;0 3 4 5 4;23 13 13 0 3,3、利用MATLAB函數(shù)建立矩陣,幾個產(chǎn)生特殊矩陣的函數(shù)：zeros、ones、 eye、rand、randn、pascal、 magic(n)、 vander

3、(V)、 hilb(n) 、 toeplitz(x,y)、compan(P) 。,這幾個函數(shù)的調(diào)用格式相似，下面以產(chǎn)生零矩陣的zeros函數(shù)為例進(jìn)行說明。其調(diào)用格式是：,zeros(m) 產(chǎn)生mm零矩陣,zeros(m,n) 產(chǎn)生mn零矩陣。,zeros(size(A) 產(chǎn)生與矩陣A同樣大小的零矩陣,相關(guān)的函數(shù)有：length(A)給出行數(shù)和列數(shù)中的較大者，即length(A)=max(size(A)；ndims(A)給出A的維數(shù)。,【例1-4】 分別建立33、32和與矩陣A同樣大小的零矩陣。,(1) 建立一個33零矩陣：zeros(3),(2) 建立一個32零矩陣：zeros(3,2),(3

4、) 建立與矩陣A同樣大小零矩陣：zeros(size(A),此外，常用的函數(shù)還有reshape(A,m,n)，它在矩陣總元素保持不變的前提下，將矩陣A重新排成mn的二維矩陣。,【例1-5】 將4階魔方陣形變?yōu)?行8列的矩陣,reshaps(magic(4),2,8),hilb(n)指令可以產(chǎn)生nn的 Hilbert 矩陣,Hilbert矩陣,Hilbert 矩陣的特性: 當(dāng)矩陣變大時(shí)，其矩陣會接近奇異（即矩陣的行列式會接近于0），Hilbert矩陣常被用來評估各種逆矩陣計(jì)算方法的穩(wěn)定性。,【例1-6】 計(jì)算6階hilbert矩陣的行列式,det(hilb(6),x1 = rand(10000,

5、 1); x2 = randn(10000, 1); subplot(2,1,1); hist(x1, 40); title(均勻分布); subplot(2,1,2); hist(x2, 40); title(高斯分布); set(findobj(gcf, type, patch),EdgeColor, w); % 改邊線為白色,【例1-7】產(chǎn)生10000個均勻均勻與正態(tài)分布的隨機(jī)數(shù),產(chǎn)生隨機(jī)數(shù)矩陣指令rand 和 randn,4、建立大矩陣,大矩陣可由方括號中的小矩陣建立起來。,A=1 2 3;4 5 6;7 8 0;,5、冒號表達(dá)式,冒號表達(dá)式的一般格式：v=s1:s2:s3,還可以用l

6、inspace函數(shù)產(chǎn)生行向量：linspace(a,b,n),linspace(a,b,n)與a:(b-a)/(n-1):b等價(jià),C=A,eye(size(A);ones(size(A),A,A=A;1 2 3;1;2;3;4;,【例1-8】使用小矩陣構(gòu)造大矩陣,【例1-9】 用不同的步距生成 (0,p) 間向量,V1=0:0.2:pi,V3=0:-1:pi,V2=0:pi,V4=pi:-1:0,V5=0:0.2:pi,pi,V7=0:pi/100:pi,V6=linspace(0,pi,100),二、對矩陣元素的操作,矩陣A中，位于第i行、第j列的元素可表示為A(i,j) i與j即是此元素的

7、下標(biāo)（Subscript）或索引（Index） MATLAB中所有矩陣的內(nèi)部表示法都是以列為主的一維向量 A(i,j)和A(i+(j-1)*m)是完全一樣的m為矩陣A的列數(shù),1、矩陣的索引或下標(biāo),我們可以使用一維或二維下標(biāo)來存取矩陣,A=4 10 1 6 2;8 2 9 4 7; 7 5 7 1 5;0 3 4 5 4;23 13 13 0 3,2 、矩陣元素,MATLAB允許用戶對一個矩陣的單個元素進(jìn)行賦值和操作。例如 A(3,2)=200,也可以采用矩陣元素的序號來引用矩陣元素。矩陣元素按列編號，先第一列，再第二列，依次類推。,以mn矩陣A為例，矩陣元素A(i,j)的序號為：(j-1) *

8、m+i。其相互轉(zhuǎn)換關(guān)系也可利用sub2ind和ind2sub函數(shù)求得。,3. 子矩陣提取,(1) 利用冒號表達(dá)式獲得子矩陣, A(:,j)表示取A矩陣的第j列全部元素；A(i,:)表示A矩陣第i行的全部元素；A(i,j)表示取A矩陣第i行、第j列的元素。, A(i:i+m,:)表示取A矩陣第ii+m行的全部元素；A(:,k:k+m)表示取A矩陣第kk+m列的全部元素，A(i:i+m,k:k+m)表示取A矩陣第ii+m行內(nèi)，并在第kk+m列中的所有元素。, 利用一般向量和end運(yùn)算符等來表示矩陣下標(biāo)，從而獲得子矩陣。end表示某一維的末尾元素下標(biāo)。,【例2-1】子矩陣提取,提取矩陣A 全部奇數(shù)行

9、，所有列,將矩陣A 左右翻轉(zhuǎn),提取矩陣A的 3,2,1 行、2,3,4 列構(gòu)成子矩陣,將矩陣A 2,4列全變?yōu)?,B1=A(1:2:end,:),B3=A(:,end:-1:1),B2=A(3,2,1,2:4),A(:,2,4)=ones(4,2),4、利用空矩陣刪除矩陣的元素,在MATLAB中，定義 為空矩陣。給變量X賦空矩陣的語句為X= 。,【注意】X= 與clear X不同，clear是將X從工作空間中刪除，而空矩陣則存在于工作空間，只是維數(shù)為0。,將某些元素從矩陣中刪除，采用將其置為空矩陣的方法就是一種有效的方法。,【例2-2】將矩陣A 2,4行刪去,A(2,3,:=,5 矩陣翻轉(zhuǎn),左

10、右翻轉(zhuǎn):,上下翻轉(zhuǎn):,旋轉(zhuǎn) 90o:,如何旋轉(zhuǎn)180o？,fliplr(A),flipud(A),rot90(A),rot90(rot90(A),A=1 2 3 4;5 6 7 8;9 10 11 12 A1=diag(A) A2=diag(A,1) A3=diag(A,-1) B1=diag(A1,1) B2=diag(A3,1),6 選取對角元素,7 選取上下）三角元陣,c1=tril(A) c2=tril(A,1) c3=tril(A,-1) d1=triu(A) d2=triu(A,1) D3=triu(A,-1),1、直接創(chuàng)建稀疏矩陣,例：,S=sparse(i,j,s,m,n),其

11、中i 和j 分別是矩陣非零元素的行和列指標(biāo)向量，s 是非零元素值向量，m，n 分別是矩陣的行數(shù)和列數(shù)。,S2=sparse(1 2 2 3 4 4,3 1 3 4 2 3,. 5 3 3 1 4 3,4,4),稀疏矩陣的輸入與生成,2、從文件中創(chuàng)建稀疏矩陣,例：設(shè)文本文件 T.txt 中有三列內(nèi)容，第一列是一些行下標(biāo)，第二列是列下標(biāo)，第三列是非零元素值。,1 3 5 2 1 3 2 3 3 3 4 1 4 2 4 4 3 3,利用load和spconvert函數(shù)可以從包含一系列下標(biāo)和非零元素的文本文件中輸入稀疏矩陣。,load T.txt S=spconvert(T),3、稀疏帶狀矩陣的創(chuàng)建,

12、例：,S=spdiags(B,d,m,n),其中m 和n 分別是矩陣的行數(shù)和列數(shù)；d是長度為p的整數(shù)向量，它指定矩陣S的對角線位置；B是全元素矩陣，用來給定S對角線位置上的元素，行數(shù)為min(m,n)，列數(shù)為p 。,B=rand(4,2); S3=spdiags(B,0 1,4,4),4、滿矩陣與稀疏矩陣之間的轉(zhuǎn)化,例：,S=sparse(A),A=full(S),A=0 0 5 0; 3 0 3 0 ; 0 0 0 1; 0 4 3 0 S=sparse(A) whos,1. 數(shù)學(xué)運(yùn)算,2. 邏輯運(yùn)算,三、 Matlab基本運(yùn)算,3. 比較運(yùn)算,4. 位運(yùn)算符,5. 集合運(yùn)算符,一、數(shù)學(xué)運(yùn)算

13、,MATLAB的基本算術(shù)運(yùn)算有： (加)、(減)、*(乘)、/(右除)、(左除)、(乘方)、共軛轉(zhuǎn)置（）正號（+）,負(fù)號（-）,【注意】運(yùn)算是在矩陣意義下進(jìn)行的，單個數(shù)據(jù)的算術(shù)運(yùn)算只是一種特例。,2. 點(diǎn)運(yùn)算,點(diǎn)運(yùn)算符有.*、./、.和.。兩矩陣進(jìn)行點(diǎn)運(yùn)算是指它們的對應(yīng)元素進(jìn)行相關(guān)運(yùn)算，要求兩矩陣的維參數(shù)相同。,1. 基本算術(shù)運(yùn)算,1.轉(zhuǎn)置（.）,冪次（.）.共軛轉(zhuǎn)置（）,矩陣冪次（） 2.正號（+）,負(fù)號（-） 3.乘法（.*）,元素右除（./）,元素左除（.）,矩陣乘法（*）,矩陣右除（/）,矩陣左除（） 4.加法（+）,減法（-） 5.冒號（:）（例如： x = 1:2:5 等）,數(shù)學(xué)運(yùn)

14、算符之優(yōu)先級:,同一類的運(yùn)算符均具有相同的優(yōu)先度(Priority),因此在計(jì)算上，是由左至右依次完成,矩陣的數(shù)學(xué)運(yùn)算,矩陣的加減與一般標(biāo)量（Scalar）的加減類似 相加或相減的矩陣必需具有相同的維度,A = 12 34 56 20; B = 1 3 2 4; C = A + B C1= A-B,例1: 加減法運(yùn)算,矩陣與純量可以直接進(jìn)行加減，MATLAB 會直接將加減應(yīng)用到每一個元素,A=1 2 3 2 1 + 5,矩陣的加減法運(yùn)算,矩陣的乘法與除法,標(biāo)量對矩陣的乘或除，可比照一般寫法,A = 1 2 3; 4 4 2; C = A/3 B = 2*A,欲進(jìn)行矩陣相乘，必需確認(rèn)第一個矩陣的

15、列數(shù)（ Column Dimension） 必需等于第二個矩陣的行數(shù)（Row Dimension）,例2:矩陣的乘法,A = 1; 2; B = 3, 4, 5; C = A*B,矩陣除法,若A為非奇異方陣，則X=A-1B,若A為非奇異方陣，則X=BA-1,矩陣左除：AX = B，求 X MATLAB 求解：X=AB,最小二乘解,矩陣右除：XA = B，求 X MATLAB求解：X=B/A最小二乘解,注：若A為非奇異方陣，則A-1在Matlab中為inv(A),矩陣的次方運(yùn)算,矩陣的次方運(yùn)算，可由“”來達(dá)成，但矩陣必需是方陣，其次方運(yùn)算才有意義,A = magic(3); B = A2,例3：

16、,復(fù)數(shù)矩陣z，其“共軛轉(zhuǎn)置”矩陣(Conjugate Transpose)可表示成矩陣z,轉(zhuǎn)置和共軛轉(zhuǎn)置矩陣,i = sqrt(-1);% 單位虛數(shù) z = 1+i, 2; 3, 1+2i; w = z % 共軛轉(zhuǎn)置(注意z后面的單引號),例4: 矩陣共軛轉(zhuǎn)置,i = sqrt(-1);% 單位虛數(shù) z = 1+i, 2; 3, 1+2i; w = z.% 單純轉(zhuǎn)置(注意z后面的句點(diǎn)及單引號),想得到任何矩陣z的轉(zhuǎn)置（Transpose），則可表示成矩陣矩陣z.,例5:矩陣的轉(zhuǎn)置,若z為實(shí)數(shù)，則z和z.的結(jié)果是一樣的,2、點(diǎn)運(yùn)算,矩陣對應(yīng)元素的直接運(yùn)算,例如,A =1 2 3；4 5 6；7

17、8 0,B = A.A,3MATLAB數(shù)組運(yùn)算函數(shù),1. 規(guī)則：設(shè)X=(xij)mn，則 f (X)=(f (xij)mn,sin cos tan cot sec csc asin acos atan acot asec acsc,2. 常用數(shù)學(xué)函數(shù),三角與反三角函數(shù),exp log log10 sqrt pow2,ceil fix floor round rem sign,abs angle conj imag real,指數(shù)與對數(shù),復(fù)變函數(shù),取整函數(shù),cart2sph cart2pol pol2cart sph2pol,坐標(biāo)轉(zhuǎn)換,注：若要對函數(shù)進(jìn)行矩陣運(yùn)算，方法為：funm(X,f),4、

18、 MATLAB矩陣運(yùn)算函數(shù),size(A) %矩陣的大小,inv(A) %矩陣的逆,det(A) %矩陣的行列式,rank(A) %矩陣的秩,rref(A) %矩陣列的極大無關(guān)組,orth(A) %矩陣化為正交陣,trace(A) %矩陣的跡,chol(A) %三角分解,A=L*LT,lu(A) %A對稱，A=L*U,U,S,V=svd(A) %A=U*S*V,qr(A) %A=q*r(q正交陣,r三角陣）,norm(A,1|2|p|inf|fro) %矩陣的模,cond(A) %矩陣的條件數(shù),4、 MATLAB矩陣運(yùn)算函數(shù),eig(A) %A*x=*x,V,D=eig(A) %A*V=V*D

19、,eig(A,B) %A*x=*B*x,V,D=eig(A,B) %A*V=B*V*D,例：設(shè)A=pascal(4),B=magic(4)，并使用下列函數(shù)進(jìn)行計(jì)算,向量的p-norm,一個向量a的p-norm可以定義為,p=2 時(shí)，此即為向量 a 的長度，或稱歐氏長度（Euclidean Length）,欲求一向量的p-norm，可使用norm指令norm(x,p),a = 3 4; x = norm(a, 1) % x = 7 y = norm(a, 2) % y = 5 z = norm(a, inf),注意：,矩陣A的p-norm可以定義如下：,例9: normMatrix,norm 指

20、令亦可用于計(jì)算矩陣的 p-norm,A = 1 2 3; 4 5 6; 7 8 9; norm(A, 2),Cholesky分解是將一個對稱的正定矩陣分解為一個上三角矩陣和其轉(zhuǎn)置的乘積。,X=pascal(5) R=chol(X) C=R*R,LU分解是將原正方 (square) 矩陣分解成一個上三角形矩陣或是排列(permuted) 的上三角形矩陣和一個 下三角形矩陣，這樣的分解法又稱為三角分解。,命令格式為：,L,U=lu(A),L,U,P=lu(A),例如：,B=1 3 2; -2 -6 1; 2 5 7 L,U=lu(B) L*U,L,U,P=lu(B) P*L*U,QR分解法是將矩陣

21、分解成一個正規(guī)正交矩陣與上三角形矩陣，矩陣A不必為正方矩陣。,A=1 3 2 5; -2 -6 1 6; 2 5 7 2 Q,R=qr(A) Q*Q Q*R,奇異值分解 (sigular value decomposition,SVD) 是另一種正交矩陣分解法；SVD是最可靠的分解法，但是它比QR 分解法要花上近十倍的計(jì)算時(shí)間。,U,S,V=svd(A)，其中U和V代表二個相互正交矩陣，而S代表一對角矩陣。 和QR分解法相同者， 原矩陣A不必為正方矩陣。,A=1,2;3,1 S,V,D=svd(A),矩陣的比較運(yùn)算, =, , =, =,=, find(), all(), any()，isem

22、pty , isnan, isfinit, islogical, isglobal, isstr 等,MATLAB的關(guān)系運(yùn)算符總是以矩陣內(nèi)元素對元素的方式作運(yùn)算,比較兩元素之間關(guān)系：若成立結(jié)論為1，否則為0,1、數(shù)量間比較,各種允許的比較關(guān)系,例1:,a=3; b=5; a5 ab a=b,關(guān)系運(yùn)算符的用法:,關(guān)系運(yùn)算符的用法,A=1 2 3; 4 5 6;7 8,0b=A5A(b),3、同維矩陣間比較；返回0-1矩陣,2、矩陣與一數(shù)比較返回0-1矩陣,例2:給出矩陣A中大與5的元素值,例3: 接上例,設(shè)B=magic(3),B=magic(3);BA,find指令 : 可傳回非零元素索引,顯

23、示一維索引或下標(biāo),元素索引,index=find(A5),可以將關(guān)系運(yùn)算符返回的0-1矩陣,來進(jìn)行矩陣的索（Indexing）求出在A矩陣中，滿足AijBij的元素,x = 0 1 2; 0 0 3 index = find(x) x(index),x(index)則是非零元素所形成的向量,要找出滿足某條件的二維索引或下標(biāo),2維索引,idx1,idx2=find(x5),idx1和idx2分別是符合條件之矩陣元素的列索引(Row Index)及行索引(Column Index),例:如何找出介于7和11的元素?,x = magic(5) x(find(7x y=sin(x); y1=(x2*p

24、i).*y; %消去負(fù)半波 q=(xpi/3 %按要求處理第(2)步,方法2： x=0:pi/100:3*pi;y=sin(x); y1=(y=0).*y; %消去負(fù)半波 p=sin(pi/3); y2=(y=p)*p+(yp).*y1; %按要求處理第(2)步,MATLAB 的位運(yùn)算符（Bitwise Operators）可對非負(fù)整數(shù)進(jìn)行位對位的邏輯運(yùn)算,位運(yùn)算符,數(shù)字10的二進(jìn)制表示法是1010，數(shù)字4的二進(jìn)制表示法是0100，因此數(shù)字10和數(shù)字4的 bitor 應(yīng)傳回1110 的十進(jìn)制數(shù)值，即數(shù)字14,數(shù)字18的二進(jìn)制表示法向左平移一個單位，其效果相當(dāng)乘2,bitor(10,4),bit

25、shift(18, 1),MATLAB 的集合運(yùn)算符（Set Operators）可將向量視為集合，來進(jìn)行各種集合運(yùn)算,集合運(yùn)算符,x = 1 2 3 4 5 6; y = 1 1 5 5 9 9 9; un= union(x,y) %并集 it= intersect(x,y) %交集 sf= setdiff(x,y) %差集 sr= setxor(x,y) %XOR 運(yùn)算 sq= unique(y) %取相異元素 im= ismember(9,y) %元素9是否屬于集合y,例：,多項(xiàng)式,1） P=a0, a1, , an-1, an表示如下n次多項(xiàng)式:,P(x)=a0 xn+ a1xn-1+

26、 + an-1x+ an,2） P=poly(A),若A為方陣，則P為特征多項(xiàng)式。,若A為一數(shù)組，則P表示以A中全部元為根多項(xiàng)式。,例：A=11 12 13; 14 15 16; 17 18 19,PA=poly(A) PPA=poly2str(PA,x),多項(xiàng)式運(yùn)算,1) P=conv(p1,p2) %多項(xiàng)式乘積,2) q,r=deconv(p1,p2) %多項(xiàng)式相除,3) P=poly(A); dp=polyder(P) %多項(xiàng)式求導(dǎo),4) pf=polyfit(X,Y,n) %多項(xiàng)式擬合,5) PA=polyval(P,s) %將值s代入多項(xiàng)式P,6) PM=polyvalm(P,s)

27、% s是矩陣，矩陣運(yùn)算,s可以是矩陣，數(shù)組運(yùn)算,字符與字符串,本節(jié)介紹MATLAB用于處理字符(Characters)與字符串(Strings)的指令及相關(guān)數(shù)據(jù)型態(tài)。MATLAB 處理字符串的相關(guān)指令大部分都放在下列目錄之中： MATLAB根目錄toolboxmatlabstrfun 其中的”strfun”就是代表”String Functions”.若要查詢與字符和字符串相關(guān)的指令，可在 MATLAB 下輸入： help strfun 或是 help strings,本節(jié)簡介,數(shù)個字符（Characters）可以構(gòu)成一個字符串（Strings） 一個字符串是被視為一個行向量（Row Vect

28、or）進(jìn)行儲存 此一字符串中的每一字符（含空格符），是以其 ASCII 碼的形式存放于此列向量中的每一個元素（Element） Matlab 用“單引號”來界定字符串變量，多個字符串變量可直接合并，以得到一個新字符串變量，例如：,1、字符與字符串的基本概念,t=Anhui University of Technology.,q=School of Mathematics and Physics,tq=q,t,欲輸入含有單引號的字符串，可重復(fù)單引號的使用 若要計(jì)算字符串變量的長度（即組成字符的個數(shù)）可用 length 指令,1、字符與字符串的基本概念（續(xù)）,ste = Ive got a dat

29、e! length(ste)% 計(jì)算字符串變量ste的長度,字符串及其ASCII code,abs和double函數(shù)都可以用來獲取字符串矩陣所對應(yīng)的ASCII碼數(shù)值矩陣。相反，char或setstr函數(shù)可以把ASCII碼矩陣轉(zhuǎn)換為字符串矩陣。例如：,t1=abs(t) %給出字符串t對應(yīng)的ASCII碼值,char(t1) %把ASCII碼矩陣t1轉(zhuǎn)換為字符串,t=Anhui University of Technology:2010,提?。?u=t(5:11) t中第5到11字符子串,倒排：,u=t(end:-1:1),u=t u變?yōu)榱邢蛄?s=find(t=a); end whos x*,x

30、3,x4,x5,x6 都是在for-loop中產(chǎn)生的變量，分別代表維度為33、44、55、66的魔方陣,class 或 ischar 指令: 判斷某一個變量是否為字符串,字符串的判斷,chinese=今日事，今日畢 out1=class(chinese) %out1的值是“char”，代表chinese是字符串變量 x = chinese+1; out2 = ischar(x) %out2的值是0，代表x不是一個字符串變量,2、一個變量來儲存多個字符串,第一種方法是使用二維字符數(shù)組（Two Dimensional Character Arrays） 必須先確認(rèn)每個字符串（即每一行）的長度一樣，

31、否則就必須在短字符串結(jié)尾補(bǔ)上空格符,departments=ee ;cs ;econ %注意空格符的使用,用char 指令儲存多個字符串,departments=char(ee,cs,econ) %注意“()”及“,”的使用,得到結(jié)果和前例一樣;,從二維字符數(shù)組抽取出字符串時(shí)，切記要使用 deblank 指令來移除尾部的空格符,departments = char(ee,cs,econ) dept1 = departments(1,:) % (1,:)代表第一行的元素 dept2 = deblank(dept1) % 使用 deblank 指令來移除尾部的空格符 len1 = length(dept1) % 顯示變量 dept1 的長度=4 len2 = length(dept2) % 顯示變量 dept2 的長度=