【求助】ContourPlot画分段函数时输出不完整
mathematica吧
全部回复
仅看楼主
level 3
一剑霜寒44 楼主
电机解析法遇到的，分段函数有点复杂，但是输出的值是对的，用ContourPlot画等A图，图片不完整
程序：
aa[x_, y_] = \!\(\*
TagBox[GridBox[{
{"\[Piecewise]", GridBox[{
{
RowBox[{
RowBox[{
SuperscriptBox["E",
RowBox[{
RowBox[{"-", "I"}], " ",
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}]}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], "-",
FractionBox[
RowBox[{"0.`", " ", "+",
RowBox[{"4.2724486700282686`*^-6", " ", "I"}]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]], "-",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.0034877
13200023075
7`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]}]}], ")"}]}], "+",
RowBox[{
SuperscriptBox["E",
RowBox[{"I", " ",
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}]}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], "+",
FractionBox[
RowBox[{"0.`", " ", "+",
RowBox[{"4.272448670028269`*^-6", " ", "I"}]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]], "+",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.0034877
13200023077
5`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]}]}], ")"}]}], "+",
RowBox[{
SuperscriptBox["E",
RowBox[{
RowBox[{"-", "2"}], " ", "I", " ",
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}]}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], "-",
FractionBox[
RowBox[{"0.`", " ", "+",
RowBox[{"1.8957
13554037898
3`*^-6", " ", "I"}]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "1.`"]], "+",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.348122308914081`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]}], "-",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"1.263282668246829`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "1.`"]}]}], ")"}]}], "+",
RowBox[{
SuperscriptBox["E",
RowBox[{"2", " ", "I", " ",
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}]}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], "+",
FractionBox[
RowBox[{"0.`", " ", "+",
RowBox[{"1.8957
13554037898
`*^-6", " ", "I"}]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "1.`"]], "-",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.348122308914081`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"1.263282668246829`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "1.`"]}]}], ")"}]}]}],
RowBox[{
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"], ">=", "0.035`"}], "&&",
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"], "<=", "0.042`"}]}]},
{
RowBox[{
RowBox[{
SuperscriptBox["E",
RowBox[{
RowBox[{"-", "I"}], " ",
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}]}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], "-",
FractionBox[
RowBox[{"0.`", " ", "+",
RowBox[{"4.312924499533802`*^-6", " ", "I"}]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]], "-",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.0034647677184439776`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]}]}], ")"}]}], "+",
RowBox[{
SuperscriptBox["E",
RowBox[{"I", " ",
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}]}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], "+",
FractionBox[
RowBox[{"0.`", " ", "+",
RowBox[{"4.312924499533802`*^-6", " ", "I"}]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]], "+",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.0034647677184439785`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"]}]}], ")"}]}], "+",
RowBox[{
SuperscriptBox["E",
RowBox[{"2", " ", "I", " ",
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}]}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], "-",
FractionBox[
RowBox[{"0.`", " ", "+",
RowBox[{"0.00001095630834195777`", " ", "I"}]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "1.`"]], "-",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"2.8951
13179815052
5`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "1.`"]}]}], ")"}]}], "+",
RowBox[{
SuperscriptBox["E",
RowBox[{
RowBox[{"-", "2"}], " ", "I", " ",
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}]}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], "+",
FractionBox[
RowBox[{"0.`", " ", "+",
RowBox[{"0.00001095630834195777`", " ", "I"}]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "1.`"]], "+",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"2.8951
13179815052
5`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "1.`"]}]}], ")"}]}]}],
RowBox[{
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"], ">", "0.042`"}], "&&",
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"], "<=", "0.044`"}]}]},
{
TagBox[GridBox[{
{"\[Piecewise]", GridBox[{
{
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "-",
RowBox[{"1.65569
18120084147
`*^-18", " ", "I"}]}], ")"}], "+",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"0.`", " ", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], "-",
FractionBox[
RowBox[{"1.
18875149331
94207`*^-27", "+",
RowBox[{"0.`", " ", "I"}]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "9.`"]], "-",
RowBox[{
RowBox[{"(",
RowBox[{"8.169038059004874`*^19", "+",
RowBox[{"0.`", " ", "I"}]}], ")"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "9.`"]}]}], ")"}], " ",
RowBox[{"Cos", "[",
RowBox[{"18", " ",
RowBox[{"(",
RowBox[{
FractionBox["\[Pi]", "36"], "+",
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}]}], ")"}]}], "]"}]}]}],
RowBox[{
RowBox[{
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}], ">=",
RowBox[{"-",
FractionBox["\[Pi]", "36"]}]}], "&&",
RowBox[{
RowBox[{"ArcTan", "[",
RowBox[{"x", ",", "y"}], "]"}], "<=",
FractionBox["\[Pi]", "36"]}]}]},
{"0",
TagBox["True",
"PiecewiseDefault",
AutoDelete->True]}
},
AllowedDimensions->{2, Automatic},
Editable->True,
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxItemSize->{
"Columns" -> {{Automatic}}, "ColumnsIndexed" -> {},
"Rows" -> {{1.}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.84]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}},
Selectable->True]}
},
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxItemSize->{
"Columns" -> {{Automatic}}, "ColumnsIndexed" -> {},
"Rows" -> {{1.}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.35]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}],
"Piecewise",
DeleteWithContents->True,
Editable->False,
SelectWithContents->True,
Selectable->False],
RowBox[{
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"], ">", "0.044`"}], "&&",
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], "0.5`"], "<=", "0.05`"}]}]},
{"0",
TagBox["True",
"PiecewiseDefault",
AutoDelete->True]}
},
AllowedDimensions->{2, Automatic},
Editable->True,
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxItemSize->{
"Columns" -> {{Automatic}}, "ColumnsIndexed" -> {},
"Rows" -> {{1.}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.84]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}},
Selectable->True]}
},
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxItemSize->{
"Columns" -> {{Automatic}}, "ColumnsIndexed" -> {},
"Rows" -> {{1.}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.35]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}],
"Piecewise",
DeleteWithContents->True,
Editable->False,
SelectWithContents->True,
Selectable->False]\);
ContourPlot[Chop@Re[aa[x, y]], {x, -0.05, 0.05}, {y, 0, 0.05},
MaxRecursion -> 0, Contours -> 30, AspectRatio -> Automatic,
PlotRange -> All]
先输出数据后画图就没有问题：
data = Table[
Chop@Re[aa[x, y]], {x, -0.05, 0.05, 0.001}, {y, 0, 0.05, 0.001}];
ListContourPlot[data, PlotRange -> All, Contours -> 30,
AspectRatio -> Automatic]
求助各位大佬
2019年04月16日 07点04分 1
吧务
level 15
xzcyr
应该是跟舍入误差有关。改成
Rationalize[Chop@Re[aa[x, y]], 0]
再把
MaxRecursion -> 0
去掉就正常了。
2019年05月04日 07点05分 2
1