1樓:魂
用 [d,v] = eig(a) 就可以瞭如:>> a=[1,2;3,4]
a =1 2
3 4
>> [d,v]=eig(a)
d =-4216/5113 -250/601671/1186 -1736/1909v =-736/1977 00 1977/368
v 中是
特徵值, d中是對應內
的特徵向量容
matlab怎麼計算矩陣的特徵值和特徵向量
2樓:天雲一號
在matlab中,可以用eig函式計算矩陣的特徵值和特徵向量。舉例如下:
>> a = [1, 2, 3; 4, 5, 6; 7, 8, 9] % 原始資料矩陣
a =1 2 3
4 5 6
7 8 9
>> [v, d] = eig(a) % 特徵值分解,其中v的每一列表示矩陣a的一個特徵向量,d是一個對角矩陣,對角線上的元素表示矩陣a的特徵值
v =-0.2320 -0.7858 0.
4082-0.5253 -0.0868 -0.
8165-0.8187 0.6123 0.
4082d =16.1168 0 00 -1.1168 00 0 -0.
0000
3樓:我行我素
類似這樣:
a=[....];
[v,d]=eig(a);%v是特徵向量組成的矩陣,d的對角線元素就是特徵值
4樓:今天
使用庫函式eig()
eig: find eigenvalues and eigenvectors(返回矩陣的特徵值和特徵向量; )
[v,d] = eig(a)
d是特徵值
v特徵向量
5樓:匿名使用者
[v, d]=eig(a)
matlab中如何求矩陣的特徵值和特徵向量
6樓:枕風宿雪流年
具體步驟分析如下:
1、第一步我們首先需要知道計算矩陣的特徵值和特徵向量要用eig函式,可以在命令列視窗中輸入help eig,檢視一下eig函式的用法,如下圖所示:
2、第二步在命令列視窗中輸入a=[1 2 3;2 4 5;7 8 9],按回車鍵之後,輸入[x,y]=eig(a),如下圖所示:
3、第三步按回車鍵之後,得到了x,y的值,其中x的每一列值表示矩陣a的一個特徵向量,這裡有3個特徵向量,y的對角元素值代表a矩陣的特徵值,如下圖所示:
4、第四步如果我們要取y的對角元素值,可以使用diag(y),如下圖所示:
5、第五步按回車鍵之後,可以看到已經取出y的對角線元素值,也就是a矩陣的特徵值,如下圖所示:
6、第六步我們也可以在命令列視窗help diag,可以看到關於diag函式的用法,如下圖所示:
7樓:子衿悠你心
可以運用eig函式求特徵值和特徵向量。
e=eig(a):求矩陣a的全部特徵值,構成向量e。
[v,d]=eig(a):求矩陣a的全部特徵值,構成對角陣d,並求a的特徵向量構成v的列向量。
[v,d]=eig(a,'nobalance'):與第2種格式類似,但第2種格式中先對a作相似變換後求矩陣a的特徵值和特徵向量,而格式3直接求矩陣a的特徵值和特徵向量。
例項:求矩陣a=[1,2;2,1]的特徵值和特徵向量。
拓展說明:
在matlab中,還有個函式eigs,可以求特徵向量和特徵值的子集。
d = eigs(a) %求稀疏矩陣a的6個絕對值最大特徵值d,d以向量形式存放。
d = eigs(a,k) %返回k個最大特徵值
8樓:百度使用者
a=[1 1/4;4 1]
a =1.0000 0.2500
4.0000 1.0000
>> [v,d]=eig(a)
v =0.2425 -0.2425
0.9701 0.9701
d =2 0
0 0
按照這道題的計算過程算就可以了,eig是求特徵值和特徵向量命令,v是特徵向量,是列向量,d是特徵值矩陣,主對角線元素就是特徵值,與特徵向量的列對應的
9樓:匿名使用者
[v.d]=eig(a) a為矩陣
在matlab中求矩陣特徵值和特徵向量的**
10樓:大野瘦子
>>clc;clear;close;
>>a=[3,-1,-2;2,0,-2;2,-1,-1];
>>[x,b]=eig(a) %求矩陣a的特徵值和特徵向量,其中b的對角線元素是特徵值,
%x的列是相應的特徵向量
最後的結果是:
x =0.7276 -0.5774 0.6230
0.4851 -0.5774 -0.2417
0.4851 -0.5774 0.7439
b =1.0000 0 0
0 0.0000 0
0 0 1.0000
特徵值和特徵向量的求解根據專案的需求或者是矩陣的具體形式,主要可以分成如下三種形式:
1、只需要獲得矩陣的最大特徵值和特徵值所對應的特徵向量。
2、需要求取矩陣的所有特徵值。
3、需要求取特徵值和特徵向量的矩陣為實對稱矩陣,則可以通過另一種方法進行求解。
這三種形式特徵值和特徵向量的求取:
1.如果自己僅僅要求最大特徵值的話肯定採用形式1的演算法,該演算法的優點是時間複雜度較低,計算量相對較小,該方法不但能夠求取特徵值和特徵向量,而且只要特徵值不全為0,該方法都能獲得想要的結果。
2.如果需要獲得一個矩陣的所有特徵值,則通過形式2可以很好的解決該問題,但是該方法的缺點是僅僅能夠獲得特徵值,獲得特徵值之後利用其它方法進行求解,這樣做自然而然計算量就大了起來。
3.如果矩陣為實對稱矩陣,那麼可以通過形式3對其進行特徵值和特徵向量的求取,該方法相對於形式2的好處就是能夠一次性將特徵值和特徵向量求取出來,缺點就是矩陣必須是實對稱矩陣,至於演算法複雜度方面我沒有進行測試。
11樓:匿名使用者
1、首先開啟自己的電腦,然後在桌面上開啟matlab軟體,進入matlab主介面。
2、然後需要知道計算矩陣的特徵值和特徵向量要用eig函式,可以在該軟體的命令列視窗中輸入help eig,檢視一下eig函式的用法。
3、在該軟體命令列視窗中輸入a=[1 2 3;2 4 5;7 8 9],你按鍵盤上的回車鍵之後,輸入[x,y]=eig(a)。
4、當你按了鍵盤上的回車鍵之後,得到了x,y的值,其中x的每一列值表示矩陣a的一個特徵向量,裡面有3個特徵向量,y的對角元素值代表a矩陣的特徵值。
12樓:
matlab具體**如下:
>>clc;
>>clear;
>>close;
>>a=[3,-1,-2;2,0,-2;2,-1,-1];
>>[x,b]=eig(a)
**中最後一行指的是求矩陣a的特徵值和特徵向量,其中b的對角線元素是特徵值,而x的列即為相應的特徵向量。
13樓:g用事實說話
不明白你說的是什麼意思,我看不懂啊,你能把那個意思詳細說一下嗎?
14樓:匿名使用者
>> a=[3 -1 -2;2 0 -2;2 -1 -1]a =3 -1 -2
2 0 -2
2 -1 -1
>> [v,d]=eig(a)
v =0.7276 -0.5774 0.
62300.4851 -0.5774 -0.
24170.4851 -0.5774 0.
7439d =1.0000 0 00 0.0000 00 0 1.
0000d為特徵值,v為每個特徵值對應
的特徵向量
15樓:匿名使用者
[d,v]=eig(a)
16樓:匿名使用者
a=[3 -1 -2 ;2 0 -2;2 -1 -1];
[u v]=eigs(a)
如何用matlab求矩陣的特徵值及其特徵向量
17樓:射手座
[v,d]=eig(a):求矩陣a的全部特徵值,構成對角陣d,並求a的特徵向量構成v的列向量
a=[1 2 3;4 5 3;5 4 2];
>> [v,d]=eig(a)
v =-0.370937915979909 -0.691545561718756 0.482529180116882
-0.703281764566580 0.083451238716327 -0.753226224556934
-0.606465021346326 0.717496081399613 0.447007656506505
d =9.696761517759944 0 0
0 -2.353923226499950 0
0 0 0.657161708740017
18樓:戊遐思衛詞
用函式[v,d]=eig(a)
矩陣d的對角元儲存的是a的所有特徵值,
而且是從小到大排列的
矩陣v的每一列儲存的是相應的特徵向量
所以應該是v的最後一個列
就是最大特徵值的特徵向量
19樓:朱培勝鈔雨
這有個我們以前的matlab冪法求特徵值和特徵響量的程式:
[maxnorm.m]
function
t=maxnorm(a)
%求數列中按模最大的分量
n=length(a);
t=0;
fori=1:n
ifabs(a(i)/max(abs(a)))>=1t=a(i);
endend
function
[mt,my]=maxtr(a,eps)
%用冪法求矩陣的主特徵值和對應的特徵向量
n=length(a);
x0=diag(ones(n));
k=1x=a*x0
while
norm(x-x0)>eps
k=k+1
q=x;
y=x/maxnorm(x)
x=a*y;
x0=q;
endmt=maxnorm(x)
my=y
[main1.m]
a=[3
2;45]
maxtr(a,0.0001)
[invmaxtr.m]
function
[mx,mt,my]=invmaxtr(a,eps)%求矩陣按模最小的特徵值和對應的特徵向量
n=length(a);
x0=diag(ones(n));
x=inv(a)*x0;
k=0while
norm(x-x0)>eps
k=k+1
q=x;
y=x/maxnorm(x)
x=inv(a)*y;
x0=q;
endmt=1/maxnorm(x)
my=y
[main.m]
a=[3
2;45]
invmaxtr(a)
eigen矩陣matlab怎麼用
專案要進行比較多的矩陣操作,特別是二維矩陣。剛開始做實驗時,使用了動態二維陣列,於 專是寫了一堆屬matrix函式,作矩陣的乘除加減求逆求行列式。實驗做完了,開始做 優化,發現matrix.h檔案裡適用性太低,而且動態二維陣列的空間分配與釋放也影響效率,於是尋找其他解決方案。首先考慮的是與matla...
matlab中如何進行矩陣的特徵分解
比如你的矩陣專 是屬a a 4 7 10 13 5 8 11 14 6 9 12 15 7 10 13 16 u,v eig a u 0.4252 0.7922 0.1848 0.2559 0.4731 0.3667 0.1379 0.0197 0.5211 0.0588 0.8302 0.8072...
這個矩陣的n次方怎麼求? 用特徵值
第一步,求特徵值 第二步,求特徵向量,對應可逆矩陣 具體請看 矩陣的n次方怎麼算?先算兩抄 次方,三次襲方,最多算到4次方,就可bai以知道n次方,du嚴格證明需要用數學zhi歸納法dao。矩陣運算在科學計算中非常重要,而矩陣的基本運算包括矩陣的加法,減法,數乘,轉置,共軛和共軛轉置。矩陣的n次方怎...