窓鹇 窓鹇
关注数: 56 粉丝数: 18 发帖数: 409 关注贴吧数: 54
关于（遵循Fokker-Planck的）波包演化的问题 在初始时刻将一个波包放在rbs处考察波包随时间的演化 r1 = 0.05; r2 = 5; xx1 = (D[\[Rho][t, rb], rb]) == -(1/TEb) (D[GLb, rb]) \[Rho][t, rb] /. rb -> r1; xx2 = (D[\[Rho][t, rb], rb]) == -(1/TEb) (D[GLb, rb]) \[Rho][t, rb] /. rb -> r2; GLb = (1 + 45 rb^8 + 24 PEb rb^10 - 64 rb^9 TEb)/( 6 rb^4); rbs = 0.6804971882808477`; rbl = 2.780560197330779`; PEb = 0.2379247924446251`; TEb = 0.5`; pdeSoln = NDSolveValue[{D[\[Rho][t, rb], t] == 500 TEb D[E^(-GLb/TEb) D[E^(GLb/TEb) \[Rho][t, rb], rb], rb], \[Rho][0, rb] == 1/(0.1 Sqrt[\[Pi]]) E^(-((rb - rbs)^2/0.01)), xx1, xx2}, \[Rho][t, rb], {t, 0, 20000}, {rb, r1, r2}] Plot3D[pdeSoln, {t, 0, 20000}, {rb, r1, r2}, PlotRange -> {0, 6}, AspectRatio -> 10/13, PlotPoints -> 150, MaxRecursion -> 2, ColorFunction -> "Rainbow"] 给出了一个有问题的结果
为什么点击login后会自动跳转到浏览器的网页上一直登录
萌新询问循环问题（整个循环有660秒） Pb = (-6 Vb^(2/3) + 1 + 8 Tb Vb)/(3 Vb^(4/3)); IPVb = \[Integral]Pb \[DifferentialD]Vb IPVb1 = IPVb /. Vb -> Vb1; IPVb2 = IPVb /. Vb -> Vb2; y1 = IPVb2 - IPVb1 == pb*(Vb2 - Vb1) /. Tb -> tb // Simplify y2 = pb == Pb /. Vb -> Vb1 /. Tb -> tb // Simplify y3 = pb == Pb /. Vb -> Vb2 /. Tb -> tb // Simplify tbList1 = Table[kk, {kk, 0.001, 0.999, (0.999 - 0.001)/200}]; tbList2 = Table[kk, {kk, 0.999, 0.99999, (0.99999 - 0.999)/100}]; tbList = Union[tbList1, tbList2]; LPTb = {{0, 0}}; LVb = {}; LTVb1 = {}; LTVb2 = {}; beginTime = SessionTime[]; For[i = 1, i <= Length[tbList], nsol = NSolve[{y1, y2, y3} /. tb -> tbList[[i]], {Vb1, Vb2, pb}, Reals]; Vb1s = Vb1 /. nsol[[1, 1]]; Vb2s = Vb2 /. nsol[[1, 2]]; pbs = pb /. nsol[[1, 3]]; LPTb = AppendTo[LPTb, {tbList[[i]], pbs}]; If[Vb1s <= Vb2s, LVb = AppendTo[LVb, {Vb1s, Vb2s}]; LTVb1 = AppendTo[LTVb1, {Vb1s, tbList[[i]]}]; LTVb2 = AppendTo[LTVb2, {Vb2s, tbList[[i]]}];, LVb = AppendTo[LVb, {Vb2s, Vb1s}]; LTVb1 = AppendTo[LTVb1, {Vb2s, tbList[[i]]}]; LTVb2 = AppendTo[LTVb2, {Vb1s, tbList[[i]]}]]; i++;] Print["Time used:", IntegerPart[SessionTime[] - beginTime], " second"]; LPTb LVb LTVb1 LTVb2
这位松见真一似乎是宝石之国的演出
1 下一页