Re: [討論]解方程式SOLVE
※ 引述《marco110232 (我要戒煙!!!)》之銘言：
: 問題如下
: 我想解一組方程式
: syms T1 M0
: [M0,T1]=solve('M0*(1-2*exp(-2600/T1)+exp(-3000/T1)*exp(137.6/T1))=873'
: ,'M0*(1-2*exp(-1600/T1)+exp(-3000/T1)*exp(137.6/T1))=555');
: 解出來的值卻是有兩種情況
: 一個是M0=0
: T1=T1
: 另外一個狀況是
: Unable to find closed form solution.
: 想請教我的程式有打錯嗎??
: 還是有高手可以推薦更好的解方程式的指令呢?
: 謝謝指教....感激不盡!!
我想你solve用錯方法了 正確方法我寫在下面
syms T1 M0
eq1 = M0*(1-2*exp(-2600/T1)+exp(-3000/T1)*exp(137.6/T1))-873;
eq2= M0*(1-2*exp(-1600/T1)+exp(-3000/T1)*exp(137.6/T1))-555;
sol = solve(eq1,eq2,T1,M0); % 解出來的sol是結構型態 包含了M0 T1
M = sol.M0;
T = sol.T1;
% M T 就是你的答案
subs(eq1,{M0,T1},{M,T}) % 帶回去驗證答案
subs(eq2,{M0,T1},{M,T}) % 帶回去驗證答案
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 111.251.194.72
推
11/30 14:32, , 1F
11/30 14:32, 1F
推
11/30 14:47, , 2F
11/30 14:47, 2F
→
11/30 14:47, , 3F
11/30 14:47, 3F
推
11/30 14:53, , 4F
11/30 14:53, 4F
→
11/30 14:54, , 5F
11/30 14:54, 5F
→
11/30 14:55, , 6F
11/30 14:55, 6F
推
11/30 15:26, , 7F
11/30 15:26, 7F
→
11/30 18:19, , 8F
11/30 18:19, 8F
討論串 (同標題文章)