Re: [問題] 解非線性聯立方程式(已爬文)已回收
我的程式
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 =((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 =((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*b*D*(lP^(b-a)-uP^(b-a)))+(1-lP)/D-K ;
[uP,lP]=solve(F,G,uP,lP)
Warning: Warning, solutions may have been lost
uP =
.83500855158033480860663378178368
lP =
1.2091851892057315981352098603514
%問題來了，因為答案不是我要的我的uP必大於1
%而且我把此值帶入
subs(F)
ans =
0.2000
subs(G)
ans =
-0.2000
%為什麼我求解不是讓這兩個為0嗎？
%怎麼我把算出來的解帶入竟然會等於0.2
%請大大幫忙看一下吧
%而我把原本的答案做顛倒，發現非常接近0，請問為什麼呢？
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.198.186
※ 編輯: hagry 來自: 140.113.198.186 (11/05 00:11)
※ 編輯: hagry 來自: 140.113.198.186 (11/05 00:13)
討論串 (同標題文章)