24小时热门版块排行榜    

查看: 594  |  回复: 1

freebjx

金虫 (小有名气)


[交流] 【求助】matlab解方程得到一空数组

请大虾帮忙看一下程序,不知道问题出在哪?先谢谢了!

clc
clear
syms wm rm ra k dt dta n0 Gm  Ga  

a11='wm*p21+p12*wm=0';
a21=' -wm*p11-rm*p21+Gm*p31+p22*wm=0';
a31='-k*p31+dt*p41+Ga*p61+p32*wm=0';
a41='Gm*p11-dt*p31-k*p41-Ga*p51+p42*wm=0';
a51='Ga*p41-ra*p51+dta*p61+p52*wm=0';
a61=' -Ga*p31-dta*p51-ra*p61+p62*wm=0';

a12='p22*wm-wm*p11-p12*rm+p13*Gm=0';
a22=' -p12*wm-2*rm*p22+Gm*p32-wm*p21+p23*Gm+n0=0';
a32='-k*p32+dt*p42+Ga*p62-p31*wm-p32*rm+Gm*p33=0';
a42=' Gm*p12-dt*p32-k*p42-Ga*p52-p41*wm-p42*rm+p43*Gm=0';
a52='Ga*p42-ra*p52+dta*p62-p51*wm-p52*rm+p53*Gm=0';
a62=' -Ga*p32-dta*p52-ra*p62-p61*wm-p62*rm+p63*Gm=0';

a13='wm*p23-p13*k+p14*dt+p16*Ga=0';
a23='-wm*p13-rm*p23+Gm*p33-p23*k+p24*dt+p26*Ga=0';
a33='-2*k*p33+dt*p43+Ga*p63+dt*p34+Ga*p36+k';
a43=' p13*Gm-dt*p33-2*k*p43-Ga*p53+dt*p44+Ga*p46=0';
a53='Ga*p43-ra*p53+dta*p63-p53*k+p54*dt+Ga*p56=0';
a63='-Ga*p33-dta*p53-ra*p63-p63*k+p64*dt+Ga*p66=0';

a14='wm*p24+Gm*p11-p13*dt-p14*k-p15*Ga=0';
a24='-wm*p14-rm*p24+Gm*p34+p21*Gm-p23*dt-p24*k-p25*Ga=0';
a34='-2*k*p34+dt*p44+Ga*p64+Gm*p31-dt*p33-Ga*p35=0';
a44='Gm*p14-dt*p34-2*k*p44-Ga*p54+p41*Gm-dt*p43-Ga*p45+k=0';
a54='Ga*p44-ra*p54+dta*p64+p51*Gm-p53*dt-p54*k-Ga*p55';
a64=' -Ga*p34-dta*p54-ra*p64+p61*Gm-p63*dt-p64*k-Ga*p65=0';

a15='wm*p25+p14*Ga-p15*ra+p16*dta=0';
a25='-wm*p15-rm*p25+Gm*p35+p24*Ga-p25*ra+p26*dta=0';
a35=' -k*p35+dt*p45+Ga*p65+Ga*p34-p35*ra+p36*dta=0';
a45=' Gm*p15-dt*p35-k*p45-Ga*p55+Ga*p44-p45*ra+p46*dta=0';
a55='Ga*p45-2*ra*p55+dta*p65+Ga*p54+dta*p56+ra=0';
a65=' -Ga*p35-dta*p55-2*ra*p65+Ga*p64+dta*p66=0';

a16='  wm*p26-p13*Ga-p15*dta-p16*ra=0';
a26=' -wm*p16-rm*p26+Gm*p36-p23*Ga-p25*dta-p26*ra=0';
a36=' -k*p36+dt*p46+Ga*p66-Ga*p33-p35*dta-p36*ra';
a46= 'Gm*p16-dt*p36-k*p46-Ga*p56-Ga*p43-p45*dta-p46*ra=0';
a56='Ga*p46-2*ra*p56+dta*p66-Ga*p53-dta*p55=0';
a66='-Ga*p36-dta*p56-2*ra*p66-Ga*p63-dta*p65+ra=0';
S=solve(a11,a21,a31,a41,a51,a61, a12,a22,a32,a42,a52,a62, a13,a23,a33,a43,a53,a63,  a14,a24,a34,a44,a54,a64, a15,a25,a35,a45,a55,a65, a16,a26,a36,a46,a56,a66, 'p11','p12', 'p13','p14','p15','p16','p21','p22','p23','p24','p25','p26','p31','p32','p33','p34','p35','p36','p41','p42','p43','p44','p45','p46','p51','p52','p53','p54','p55','p56','p61','p62','p63','p64','p65','p66')

运行结果为


S =

[ empty sym ]
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

freebjx

金虫 (小有名气)


robert2020:建议虫友分享解决方法。PS:求解这种符号方程组看着都有点晕呢,建议虫友编程时选取符号变量名称与求解量联系在一起,这样方便自己也方便虫友们查看程序。 2010-12-06 09:35:52
[quote]Originally posted by freebjx at 2010-12-05 17:35:17:
已经解决
2楼2010-12-05 21:15:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 freebjx 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +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
[基金申请] 关于代码变化问题,想知道的进来 +9 且听虎啸 2026-08-07 12/600 2026-08-07 21:53 by 理想三旬2025
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见