3

アイゼンシュタイン級数の収束が遅すぎる!

379
0
$$$$

んで高速化しようというこのプロジェクトの結論。
まず定義。

アイゼンシュタイン級数

$$G_{2k}(\tau) := \sum_{(m,n)\in \mathbb{Z}^2 \backslash (0,0)} \frac{1}{(m\tau + n)^{2k}}$$

で、高速版。(ついでに上記は上半平面上でしか収束しないのにこっち実軸以外全平面で収束するバグ)

$$G_{2k}(\tau)=2\zeta(2k)-\frac{2\pi^{2k}}{(2k-1)!}\sum^{\infty}_{n=1} \left( \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z)\right|_{z=\pi n \tau} \right)$$

よーし頑張ってこー

$$ \begin{align} &G_{2k}(\tau) := \sum_{(m,n)\in \mathbb{Z}^2 \backslash (0,0)} \frac{1}{(m\tau + n)^{2k}} \\ &\qquad\quad \ = 2\sum^{\infty}_{m=1} \sum^{\infty}_{n=1} \left( \frac{1}{(m\tau + n)^{2k}} + \frac{1}{(m\tau - n)^{2k}} \right)+2\sum^{\infty}_{m=1}\frac{1}{(m\tau)^{2k}}+2\sum^{\infty}_{n=1}\frac{1}{n^{2k}} \\ &\qquad\quad \ = 2\frac{1}{(2k-1)!} \sum^{\infty}_{m=1} \left( \psi^{(2k-1)}(1+m\tau) + \psi^{(2k-1)}(1-m\tau) \right)+2\frac{1}{\tau^{2k}}\zeta(2k)+2\zeta(2k) \\ &\qquad\quad \ = 2\frac{1}{(2k-1)!} \sum^{\infty}_{m=1} \left( \psi^{(2k-1)}(1+m\tau) - \psi^{(2k-1)}(m\tau) - \pi^{2k} \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z) \right|_{z=\pi m \tau} \right)+2\frac{1}{\tau^{2k}}\zeta(2k)+2\zeta(2k) \\ &\qquad\quad \ = 2\frac{1}{(2k-1)!} \sum^{\infty}_{m=1} \left( -\frac{(2k-1)!}{(\tau m)^{2k}} - \pi^{2k} \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z) \right|_{z=\pi m \tau} \right)+2\frac{1}{\tau^{2k}}\zeta(2k)+2\zeta(2k) \\ &\qquad\quad \ = - 2\frac{1}{\tau^{2k}}\zeta(2k)- \frac{2\pi^{2k}}{(2k-1)!} \sum^{\infty}_{m=1} \left( \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z) \right|_{z=\pi m \tau} \right)+2\frac{1}{\tau^{2k}}\zeta(2k)+2\zeta(2k) \\ &\qquad\quad \ = 2\zeta(2k) - \frac{2\pi^{2k}}{(2k-1)!} \sum^{\infty}_{m=1} \left( \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z) \right|_{z=\pi m \tau} \right) \end{align} $$

ほぼゼータの定義とディガンマ関数の級数展開、相反公式をただぶっこむだけだね。
あってるよな...?(自分を信じたら負けだと思ってる)
DomeainColoringでもするのだ。とりあえず$k=2,3$でやろうではないか。
$$G_4(\tau) = \frac{\pi^4}{45}-\frac{2\pi^4}{3}\sum^{\infty}_{n=1}\left( \frac{2}{\sin^2(\pi n\tau)} - \frac{3}{\sin^4(\pi n\tau)} \right)$$
$$G_6(\tau) = \frac{2\pi^6}{945}+\frac{2\pi^6}{15}\sum^{\infty}_{n=1}\left( \frac{2}{\sin^2(\pi n\tau)} - \frac{15}{\sin^4(\pi n\tau)} + \frac{15}{\sin^6(\pi n\tau)} \right)$$
計算ミスが怖いけど、まぁグラフ書かせるかぁー
G4の画像 G4の画像

お、いい感じ。

G4の拡大図 G4の拡大図

思ったより緻密だねぇ?

G6の画像 G6の画像
うむうむ、よろしい。

ちなみに、一緒にこうやって遊んでたやつが、「微分のところ一般化できないかな」なんていってこんな式を得ていた。

$$ \begin{align} &G_{2k}(\tau)=(-1)^{k+1}\frac{(2\pi)^{2k}}{(2k)!}\left( B_{2k} -2k\sum^{\infty}_{n=1}\frac{A_{2k-1}(q^n)}{(1-q^n)^{2k}} \right) \\ &ただしq=e^{2\pi i \tau} \end{align} $$

なお、$A_{n}(x)$はオイラー多項式、$B_{n}$はおなじみベルヌーイ数ですぞ。
要望があれば導出かくかも。でも収束遅いしまた下半平面で収束してないし、作った本人もじゃあ何ってところまで行ってなかったから多分そんな要望は来ない。たぶん。

ではまた。

投稿日:202354
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中