資源描述:
《第3章 MATLAB矩陣分析與處理.ppt》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、第3章MATLAB矩陣分析與處理3.1特殊矩陣3.2矩陣結(jié)構(gòu)變換3.3矩陣求逆與線性方程組求解3.4矩陣求值3.5矩陣的特征值與特征向量3.6矩陣的超越函數(shù)特殊矩陣通用的特殊矩陣常用的產(chǎn)生通用特殊矩陣的函數(shù)有:zeros:產(chǎn)生全0矩陣(零矩陣)。ones:產(chǎn)生全1矩陣(幺矩陣)。eye:產(chǎn)生單位矩陣。rand:產(chǎn)生0~1間均勻分布的隨機(jī)矩陣。randn:產(chǎn)生均值為0,方差為1的標(biāo)準(zhǔn)正態(tài)分布隨機(jī)矩陣。例3.1分別建立3×3、3×2和與矩陣A同樣大小的零矩陣。(1)建立一個(gè)3×3零矩陣。zeros(3)(2)建立一個(gè)3×2零矩陣。ze
2、ros(3,2)(3)設(shè)A為2×3矩陣,則可以用zeros(size(A))建立一個(gè)與矩陣A同樣大小零矩陣。A=[123;456];%產(chǎn)生一個(gè)2×3階矩陣Azeros(size(A))%產(chǎn)生一個(gè)與矩陣A同樣大小的零矩陣?yán)?.2建立隨機(jī)矩陣:(1)在區(qū)間[20,50]內(nèi)均勻分布的5階隨機(jī)矩陣。(2)均值為0.6、方差為0.1的5階正態(tài)分布隨機(jī)矩陣。命令如下:x=20+(50-20)*rand(5)%yi=a+(b-a)xiy=0.6+sqrt(0.1)*randn(5)%yi=μ+σxi此外,常用的函數(shù)還有reshape(A,m,n
3、),它在矩陣總元素保持不變的前提下,將矩陣A重新排成m×n的二維矩陣。例3.2建立隨機(jī)矩陣:(1)在區(qū)間[20,50]內(nèi)均勻分布的5階隨機(jī)矩陣。(2)均值為0.6、方差為0.1的5階正態(tài)分布隨機(jī)矩陣。命令如下:x=20+(50-20)*rand(5)y=0.6+sqrt(0.1)*randn(5)此外,常用的函數(shù)還有reshape(A,m,n),它在矩陣總元素保持不變的前提下,將矩陣A重新排成m×n的二維矩陣。用于專門學(xué)科的特殊矩陣(1)魔方矩陣魔方矩陣有一個(gè)有趣的性質(zhì),其每行、每列及兩條對(duì)角線上的元素和都相等。對(duì)于n階魔方陣,
4、其元素由1,2,3,…,n2共n2個(gè)整數(shù)組成。MATLAB提供了求魔方矩陣的函數(shù)magic(n),其功能是生成一個(gè)n階魔方陣。例3.3將101~125等25個(gè)數(shù)填入一個(gè)5行5列的表格中,使其每行每列及對(duì)角線的和均為565。M=100+magic(5)%5階魔方矩陣每行、每列及對(duì)角線的和均為65范得蒙矩陣范得蒙(Vandermonde)矩陣最后一列全為1,倒數(shù)第二列為一個(gè)指定的向量,其他各列是其后列與倒數(shù)第二列的點(diǎn)乘積??梢杂靡粋€(gè)指定向量生成一個(gè)范得蒙矩陣。在MATLAB中,函數(shù)vander(V)生成以向量V為基礎(chǔ)向量的范得蒙矩陣。例
5、如,A=vander([1;2;3;5])即可得到上述范得蒙矩陣。A=[1111;8421;27931;1252551]希爾伯特矩陣在MATLAB中,生成希爾伯特矩陣的函數(shù)是hilb(n),hij=1/(i+j-1)。使用一般方法求逆會(huì)因?yàn)樵紨?shù)據(jù)的微小擾動(dòng)而產(chǎn)生不可靠的計(jì)算結(jié)果。MATLAB中,有一個(gè)專門求希爾伯特矩陣的逆的函數(shù)invhilb(n),其功能是求n階的希爾伯特矩陣的逆矩陣。例3.4求4階希爾伯特矩陣及其逆矩陣。命令如下:formatrat%以有理形式輸出H=hilb(4)H=invhilb(4)托普利茲矩陣托
6、普利茲(Toeplitz)矩陣除第一行第一列外,其他每個(gè)元素都與左上角的元素相同。生成托普利茲矩陣的函數(shù)是toeplitz(x,y),它生成一個(gè)以x為第一列,y為第一行的托普利茲矩陣。這里x,y均為向量,兩者不必等長(zhǎng)。toeplitz(x)用向量x生成一個(gè)對(duì)稱的托普利茲矩陣。例如T=toeplitz(1:6)伴隨矩陣MATLAB生成伴隨矩陣的函數(shù)是compan(p),其中p是一個(gè)多項(xiàng)式的系數(shù)向量,高次冪系數(shù)排在前,低次冪排在后。例如,為了求多項(xiàng)式的x3-7x+6的伴隨矩陣,可使用命令:p=[1,0,-7,6];compan(p)ans
7、=07-6100010帕斯卡矩陣二次項(xiàng)(x+y)n展開后的系數(shù)隨n的增大組成一個(gè)三角形表,稱為楊輝三角形。由楊輝三角形表組成的矩陣稱為帕斯卡(Pascal)矩陣。函數(shù)pascal(n)生成一個(gè)n階帕斯卡矩陣。例3.5求(x+y)5的展開式。在MATLAB命令窗口,輸入命令:pascal(6)矩陣次對(duì)角線上的元素1,5,10,10,5,1即為展開式的系數(shù)。矩陣結(jié)構(gòu)調(diào)整變換對(duì)角陣1.對(duì)角陣只有對(duì)角線上有非0元素的矩陣稱為對(duì)角矩陣,對(duì)角線上的元素相等的對(duì)角矩陣稱為數(shù)量矩陣,對(duì)角線上的元素都為1的對(duì)角矩陣稱為單位矩陣。提取矩陣的對(duì)角線元素
8、設(shè)A為m×n矩陣,diag(A)函數(shù)用于提取矩陣A主對(duì)角線元素,產(chǎn)生一個(gè)具有min(m,n)個(gè)元素的列向量。diag(A)函數(shù)還有一種形式diag(A,k),其功能是提取第k條對(duì)角線的元素。>>A=[123;456;7