能不能问一下这FindRoot的警告是什么意思
mathematica吧
全部回复
仅看楼主
level 4
Biasº 楼主
为什么会出现这种情况?这怎么办?就是这个红色的警告
代码
\[CurlyEpsilon]0 = 8.85*10^-12; e = 1.6*10^-19;
kb = 1.380649*10^-23;
\[Mu]c = 1.6*10^-20;
\[CapitalGamma] = 2.857143*10^12;
\[Gamma] = 1.519568*10^13;
T = 300;
h = 1.05*10^-34;
\[CurlyEpsilon]0 = 8.85*10^-12;
fd[\[CurlyEpsilon]_] := (Exp[(\[CurlyEpsilon] - \[Mu]c)/(kb T)] +
1)^-1;
\[Sigma] = (
I e^2 kb T)/(\[Pi] h^2 (\[Omega] + I \[CapitalGamma])) (\[Mu]c/(
kb T) - 2 Log[Exp[-(\[Mu]c/(kb T))] + 1]) + (
I (e^2) (\[Omega] + I \[Gamma]) )/(\[Pi] h^2) \!\(
\*SubsuperscriptBox[\(\[Integral]\), \(0\), \(\[Infinity]\)]\(
\*FractionBox[\(fd[\(-\[CurlyEpsilon]\)] - fd[\[CurlyEpsilon]]\), \(
\*SuperscriptBox[\((\[Omega] + I\ \[Gamma])\), \(2\)] - 4
\*SuperscriptBox[\((\[CurlyEpsilon]/
h)\), \(2\)]\)] \[DifferentialD]\[CurlyEpsilon]\)\);
\[Omega] = 2 Pi f;
\[Kappa]g = (I \[CurlyEpsilon]0 \[Omega] )/\[Sigma];
k = 1.5*10^8; \[CurlyEpsilon]g = 12; \[CurlyEpsilon]s = 4;
FindRoot[k == (\[CurlyEpsilon]g + \[CurlyEpsilon]s) \[Kappa]g, {f,
2.6*10^13}]
2019年07月04日 08点07分 1
吧务
level 15
警告的意思就是字面意思,要是感觉理解不了就试着站在FindRoot的角度上去理解一下——要是这对你来说还太难了的话,那么我可以告诉你,可能导致这一错误的“人类一侧”的原因有很多,其中之一就是这个方程(至少在给定的初值附近)根本没有解。你画画图就能明白了。
2019年07月06日 05点07分 2
1