スキップしてメイン コンテンツに移動

オイラーの式の証明その後、あるいは年初の目標進捗

覚え書きもあり、今週末の感想を。

年初の目標の一つ、自然数の和


1+2+3+4+...=-1/12


の証明ができるようになるため、年末に「リーマン予想のこれまでとこれから」を購入した。



ざっと一回読んでみて、なるほどリーマンのゼータ関数は広がりと深みがあることがわかった。これはこれで面白い。第I部「リーマン予想への助走」の中の「有限ゼータ関数」と「無限への接近」、第II部の「ピタゴラスからオイラーまで」で、その式はオイラーにより導出されている。一つ一つ手で確かめながら進めており、解るたびに自分が好きになる。しかし、この本はそんな簡単な式に拘泥すること無く先に進む。そこで、同じ著者の黒川信重が書いた「オイラー探検 無限大の滝と12連峰」を買い求めてこちらも読んでみた。



その式は12章の「自然数全体の和:オイラー瀑布」に取り扱われている。


  • 1-2+3-4+5-6+7-8+9-10...=1/4


で、自然数の全体の和から偶数の和の2倍を引いた式になっており、その結果自然数の和が-1/12と出てくる。

つまり

1-2+3-4+5-6+7-8+9-10...
=(1+2+3...)-2(2+4+6...)
=(1+2+3...)-4(1+2+3...)
=-3(1+2+3...)

より
1+2+3+...=-1/12

オイラーは冪級数
1+x+x^2+x^3+... = 1/(1-x)
を用いる。これは1/(1-x)のテーラー展開したもので、これを二乗し
(1+x+x^2+x^3...)(1+x+x^2+x^3...)=1/(1-x)^2
の左辺は、
(1+2x+3x^2+4x^3...)であるからx=-1を代入すると
左辺=1-2+3-4+5-6...
右辺=1/(1+1)^2=1/4
である。ここまでは収束性を考慮に入れていないが、ゼータ関数のその後の研究でこの考え方は基本的に正しいのだそうだ。この正しいそうだを確信にすることが次の目標。問題をブレークダウンすることができた、というのが1月の成果。

今更ながら無限とか収束とかキチンと理解しておらず、そういえばテーラー展開も証明はできないことに気がつき、書斎をごそごそ探すと、「改訂微分積分学要論(青木利夫・吉原健一共著)」が出てきた。これは古本屋で購入したもの。自分は本当は理解していなかったことを改めて感ずる。



この本が意外と良い。解りやすいのだ。大学で微分積分を学ぶためのテキストで大学卒業後20年を超えてもまだ読める。たしかに1/(1-x)のテーラ展開が可能なのは|x|<1の区間であるとしている。

次に攻めるべきは、この微積の演習をいくつか解いてみてまた黒川の2番目の本に戻ってこよう。

コメント

このブログの人気の投稿

瞬間日記のデータをday oneに移動させるのだ

瞬間日記をiPodを購入した2010年から使ってきたが、day oneが人気なので使ってみた。なかなかよかった。コマンドライン(CLI:Command Line Interface)から入力できるし、クラウドで同期できるのもいい。見た目もきれい。安定しているのもいい。 という訳で、瞬間日記からday oneにデータを引っ越しすることにした。 (このページで半角の &gt; と表示されたら半角の”>”記号だと思ってください。HTMLの仕様のインプリバグでしょう。) 特定の日記アプリ間のデータの移動とは人生で一回しかないだろうから汎用性など考えずに間違えないように慎重にすることが大切だ。確かめて確かめて。このページはそんな備忘録である。また、瞬間日記のデータをday oneに移動させたという記事も見つけられなかったので誰かのお役に立てるかもしれない。でも自己責任でね。 瞬間日記は独自フォーマット(.dat3)や標準的な(.sqlite)、(.csv)でもMacにバックアップできる。使い勝手のいいcsvファイルは本文だけで写真情報が欠落している。瞬間日記のデータを引き上げるのは色々試したが、テキストばかりでなく写真も移動させたいの結局メールでバックアップすることにした。瞬間日記側で自分にひと月毎のバックアップをメールする。3年半分。ソチオリンピック開会式を見ながら作業でもれなく終了。 これは瞬間日記が30枚しか1つのメールに写真を貼付できないので、ひと月分だけせっせとメールする。 届いたメールはOSX側のmail.appで開くが、~/Libraryに保存されているelmsでなく、単にテキストファイルでmail.app側の機能を使って出力することにする。 瞬間日記のデータを取り出し、day oneにエントリできるくらいの粒度のファイルにするまでの流れ: iPodの瞬間日記のひと月分に区切ったバックアップデータを自分のメールアドレスに送信する OSXのmail.appで受信したメール全部(件名"MomentDialy"で始まるメールを一つのメールボックスに束ねる。例えばMomentDiaryというメールボックスに全部入れる mail.app上のメールボックスMomentDiary内の全部のメールを選択する。つまり⌘+a。 ...

電解コンデンサの容量抜けを測ってみた

中古オーディオ屋でスーパーウーハーを買ってから7〜8年経つだろうか。 1989年に発売されたONKYOの SL-10 という機種で、現役で商品が店頭に並んでいた時にピアノ曲に深みが増すことにちょっとした感動を覚えて中古屋で見つけた時2万円で即購入したのだ。 スピーカのエッジのクズがエンクロージャのダクトから出てきて開けてみたら見事エッジ全体が下手っていた。これを奇麗にして、ウレタンエッジに交換することにした。写真1、と写真2がそのBefore/Afterだ。 写真1 SL-10のスピーカのエッジを交換してみた (エッジを取り除いた状態) 写真2 SL-10のエッジを交換したスピーカ また、SL-10が勝手に自己発振する様になっていたのでコンデンサの容量が抜けたと想像してこれらを交換することにした。これは電解コンデンサの容量が抜ける事でフィードバック回路の時定数が変わってしまいネガティブフィードバックがポジティブフィードバックに位相が回って知ったのでは予想した。 電解コンデンサを交換したSL-10のコンデンサ 交換した電解コンデンサ達(容量は抜けていなかった)  マルツ電波で発注しておいたオーディオ用のコンデンサを付け替えた。ただ、それだけでは詰まらないので交換したコンデンサを測ってみた。使った測定器は卓上テスタ(GBW 9000A)の付属機能。最大20μFまでしか測れなかったが、結論からすると交換した電解コンデンサの容量抜けは無かった。 測定例 定格(μF)   実測(μF) 10       10.19 10       10.00 4.7       4.78 3.3       3.27 1.0       1.01 0.47       0.46 というわけで直しても自己発振は直らないかも。明日組立ててみよう。

声紋・ソノグラム・スペクトルグラム

声紋(sonogram, spectrogram)は野鳥の声の分析に不可欠なツールです。横軸に時間、縦軸に周波数を取って、音の強さを色で表示します。 ここで用いられている処理は、音声である波形を周波数の信号に変換する時間−周波数変換(フーリエ変換)です。フーリエ変換するときに幾つかの調整できるパラメータがありそれぞれがどう関係しているのかをまとめてみました。案外、まとめて書かれているものがありません。そのため自分で調べました。厳密に定義しているのではなく備忘録です。初心者向けにコ ーネル大学の作成したCanary/Ravenのマニュアルの付録に鳥の声のスペクトル解析の章 がありこれを参考にしています[1]。 音声データのパラメータ 波形とは時間に対する強度が連続的に変化する信号を表します。音声は波形であり、ICレコーダは連続の波形を離散化(デジタル化)されて記録されます。離散化された信号も波形と呼びまし、今はディジタル信号を意味する事が多いです。 連続な信号をとびとびの信号で表すので失われる情報があります。しかし、必要とする周波数帯域の2倍以上でサンプリングレート(サンプリング周波数)で離散化すれば完全に復元できるとシャノンさんが証明しました。 この関係は、次の表のようになります。 サンプリングレート(SR) 標本点の周期(ΔT) 192KHz 5.2μs 96KHz 10.42μs 48KHz 20.83μs 44.1KHz 22.67μs 音を声紋(スペクトログラム)にするには短い一定期間の波形データを切り出して(スライス)これを周波数変換し並べます。この短時間の周波数変換を短時間フーリエ変換:STFT(Short Time Fourier Transform)と呼びます。問題はこの時間の選び方です。時間と言っても離散化されているのでその切り出す波形の個数です。切り出された波形をフレームと呼んだりチャンクと呼んだりします。ここではフレームと呼びます。1フレームの波形の個数をフレーム長(Flame Length)と呼びます。各アプリではもっとわかりやすい言葉を使っていますが同じ意味です。 Amadeus Pro : FFT Size Audacity : ウィンドウサイズ 鳥の...