怎么用Matlab求矩阵的特征值和特征向量 怎么用Matlab求矩阵的特征值和特征向量
\u5728MATLAB\u4e2d\u6c42\u77e9\u9635\u7279\u5f81\u503c\u548c\u7279\u5f81\u5411\u91cf\u7684\u4ee3\u7801\u6c42\u77e9\u9635\u7279\u5f81\u503c\u548c\u7279\u5f81\u5411\u91cf\u7684\u65b9\u6cd5\u7279\u5f81\u591a\u9879\u5f0f
eig\u51fd\u6570\u76f4\u63a5\u53ef\u4ee5\u6c42\u7279\u5f81\u503c\u548c\u7279\u5f81\u5411\u91cf
\u5728matlab\u4e2d\uff0c\u8ba1\u7b97\u77e9\u9635a\u7684\u7279\u5f81\u503c\u548c\u7279\u5f81\u5411\u91cf\u7684\u51fd\u6570\u662feig(a)\uff0c\u5e38\u7528\u7684\u8c03\u7528\u683c\u5f0f\u67095\u79cd\uff1a
e=eig(a)\uff1a\u6c42\u77e9\u9635a\u7684\u5168\u90e8\u7279\u5f81\u503c\uff0c\u6784\u6210\u5411\u91cfe\u3002
[v,d]=eig(a)\uff1a\u6c42\u77e9\u9635a\u7684\u5168\u90e8\u7279\u5f81\u503c\uff0c\u6784\u6210\u5bf9\u89d2\u9635d\uff0c\u5e76\u6c42a\u7684\u7279\u5f81\u5411\u91cf\u6784\u6210v\u7684\u5217\u5411\u91cf\u3002
[v,d]=eig(a,'nobalance')\uff1a\u4e0e\u7b2c2\u79cd\u683c\u5f0f\u7c7b\u4f3c\uff0c\u4f46\u7b2c2\u79cd\u683c\u5f0f\u4e2d\u5148\u5bf9a\u4f5c\u76f8\u4f3c\u53d8\u6362\u540e\u6c42\u77e9\u9635a\u7684\u7279\u5f81\u503c\u548c\u7279\u5f81\u5411\u91cf\uff0c\u800c\u683c\u5f0f3\u76f4\u63a5\u6c42\u77e9\u9635a\u7684\u7279\u5f81\u503c\u548c\u7279\u5f81\u5411\u91cf\u3002
e=eig(a,b)\uff1a\u7531eig(a,b)\u8fd4\u56den\u00d7n\u9636\u65b9\u9635a\u548cb\u7684n\u4e2a\u5e7f\u4e49\u7279\u5f81\u503c\uff0c\u6784\u6210\u5411\u91cfe\u3002
[v,d]=eig(a,b)\uff1a\u7531eig(a,b)\u8fd4\u56de\u65b9\u9635a\u548cb\u7684n\u4e2a\u5e7f\u4e49\u7279\u5f81\u503c\uff0c\u6784\u6210n\u00d7n\u9636\u5bf9\u89d2\u9635d\uff0c\u5176\u5bf9\u89d2\u7ebf\u4e0a\u7684n\u4e2a\u5143\u7d20\u5373\u4e3a\u76f8\u5e94\u7684\u5e7f\u4e49\u7279\u5f81\u503c\uff0c\u540c\u65f6\u5c06\u8fd4\u56de\u76f8\u5e94\u7684\u7279\u5f81\u5411\u91cf\u6784\u6210n\u00d7n\u9636\u6ee1\u79e9\u77e9\u9635\uff0c\u4e14\u6ee1\u8db3av=bvd\u3002
具体步骤分析如下：
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函数的用法，如下图所示：
扩展资料：
MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。
MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。
参考资料：百度百科——MATLAB
eig函数直接可以求特征值和特征向量
在MATLAB中，计算矩阵A的特征值和特征向量的函数是eig(A)，常用的调用格式有5种：
E=eig(A)：求矩阵A的全部特征值，构成向量E。
[V,D]=eig(A)：求矩阵A的全部特征值，构成对角阵D，并求A的特征向量构成V的列向量。
[V,D]=eig(A,'nobalance')：与第2种格式类似，但第2种格式中先对A作相似变换后求矩阵A的特征值和特征向量，而格式3直接求矩阵A的特征值和特征向量。
E=eig(A,B)：由eig(A,B)返回N×N阶方阵A和B的N个广义特征值，构成向量E。
[V,D]=eig(A,B)：由eig(A,B)返回方阵A和B的N个广义特征值，构成N×N阶对角阵D，其对角线上的N个元素即为相应的广义特征值，同时将返回相应的特征向量构成N×N阶满秩矩阵，且满足AV=BVD。
绛旓細1..鎵撳紑matlab鐣岄潰锛岃緭鍏鐭╅樀A锛屽A=[1 2 1;1/2 1 1;1 1 1],鐒跺悗鐐笶nter杩涘叆銆2..鍐嶇户缁緭鍏V,D]=eig(A),鐐笶nter杩涘叆 .杩欐椂鍊欏嚭鐜癡=... D=... 寮濮嬩粠涓鎵鐗瑰緛鍚戦噺杩樻湁鏈澶х壒寰佸笺傚浜庤繖涓煩闃 V = 0.7024 0.7024 0.7024 0.4425 -0.2212 + 0.3832i ...
绛旓細eig鍑芥暟鐩存帴鍙互姹傜壒寰鍊煎拰鐗瑰緛鍚戦噺 鍦∕ATLAB涓紝璁＄畻鐭╅樀A鐨勭壒寰鍊煎拰鐗瑰緛鍚戦噺鐨勫嚱鏁版槸eig(A)锛屽父鐢ㄧ殑璋冪敤鏍煎紡鏈5绉嶏細E=eig(A)锛姹傜煩闃A鐨勫叏閮ㄧ壒寰佸硷紝鏋勬垚鍚戦噺E銆俒V,D]=eig(A)锛氭眰鐭╅樀A鐨勫叏閮ㄧ壒寰佸硷紝鏋勬垚瀵硅闃礑锛屽苟姹侫鐨勭壒寰佸悜閲忔瀯鎴怴鐨勫垪鍚戦噺銆俒V,D]=eig(A,'nobalance')锛氫笌绗2绉...
绛旓細5;1/5 1/5 1/3 1 3 2;1/7 1/7 1/5 1/3 1 1/2;1/7 1/7 1/3 1/2 2 1锛鍒╃敤MATLAB杈撳叆 [x,lumda]=eig(A);r=abs(sum(lumda));n=find(r==max(r));max_lumda_A=lumda(n,n)max_x_A=x(:,n)鍗冲彲寰楀埌鐗瑰緛鍊煎拰鐗瑰緛鍚戦噺 鍐嶈緭鍏 w=A/sum(A)鍗冲緱鍒版潈鍚戦噺 ...
绛旓細>>clc;clear;close;>>A=[3,-1,-2;2,0,-2;2,-1,-1];>>[X,B]=eig(A)%姹傜煩闃A鐨勭壒寰鍊煎拰鐗瑰緛鍚戦噺,鍏朵腑B鐨勫瑙掔嚎鍏冪礌鏄壒寰佸,%X鐨勫垪鏄浉搴旂殑鐗瑰緛鍚戦噺鏈鍚庣殑缁撴灉鏄:X=0.7276-0.57740.6230 0.4851-0.5774-0.2417 0.4851-0.57740.7439B=1.000000 00.00000 001.0000 ...
绛旓細鍦澶╀綋鐗╃悊銆侀噺瀛愬姏瀛︾瓑棰嗗煙锛屼篃浼氬嚭鐜版棤绌风淮鐨勭煩闃碉紝鏄煩闃电殑涓绉嶆帹骞裤鐭╅樀鐨勭壒寰鍊间笌鐗瑰緛鍚戦噺 n脳n鐨勬柟鍧楃煩闃礎鐨勪竴涓壒寰佸煎拰瀵瑰簲鐗瑰緛鍚戦噺鏄弧瓒崇殑鏍囬噺浠ュ強闈為浂鍚戦噺 銆傚叾涓璿涓虹壒寰佸悜閲忥紝涓虹壒寰佸笺侫鐨勬墍鏈夌壒寰佸肩殑鍏ㄤ綋锛屽彨鍋欰鐨勮氨 锛岃涓 銆傜煩闃电殑鐗瑰緛鍊煎拰鐗瑰緛鍚戦噺鍙互鎻ず绾挎у彉鎹㈢殑娣卞眰鐗规с
绛旓細MATALB涓殑eig锛圓锛夊嚱鏁扮敤浜姹傜煩闃电殑鐗瑰緛鍊煎拰鐗瑰緛鍚戦噺锛屽父鐢ㄧ殑璋冪敤鏍煎紡鏈変笁绉嶏細1銆丒=eig(A)锛氭眰鐭╅樀A鐨勫叏閮ㄧ壒寰佸硷紝鏋勬垚鍚戦噺E銆2銆乕V,D]=eig(A)锛氭眰鐭╅樀A鐨勫叏閮ㄧ壒寰佸硷紝鏋勬垚瀵硅闃礑锛屽苟姹侫鐨勭壒寰佸悜閲忔瀯鎴怴鐨勫垪鍚戦噺銆3銆乕V,D]=eig(A,'nobalance')锛氫笌绗簩绉嶇被浼硷紝浣嗙浜岀鏍煎紡涓厛瀵笰...
绛旓細[V,D]=eig(A);V鐗瑰緛鍊硷紝D鐗瑰緛鍚戦噺锛泃z=max(D);%鏈澶х壒寰佸 [max_column, index_row] = max(D);%鏈澶х壒寰佸兼墍鍦ㄤ綅缃 a=V(:,index_row(2))瀵瑰簲鐗瑰緛鍚戦噺
绛旓細1.棣栧厛锛屾垜浠缓绔嬩竴涓垜浠渶瑕璁＄畻鐗瑰緛鍊煎拰鐗瑰緛鍚戦噺鐨勬柟闃点2.鐒跺悗灏遍渶瑕佺敤鍒matlab鑷甫鐨勫嚱鏁拌〃杈惧紡鏉ュ疄鐜版柟闃电殑鐗瑰緛鍊煎拰鐗瑰緛鍚戦噺鐨勮绠椾簡銆傛牸寮忓涓嬶細[V,D]=eig(a)銆3.鐒跺悗鎸夊洖杞﹂敭灏卞彲浠ュ緱鍒版垜浠渶瑕佹眰寰楃殑鐭╅樀鐨勭壒寰鍊煎拰鐗瑰緛鍚戦噺浜嗐4.杩欎釜鎴戜滑鍙互鐢ㄧ櫨搴︽悳绱㈠緱鍒帮紝褰撶劧matlab涔熷唴缃湁鍑芥暟鐨勭敤娉曪紝...
绛旓細鍏朵粬鍥炵瓟 鐪媘atlab甯姪鏂囦欢,閲岃竟閮芥湁渚嬬▼,鐢佃剳娌℃湁瀹夎matlab,鎳掑緱瀹夎浜,浣嗘槸鎴戠煡閬撳府鍔╂枃浠堕噷鏈夎繖涓緥瀛 godvmxi | 鍙戝竷浜2010-12-11 涓炬姤| 璇勮 0 0 涓烘偍鎺ㄨ崘: 鐭╅樀鐨勭壒寰鍊 MATLAB绉垎 MATLAB 鈭 鐗瑰緛鍊间笌鐗瑰緛鍚戦噺 鐗瑰緛鏂归樀闃 姹傜壒寰鍊肩殑MATLAB MATLAB姹傚 MATLAB姹A*x=C 鐭╅樀鐗瑰緛鍊肩殑...
绛旓細鐢 [d,v] = eig(A) 灏卞彲浠ヤ簡 濡:>> A=[1,2;3,4]A = 1 2 3 4 >> [d,v]=eig(A)d = -4216/5113 -250/601 671/1186 -1736/1909 v = -736/1977 0 0 1977/368 v 涓槸鐗瑰緛鍊, d涓槸瀵瑰簲鐨勭壒寰鍚戦噺 婊℃剰璇烽噰绾砠_^ ...
