24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

啦啦的哇咔咔

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 591.1
散金: 93
帖子: 30
在线: 3.9小时
虫号: 4422914
注册: 2016-02-21
性别: MM
专业: 通信理论与系统

[求助] matlab 矩阵方面的编程求帮忙已有3人参与

a=
   [1    0    0    0    0    1    0    0    0    0;
   0    1    0    0    0    0    1    0    0    0;
   0    0    1    0    0    0    0    1    0    0;
   0    0    0    1    0    0    0    0    1    0;
   0    0    0    0    1    0    0    0    0    1;
   0    0    0    0    0    1    0    0    0    0;
   0    0    0    0    0    0    1    0    0    0;
   0    0    0    0    0    0    0    1    0    0;
   0    0    0    0    0    0    0    0    1    0;
   0    0    0    0    0    0    0    0    0    1];
怎样将这个矩阵变成对角矩阵，不用函数eig(a)，求大家指点帮忙怎样写这个程序，因为我需要知道变成对角矩阵的过程，可能用到过程中的一些计算结果。 @laosam280

回复此楼

» 猜你喜欢

1楼 2016-09-11 16:04:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

guanmm

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 1724
帖子: 152
在线: 40.8小时
虫号: 450384
注册: 2007-11-04
性别: GG
专业: 其他无机非金属材料

不懂，学习了

发自小木虫Android客户端

回复此楼

5楼2016-09-12 16:13:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 10 个回答

FMStation

至尊木虫 (知名作家)

应助: 591 (博士)
贵宾: 0.03
金币: 18252.2
红花: 97
帖子: 8891
在线: 1487.9小时
虫号: 2400059
注册: 2013-04-04
专业: 计算机应用技术

【答案】应助回帖

感谢参与，应助指数 +1

http://www.ece.utah.edu/eceCTools/LinearAlgebra/EigenValVec/Diagonalization/LinAlgEigenDiagMatlabEx.pdf

CODE:

syms a b c d

syms A S D

A = [a, b; c, d]

A =

[ a, b]

[ c, d]

[V,D] = eig(A)

V =

[ -(-1/2*a+1/2*d-1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2))/c, -(- 1/2*a+1/2*d+1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2))/c]

[ 1, 1]

D =

[ 1/2*a+1/2*d+1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2), 0]

[ 0, 1/2*a+1/2*d-1/2*(a^2-2*a*d+d^2+4*b*c)^(1/2)]

Snum = [-1,1;1,0]'

Snum =

 -1 1

 1 0

Dnum = [2, 0; 0, 3]

Dnum =

 2 0

 0 3

Anum = Snum * Dnum * inv(Snum)

Anum =

 3 1

 0 2

[Vnum,Dnum] = eig(Anum)

Vnum =

 1 -0.70711

 0 0.70711

Dnum =

 3 0

 0 2