blog

特性方程式と重根

今日の講義で微分方程式の特性根と一般解の話をしたら,重根をもつ場合について,高校で習う隣接3項漸化式で重根を持つ場合の解き方と同じ考え方で微分方程式も説明できることが分かったと講義の後で教えてくれました.以下にまとめてみます.

隣接3項漸化式で特性方程式が重根をもつ場合

[latex]y(k+2)=2ay(k+1)-a^2y(k)[/latex]
を変形して,差が等比数列の形になるようにする.

[latex]y(k+2)-ay(k+1)=a(y(k+1)-ay(k))[/latex]

右辺は
[latex]a^{k+1}(y(1)-ay(0))[/latex]
となる.

両辺を[latex]a^{k+1}[/latex]で割る.

[latex]
\displaystyle\frac{y(k+2)}{a^{k+1}}-\frac{y(k+1)}{a^k}=y(1)-ay(0)
[/latex]
となり,[latex]\displaystyle\frac{y(k+1)}{a^k}[/latex]は,初項[latex]y(0)/a^{-1}[/latex],公差[latex]y(1)-ay(0)[/latex]の等差数列になる.

したがって
[latex]
\displaystyle\frac{y(k)}{a^{k-1}}=\frac{y(0)}{a^{-1}} + (k-1)(y(1)-ay(0))
[/latex]
となる.

両辺を[latex]a^{k-1}[/latex]倍して
[latex]
y(k)=a^ky(0)+a^{k-1}(k-1)(y(1)-ay(0))
[/latex]
になる.

隣接3項漸化式では,[latex]y(k)[/latex]を求めることになるが,線形差分方程式とみた場合には,基底に[latex]a^k[/latex]と[latex]ka^k[/latex]を用いればよいことが分かる.

微分方程式

微分方程式の特性根が重根の場合,
[latex]
\displaystyle\frac{d^2y(t)}{dt^2} -2a\frac{dy(t)}{dt} + a^2y(t)=0
[/latex]
を変形して
[latex]
\displaystyle\frac{d^2y(t)}{dt^2} – a\frac{dy(t)}{dt} = a\left(\frac{dy(t)}{dt} – ay(t)\right)
[/latex]
となる.
これは,
[latex]
\displaystyle\frac{d}{dt}\left(\frac{dy(t)}{dt}-ay(t)\right)= a\left(\frac{dy(t)}{dt} – ay(t)\right)
[/latex]
なので,
[latex]
\displaystyle\frac{dy(t)}{dt}-ay(t) = e^{at}\left(\frac{dy(0)}{dt} – ay(0)\right)
[/latex]
になる.
[latex]e^{at}[/latex]で両辺を割ると,
[latex]
e^{-at}\displaystyle\frac{dy(t)}{dt}-ae^{-at}y(t) = \frac{dy(0)}{dt} -ay(0)
[/latex]
なので,左辺を次のように書き換えて
[latex]
\displaystyle\frac{d}{dt}\left(e^{at}y(t)\right) = \frac{dy(0)}{dt} -ay(0)
[/latex]
になる.
積分して,
[latex]
e^{-at}y(t) = \left(\displaystyle\frac{dy(0)}{dt} -ay(0)\right)t + c
[/latex]
になり([latex]c[/latex]は積分定数),
[latex]
y(t)=\left(\displaystyle\frac{dy(0)}{dt} -ay(0)\right)te^{at} + ce^{at}
[/latex]
になる,
つまり,基底は[latex]e^{at}[/latex]と[latex]te^{at}[/latex]になる.

実際問題として,重根のときの基底を導出する時間はないでしょうから,重根のときの基底を覚えておくとよいでしょう.また,講義で説明した定数変化法を理解しておくとよいと思います.

線形と非線形

月曜日は2時限に物理情報数学B,4時限に生体制御を担当しています.

物理情報数学Bは線形代数です.先週提出の宿題,ランチェスターのモデルの問題は,講義中に例を板書した,α=0.1,β=0.1では,さすがにつまらないです.ノートを右から左に写すだけになってしまいます.自分でαとβの値を設定して,つまり問題を自ら設定できないと評価されません.ランチェスターのモデルでは固有値が簡単な値になるようにαとβを定めることは簡単なことです.

生体制御は,非線形振動を扱いました.BVPモデル,好中球の密度の振動,心循環系の圧反射のモデル,サーカディアンリズムです.私が大学院生のとき,非線形振動だけで1つの講義があり,van der Polの式を数値的に解いて,リミットサイクルを描き,また位相面にメッシュを切って各格子点の速度ベクトルを描くことが課題だったように思います.当時のPCは,8086 8 MHzやV30 10 MHzの16ビットでMS-DOSで動いていました.XYプロッターに出力するプログラムを書いた記憶があります.

7セグメント表示器

デジタル基礎の講義で7セグメント表示器の任意のセグメントを光らせる組み合わせ回路を設計する問題でした.Don’t careがあるカルノーマップの問題です.7セグメント表示器では,10〜15はdont’ careになります.

 

丁度,ネットで
1=2, 2=5, 3=5, 4=4…
という問題がアナログな頭では解けないというのを見かけたところでした.7セグメント表示の文字はやはりデジタルの代名詞なのかと思ったところです.

 

World Congress 2012

World Congress 2012 Medical Physics and Biomedical Engineeringで”Independent Component Analysis of Mechanomyogram Detected with an Acceleration Sensor in Motion”で発表してきました.

会場は,Beijing International Convention Centerです.道路向かいには鳥の巣がありました.

金環日食

自宅で撮影しました.EOS 60D,EF 70-200 mm F4 IS,  1/400の減光フィルタです.ISO 200, F 10です.この日は,雲で太陽が隠れたり雲がきれたりと光量の変化が激しかったので,シャッター速度を大きく変えて撮影しました.上から,1/60,1/4000,1/2です.雲が切れたときには1/400の減光フィルタを2枚使いたかったところです.写真はドットバイドットでトリミングしたものを1/2に縮小しています.

 

第51回生体医工学会大会

第51回生体医工学会大会が5月10日から5月12日まで,福岡国際会議場で開催されました.本研究室では,以下の3件の発表を行いました.

  • 有岡:機械刺激による筋・脂肪の力学モデルのシステム同定
  • 田村:筋音図のシステム同定による単一運動単位収縮モデルの検討
  • 酒井:短母指外転筋における誘発筋音図のシステム同定

スーパームーン

金環日食の写真を載せたら,スーパームーンも載せておいた方がという気になったので載せます.EOS 60D EF 70-200mm, ISO 200, F 10, 1/200です.ドットバイドットに拡大してトリミングしています.

月曜日代替日

今日は火曜日ですが,月曜日の授業が行われます.月曜日は祝日と振替休日で休みになることが多いので,代替日が設けられることが多くなります.

開校記念日

4月23日は開校記念日です.一昨年までは休講でしたが,今年度はカレンダーの都合上,授業日です.