日曜開発者のブログ -12ページ目

日曜開発者のブログ

日曜プログラマーが自作ソフトを細々と公開していくブログです。一般化調和解析(GHA)を応用したノイズ除去ソフトを公開しました。

『最大エントロピー法による時系列解析 ☆〔第2版〕☆ ― MemCalcの理論と実際』 に掲載されている計算例題を使って、新しいGHA実装の精度検証をしています。

今回の例題である計算例17番の入力波形は次式で与えられます。


式


標本化周波数は1秒です。
データの入力画面は下のようになります。
以下の例では解析時間を1300秒に設定しています。




解析区間1300秒の場合は次の解析結果になります。



解析区間650秒の場合は次の解析結果になります。


解析区間325秒の場合は次の解析結果になります。


T=325


解析区間169秒の場合は次の解析結果になります。

T=169


解析区間78秒の場合は次の解析結果になります。

T=78


解析区間50秒の場合は次の解析結果になります。


T=50

解析区間25秒の場合は次の解析結果になります。
波を構成する成分の周期の半分程度の時間の情報でもほぼ正確に推定が出来ています。

T=25