求问为什么代码跑不出结果
mathematica吧
全部回复
仅看楼主
level 1
qdcvhjm 楼主
a = 9.6
b = 0.28
epsln = 20
epsl = 1
c = 3
fmhz = 900000000
R = 5
h = 200
u = 1
Dm = 2
T = 5
\[Lambda] = 1
l[r_] := ((epsln*4*\[Pi]*fmhz/c)*(1 -
1/(1 + a*
E^(a*b - b*ArcTan[(h/r)]*180/\[Pi]))) + (epsl*4*\[Pi]*
fmhz/c)*(1/(1 +
a*E^(a*b - b*ArcTan[(h/r)]*180/\[Pi]))))^(-1)*((r^2 +
h^2)^(-1))
p[rr_] := (1/(1 + (T *l[rr])/l[rm])) (1/(1 + (T *l[rr + v])/l[rm + v]))
f[rrr_] := If[0 < rrr < R, 1/(\[Pi] *rrr^2), 0]
s[v_] := 2/Dm^2 \!\(
\*SubsuperscriptBox[\(\[Integral]\), \(0\), \(Dm\)]\(
\*SuperscriptBox[\(E\), \(\((\(-\[Lambda]\))\)*2\ \[Pi] \(
\*SubsuperscriptBox[\(\[Integral]\), \(-10000\), \(10000\)]\((1 -
\*SuperscriptBox[\(E\), \(\(-c\)*\((1 -
\*SubsuperscriptBox[\(\[Integral]\), \(-10000\), \(10000\)]2*\[Pi]*\
p[x + y]\ \ *f[y]*y \[DifferentialD]y)\)\)])\)*
x \[DifferentialD]x\)\)] rm \[DifferentialD]rm\)\)
Plot[s[v], {v, 0, 20}]
2019年04月23日 12点04分 1
level 1
qdcvhjm 楼主
代码一直运行,就是不出图
2019年04月23日 12点04分 2
吧务
level 15
在不确定积分积不积得出来的情况下直接往下编程就已经是自杀行为了,你居然还延迟赋值后扔进了Plot里?这意味着你这个积分就算积得出来(这积分套积分套积分的玩意你确定真有符号解?)Mathematica也会在
每个点处都把它符号计算一遍,不慢才怪。
总之先试试NIntegrate有没有用。
2019年05月04日 06点05分 3
1