24小时热门版块排行榜

返回列表

当前主题已经存档。

string121

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 809.5
红花: 3
帖子: 117
在线: 173.2小时
虫号: 810697
注册: 2009-07-17
性别: GG
专业: 计算数学与科学工程计算

的地方是: ) 冒号+右括号

赞一下

回复此楼

11楼2009-08-04 20:51:39

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

string121

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 809.5
红花: 3
帖子: 117
在线: 173.2小时
虫号: 810697
注册: 2009-07-17
性别: GG
专业: 计算数学与科学工程计算

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

的地方是: ) 冒号+右括号

赞一下(1人)

回复此楼

12楼2009-08-04 20:52:34

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

limengiiii

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 465.2
帖子: 101
在线: 6.8小时
虫号: 741527
注册: 2009-04-06
性别: GG
专业: 控制理论与方法

不好意思，好久没看这个帖子，以为不会有结果了，今天试了一下，果然可以，厉害+佩服

赞一下

回复此楼

13楼2009-10-09 15:00:04

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

mao_kun

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 25.3
帖子: 12
在线: 10.9小时
虫号: 829102
注册: 2009-08-16
专业: 控制理论与控制工程

求解

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
suton(金币+1,VIP+0):谢谢交流! 10-17 15:26

这个矩阵方程应不难求解，因为矩阵方程Kronecker乘积后，所求的解应是对称的，求解的MATLAB程序如下：
function [a,Z]=stomatequ(S1,S2)
%(A+B*K)*Z*(A+B*K)'+(C+D*K)*Z*(C+D*K)'=a*Z
% 其中，A,B,C,D,K均为已知矩阵，Z为对称矩阵，a为一标量，目的是要求a的值
% A=rand(dim);B=rand(dim);C=rand(dim);D=rand(dim);
% K=rand(dim);S1=A+B*K;S2=C+D*K;
m=size(S1,2);
A1=kron(S1,S1);
A2=kron(S2,S2);
A3=A1+A2;
[lamdvec,lamd]=eig(A3);
lamd=diag(lamd);
z=lamdvec(:,1);
Z=reshape(z,m,m);
a=lamd(1);
end

其中a和Z是相对应的，这样的解应有dim组（dim为系统的维数）。

赞一下(11人)

回复此楼

14楼2009-10-11 23:19:53

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