1/fゆらぎ
非整数ブラウン運動の説明用のアニメを作ってみました． ブラウン運動 (H = 0.5)のサンプルパスの拡大非整数ブラウン運動 (H = 0.8)のサンプルパス非整数ブラウン運動 (H = 0.3)のサンプルパス高解像度な，Youtube版は以下です． www.youtube.com www.youtub…
前回に続き，ARFIMA(0, d, 0)を，無限次の自己回帰過程 (AR過程: autoregressive process)で表現する話です．今回は数値実験です． 1. 基本事項：long AR過程表現 2. 数値計算のポイント 3. 数値計算結果 Rスクリプト 前回の記事は，これです． chaos-kiyono…
今回は，1/fノイズ (ピンクノイズ)の自己共分散 (自己相関)関数についての話です．つまり，パワースペクトルが となる時系列の自己共分散 (自己相関)関数は， ということを説明します．解析的に導出しようかと思いましたが，この論文 doi.org の流れをなぞる…
これまでに，ピンクノイズ (1/fゆらぎ)のサンプル時系列の生成方法について，私が知っているものは，ほぼすべて紹介しました (自己相関関数を指定するものだけ説明してません)．今回は，ピンクノイズ (1/fゆらぎ)の生成方法をまとめておきます． Rのパッケー…
ピンクノイズ (1/fノイズ)を聴く効果については，いろいろな説があるようです．例えば，以下のページ muysalud.com では， 眠りにつくのを助ける 深い眠りにいざなう 記憶力を改善させる 集中力を高める という効果があると説明されています．私は，ピンクノ…
今回は，中点変位法 (midpoint displacement algorithm)で，1/fノイズ (ピンクノイズ)のサンプル時系列を生成します．元々，中点変位法は，非整数ブラウン運動のサンプルパス (サンプル時系列)を作る方法ですが，非整数ブラウン運動をわずかにはみ出た領域を…
前に，階積分 (和分)過程のパワースペクトルの記事を書いたときに導いた結果，分散の白色ノイズの階差分過程の時系列のパワースペクトルは， 階積分 (和分)過程の時系列のパワースペクトルは， について，数値実験で確認していなかったので，やっておきます…
今回はパワースペクトルが，を定数として， の形になる，差分方程式で記述される確率過程を考えます．前の記事 chaos-kiyono.hatenablog.com で，m階積分 (和分)過程のパワースペクトル（もどき）を計算しました．つまり， 分散の白色ノイズの階積分 (和分)…
今回は，Autoregressive fractionally integrated moving average process (ARFIMA)への導入の手前です．ARFIMAの日本語訳は，私が訳せば「自己回帰非整数積分移動平均過程」，インターネット検索でよく見る訳は「自己回帰実数和分移動平均過程」です．MA (m…
時系列の長期記憶 (long memory)，長時間相関 (long-range correlation)，フラクタル性 (fractality)について整理しておきます． 長期記憶過程 (long memory process) 長期記憶過程は自己共分散関数が， となる確率過程です．何で全部足すの？という疑問がわ…
「ランダムウォーク解析」の話の続きです．前回の記事は，これです． chaos-kiyono.hatenablog.com ランダムウォーク解析では，(1) 観測時系列の平均を０にして，(2) 積分して，(3) その増分の2乗平均の平方根をいろんなスケールで計算して，(4) 対のプロッ…
今回も，ランダムウォーク解析についてです．ランダムウォーク解析 (fluctuation analysisと呼ばれることも)は，Natureに出版された論文でも使われています (以下のリンク)． www.nature.com とはいえ，ランダムウォーク解析は，今では実用的な時系列解析法…
ランダムウォーク解析について前回に引き続き考えていきます．今回は自己相関（自己共分散）関数との関係です．前回の記事は以下です．chaos-kiyono.hatenablog.com ゆらぎ関数と自己相関（自己共分散）関数との関係 ここでは．平均0，分散の弱定常過程の時…
今回は，ランダムウォーク解析の導入です．detrended fluctuation analysis (DFA)とか，detrending moving average analysis (DMA)とかを使ってる研究者の方が，そこそこいると思います．ランダムウォーク解析は，DFA，DMAの原型となる，素朴な方法です．ラ…
私の記事で扱いたいテーマの一つに，1/fノイズ（1/fゆらぎ，ピンクノイズ）があります．私は，1/fゆらぎの科学的意義を伝えたいのですが，世の中では「1/fゆらぎには，リラックス作用がある」と認識されていて，根強い人気があるようです． 私のYoutubeチャ…
1/fノイズやホワイトノイズを音として体験する動画資料です．視聴する際は音量を低めに設定してください．1/fノイズ（ピンクノイズ）の音です． www.youtube.comホワイトノイズの音です． www.youtube.comホワイトノイズのハイパス信号です．カットオフ周波…
前回は，時系列の積分がパワースペクトルにもたらす変化「ゆらぎを積分すると，ゆらぎになる」を解析的に示しました．chaos-kiyono.hatenablog.com今回は，数値実験です． 非整数ガウスノイズ（左）とその積分（右），下は対応するパワースペクトルのlog-log…
今回は，「ゆらぎを積分すると，ゆらぎになる」ことの説明です．つまり，時系列のパワースペクトルが，低周波数側で に従うとき，を積分した時系列 ] ()のパワースペクトルが，低周波数側で になることの説明です（下図参照）．逆に，差分をとるとになります…
パワースペクトル解析，長時間相関解析などについて，私の古い解説動画のリンクをまとめておきます．興味があれば見てみてください．中途半端なところで解説動画が止まっているので，このブログで文章での説明がある程度まとまれば，新しく解説動画を作って…
今回は，周波数がの領域で1/fゆらぎ（ピンクノイズ）の特性を持つサンプル時系列の生成をやります．詳しい説明はしませんが，1/fゆらぎの生成メカニズムとして，自己共分散関数が指数関数的な減衰を示すサブシステムの和を考えるものが提案されています（下…
今回は，狭い意味での1/fゆらぎを生成します．ピンクノイズとも呼ばれます．なぜピンクかといえば，色とのアナロジーです．可視光の低周波数側の端が赤色なので，低周波側，つまり，白色よりも赤色が強調されていると考えられるので，色で言えば，どちらかと…
今回は，何で「1/f（えふぶんのいち）ゆらぎ」とか，1/fノイズに注目するのか，何が不思議なのか，何がおもしろいのか，というお話です．1/fゆらぎについて，歴史と最新の成果を正確に記述することは，私には難しいので，まずは，現状の私の知識，印象，思い…
