[問題] 解非線性聯立方程式(已爬文)已回收
我的問題是要解聯立方程式
在前面解的時候都不會有問題,
但是用solve解時,卻出現此方程式非有效式子,
請各位大大幫我看ㄧ下,給我ㄧ下意見,
假如有數值解和逼近解怎麼寫??
我的程式
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;
%uP = sym('uP','real');
%uP = sym('lP','real');
F ='0=((b*uP^a*(uP^(1-a)*lP^(b-a)-lP^(1-a)*uP^(b-a)))-(a*uP^b*(uP^(1-a)-lP^(1-a))))/(a*b*D*(lP^(b-a)-uP^(b-a))+(1-uP )/D+K';
G ='0=((b*lP^a*(uP^(1-a)*lP^(b-a)-lP^(1-a)*uP^(b-a)))-(a*lP^b*(uP^(1-a)-lP^(1-a))))/(a*D*(lP^(b-a)-uP^(b-a))+(1-lP) /D-K ';
%以上在跑的時候都沒問題
[uP,lP]=solve(F,G,'uP,lP')
%但是在解這步的時候卻出現is not a valid expression or equation
%我就換一個fsolve去解,結果他說FSOLVE only accepts inputs of data type double.
我要探討的是在K不同變化下,uP和lP的變化,
前提是要先解出來這個問題,我才能繼續下去,謝謝各位,感謝。
還是說我要去學其他程式??
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.198.186
→
11/04 19:55, , 1F
11/04 19:55, 1F
討論串 (同標題文章)
以下文章回應了本文:
完整討論串 (本文為第 1 之 3 篇):