matlab 解矩阵方程
答案:3 悬赏:50 手机版
解决时间 2021-03-16 08:11
- 提问者网友:感性作祟
- 2021-03-16 03:42
matlab 解矩阵方程
最佳答案
- 五星知识达人网友:时间的尘埃
- 2021-03-16 05:02
这个问题俺也感兴趣,提供一种方法,仅供参考。
clc;clear
P=sym('[p11 p12 p13;p21 p22 p23;p31 p32 p33]')
A=rand(3,3)
At=A'
Q=diag(diag(A))
f=At*P+P*A+Q
结果:
P =
[ p11, p12, p13]
[ p21, p22, p23]
[ p31, p32, p33]
A =
0.3311 0.5676 0.1526
0.6586 0.9805 0.8330
0.8636 0.7918 0.1919
At =
0.3311 0.6586 0.8636
0.5676 0.9805 0.7918
0.1526 0.8330 0.1919
Q =
0.3311 0 0
0 0.9805 0
0 0 0.1919
f =
[ .6622*p11+.6586*p21+.8636*p31+.6586*p12+.8636*p13+.3311, 1.312*p12+.6586*p22+.8636*p32+.5676*p11+.7918*p13, .5230*p13+.6586*p23+.8636*p33+.1526*p11+.8330*p12]
[ .5676*p11+1.312*p21+.7918*p31+.6586*p22+.8636*p23, .5676*p12+1.961*p22+.7918*p32+.5676*p21+.7918*p23+.9805, .5676*p13+1.172*p23+.7918*p33+.1526*p21+.8330*p22]
[ .1526*p11+.8330*p21+.5230*p31+.6586*p32+.8636*p33, .1526*p12+.8330*p22+1.172*p32+.5676*p31+.7918*p33, .1526*p13+.8330*p23+.3837*p33+.1526*p31+.8330*p32+.1919]
然后用solve解9个方程组,求得P,不过,俺觉着,这方法太笨。
clc;clear
P=sym('[p11 p12 p13;p21 p22 p23;p31 p32 p33]')
A=rand(3,3)
At=A'
Q=diag(diag(A))
f=At*P+P*A+Q
结果:
P =
[ p11, p12, p13]
[ p21, p22, p23]
[ p31, p32, p33]
A =
0.3311 0.5676 0.1526
0.6586 0.9805 0.8330
0.8636 0.7918 0.1919
At =
0.3311 0.6586 0.8636
0.5676 0.9805 0.7918
0.1526 0.8330 0.1919
Q =
0.3311 0 0
0 0.9805 0
0 0 0.1919
f =
[ .6622*p11+.6586*p21+.8636*p31+.6586*p12+.8636*p13+.3311, 1.312*p12+.6586*p22+.8636*p32+.5676*p11+.7918*p13, .5230*p13+.6586*p23+.8636*p33+.1526*p11+.8330*p12]
[ .5676*p11+1.312*p21+.7918*p31+.6586*p22+.8636*p23, .5676*p12+1.961*p22+.7918*p32+.5676*p21+.7918*p23+.9805, .5676*p13+1.172*p23+.7918*p33+.1526*p21+.8330*p22]
[ .1526*p11+.8330*p21+.5230*p31+.6586*p32+.8636*p33, .1526*p12+.8330*p22+1.172*p32+.5676*p31+.7918*p33, .1526*p13+.8330*p23+.3837*p33+.1526*p31+.8330*p32+.1919]
然后用solve解9个方程组,求得P,不过,俺觉着,这方法太笨。
全部回答
- 1楼网友:玩世
- 2021-03-16 05:57
我忘了是哪个函数了
不过你自己看看这个
demo
试试
不过你自己看看这个
demo
试试
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
推荐资讯