[問題] 解方程式的問題--次方向非整數已回收
我是matlab新手，
我想問各位大大，最近在研究paper裡面的圖，
要試著模擬出來跟上面相似的圖，
但我現在卡在程式的問題，
我想解出兩個滿醜的聯立方程式，但是有時候會出現：
此式子非有效式子
詳細的程式碼如下：
K=0.1;
D=0.09;
S=0.2;
a=0.5+sqrt(0.25+2*D/S^2);
b=0.5-sqrt(0.25+2*D/S^2);
syms uP lP;
F =
'0=(uP^a*(uP^(1-a)*lP^(b-a)-lP^(1-a)*uP^(b-a)))/(a*D*(lP^(b-a)-uP^(b-a))+
1/D+K-(uP^b*(uP^(1-a)-lP^(1-a)))/(b*D*(lP^(b-a)-uP^(b-a))-uP/D';
G =
'0=(lP^a*(uP^(1-a)*lP^(b-a)-lP^(1-a)*uP^(b-a)))/(a*D*(lP^(b-a)-uP^(b-a))+
1/D-K-(lP^b*(uP^(1-a)-lP^(1-a)))/(b*D*(lP^(b-a)-uP^(b-a))-lP/D';
[uP,lP]=solve(F,G,'uP,lP')
希望各位大大能幫忙解釋一下，還是matlab沒辦法解，
我主要是要畫圖，連式子都沒辦法解，我就無法探討上面變數那些變數變動下的圖形了，
拜託各位了。
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.4.234
推
11/03 20:42, , 1F
11/03 20:42, 1F
※ 編輯: hagry 來自: 123.194.80.243 (11/03 22:32)
→
11/03 22:35, , 2F
11/03 22:35, 2F
→
11/03 22:39, , 3F
11/03 22:39, 3F
→
11/04 00:02, , 4F
11/04 00:02, 4F