作者bennygameii (硷性椅子水)
看板MATLAB
标题[问题] 方程式移项问题
时间Wed Aug 1 16:06:21 2012
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)
1F：→ doom8199:原po你不觉得 f 很多 term 都 repeat 吗 08/02 13:37
2F：→ doom8199:像 f 里头还出现 exp(log( , 这些都是可以在化简的 08/02 13:41
3F：→ doom8199:假设好变数, f 是可以写得很好看的,对 matlab 负担也较低 08/02 13:42
4F：→ bennygameii:恩...我知道有很多项都重复了,可是这个函式其实也是用 08/03 17:02
5F：→ bennygameii:MATLAB跑出来的,所以我要把它写出来化简吗@@? 08/03 17:03
6F：→ bennygameii:还是里面有什麽function可以帮忙的?? 08/03 17:04
7F：→ bennygameii:感谢你~ 08/03 17:04
8F：推 doom8199:程式解不出来通常有两种, 一种是对应不到内建的 pattern 08/05 21:09
9F：→ doom8199:一种是 stack overflow 08/05 21:09
10F：→ doom8199:若原po只是单纯想跑给定 (v_d, m, n), 解 f(v_x_1)=0 08/05 21:11
11F：→ doom8199:可以用 interpolation 近似就好, 不一定非得得到解析解 08/05 21:12
12F：→ bennygameii:我是要跑v_x_1的参数解@@"也就是v_x_1=g(m,n,v_d) 08/08 13:49
13F：→ bennygameii:谢谢喔^^ 我在试看看 08/08 13:49