hirax.net::Keywords::「いわき」のブログ

■渡り鳥の秘密　

3000kmの彼方へ

　先日、福島のいわき市の近くに友人を訪ねた。その近くで、はるか海を渡ってきた白鳥の群れを見た。そこでどこから来るのかとか、どうやって来るのかなどの話をした。大体、以下のような疑問である。

• Q.1:海を渡るのにどの位の時間がかかるのか?
• Q.2:位置(あるいは方向)はどのように認識するのか?
• Q.3:どこから来るのか?
• Q.4:なぜ渡るのか?

Q.3のどこから来るのかについては、仙台市科学館の中のGooseSchool(http://www.smus.city.sendai.jp/goose/home.html)に簡単な説明があった。そこにガンの「渡り」の地図が示してあった。その図を引用させて頂く。いや、もちろん「白鳥でなくガンであるし、そもそも飛んでくる場所が違う」のはわかっている。しかし、先の疑問の「どこから来るのか」はべつにいわき周辺の白鳥にこだわった疑問ではないのだから、他の場所の他の動物でもいいのだ。

(画像リンク先は仙台市科学館　http://www.smus.city.sendai.jp/)

• Q.1: 海を渡るのにどの位の時間がかかるのか?
• Q.2: 位置(あるいは方向)はどのように認識するのか?

• A.1: 人工衛星用位置送信機を取り付け経路を追跡した結果では2日間で2000ｋｍをノンストップで渡り続ける。
• A.2: 時間は体内時計で知り，方向は星や太陽で分かる。鳥はコンパスや六分儀のような能力を両方持っている。

　そこで、簡単に手に入った

1. 「鳥の渡りの謎」 (ロビン・ベーカー著　網野ゆき子訳　平凡社発行　ISBN4-582-52717-5)
2. 「鳥学の世界へようこそ」　(ディヴィッド・ラック著　蓮尾純子　訳　平河出版社ISBN4-89203-192-5)

　しかし、今回は自分で実験も解析もしていないので、簡単なまとめだけにしておく。

• 鳥(種類、地域、条件)により「渡り」の際に使用する航法(navigation)は異なる。
• 方位を頼りにするやり方も、位置を頼りにするやり方もある。使用されるのはどちらか片方のみというわけでもない。
• 位置のセンサーとしては、目印、音、におい、磁場*、重力*
• 方位のセンサーとしては、太陽の位置、光の偏光パターン、星、月などの位置、地磁気*

　さて、今回は下調べばかりであった。これでは、いけない。いけない。自分で実験もしなければ意味が無い。というわけで、近日中にある実験をする予定である。

■分数階微分の謎　

線形代数、分数階微分、シュレディンガー方程式の三題話

分数階微分？

InterLabの榊原教授の記事を引用すると、-通常微分・積分は整数回実行できるが、分数階微分はこれを分数に一般化したものである。さまざまな物理や工学の現象の記述に使われるようになった-とある。一階微分とか二階微分というものはよく使うが、0.5階微分などというものは使ったことがない。どのようなモノなのかさえよくわからない。

参考：

• InterLab (http://www.inter-lab.gr.jp/)

いわき明星大学の清水・榊原研究室の「粘弾性動モデル」が引っ掛かる。

参考：

• いわき明星大学清水研究室 (http://www.iwakimu.ac.jp/lab/nshimlab/siricon.html)

もう少し調べてみると「バナッハ空間バナッハスケールにおける分数階積分作用素」というようなキーワードも引っ掛かる。

そこで、まずは勝手に分数階微分について考えてみた。

分数階微分・積分の勝手な想像図

まずは、イメージを考えるためにグラフを作成してみる。x^2の関数、および、それを微分・積分した関数である。微分は3階まで、積分は2階まで行っている。

図.1:x^2を微分(3階まで)したものと、2階まで積分したもの

このグラフ形式の表示をちょっとだけ変えてみる。

図.2:x^2を微分(3階まで)したものと、2階まで積分したもの

ここまでくると、平面グラフにしてみたくなる。つまり、微分・積分の階数を離散的な整数値でなく、連続的な値としてのイメージに変えたくなる。

図.3:x^2を微分(3階まで)したものと、2階まで積分したもの

これで、微分・積分が整数階でない場合のイメージ(勝手な)ができた。微分・積分が離散的なものではなくスムーズにつながっているものであるというイメージである。図.2から図.3への変化をよく覚えていてほしい。

といっても、これは数学的なイメージのみで物理的なイメージはまだここでは持っていない。位置、速度、加速度などの微分・積分で選られるものに対して同じようなイメージを適用すると、位置なんだけれどちょっと加速度っぽいもの、とか、速度と加速度の「合いの子」みたいなものというような感じだろうか？

さらに、これから先は、f(x)という関数が示す無限個の値を位置ベクトルと考えて、f(x)というのは無限次元空間の一つの点だというイメージを持つことにする。線形代数を考えるならそれが一番わかりやすいだろう。任意の階で微分された関数群が集まって、さらに高次元の空間をなしているというイメージである。

分数階微分を調べる

その中の言葉を少し引くと、
フーリエ変換は等距離作用素である、関数空間L^2(R)における回転といえる。結局、

ここで、fは元の関数であり、Fはフーリエ変換

さて、この式自体は非常に簡単である。それにイメージも湧きやすい。
i を掛ける演算、私のイメージでは複素数空間の中で90度回転をする(言い換えれば、位相が90度ずれる)演算、が微分・積分であるというイメージはスムーズに受け入れやすい(それが正しいかどうかは知らないが)。なぜなら、微分が空間の中での回転であるとすると、三角関数の微分・積分に関する性質(例えば、Sinを微分するとCosに、Sinを2階微分すると-Sinになる、すなわち、一回の微分につき位相が90°ずつ回転する(位相がずれる)というような性質)が納得でき、それがフーリエ変換という形で登場してくることがスムーズに受け入れられるのである。また、微分といえばとりあえず三角関数の登場というイメージもある。

　もう少しわかりやすく書くと、

• 三角関数では一階微分の結果は90度位相がずれる(回転する)。
• ならば、(例えば)0.5階微分は45度位相をずらせば良い。
• 任意の関数もフーリエ変換により、三角関数に分解される。
• ならば、任意の関数に任意の実数値の微分が成立する。

　任意の関数をフーリエ変換し三角関数に分解した時の位相、言い換えれば、周波数領域での位相ずらし、で分数階微分が定義されるということは、物理的実用的に大きな意味を持つ。例えば、電磁波、弾塑性運動などの物理現象の中での位相変化を分数階微分で解けることになる。例えば、複素貯蔵弾性率などについて分数階微分との関係は深そうである。あるいは、媒体中の電磁波の位相などについて適用するのも面白そうである。

分数階微分を使ってみる

よく分からないところも多いが、とりあえず、

ガウス分布(左)とその一階微分の解析解(右)

それでは、今回の方法による一階微分の結果と、それと解析解との比較を示す。なお、本来無限領域のフーリエ変換を有限の領域で行っているため、端部近くで変なことが生じるのはしかたがないだろう。また、色々な事情により係数の違いは無視して欲しい。

フーリエ変換を用いた方法(左)と解析解(右)の比較

ちょっとずれが生じているが、こんなものだろう。しかし、これだけでは今回のフーリエ変換を用いた微分の面白さはでてこないので、0から2の範囲で連続的に分数階微分をしてみる。

ガウス分布の0から2の範囲における連続的な分数階微分

1/10 (=0.1)階微分	1/2 (=0.5)階微分	7/10 (=0.7)階微分	1階微分
13/10 (=1.3)階微分	15/10 (=1.5)階微分	17/10 (=1.7)階微分	2階微分

モーフィングのようで面白い。

さて、今回は分数階微分を勉強してみる所までで、これの応用は別に行ってみたい。もちろん、言うまでもないと思うが、間違いは多々あると思う。いや、田舎に住んでいるもので資料がないんですよ。

■いわき市三崎　

　いわき市三崎で。

■いわき　

　茶碗に晩春の草木を入れて、飲んでみる。

■いわき