公告:本站正式转型为非交互式静态网站!
转型:本站将通过笔记和博客的形式继续为大家服务,关于 Mathematica 问答服务请移步至QQ群:365716997
联系:如有问题请联系QQ群管理员,或发送邮件至:lixuan.xyz@qq.com。
感谢:最后非常感谢大家多年来的支持与帮助!
参考《互联网跟帖评论服务管理规定》《中华人民共和国网络安全法》《网络信息内容生态治理规定》《互联网用户账号信息管理规定》

—— 2022-11-27

欢迎来到 Mathematica 问答社区

提问时请贴上文本代码

语法高亮:在编辑器中点击

被禁止的话题:广告破解

请阅读:《提问的智慧》

备用域名:mma.ooo

支持LaTex数学公式
行内公式标识符:\$ 或“$\backslash ($”+“$\backslash )$”,
行间公式标识符:\$\$ 或 “$\backslash [$”+“$\backslash ]$”

社区建议QQ群:365716997

分类

0 投票
2.1k 浏览
g[w_](*电子电场漂移作图*):= 
 Module[{phy = w(*初相位*), \[Omega]t, tra, x, p1, p2, u, result}, 
   trajectory[\[Omega]t_, \[Omega]t0_] := (x[\[Omega]t] /. 
      First@DSolve[
        {x''[\[Omega]t] == -Cos[\[Omega]t], x'[\[Omega]t0] == 0, 
         x[\[Omega]t0] == 0}
        , x, \[Omega]t
        ]);
   result = \[Omega]t /. 
     NSolve[{trajectory[\[Omega]t, phy] == 0, \[Pi]/2 <= \[Omega]t < 
         2 \[Pi]}, \[Omega]t][[1]];
   u = Evaluate@2 (Sin[phy] - Sin[result])^2;
   (*Print[u];*)tra = trajectory[\[Omega]t, phy];
   result = \[Omega]t /. FindRoot[tra == 0, {\[Omega]t, #}] & /@ 
     Range[0, 2 Pi];
   p1 = Plot[-tra, {\[Omega]t, Min@result, Max@result}, 
     PlotStyle -> Blend[{Red, Green, Blue}, RandomReal[]], 
     ImageSize -> 600, PlotRange -> All, 
     Frame -> {{True, False}, {True, False}}, 
     FrameLabel -> {Style["\[Omega]t", FontSize -> 20], 
       Style["x(a.u)", FontSize -> 20]}, 
     LabelStyle -> 
      Directive[Red, Bold, FontSize -> 20, Arrowheads[{0.03}]](*,
     PlotLegends\[Rule]Placed[LineLegend[{"E="<>ToString[u]<>
     "  \[Phi]="<>ToString[phy]},
     LabelStyle\[Rule]{FontSize\[Rule]20}],{Scaled[{1,1}], {1, 
     1}}]*)]] // Quiet
     
     
     
     Table[g[i], {i, 0.1, Pi/2, 0.1}] // Show

我想给边框加上箭头,向坐标轴那样

 

 

或者说怎么将坐标轴的标签放到下面

用户: 落雨流觞 (1.1k 分)
修改于 用户:落雨流觞

1个回答

0 投票
 
已采纳

貌似两个不可兼得,需要变通一下:用Axes画坐标,用Inset显示坐标轴标签。

g[w_](*电子电场漂移作图*):= 
 Module[{phy = w(*初相位*), \[Omega]t, tra, x, p1, p2, u, result}, 
   trajectory[\[Omega]t_, \[Omega]t0_] := (x[\[Omega]t] /. 
      First@DSolve[{x''[\[Omega]t] == -Cos[\[Omega]t], 
         x'[\[Omega]t0] == 0, x[\[Omega]t0] == 0}, x, \[Omega]t]);
   result = \[Omega]t /. 
     NSolve[{trajectory[\[Omega]t, phy] == 0, \[Pi]/2 <= \[Omega]t < 
         2 \[Pi]}, \[Omega]t][[1]];
   u = Evaluate@2 (Sin[phy] - Sin[result])^2;
   (*Print[u];*)tra = trajectory[\[Omega]t, phy];
   result = \[Omega]t /. FindRoot[tra == 0, {\[Omega]t, #}] & /@ 
     Range[0, 2 Pi];
   p1 = Plot[-tra, {\[Omega]t, Min@result, Max@result}, 
     PlotStyle -> Blend[{Red, Green, Blue}, RandomReal[]], 
     ImageSize -> 600, PlotRange -> All,
     AxesStyle -> Directive[Arrowheads[0.05], Red],
     Epilog -> {Inset[
        Style["\[Omega]t", Red, FontSize -> 20], {Center, -0.2}], 
       Inset[Rotate[Style["x(a.u)", Red, FontSize -> 20], 
         90 Degree], {-0.4, Center}]},
     LabelStyle -> Directive[Red, Bold, FontSize -> 20]
     ]
   ] // Quiet

Show[Table[g[i], {i, 0.1, Pi/2, 0.1}], 
 PlotRange -> {{-0.5, 5}, {-0.2, 2}}]

 

用户: 野鹤 (5.1k 分)
采纳于 用户:落雨流觞
谢谢老大,但是有个问题,怎么将横轴稍微延长一些?
可以将PlotRange的范围设置得更宽些
明白了 谢谢坛主
...