Question

求助mathematica问题

f = 1/8 E^(-2 r) (-2 + (-1 + E^(4 r))/(I - 2 l)^2 - (
4 E^(4 r))/(-I + l)^2 - (
4 E^(4 r))/(I + l)^2 + (-1 + E^(4 r))/(I + 2 l)^2 + (
2 T^2 HurwitzZeta[2, T/c])/c^2 +
6 E^(4 r) (-1 + (T^2 HurwitzZeta[2, T/c])/c^2) + (
2 (1 + 3 E^(4 r)) T^2 PolyGamma[1, T/c])/c^2 - (
8 E^(4 r) T^2 PolyGamma[1, (T - I l T)/c])/c^2 - (
8 E^(4 r) T^2 PolyGamma[1, (T + I l T)/c])/c^2 + (
2 (-1 + E^(4 r)) T^2 PolyGamma[1, (T - 2 I l T)/c])/c^2 + (
2 (-1 + E^(4 r)) T^2 PolyGamma[1, (T + 2 I l T)/c])/c^2)
T = 0.8 c
ct = l
a = 0.01
r = 0.31
Plot3D[2 E^(-4 a*f)*Abs[\[Alpha]*Sqrt[1 - \[Alpha]^2]], {\[Alpha], 0,
1}, {l, 0, 20},
AxesLabel -> {\[Alpha]^2,
"\!\(\*SubscriptBox[\"\[Omega]\", \"c\"]\)t",
"\!\(\*SubscriptBox[\"C\", \"\[CapitalPsi]\"]\)(t)"},
LabelStyle -> Directive[Bold], AxesStyle -> Directive[14]]为什么这个有虚数就能出图？那么下面的就出不了？r = 0.31
\[Omega]*t*g = l

h = (2 g^2)/\[Omega]^2 (1 -
Cos[\[Omega]*t])*((-3 E^(-I*\[Omega]*t) + Cos[\[Omega]*t]) Sinh[
r] Cosh[r] + Sinh[r] Sinh[r] + Cosh[r] Cosh[r])
Plot3D[2 E^h*Abs[\[Alpha]*Sqrt[1 - \[Alpha]^2]], {\[Alpha], 0, 1}, {l,
0, 2000},
AxesLabel -> {\[Alpha]^2, "\!\(\*SubscriptBox[\"g\", \"0\"]\)t",
"\!\(\*SubscriptBox[\"C\", \"\[CapitalPsi]\"]\)(t)"},
LabelStyle -> Directive[Bold], AxesStyle -> Directive[14]]

Answer 1

你第二个方程 w×t×g=l的赋值语句不对。

估计你想问的是这个，Plot[Sin[x] + 0.1^20 I, {x, 0, 6 Pi}]，有虚数但是能画出图，

Plot[Sin[x] + 0.1 I, {x, 0, 6 Pi}]，但是这个就不可以。

这是由于，在计算时候，有些小虚数是多次计算产生的误差，在画图时候软件自动处理为实数。而0.1 I软件并不认为是误差，一个虚数是不能简单画在二维平面图中的。

追问

请问，那系统是自动的不考虑虚数那部分了还是用一些其他的东西代替呢？如果是不考虑虚数了的情况，那是不是就是前面直接取一个Re就行了呢？

追答

可以试试，在实数域上求解。或者，能确定所求是实数的话，就re个把。。