Wolfram mathematica 利用条件求虚部和实部
我用Solve在Mathematica中找到了一个方程的解（我在这里发布的原因是没有人能在Mathematica stack中回答我的问题）。这个解叫做s，它是两个变量v和ro的函数。我想找到s的虚部和实部，我想利用v和ro是实部的信息，它们在下面的区间内： $0.02Wolfram mathematica 利用条件求虚部和实部,wolfram-mathematica,mathematica-8,Wolfram Mathematica,Mathematica 8,我用Solve在Mathematica中找到了一个方程的解（我在这里发布的原因是没有人能在Mathematica stack中回答我的问题）。这个解叫做s，它是两个变量v和ro的函数。我想找到s的虚部和实部，我想利用v和ro是实部的信息，它们在下面的区间内： $0.02
ClearAll["Global`*"]
d = 1; l = 100; k = 0.001; kk = 0.001;ke = 0.0014;dd = 0.5 ; dr = 0.06; dc = 1000; p = Sqrt[8 (ro l /2 - 1)]/l^2;
m = (4 dr + ke^2 (d + dd)/2) (-k^2 + kk^2) (1 - l ro/2) (d - dd)/4 -
I v p k l (4 dr + ke^2 (d + dd)/2)/4 - v^2 ke^2/4 + I v k dr l p/4;
xr = 0.06/n;
tr = d/n;
dp = (x (v I kk/2 (4 dr + ke^2 (d + dd)/2) - I v kk ke^2 (d - dd)/8 - dr l p k kk (d - dd)/4) + y ((xr I kk (ro - 1/l) (4 dr + ke^2 (d + dd)/2)) - I v kk tr ke^2 (1/l - ro/2) + I dr xr 4 kk (1/l - ro/2)))/m;
a = -I v k dp/4 - I xr y kk p/2 + l ke^2 dp p (d + dd)/8 + (-d + dd)/4 k kk x + dr l p dp;
aa = -v I kk dp/4 + xr I y k p/2 - tr y ke^2 (1/l - ro/2) - (d - dd) x kk^2/4 + ke^2 x (d - dd)/8;
ca = CoefficientArrays[{x (s + ke^2 (d + dd)/2) +
dp (v I kk - l (d - dd) k p kk/2) + y (tr ro ke^2) - (d -
dd) ((-kk^2 + k^2) aa - 2 k kk a)/(4 dr + ke^2 (d + dd)/2) == 0, y (s + dc ke^2) + n x == 0}, {x, y}];
mat = Normal[ca];
matt = Last@mat;
sha = Solve[Det[matt] == 0, s];
shaa = Assuming[v < 100 && v > 40 && ro < 1 && ro > 0.03,Simplify[%]];
reals = Re[shaa];
ims = Im[shaa];
Solve[reals == 0, ro]
ClearAll[“全局”*]
d=1；l=100；k=0.001；kk=0.001；ke=0.0014；dd=0.5；dr=0.06；dc=1000；p=Sqrt[8（ro l/2-1）]/l^2；
m=（4dr+ke^2（d+dd）/2（-k^2+kk^2）（1-lro/2）（d-dd）/4-
I v p k l（4 dr+ke^2（d+dd）/2）/4-v^2 ke^2/4+I v k dr l p/4；
xr=0.06/n；
tr=d/n；
dp=（x（v I kk/2（4 dr+ke^2（d+dd）/2）-I v kk^2（d-dd）/8-dr l p k kk（d-dd）/4）+y（（xr I kk（ro-1/l）（4 dr+ke^2（d+dd）/2））-I v kk tr ke^2（1/l-ro/2）+I dr xr 4 kk（1/l-ro/2））/m；
a=-I v k dp/4-I xr y kk p/2+l ke^2 dp p（d+dd）/8+（-d+dd）/4 k kk x+dr l p dp；
aa=-vikkdp/4+xriykp/2-tryke^2（1/l-ro/2）-（d-dd）xkk^2/4+ke^2x（d-dd）/8；
ca=系数[{x（s+ke^2（d+dd）/2）+
dp（v I kk-l（d-dd）k p kk/2）+y（tr ro ke^2）-（d-
dd）（（-kk^2+k^2）aa-2kkka）/（4dr+ke^2（d+dd）/2）==0，y（s+dc-ke^2）+nx==0}，{x，y}）；
mat=正常[ca]；
马特=Last@mat;
sha=求解[Det[matt]==0，s]；
shaa=假设[v<100&&v>40&&ro<1&&ro>0.03，简化[%]；
reals=Re[shaa]；
ims=Im[shaa]；
求解[reals==0，ro]
但它没有给出答案。有人能帮忙吗？我真的很感激任何解决这个问题的方法。我一直运行您的代码
mat = Normal[ca]
sha = Solve[Det[matt] == 0, s]
mat = Normal[ca]
sha = Solve[Det[matt] == 0, s]
sha = s/.Solve[Det[matt] == 0, s]
shaa = Assuming[40<v<100 && .03<ro<1, Simplify[sha]];
reals = Re[shaa]
Arg[-1. + 50. ro]
reals=reals/.Arg[-1. + 50. ro]->0
Solve[reals == 0, ro]