24小时热门版块排行榜    

查看: 1150  |  回复: 9

啦啦的哇咔咔

新虫 (初入文坛)

[求助] 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yogort

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
rref或rrefmovie

» 本帖已获得的红花(最新10朵)

2楼2016-09-11 16:20:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +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

» 本帖已获得的红花(最新10朵)

3楼2016-09-12 10:25:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我怀念的TQ

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
将一个矩阵变成对角阵???是P‘AP=B,B是对角阵,这个意思么?这个矩阵论或线性代数都有啊,你可以查一下

» 本帖已获得的红花(最新10朵)

MindandHand
4楼2016-09-12 15:52:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guanmm

木虫 (小有名气)

5楼2016-09-12 16:13:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

啦啦的哇咔咔

新虫 (初入文坛)

送红花一朵
引用回帖:
4楼: Originally posted by 我怀念的TQ at 2016-09-12 15:52:58
将一个矩阵变成对角阵???是P‘AP=B,B是对角阵,这个意思么?这个矩阵论或线性代数都有啊,你可以查一下

谢谢了
6楼2016-09-19 16:24:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

啦啦的哇咔咔

新虫 (初入文坛)

送红花一朵
引用回帖:
3楼: Originally posted by FMStation at 2016-09-12 10:25:48
http://www.ece.utah.edu/eceCTools/LinearAlgebra/EigenValVec/Diagonalization/LinAlgEigenDiagMatlabEx.pdf


syms a b c d
syms A S D

A =
A =



= eig(A)
V =


D =



Snum = '
Snum ...

谢谢帮忙  我的问题已经解决了
7楼2016-09-19 16:25:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

啦啦的哇咔咔

新虫 (初入文坛)

送红花一朵
引用回帖:
2楼: Originally posted by yogort at 2016-09-11 16:20:36
rref或rrefmovie

谢谢回复  问题已经解决
8楼2016-09-19 16:25:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

啦啦的哇咔咔

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by guanmm at 2016-09-12 16:13:45
不懂,学习了

一起学习
9楼2016-09-19 16:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茶水人生

新虫 (知名作家)

10楼2018-03-06 15:52:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 啦啦的哇咔咔 的主题更新
信息提示
请填处理意见