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

回复此楼

» 猜你喜欢

最近几年招的学生写论文不引自己组发的文章已经有11人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有4人回复
需要合成515-64-0，50g，能接单的留言已经有3人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复
临港实验室与上科大联培博士招生1名已经有8人回复
想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有4人回复
带资进组求博导收留已经有9人回复

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
金币: 18130.2
红花: 97
帖子: 8891
在线: 1485.3小时
虫号: 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