应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 使用contourplot画图，当用Exclusions设置间断点时两侧的等高线错位严重
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
51/1返回列表
查看: 1933  |  回复: 11
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kabaisun

专家顾问 (正式写手)

[求助] 使用contourplot画图，当用Exclusions设置间断点时两侧的等高线错位严重已有1人参与

想画出一隐函数的等值线图，代码如下：
CODE:
a = 0; b = 10; ContourPlot[{y + 3 ArcTan[y/(x - a)]}/20, {x, -b,
  b}, {y, -b, b}, Exclusions -> {x == a}, PlotPoints -> 50,
ContourStyle -> {Red}, Frame -> True, Axes -> True,
ContourShading -> None, PlotLegends -> {"\[Psi](x,y)"}]

由于有间断点，设置Exclusions 指定，图中在此间断点两侧（图中为Y轴）的等高线应该相连才正确，这个错位怎么调整？一直没想出来。。

使用contourplot画图，当用Exclusions设置间断点时两侧的等高线错位严重
平面无旋流叠加图.png
回复此楼

» 猜你喜欢
严实诚勇
1楼2015-01-07 22:00:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xzczd

木虫 (小有名气)
引用回帖:
8楼: Originally posted by kabaisun at 2015-01-08 13:49:46
差别可能只在原点处，arctan(0)是间断点，而arg不是。...

……不，ArcTan值域是-Pi/2到Pi/2，所以你转了半圈就又循环回去了，导致图像错位。Arg不一样，可以转够一圈。所以说你画下图就都明白了：
CODE:
(* 注意值域 *)
Plot3D[ArcTan[y/x], {x, -10, 10}, {y, -10, 10}]

使用contourplot画图，当用Exclusions设置间断点时两侧的等高线错位严重-1
CODE:
(* 注意值域 *)
Plot3D[Arg[x + I y], {x, -10, 10}, {y, -10, 10}]

使用contourplot画图，当用Exclusions设置间断点时两侧的等高线错位严重-2
赞一下 回复此楼
小木虫Mathematica版块已毁（当然原本也不咋的），建议大家前往百度贴吧或Stackexchange。
9楼2015-01-08 14:05:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

xzczd

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
kabaisun: 金币+15, ★★★★★最佳答案 2015-01-08 13:50:39
很遗憾，如果你坚持要用这个方程，那么这就是个不可能的任务，因为ContourPlot忠实地画出了你输入的东西：
CODE:
a = 0;
Limit[y + 3 ArcTan[y/(x - a)], x -> 0, Direction -> #] & /@ {-1, 1}
(* {y + (3 \[Pi] Sqrt[y^2])/(2 y), y - (3 \[Pi] Sqrt[y^2])/(2 y)} *)

如你所见，这里左右极限不等，本来就是间断的。好了我睡觉去了——正当我这么想的时候我突然反应过来了。说的具体一点，就是我明白你其实是想要做什么了。你所需要的其实不是ArcTan，而是Arg（复数幅角）：
CODE:
a = 0; b = 10;
ContourPlot[{y + 3 Arg[(x - a) + I y]}/20, {x, -b, b}, {y, -b, b},
  ContourStyle -> Red, Frame -> True, Axes -> True,
  ContourShading -> None] /.
Graphics[a__] :>
  Legended[Graphics[a],
   Placed[LineLegend[{Red}, {"\[Psi](x,y)"}], After]]

使用contourplot画图，当用Exclusions设置间断点时两侧的等高线错位严重-3

注意现在Exclusions和PlotPoints也用不着了。因为奇点消失了。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

kabaisun
小木虫Mathematica版块已毁（当然原本也不咋的），建议大家前往百度贴吧或Stackexchange。
2楼2015-01-07 22:28:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kabaisun

专家顾问 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by xzczd at 2015-01-07 22:28:07
很遗憾，如果你坚持要用这个方程，那么这就是个不可能的任务，因为ContourPlot忠实地画出了你输入的东西：

a = 0;
Limit & /@ {-1, 1}
(* {y + (3 \ Sqrt)/(2 y), y - (3 \ Sqrt)/(2 y)} *)


如你所见 ...

错位纠正了，还有一个疑问：ArcTan求角度和Arg求的不一样吗？
赞一下 回复此楼
严实诚勇
3楼2015-01-07 23:00:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kabaisun

专家顾问 (正式写手)
引用回帖:
2楼: Originally posted by xzczd at 2015-01-07 22:28:07
很遗憾，如果你坚持要用这个方程，那么这就是个不可能的任务，因为ContourPlot忠实地画出了你输入的东西：

a = 0;
Limit & /@ {-1, 1}
(* {y + (3 \ Sqrt)/(2 y), y - (3 \ Sqrt)/(2 y)} *)


如你所见 ...

如果方程中不止这个幅角函数有间断点，还有对数函数，如下面代码中的方程，这个错位可以调整吗？
CODE:
z1 = -(0.5*Log[x^2 + y^2] - Arg[x + y I])/(2 Pi);
ContourPlot[{z1}/300, {x, -10, 10}, {y, -10, 10},
ContourStyle -> {Red}, Frame -> True, Axes -> False,
ContourShading -> None,
Exclusions -> {x == 0, y == 0}, PlotPoints -> 50]

使用contourplot画图，当用Exclusions设置间断点时两侧的等高线错位严重-4
平面无旋流叠加图.png
赞一下 回复此楼
严实诚勇
4楼2015-01-07 23:31:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭