[問題] 方程式移項問題

[bennygameii (鹼性椅子水)]

GOOGLE之後發現可以用solve做移項的動作



現在有一個頗為複雜的方程式在下一頁



請各位做好眼花的準備







方程式如下:

f = exp((3012*v_d)/125 - (75917*v_x_1)/2500) + m*(1/exp((3012*v_d)/125 +
(75917*v_x_1)/2500) + exp((3012*v_d)/125 - (75917*v_x_1)/2500) -
exp(log((m*exp((3012*v_d)/125 + (75917*v_x_1)/2500) + m*exp((3012*v_d)/125 +
75917/25000) + m*exp((3012*v_d)/125 - (75917*v_x_1)/2500)*exp((3012*v_d)/125
+ (75917*v_x_1)/2500)*exp((3012*v_d)/125 + 75917/25000) +
m*exp((3012*v_d)/125 + (75917*v_x_1)/2500)*exp((3012*v_d)/125 -
75917/25000)*exp((3012*v_d)/125 +
75917/25000))/(2*m*exp(75917/25000)*exp((3012*v_d)/125)*exp((75917*v_x_1)/2500)
+ 2*m*exp(75917/25000)*exp((9036*v_d)/125)*exp((75917*v_x_1)/2500) +
n*exp(75917/25000)*exp((9036*v_d)/125)*exp((75917*v_x_1)/2500))) -
(3012*v_d)/125) - exp((3012*v_d)/125 + log((m*exp((3012*v_d)/125 +
(75917*v_x_1)/2500) + m*exp((3012*v_d)/125 + 75917/25000) +
m*exp((3012*v_d)/125 - (75917*v_x_1)/2500)*exp((3012*v_d)/125 +
(75917*v_x_1)/2500)*exp((3012*v_d)/125 + 75917/25000) + m*exp((3012*v_d)/125
+ (75917*v_x_1)/2500)*exp((3012*v_d)/125 - 75917/25000)*exp((3012*v_d)/125 +
75917/25000))/(2*m*exp(75917/25000)*exp((3012*v_d)/125)*exp((75917*v_x_1)/2500)
+ 2*m*exp(75917/25000)*exp((9036*v_d)/125)*exp((75917*v_x_1)/2500) +
n*exp(75917/25000)*exp((9036*v_d)/125)*exp((75917*v_x_1)/2500)))))




簡單來說就是一個有( v_x_1, v_d, m, n )四個變數的方程式 f

現在想要找出v_x_1的表示式,故使用solve函式


v_x_1 = solve(f,v_x_1);


卻出現了無解的warning如下:

Warning: Explicit solution could not be found.
> In solve at 160


之前使用都沒有問題


想請問各位有沒有解決的辦法??

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.118.106.75
※ 編輯: bennygameii     來自: 140.118.106.75       (08/01 16:07)
※ 編輯: bennygameii     來自: 140.118.106.75       (08/01 16:07)
→ doom8199:原po你不覺得 f 很多 term 都 repeat 嗎
→ doom8199:像 f 裏頭還出現 exp(log( , 這些都是可以在化簡的
→ doom8199:假設好變數, f 是可以寫得很好看的,對 matlab 負擔也較低
→ bennygameii:恩...我知道有很多項都重複了,可是這個函式其實也是用
→ bennygameii:MATLAB跑出來的,所以我要把它寫出來化簡嗎@@?
→ bennygameii:還是裡面有什麼function可以幫忙的??
→ bennygameii:感謝你~
推 doom8199:程式解不出來通常有兩種, 一種是對應不到內建的 pattern
→ doom8199:一種是 stack overflow
→ doom8199:若原po只是單純想跑給定 (v_d, m, n), 解 f(v_x_1)=0
→ doom8199:可以用 interpolation 近似就好, 不一定非得得到解析解
→ bennygameii:我是要跑v_x_1的參數解@@"也就是v_x_1=g(m,n,v_d)
→ bennygameii:謝謝喔^^ 我在試看看