| 查看: 608 | 回复: 1 | |||
[交流]
【求助】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 ] |
回复此楼» 猜你喜欢
基金系统什么内容也没有
已经有10人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有15人回复
-
我不理解！
已经有26人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有8人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
云南大学丁军桥团队2027年博士招生
+1/810
-
【征稿中！】2026年机械、制造技术与材料国际会议 (MEMTM 2026)
+1/789
-
[求助]请问美国有哪些PC品牌厂商？
+1/193
-
哈尔滨工业大学（校本部）电信学院2027级博士招生
+1/186
-
2027级博士研究生——质谱分析
+5/160
-
评职称需要，卖发明专利了:D有技术储备需求的设备厂商可联系
+1/19
-
山东大学-山东字节信息智能技术实验室2026年招聘博士后
+2/16
-
哈工大马樱招收2027级计算机类/集成电路类博士生
+1/10
-
2027硕士研究生推免-军事医学研究院生物技术药物代谢研究方向
+1/10
-
北京理工大学机械与车辆学院李伟老师课题组招收2027级硕士研究生（推免生）
+1/6
-
吉林大学高压与超硬材料全国重点实验室李达教授课题组招收2027申请考核制博士
+1/6
-
科研助理招聘
+1/5
-
义翘神州推出全长p-Tau217蛋白，助力阿尔茨海默病下一代生物标志物检测开发
+1/4
-
哈尔滨工程大学动力学院赵建辉团队招收2027级博士研究生
+1/3
-
哈尔滨工程大学动力学院赵建辉团队招收2027级博士研究生
+1/3
-
结肠炎多维度分子解析：炎症‑免疫‑肠道屏障一体化检测策略
+1/3
-
哈尔滨工程大学动力学院赵建辉团队招收2027级博士研究生
+1/2
-
【即将截止】东北大学理学院 力学 李健组招2027级保研
+1/2
-
项目助理招聘启事
+1/2
-
哈尔滨工程大学动力学院赵建辉团队招收2027级博士研究生
+1/2
2楼2010-12-05 21:15:45
