| 疲れた | 2005年9月20日(火) | |
| Infoseek | ||
| 今日は曇。疲れた。 | ||
| 采女祭 | 2005年9月19日(月) | |
| Infoseek | ||
| 今日は晴後曇。下の写真は昨日行われた采女祭の最後、花扇を猿沢池へ沈めるシーン。 | ||
| ||
| JWordは有害サイトではないのか | 2005年9月17日(土) | |
| Infoseek | ||
| 今日は晴後雨。白い彼岸花が咲いている。 | ||
| ||
さて最近、あの悪質スパイウエアJWordが、悪質·犯罪関連サイトの被害抑止·軽減に向け警視庁と協力し始めた、とある。
| ||
| 日経夕刊より | 2005年9月15日(木) | |
| Infoseek | ||
今日は晴時々曇。朝晩、涼しい。さて久し振りに心理学寄りの話。日経夕刊表紙の左下の囲み記事のタイトルは「脳」整理法。
| ||
| 疲れた | 2005年9月14日(水) | |
| Infoseek | ||
| 今日は晴時々曇。疲れた。下の写真は苦瓜と豚の炒め物。 | ||
| ||
| 共役勾配法への疑問 | 2005年9月13日(火) | |
| Infoseek | ||
| 今日は晴。相変わらず暑い。下の写真は数日前に行ったしゃぶしゃぶ。最近、「バサラ邪馬台国」や「愛元麻樹」をキーワードとする検索でこの日記に来る事がある様だ。 | ||
| ||
| さてこの所、数値計算技法の話を出しているが、連立1次方程式の非定常反復解法として"共役勾配法"(conjugate gradient method)は非常に有名だ。本来、係数行列が対称(複素Hermit)·正定値でなければならない制約があったが、これに縛られない進化した方法も提案されている。共役勾配法のアルゴリズムは以下の通り。先ず | ||
|
r0=b0-Ax0 r0'=Mr0 p0=r0' | ||
| が初期化の段階であり、次の一連の手続きを残差ri+1が一定以下になる迄続ける。なおAが係数行列、Mが前処理行列だ。 | ||
|
qi=Api αi=<pi|ri>/<pi|qi> xi+1=xi+αipi ri+1=ri-αiqi ri+1'=Mri+1 βi=<ri+1'|qi>/<pi|qi> pi+1=ri+1'+βipi | ||
| このやり方に、何か疑問はないだろうか。そう、初期化の段階で残差がr0=b0-Ax0の様に定義されている。しかし、途中の計算でAが作用する相手は中間変数piであり、残差は漸化式ri+1=ri-αqiで更新されている。従って数値的に本当にri=b0-Axiが保証されているか否か疑念が涌くだけでなく、非線形方程式への拡張も厄介だ。厳密にri=b0-Axiを満たす事は面倒であるにしても、予測子·修正子型の手続き | ||
|
xi+1*=xi+pi ri+1*=b0-Axi+1* qi=-(ri+1*-ri) αi=<pi|ri>/<pi|qi>-1 xi+1=xi+1*+αipi ri+1=ri+1*-αiqi ri+1'=Mri+1 βi=<ri+1'|qi>/<pi|qi> pi+1=ri+1'+βipi | ||
| に改善すれば、何をやっているか、多少は分かり易くなる。 | ||
| アルコール各種 | 2005年9月12日(月) | |
| Infoseek | ||
| 今日は晴。左から泡盛、薬用酒、そしてラムが2種類。 | ||
| ||
| ゴロゴロ | 2005年9月11日(日) | |
| Infoseek | ||
| 今日は曇時々雨。ゴロゴロ。 | ||
| 表紙へ戻る | ||