接主题所指的图。试着在Mathematics6.0上做的两个投影(主要靠用FindRoot在内部不停的解a,费时大)椭圆经线等面积伪圆柱投影(P158):R = 1ParametricPlot[{2 Sqrt[2]*\[Lambda]*R/Pi* Cos[a /. FindRoot[2*a + Sin[2*a] == Pi*Sin[\[CurlyPhi]], {a, -7}]], Sqrt[2]*R* Sin[a /. FindRoot[ 2*a + Sin[2*a] == Pi*Sin[\[CurlyPhi]], {a, -7}]]}, {\[Lambda], -Pi, Pi}, {\[CurlyPhi], -1/2*Pi, 1/2*Pi}]极点投影成线的等面积伪圆柱投影(P157):R = 1ParametricPlot[{2*R*\[Lambda]/Sqrt[Pi + 2]* Cos[1/2*a /. FindRoot[Sin[a] + a == (Pi + 2)/2*Sin[\[CurlyPhi]], {a, -4}]]^2, 2*R/Sqrt[Pi + 2]*a /. FindRoot[ Sin[a] + a == (Pi + 2)/2* Sin[\[CurlyPhi]], {a, -4}]}, {\[Lambda], -Pi, Pi}, {\[CurlyPhi], -Pi/2, Pi/2}]
findroot 是比较慢。等我有空也试验下。
发表评论
2 条评论:
接主题所指的图。
试着在Mathematics6.0上做的两个投影
(主要靠用FindRoot在内部不停的解a,费时大)
椭圆经线等面积伪圆柱投影(P158):
R = 1
ParametricPlot[{2 Sqrt[2]*\[Lambda]*R/Pi*
Cos[a /. FindRoot[2*a + Sin[2*a] == Pi*Sin[\[CurlyPhi]], {a, -7}]],
Sqrt[2]*R*
Sin[a /.
FindRoot[
2*a + Sin[2*a] ==
Pi*Sin[\[CurlyPhi]], {a, -7}]]}, {\[Lambda], -Pi,
Pi}, {\[CurlyPhi], -1/2*Pi, 1/2*Pi}]
极点投影成线的等面积伪圆柱投影(P157):
R = 1
ParametricPlot[{2*R*\[Lambda]/Sqrt[Pi + 2]*
Cos[1/2*a /.
FindRoot[Sin[a] + a == (Pi + 2)/2*Sin[\[CurlyPhi]], {a, -4}]]^2,
2*R/Sqrt[Pi + 2]*a /.
FindRoot[
Sin[a] + a == (Pi + 2)/2*
Sin[\[CurlyPhi]], {a, -4}]}, {\[Lambda], -Pi,
Pi}, {\[CurlyPhi], -Pi/2, Pi/2}]
findroot 是比较慢。
等我有空也试验下。
发表评论