应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 信息科学 » 程序算法 » 【求助】MATLAB解矩阵方程
171/2返回列表12下一页
查看: 1363  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

limengiiii

银虫 (小有名气)

[交流] 【求助】MATLAB解矩阵方程

请教一个Matlab关于矩阵方程的问题
小弟近日一直在研究一个矩阵方程的MATLAB求解问题，百思不得其解，还请给位给点指示，问题描述如下：
方程形式为：（A+B*K）Z（A+B*K）'+（C+D*K）Z（C+D*K）'=a*Z
其中，A,B,C,D,K均为已知矩阵，Z为对称矩阵，a为一标量，目的是要求a的值
曾尝试用克罗内克积的形式来变换，但是处理不了对称矩阵的问题，还请各位多多指点！！！

[ Last edited by bslt on 2009-6-23 at 14:33 ]
回复此楼

» 猜你喜欢
1楼 2009-06-22 18:01:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vflag

木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bslt(金币+1,VIP+0):感谢交流。。。 6-29 19:52
Z如果也是变量的话可能不止一组解，设出来Z的上对角线元素，用循环一个一个套吧。
赞一下(11人) 回复此楼
2楼2009-06-22 22:21:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vflag

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
搞定没？搞定再说下哈。呵呵
赞一下(1人) 回复此楼
3楼2009-06-23 09:37:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limengiiii

银虫 (小有名气)
没搞定，的确Z是可以任意解的，关键是求标量a的值，就像解普通矩阵方程似的，AX=b，要求X为非零解，必须要使得A的行列式为零，这个维数较低时可以用手算出来，但维数太高时必须用MATLAB来解答，可惜...
赞一下 回复此楼
4楼2009-06-23 10:01:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limengiiii

银虫 (小有名气)
虫子们来帮个忙啊
赞一下 回复此楼
5楼2009-06-24 20:39:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aimit

木虫 (职业作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bslt(金币+1,VIP+0):感谢交流。。。 6-29 19:52
最后能写成Mz = a*z，其中z为Z的向量化，即z = vec(Z)
因此，a为M的特征值，z为对应的特征向量。
赞一下(11人) 回复此楼
研究方向机器学习及其应用，欢迎交流讨论。
6楼2009-06-27 19:40:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limengiiii

银虫 (小有名气)
是的，但是关键的一步是Z是对称的，所以这个vec（Z）是一个特殊的情况，比如说Z=[1，2；1，1]，vec（Z）=（1，2，1，1）'，但是当Z=[1,2;2,1]时，我们只需令VEC（Z）=[1，2，1]即能表示Z的信息，如何在MATLAB中提取这个信息是我没想明白的，还望赐教！
赞一下 回复此楼
7楼2009-06-28 15:50:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aimit

木虫 (职业作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bslt(金币+1,VIP+0):感谢交流。。。 7-5 09:59
M是已知的矩阵，而且还必须满足Mz = a*z，因此在M的特征向量中，如果由特征向量恢复的矩阵不是对称的，则无解。
赞一下(11人) 回复此楼
研究方向机器学习及其应用，欢迎交流讨论。
8楼2009-06-29 21:48:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limengiiii

银虫 (小有名气)
你说的的确正确，但我想你对我问题的理解有点儿小的偏差，要是不嫌麻烦的话，您可以用手算一下我提出的那个方程，2维的即可，（Z为对称，Z=[Z11,Z12;Z12,Z22]）
赞一下 回复此楼
9楼2009-06-30 09:15:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

string121

金虫 (小有名气)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bslt(金币+2,VIP+0):感谢交流，欢迎常来。 8-7 19:12
bslt(金币+1,VIP+0):赞一个！ 8-7 19:13
矩阵都是dim*dim方阵时Matlab code:

clc;
clear;
dim=10;
A=rand(dim);
B=rand(dim);
C=rand(dim);
D=rand(dim);
K=rand(dim);
S1=A+B*K;
S2=C+D*K;
%z不对称  
for j=1:dim
    for i=1:dim
        for j1=1:dim
            for i1=1:dim
                I=(j-1)*dim+i;
                J=(j1-1)*dim+i1;
                M1(I,J)=S1(i,i1)*S1(j,j1)+ S2(i,i1)*S2(j,j1);
            end
        end
    end
end
%z对称
M2=M1;
col=[];
for j=1:dim
    for i=1:j-1
        M2(:,(i-1)*dim+j)=M1(:,(i-1)*dim+j)+M1(:,(j-1)*dim+i);
        col=[col (j-1)*dim+i];
    end
end
M2(:,col)=[];
M3=M2;
row=[];
for j=1:dim
    for i=1:j-1
        M3((i-1)*dim+j,=(M2((i-1)*dim+j,+M2((j-1)*dim+i,)/2;
        row=[row (j-1)*dim+i];
    end
end
M3(row,=[];
a=eig(M3);
disp('a=')
disp(a)


M3 就是在z对称时7楼所说的的M. a是你要的结果.
没推出向量表示,所以当dim很大时比较慢, 如要做大型问题,
将所有for循环部分改成mex文件.  你算几个特殊的试试对不对
赞一下(21人) 回复此楼
10楼2009-08-04 20:47:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 limengiiii 的主题更新
171/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[基金申请] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 十年后又回来了，论文投稿求助 +5 哈哈114477 2026-08-01 5/250 2026-08-08 15:14 by 仁砚薪传
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +5 a089 2026-08-07 5/250 2026-08-08 14:04 by alian_214
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭