JacobiCN[u,m]
ヤコビ(Jacobi)の楕円関数
を与える.
JacobiCN
JacobiCN[u,m]
ヤコビ(Jacobi)の楕円関数
を与える.
例題
すべて開く すべて閉じる例 (5)
JacobiCN[2.5, 0.3]Plot[JacobiCN[x, 1 / 3], {x, -10, 10}]ComplexPlot3D[JacobiCN[z, 1 / 3], {z, -2 - 2I, 2 + 2I}, PlotLegends -> Automatic]Series[JacobiCN[z, m], {z, 0, 7}]Series[JacobiCN[z, m], {m, 0, 1}]Series[JacobiCN[x, 1], {x, (I π/2), 3}]スコープ (33)
数値評価 (5)
N[JacobiCN[2, 1 / 3], 50]JacobiCN[2, 0.333333333333333333333333333333333333333333]JacobiCN[2.5 + I, 1 - I]JacobiCNを高精度で効率よく評価する:
JacobiCN[2, 0.3`500]//TimingJacobiCN[2, 0.3`10000];//TimingAroundを使って平均的な場合の統計区間を計算する:
JacobiCN[2, Around[1 / 2, 0.01]]JacobiCN[-(I π/3), {{1, 0}, {0, 1}}]MatrixFunctionを使って行列のJacobiCN関数を計算することもできる:
MatrixFunction[JacobiCN[-(I π/3), #]&, {{1, 0}, {0, 1}}]特定の値 (3)
{JacobiCN[z, 0], JacobiCN[z, 1]}{JacobiCN[0, m], JacobiCN[EllipticK[m] / 2, m], JacobiCN[EllipticK[m], m]}JacobiCNのいくつかの極:
{JacobiCN[-3 I EllipticK[1 - m], m], JacobiCN[-I EllipticK[1 - m], m], JacobiCN[I EllipticK[1 - m], m], JacobiCN[3I EllipticK[1 - m], m]}xzero = N@Solve[JacobiCN[x, 1 / 3] == 0 && 0 < x < π, x][[1, 1, 2]]//QuietPlot[JacobiCN[x, 1 / 3], {x, -π, 2π}, Epilog -> Style[Point[{xzero, JacobiCN[xzero, 1 / 3]}], PointSize[Large], Red]]可視化 (3)
JacobiCN関数をパラメータのさまざまな値についてプロットする:
Plot[{JacobiCN[x, 0], JacobiCN[x, 0.5], JacobiCN[x, 1], JacobiCN[x, 1.5], JacobiCN[x, 2]}, {x, -2π, 2π}]JacobiCNをパラメータ
の関数としてプロットする:
Plot[{JacobiCN[1, m], JacobiCN[2, m], JacobiCN[3, m]}, {m, -7, 7}]ComplexContourPlot[Re[JacobiCN[z, 1 / 2]], {z, -4 - 4 I, 4 + 4 I}, IconizedObject[«PlotOptions»]]ComplexContourPlot[Im[JacobiCN[z, 1 / 2]], {z, -4 - 4 I, 4 + 4 I}, IconizedObject[«PlotOptions»]]関数の特性 (8)
JacobiCNは実軸に沿って
の周期を持つ:
JacobiCN[x, m] == JacobiCN[x + 4EllipticK[m], m]JacobiCNは虚軸に沿って
の周期を持つ:
JacobiCN[x, m] == JacobiCN[x + 4I EllipticK[1 - m], m]JacobiCNはその第1引数において偶関数である:
JacobiCN[-x, m] == JacobiCN[x, m]JacobiCNは x の解析関数である:
FunctionAnalytic[JacobiCN[x, m], {x, m}]FunctionSingularities[JacobiCN[x, m], {x, m}]FunctionDiscontinuities[JacobiCN[x, m], {x, m}]FunctionMonotonicity[JacobiCN[x, 1 / 3], x]FunctionInjective[JacobiCN[x, m], x, Assumptions -> m > 1]FunctionInjective[JacobiCN[x, m], x, Assumptions -> m < 1]FunctionInjective[JacobiCN[x, 1], x]Plot[{JacobiCN[x, 1 / 3], .5}, {x, -10, 10}]FunctionSurjective[JacobiCN[x, m], x, Assumptions -> m∈Reals]Plot[{JacobiCN[x, 3], .5}, {x, -10, 10}, PlotRange -> {0, 1.5}]FunctionSign[JacobiCN[x, m], x, Assumptions -> m∈PositiveIntegers]FunctionSign[JacobiCN[x, m], {x, m}]JacobiCNは凸でも凹でもない:
FunctionConvexity[JacobiCN[x, m], {x, m}]微分 (3)
D[JacobiCN[x, m], x]derivs = Table[D[JacobiCN[x, m], {x, n}], {n, 1, 4}]//FullSimplifyPlot[Evaluate[derivs /. m -> 1 / 3], {x, -6, 6}, PlotLegends -> {"First Derivative", "Second Derivative", "Third Derivative", "Fourth Derivative"}]D[JacobiCN[x, m], m]積分 (3)
JacobiCNの不定積分:
Integrate[JacobiCN[z, m], z]Integrate[JacobiCN[z, 1 / 3], {z, -1, 1}]Integrate[JacobiCN[z, 1 / 3], {z, 0, 1}]Integrate[JacobiCN[z, m]JacobiDN[z, m], z]Integrate[(JacobiDN[t, m]/1 + JacobiCN[t, m]), t]級数展開 (3)
Series[JacobiCN[x, 1 / 3], {x, 0, 7}]terms = Normal@Table[Series[JacobiCN[x, 1 / 3], {x, 0, n}], {n, 2, 7, 2}];
Plot[{JacobiCN[x, 1 / 3], terms}, {x, -2, 2}]Series[JacobiCN[1, m], {m, 0, 3}]terms = Normal@Table[Series[JacobiCN[1, m], {m, 0, n}], {n, 1, 3}];
Plot[{JacobiCN[1, m], terms}, {m, -5, 5}]JacobiCNはベキ級数に適用できる:
JacobiCN[ArcTanh[x] + O[x] ^ 8, m]関数の恒等式と簡約 (2)
JacobiCN[-u, m]JacobiCN[u + 4 EllipticK[m] + 4I EllipticK[1 - m], m]JacobiSNを含む恒等式:
JacobiCN[z, m] ^ 2 + JacobiSN[z, m] ^ 2//FullSimplify関数表現 (3)
JacobiAmplitudeのCosによる表現:
Cos[JacobiAmplitude[z, m]]((m - 1) JacobiCD[z, m]^2/m JacobiCD[z, m]^2 - 1)//FullSimplifyJacobiCN[z, m]^2 == (JacobiCS[z, m]^2/1 + JacobiCS[z, m]^2)//FullSimplifyTraditionalFormによる表示:
JacobiCN[z, m]//TraditionalFormアプリケーション (9)
f[x_, t_] := Sqrt[(a/m)](1/m EllipticK[m]) (12 a t EllipticE[m] + (4 a (-2 + m) t + m (x - 6 t β)) EllipticK[m]);u[x_, t_] = 2(a/m)(1 - m + (m β/2a) - ( EllipticE[m]/EllipticK[m]) + m JacobiCN[f[x, t], m] ^ 2);D[u[x, t], t] + 6 u[x, t] D[u[x, t], x] + D[u[x, t], {x, 3}]//FullSimplifyPlot3D[Evaluate[u[x, t] /. {β -> 1, a -> 3, m -> 3 / 4}], {x, -2, 2}, {t, 0, 3}]w[z_] := With[{ζ = (z Gamma[1 / 3] ^ 3/2^1 / 33^1 / 4π), m = Sin[π / 12] ^ 2}, 2^11 / 6Cos[π / 12](1 - JacobiCN[ζ, m])(1 + Tan[π / 12]JacobiCN[ζ, m]) /
(2 3^1 / 4JacobiSN[ζ, m]JacobiDN[ζ, m] + (1 - JacobiCN[ζ, m]) ^ 2)]trianglePoints = With[{t = E ^ (2Pi I / 3)}, Flatten[Table[s / 8. Table[#1 + j / 16(#2 - #1), {j, 0, 16}]&@@@({{1, t}, {t, 1 / t}, {1 / t, 1}} / Sqrt[3]), {s, 8}]]];Row[{Graphics[Point[{Re[#], Im[#]}]& /@ trianglePoints, Frame -> True, ImageSize -> Small],
Graphics[Point[{Re[#], Im[#]}]& /@ w[trianglePoints], Frame -> True, ImageSize -> Small]}, Spacer[20]]ϕ[k_, φ_] = Sqrt[(2k/1 - 2k)]JacobiCN[φ / Sqrt[1 - 2k], k];D[ϕ[k, φ], {φ, 2}] == -ϕ[k, φ] ^ 3 - ϕ[k, φ]//FullSimplifyPlot[Evaluate[Table[ϕ[k, ϕ], {k, 0.1, 0.4, 0.1}]], {ϕ, 0, 6Pi}]ellipse[ϕ_, {a_, b_}] := With[{u = a EllipticE[ϕ, 1 - (b^2/a^2)]}, {a JacobiCN[u, Sqrt[1 - (b^2/a^2)]], b JacobiSN[u, Sqrt[1 - (b^2/a^2)]]}]楕円形のパラメータ化と円形のパラメータ化を使ってプロットする:
With[{b = 1, a = 2}, ParametricPlot[Evaluate[{{a Cos[ϕ], b Sin[ϕ]}, ellipse[ϕ, {a, b}]}], {ϕ, -Pi, Pi}, PlotStyle -> {{StandardGray, Thick}, {White, Dashing[0.1]}}]]Subscript[T, 1][z_] = -Sqrt[m ]JacobiSN[z, m] {{0, 1, 0}, {1, 0, 1}, {0, 1, 0}} / Sqrt[2];
Subscript[T, 2] [z_] = I Sqrt[m] JacobiCN[z, m] {{0, I, 0}, {-I, 0, I}, {0, -I, 0}} / Sqrt[2];
Subscript[T, 3][z_] = I JacobiDN[z, m] {{1, 0, 0}, {0, 0, 0}, {0, 0, -1}};Subsuperscript[T, 1, ′][z] - I (Subscript[T, 2][z]. Subscript[T, 3][z] - Subscript[T, 3][z].Subscript[T, 2][z])Subsuperscript[T, 2, ′][z] - I (Subscript[T, 3][z].Subscript[T, 1][z] - Subscript[T, 1][z]. Subscript[T, 3][z])Subsuperscript[T, 3, ′][z] - I (Subscript[T, 1][z]. Subscript[T, 2][z] - Subscript[T, 2][z].Subscript[T, 1][z])マイラー樹脂製風船(2つの平らなプラスチックシートを外周に沿って縫い合わせ,膨らませたもの)のパラメータ化:
x[u_] := JacobiCN[u, 1 / Sqrt[2]];
y[u_] := Sqrt[2](EllipticE[JacobiSN[u, 1 / Sqrt[2]], 1 / Sqrt[2]] - EllipticF[JacobiSN[u, 1 / Sqrt[2]], 1 / Sqrt[2]] / 2)ParametricPlot3D[{x[u]Cos[v], x[u]Sin[v], y[u]}, {v, 0, 2Pi}, {u, -EllipticK[1 / Sqrt[2]], EllipticK[1 / Sqrt[2]]}, MaxRecursion -> 1]JacobiCNを使って2D単体から3X3の相関行列に点をマッピングする:
Clear[corMat];corMat[{p1_, p2_}, m_] := With[{T = 2EllipticK[m]}, With[{u1 = T p1, u2 = T p2, u3 = T(1 - p1 - p2)},
(| | | |
| --------------- | --------------- | --------------- |
| 1 | JacobiCN[u1, m] | JacobiCN[u2, m] |
| JacobiCN[u1, m] | 1 | JacobiCN[u3, m] |
| JacobiCN[u2, m] | JacobiCN[u3, m] | 1 |)]]Plot3D[Det[corMat[{p1, p2}, 3 / 4]]//Evaluate, {p1, p2}∈Simplex[2], MaxRecursion -> 5]FullSimplify[D[Det[corMat[{Indexed[p, {1}], Indexed[p, {2}]}, m]], {{Indexed[p, {1}], Indexed[p, {2}]}}] /. p -> RegionCentroid[Simplex[2]]]Plot[Det[corMat[{1 / 3, 1 / 3}, m]], {m, 0, 1}]𝒟 = CopulaDistribution[{"MultivariateT", corMat[{0.2, 0.4}, 0.1], 7}, {UniformDistribution[], UniformDistribution[], UniformDistribution[]}]{PointSize[Tiny], Point[RandomVariate[𝒟, 10 ^ 4]]}//Graphics3Dp[t_, m_] := {(JacobiSN[t, m]/1 + JacobiCN[t, m] ^ 2), (JacobiSN[t, m]JacobiCN[t, m]/1 + JacobiCN[t, m] ^ 2)};パラメトリック関数がレムニスケートについて代数方程式を解くことを確認する:
Entity["PlaneCurve", "Lemniscate"][EntityProperty["PlaneCurve", "AlgebraicEquation"]][a][, ]% /. Thread[{, } -> a p[t, m]]//FullSimplify3つの物体の周期的な動き.物体は共通軌道上のレムニスケートで等しい時間間隔で互いを追いかける:
p3[t_, m_] := With[{τ = 4 EllipticK[m]}, {p[t, m], p[t + τ / 3, m], p[t - τ / 3, m]}]楕円パラメータは,3体系の重心が原点に固定されたままであることを要求することで固定される:
m3 = m /. FindRoot[With[{t = 3 / 5EllipticK[m]}, First[Mean[p3[t, m]]]] == 0, {m, 1 / 2}, WorkingPrecision -> 50]Animate[
ParametricPlot[p[t, m3]//Evaluate, {t, 0, 4EllipticK[m3]}, Epilog -> MapIndexed[{ColorData[1]@@#2, PointSize[Large], Point[#1]}&, p3[u, m3]]], {u, 0, 12EllipticK[m3] / 3}, SaveDefinitions -> True]p5[t_, m_] := With[{τ = 4 EllipticK[m]}, Table[p[t + (k/5) τ, m], {k, -2, 2}]]m5 = m /. FindRoot[With[{t = 3 / 5EllipticK[m]}, First[Mean[p5[t, m]]]] == 0, {m, 1 / 2}, WorkingPrecision -> 50]Animate[
ParametricPlot[p[t, m5]//Evaluate, {t, 0, 4 EllipticK[m5]}, Epilog -> MapIndexed[{ColorData[1]@@#2, PointSize[Large], Point[#1]}&, p5[u, m5]]], {u, 0, 4EllipticK[m5]}, SaveDefinitions -> True]5体構成では質量の中心を固定する2つの楕円パラメータが使用できる:
m5a = m /. FindRoot[With[{t = 3 / 5EllipticK[m]}, First[Mean[p5[t, m]]]] == 0, {m, 999 / 1000}, WorkingPrecision -> 50]Animate[
ParametricPlot[p[t, m5a]//Evaluate, {t, 0, 4 EllipticK[m5a]}, Epilog -> MapIndexed[{ColorData[1]@@#2, PointSize[Large], Point[#1]}&, p5[u, m5]]], {u, 0, 8 / m5a ^ 2}, SaveDefinitions -> True]コスタ(Costa)の極小曲面のパラメータ化[MathWorld]:
x[w_, m_] := With[{c = 2EllipticK[m]}, Re[((π/2) + (c^2/4))w - (c/2)(JacobiEpsilon[c w, m] + JacobiCN[c w, m] JacobiDS[c w, m]) - (π/4c)JacobiSD[c w, m](JacobiNC[c w, m] + JacobiCN[c w, m])]];y[w_, m_] := With[{c = 2EllipticK[m]}, Im[((π/2) - (c^2/4))w + (c/2)(JacobiEpsilon[c w, m] + JacobiCN[c w, m] JacobiDS[c w, m]) - (π/4c)JacobiSD[c w, m](JacobiNC[c w, m] + JacobiCN[c w, m])]];z[w_, m_] := With[{c = 2EllipticK[m]}, Sqrt[(π/8)]Log[Abs[JacobiCN[c w, m]^2]]];ParametricPlot3D[With[{w = u + I v}, {x[w, 1 / 2], y[w, 1 / 2], z[w, 1 / 2]}], {u, 0, 1}, {v, 0, 1}, Mesh -> False]特性と関係 (4)
{JacobiCN[InverseJacobiCN[z, m], m], InverseJacobiCN[JacobiCN[z, m], m]}PowerExpandを使って逆関数の多価性を無視する:
PowerExpand[%]CosをJacobiAmplitudeに適用した結果として評価する:
Cos[JacobiAmplitude[z, m]]Solve[JacobiCN[x, m]^2 + 2JacobiCN[x, m] == a, x]//QuietIntegrate[JacobiCN[z, m] ^ 3, z]考えられる問題 (2)
テクニカルノート
関連するガイド
-
▪
- 楕円関数
関連リンク
履歴
1988 で導入 (1.0)
テキスト
Wolfram Research (1988), JacobiCN, Wolfram言語関数, https://reference.wolfram.com/language/ref/JacobiCN.html.
CMS
Wolfram Language. 1988. "JacobiCN." Wolfram Language & System Documentation Center. Wolfram Research. https://reference.wolfram.com/language/ref/JacobiCN.html.
APA
Wolfram Language. (1988). JacobiCN. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/JacobiCN.html
BibTeX
@misc{reference.wolfram_2026_jacobicn, author="Wolfram Research", title="{JacobiCN}", year="1988", howpublished="\url{https://reference.wolfram.com/language/ref/JacobiCN.html}", note=[Accessed: 03-November-2026]}
BibLaTeX
@online{reference.wolfram_2026_jacobicn, organization={Wolfram Research}, title={JacobiCN}, year={1988}, url={https://reference.wolfram.com/language/ref/JacobiCN.html}, note=[Accessed: 03-November-2026]}