level 1
无畏的探险者
楼主
s11 = 0.008245;
s22 = 0.010116;
s33 = 0.007031;
s12 = 0.00000808763;
s13 = -0.00119;
s23 = 0.000464;
s44 = 0.044509;
s55 = 0.051116;
s66 = 0.021527;
Clear[l1, l2, l3, x, y]
l1 = Sin[x]*Cos[y];
l2 = Sin[x]*Sin[y];
l3 = Cos[x];
B = (s11 + s12 + s13)*l1^2 + (s12 + s22 + s23)*
l2^2 + (s13 + s23 + s33)*l3^2;
e = (l1^4*s11 + 2*l1^2*l2^2*s12 + 2*l1^2*l3^2*s13 + l2^4*s22 +
2*l2^2*l3^2*s23 + l3^4*s33 + l2^2*l3^2*s44 + l1^2*l3^2*s55 +
l1^2*l2^2*s66)^(-1);
SphericalPlot3D[B, {x, 0, Pi}, {y, 0, 2*Pi}, PlotRange -> All]
SphericalPlot3D[e, {x, 0, Pi}, {y, 0, 2*Pi}, PlotRange -> All]
L1 = Cos[y];
L2 = Sin[y];
L3 = 0;
eZ = (L1^4*s11 + 2*L1^2*L2^2*s12 + 2*L1^2*L3^2*s13 + L2^4*s22 +
2*L2^2*L3^2*s23 + L3^4*s33 + L2^2*L3^2*s44 + L1^2*L3^2*s55 +
L1^2*L2^2*s66)^(-1);
L1 = Sin[x];
L2 = 0;
L3 = Cos[x];
eY = (L1^4*s11 + 2*L1^2*L2^2*s12 + 2*L1^2*L3^2*s13 + L2^4*s22 +
2*L2^2*L3^2*s23 + L3^4*s33 + L2^2*L3^2*s44 + L1^2*L3^2*s55 +
L1^2*L2^2*s66)^(-1);
L1 = 0;
L2 = Sin[x];
L3 = Cos[x]; eX = (L1^4*s11 + 2*L1^2*L2^2*s12 + 2*L1^2*L3^2*s13 +
L2^4*s22 + 2*L2^2*L3^2*s23 + L3^4*s33 + L2^2*L3^2*s44 +
L1^2*L3^2*s55 + L1^2*L2^2*s66)^(-1); Show[
PolarPlot[eZ, {y, 0, 2 Pi}, PlotStyle -> AbsoluteThickness[2]],
PolarPlot[eY, {x, 0, 2 Pi},
PlotStyle -> {AbsoluteThickness[2], Dashed}],
LabelStyle -> FontSize -> 14, AxesStyle -> Thickness[0.01]]
Show[PolarPlot[eZ, {y, 0, 2 Pi}, PlotStyle -> AbsoluteThickness[2]],
PolarPlot[eX, {x, 0, 2 Pi},
PlotStyle -> {AbsoluteThickness[2], Dashed}],
LabelStyle -> FontSize -> 14, AxesStyle -> Thickness[0.01]]
Show[PolarPlot[eY, {x, 0, 2 Pi}, PlotStyle -> AbsoluteThickness[2]],
PolarPlot[eX, {x, 0, 2 Pi},
PlotStyle -> {AbsoluteThickness[2], Dashed}],
LabelStyle -> FontSize -> 14, AxesStyle -> Thickness[0.01]]

2016年12月27日 13点12分
1
s22 = 0.010116;
s33 = 0.007031;
s12 = 0.00000808763;
s13 = -0.00119;
s23 = 0.000464;
s44 = 0.044509;
s55 = 0.051116;
s66 = 0.021527;
Clear[l1, l2, l3, x, y]
l1 = Sin[x]*Cos[y];
l2 = Sin[x]*Sin[y];
l3 = Cos[x];
B = (s11 + s12 + s13)*l1^2 + (s12 + s22 + s23)*
l2^2 + (s13 + s23 + s33)*l3^2;
e = (l1^4*s11 + 2*l1^2*l2^2*s12 + 2*l1^2*l3^2*s13 + l2^4*s22 +
2*l2^2*l3^2*s23 + l3^4*s33 + l2^2*l3^2*s44 + l1^2*l3^2*s55 +
l1^2*l2^2*s66)^(-1);
SphericalPlot3D[B, {x, 0, Pi}, {y, 0, 2*Pi}, PlotRange -> All]
SphericalPlot3D[e, {x, 0, Pi}, {y, 0, 2*Pi}, PlotRange -> All]
L1 = Cos[y];
L2 = Sin[y];
L3 = 0;
eZ = (L1^4*s11 + 2*L1^2*L2^2*s12 + 2*L1^2*L3^2*s13 + L2^4*s22 +
2*L2^2*L3^2*s23 + L3^4*s33 + L2^2*L3^2*s44 + L1^2*L3^2*s55 +
L1^2*L2^2*s66)^(-1);
L1 = Sin[x];
L2 = 0;
L3 = Cos[x];
eY = (L1^4*s11 + 2*L1^2*L2^2*s12 + 2*L1^2*L3^2*s13 + L2^4*s22 +
2*L2^2*L3^2*s23 + L3^4*s33 + L2^2*L3^2*s44 + L1^2*L3^2*s55 +
L1^2*L2^2*s66)^(-1);
L1 = 0;
L2 = Sin[x];
L3 = Cos[x]; eX = (L1^4*s11 + 2*L1^2*L2^2*s12 + 2*L1^2*L3^2*s13 +
L2^4*s22 + 2*L2^2*L3^2*s23 + L3^4*s33 + L2^2*L3^2*s44 +
L1^2*L3^2*s55 + L1^2*L2^2*s66)^(-1); Show[
PolarPlot[eZ, {y, 0, 2 Pi}, PlotStyle -> AbsoluteThickness[2]],
PolarPlot[eY, {x, 0, 2 Pi},
PlotStyle -> {AbsoluteThickness[2], Dashed}],
LabelStyle -> FontSize -> 14, AxesStyle -> Thickness[0.01]]
Show[PolarPlot[eZ, {y, 0, 2 Pi}, PlotStyle -> AbsoluteThickness[2]],
PolarPlot[eX, {x, 0, 2 Pi},
PlotStyle -> {AbsoluteThickness[2], Dashed}],
LabelStyle -> FontSize -> 14, AxesStyle -> Thickness[0.01]]
Show[PolarPlot[eY, {x, 0, 2 Pi}, PlotStyle -> AbsoluteThickness[2]],
PolarPlot[eX, {x, 0, 2 Pi},
PlotStyle -> {AbsoluteThickness[2], Dashed}],
LabelStyle -> FontSize -> 14, AxesStyle -> Thickness[0.01]]
