2006年09月11日

2次のburg法で周波数推定(2)

2次のburg法で周波数推定の続きです。

2つの係数が与えられたときの周波数が分かったので、Rでこれを求める関数を作ってみましょう。後で窓関数を入れる事が出来るようにしておきます。

frq.burg <- function(x){
len <- length(x)
av_av <- mean(x)
x <- x - av_av
si <- 0.0
bo <- 0.0
for(k in 1:(len-2)){
# Normal Burg
W <- 1/(len-2)
# main
si <- si + (x[k] + x[k+2]) * x[k+1] *W
bo <- bo + x[k+1]*x[k+1] * W
}
a1 <- -si/(bo *2)
si <- 0.0
bo <- 0.0

for(k in 1:(len-2)){
# Normal Burg
W <- 1/(len-2)
# main
si <- si+(x[k]+a1*x[k+1])*(x[k+2]+a1*x[k+1]) *W
bo <- bo+((x[k]+a1*x[k+1])*(x[k]+a1*x[k+1])+(x[k+2]+a1*x[k+1])*(x[k+2]+a1*x[k+1]))*W
}
a2 <- - 2*si/bo
a1 <- (1 + a2 )* a1
peak <- (pi-acos(a1/(4*a2)+a1/4))/(2*pi)
return(peak)
}

ちょっと無駄がありますが..

使い方ですが、たとえば時系列 ts があったときに
freq <- ar.otburg(ts)*frequency(ts)

のようにして使います。

さて3次の場合を求めることができますか?
(4次まではなんとか出来ますが、やらないほうが身のためです。)
posted by metola at 13:28| Comment(17) | TrackBack(0) | 統計解析システムR | このブログの読者になる | 更新情報をチェックする
この記事へのコメント
出合いを求めているならココ!無料出逢いの決定版!退屈な毎日から抜け出せるチャンスがここに!日頃出合いが無いとお嘆きの方、お困りの方は当サイトにぜひ!ご利用するだけで素敵な出合いが体験できること間違いなし!
Posted by 無料出逢い at 2011年04月23日 14:56
ヴィクトリアマイルのデータ満載!オッズ・枠順・出走馬情報!弊社が収集した極上データを惜しみなく公開!最新の情報を現地からお届けします!勝ちにいく予想はここから
Posted by ヴィクトリアマイル at 2011年05月04日 00:13
出会いは京の町から♪京都の女の子はいつも人気ランキング上位!!そんな奥ゆかしい京ガール達と気軽に出会えるサイトがついに誕生☆出会いの都、京都でステキな京ガールと出会っちゃおう!
Posted by 京都 出会い at 2011年05月27日 04:39
歴史マニア必見!歴史が好きな歴女が集まる歴女専門サイトの誕生です!歴史好きな女の子と楽しく歴史トークをしてみませんか?同じ歴史好きならすぐに仲良くなれるはず!
Posted by 歴女 at 2011年07月02日 19:19

尺八してもらったけどありゃたまんねーな!!
ジュパジュパ凄い音させながら吸いつかれて、30秒で発射しちまった(笑)
しかもオレ、女にお任せして寝てただけなのに5マソも貰った件wwww
http://jmd8m64.jp.takaoka.mobi/
Posted by 吸引力の変わらないただひとつの… at 2011年07月31日 19:16
女性にいじめられるのが大好きなM男のそこのアナタ、アナタのM度を測れるサイトが誕生しました!!自分のM度はどれくらいなのか?我こそは最強の変態だというアナタの参加をお待ちしております!
Posted by M男 at 2011年08月08日 04:47
早く童貞を捨てたいと思っているあなたにピッタリのサイトが誕生!このサイトで長年守っていたあなたの童貞がついに奪われる!?あなたの童貞卒業を優しくサポートしちゃいます☆
Posted by 童貞 at 2011年08月22日 06:42
空いた時間にお小遣い稼ぎをしませんか?ライター募集のご案内です。簡単な在宅ワークで安定した副収入を得ることができます。努力しだいで高額の副収入を得ることも可能です。
Posted by ライター募集 at 2011年09月12日 07:30
ニューハーフとの相性占いしてみませんか?あなたがニューハーフと付き合える可能性がどのくらいあるのかがわかります!今すぐチェックしてみてください
Posted by ニューハーフ at 2011年09月20日 21:55
主婦の方にぴったりの副業を初めてみませんか?在宅で出来るので時間にとらわれずにお仕事していただけます。真面目に取り組んで頂ける方からの募集をお待ちしております。
Posted by 主婦 副業 at 2011年09月21日 17:53
おっぱいマニアに大人気!あなたのおっぱいへのこだわりがズバリ暴かれる!巨乳や貧乳、乳首の色や形・大きさまであなたが望む最高のおっぱいを見てみませんか?
Posted by おっぱい at 2011年09月22日 11:58
三冠牝馬誕生なるか!第15回 秋華賞!!コースが複雑なため、たとえ強豪馬であっても逆転できない状況がしばしば見られる秋華賞。予想は困難を極めますが、過去のデータを分析していくと面白いことになるかも?
Posted by 秋華賞 at 2011年10月05日 04:54
第62回朝日杯フューチュリティステークス オッズ・枠順からレース展開が見える!人気に騙されない買い方伝授します!人気馬で決めて本当にイイのか!?
Posted by 朝日杯フューチュリティステークス at 2011年12月07日 06:54
中山大障害の攻略方をこっそり教えちゃいます。JRAの裏話も混ぜつつ在り来たりな競馬予想サイトではない個性溢れるコンテツが満載です。
Posted by 中山大障害 2011 at 2011年12月14日 04:56
まだそこにハマってるの?って今鼻で笑ったそこの君!!ばかにしたらいけないよ!W大晦日の夜何かが始まるんだって!何が始まるんだろ?ワクワクだよね
Posted by モバゲー at 2011年12月30日 03:39
まだまだ呟いて火傷しちゃった有名人が後を絶たない。次誰が炎上するかリスト作ってみたんで一緒に共感してくれ
Posted by ツイッター at 2012年01月21日 03:15
グリーングリーン!!きゃほきゃほわはははは。CMの長●大好きなんだけど!!明日無料でマル秘画像が無料でダウンロードできるらしい!!どうしようどうしよう★そっこうアクセスしなきゃファン失格っしょ
Posted by グリー at 2012年01月22日 01:54
コメントを書く
お名前: [必須入力]

メールアドレス:

ホームページアドレス:

コメント: [必須入力]

認証コード: [必須入力]


※画像の中の文字を半角で入力してください。
この記事へのトラックバックURL
http://blog.seesaa.jp/tb/23588668

この記事へのトラックバック

Physical Review

Recent stories in Physical Review Focus

×

この広告は1年以上新しい記事の投稿がないブログに表示されております。