向量和矩陣運(yùn)算實(shí)驗(yàn)

向量和矩陣運(yùn)算實(shí)驗(yàn)

向量與矩陣的生成向量與矩陣運(yùn)算

向量的生成

直接輸入:a=[1,2,3,4]

冒號(hào)運(yùn)算符a=[1:4]

==>a=[1,2,3,4]b=[0:pi/3:pi]

==>

b=[0,1.0472,2.0944,3.1416]c=[6:-2:0]

==>c=[6,4,2,0]例：

從矩陣中抽取行或列向量和矩陣運(yùn)算實(shí)驗(yàn)

向量與矩陣的生成（續(xù)）向量與矩陣運(yùn)算

矩陣的生成

直接輸入:A=[1,2,3;4,5,6;7,8,9]

由向量生成

由函數(shù)生成

通過編寫m文件生成例：>>x=[1,2,3];y=[2,3,4];>>A=[x,y],B=[x;y]例：>>C=magic(3)向量和矩陣運(yùn)算實(shí)驗(yàn)常見矩陣生成函數(shù)zeros(m,n)生成一個(gè)m

行n

列的零矩陣，m=n

時(shí)可簡(jiǎn)寫為zeros(n)ones(m,n)生成一個(gè)m行n列的元素全為1的矩陣,

m=n

時(shí)可寫為ones(n)eye(m,n)生成一個(gè)主對(duì)角線全為1的m

行n

列矩陣,

m=n

時(shí)可簡(jiǎn)寫為eye(n)，即為n

維單位矩陣diag(X)若X是矩陣，則diag(X)為X的主對(duì)角線向量若X是向量，diag(X)產(chǎn)生以X為主對(duì)角線的對(duì)角矩陣tril(A)提取一個(gè)矩陣的下三角部分triu(A)提取一個(gè)矩陣的上三角部分rand(m,n)產(chǎn)生0～1間均勻分布的隨機(jī)矩陣m=n

時(shí)簡(jiǎn)寫為rand(n)randn(m,n)產(chǎn)生均值為0，方差為1的標(biāo)準(zhǔn)正態(tài)分布隨機(jī)矩陣

m=n

時(shí)簡(jiǎn)寫為randn(n)向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣操作提取矩陣的部分元素：冒號(hào)運(yùn)算符

A(:)A的所有元素

A(:,:)

二維矩陣A的所有元素

A(:,k)A的第k列，A(k,:)A的第k行

A(k:m)A的第k到第m個(gè)元素

A(:,k:m)A的第k到第m列組成的子矩陣A(:)與A(:,:)的區(qū)別?如何獲得由A的第一、三行和第一、二列組成的子矩陣？自己動(dòng)手向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣操作矩陣的旋轉(zhuǎn)

fliplr(A)

左右旋轉(zhuǎn)

flipud(A)

上下旋轉(zhuǎn)

rot90(A)

逆時(shí)針旋轉(zhuǎn)90度；

rot90(A,k)逆時(shí)針旋轉(zhuǎn)k×90

度例：>>

A=[123;456]>>

B=fliplr(A)>>

C=flipud(A)>>

D=rot90(A),E=rot90(A,-1)向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣操作矩陣的轉(zhuǎn)置與共軛轉(zhuǎn)置

’

共軛轉(zhuǎn)置

.’

轉(zhuǎn)置，矩陣元素不取共軛例：>>

A=[12;2i3i]>>

B=A’>>

C=A.’點(diǎn)與單引號(hào)之間不能有空格！向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣操作改變矩陣的形狀：reshapereshape(A,m,n):將矩陣元素按

列方向

進(jìn)行重組重組后得到的新矩陣的元素個(gè)數(shù)必須與原矩陣元素個(gè)數(shù)相等！向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣操作查看矩陣的大小：size

size(A)

列出矩陣A的行數(shù)和列數(shù)

size(A,1)返回矩陣A的行數(shù)

size(A,2)返回矩陣A的列數(shù)例：>>

A=[123;456]>>

size(A)>>

size(A,1)>>

size(A,2)

length(x)返回向量X的長(zhǎng)度

length(A)等價(jià)于max(size(A))向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣基本運(yùn)算

矩陣的加減：對(duì)應(yīng)分量進(jìn)行運(yùn)算要求參與加減運(yùn)算的矩陣具有相同的維數(shù)例：>>

A=[123;456];B=[321;654]>>

C=A+B;D=A-B;矩陣的普通乘法要求參與運(yùn)算的矩陣滿足線性代數(shù)中矩陣相乘的原則例：>>

A=[123;456];B=[21;34];>>

C=A*B向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣基本運(yùn)算

矩陣的除法：/、\右除和左除

若A可逆方陣，則A\B

<==>A的逆左乘B<==>inv(A)*BB/A

<==>A的逆右乘B<==>B*inv(A)X=A\B

<==>A*X=BX=B/A

<==>X*A=B

通常，矩陣除法可以理解為當(dāng)A和B行數(shù)相等時(shí)即可進(jìn)行左除當(dāng)A和B列數(shù)相等時(shí)即可進(jìn)行右除向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣的乘方

A是方陣，p是正整數(shù)A^p

表示A

的p

次冪，即p

個(gè)A

相乘。

若A是方陣，p不是正整數(shù)

A^p

的計(jì)算涉及到A的特征值分解，即若

A=V*D*V-1

則

A^p=V*(D.^p)/V向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣的乘方若a是標(biāo)量，A是方陣，且[V,D]=eig(A)，則

a^A＝V*(a^D)/V若A,P均是矩陣，則A^P

無定義若a是標(biāo)量，則向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣的Kronecker乘積矩陣Kronecker乘積的定義設(shè)A是n×m矩陣，B是p×q矩陣，則A與B的kronecker乘積為：Kronecker乘積的性質(zhì)是np×mq矩陣；通常任何兩個(gè)矩陣都有Kronecker乘積

Matlab中實(shí)現(xiàn)兩個(gè)矩陣Kronecker相乘的函數(shù)為kron(A,B)Kronecker乘積有時(shí)也稱張量積向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣的數(shù)組運(yùn)算數(shù)組運(yùn)算：對(duì)應(yīng)元素進(jìn)行運(yùn)算點(diǎn)與算術(shù)運(yùn)算符之間不能有空格！數(shù)組運(yùn)算包括：點(diǎn)乘、點(diǎn)除、點(diǎn)冪相應(yīng)的數(shù)組運(yùn)算符為：“.*”，“./”，“.\”和“.^”參與運(yùn)算的對(duì)象必須具有相同的形狀！例：>>

A=[123;456];B=[321;654];>>

C=A.*B;D=A./B;E=A.\B;F=A.^B;向量和矩陣運(yùn)算實(shí)驗(yàn)函數(shù)取值設(shè)x是變量，f

是一個(gè)函數(shù)

當(dāng)x=a是標(biāo)量時(shí)，f(x)=f(a)也是一個(gè)標(biāo)量當(dāng)x=[a,b,…,c]是向量時(shí)，f(x)=[f(a),f(b),…,f(c)]函數(shù)作用在矩陣上的取值若A是矩陣，則f(A)是一個(gè)與A同形狀的矩陣

f

作用在x的每個(gè)分量上向量和矩陣運(yùn)算實(shí)驗(yàn)函數(shù)取值怎樣計(jì)算eA?例：>>

x=[0:pi/4:pi];A=[123;456];>>

y1=sin(x);y2=exp(A);y3=sqrt(A);例：向量和矩陣運(yùn)算實(shí)驗(yàn)矩陣的超越函數(shù)

Matlab提供了三種矩陣函數(shù)：expm、sqrtm、logm詳情參見聯(lián)機(jī)幫助（helpexpm/sqrtm/logm）

更一般的矩陣函數(shù)：funm

funm(A,@fun)參數(shù)fun

的可以是exp,，log，cos，sin，cosh，sinh

向量和矩陣運(yùn)算實(shí)驗(yàn)Matlab中常見數(shù)學(xué)函數(shù)sin、cos、tan、cot、sec、csc、…asin、acos、atan、acot、asec、acsc、…exp、log、log2、log10、sqrtabs、conj、real、imag、signfix、f