iewoug 发表于 2008-3-21 09:25

matlab中矩阵方程的求解问题

遇到一个大型矩阵方程求解的问题，请教大家。描述如下：
A(32*32)*U(32*32)+U(32*32)*A(32*32)'=B(32*32),其中U(32*32)中元素为未知量，在我的问题中是协方差矩阵。
如何用matlab求解该协方差矩阵U呢？
谢谢

[ 本帖最后由 eight 于 2008-3-21 10:34 编辑 ]

eight 发表于 2008-3-21 10:34

原帖由 iewoug 于 2008-3-21 09:25 发表 http://www.chinavib.com/forum/images/common/back.gif
遇到一个大型矩阵方程求解的问题，请教大家。描述如下：
A(32*32)*U(32*32)+U(32*32)*A(32*32)'=B(32*32),其中U(32*32)中元素为未知量，在我的问题中是协方差矩阵。
如何用matlab求解该协方差矩阵U呢？
谢谢 搜索一下版面，关于求解方程的问题我不熟悉，不过觉得只要知道 U 是未知量就可以求解了，跟它的物理意义无关

iewoug 发表于 2008-3-21 10:48

回复 2楼 的帖子

版块搜索过了
基本上不是类似问题。或者类似问题，没有下文了 － －

咕噜噜 发表于 2008-3-21 11:01

把A对角化然后计算是不是好一点

iewoug 发表于 2008-3-21 11:08

回复 4楼 的帖子

对角化之后还是需要求解32×32方程。。
用solve命令，需要把方程一个个的列出来。。。不太可能
我之前做个一个维数较少的，可以solve通过符号计算来解决
遇到整个就无策了。。
困扰了 － －

hehy350901 发表于 2008-7-5 15:17

张量分解

可以化成（1024×1）的线性方程组

rushuidetian 发表于 2008-7-7 11:50

化成（1024×1）的线性方程组后，A怎么改写呢？

iewoug 发表于 2008-7-19 17:02

回复 6楼 的帖子

改成1024×1，组成大型的稀疏矩阵，而且对角有零元
稀疏矩阵方程如何处理 0   0
求逆不可行了

Happy99 发表于 2010-11-23 22:49

LZ解决否? 怎么解决?

zhouyang664 发表于 2010-11-23 23:29

啊哈，matlab里有解这类特殊方程的函数：lyap
U = lyap(A, -B) 即可！
GoodLuck！
页: [1]
查看完整版本: matlab中矩阵方程的求解问题