スペクトル解析パラメトリックと非パラメトリック手法

G

Guest

Guest
こんにちは
、 すべて今Matlabの初心者で
、 ちょっと助けが
必要私パラメトリックと非パラメトリック手法を比較する
iとし
、 コードはパラメトリックまたは非パラメトリック法..とこの信号のにPSDを取ってより良い方法を参照してください...私は単なる例です..希望私は自分自身を明確に必要な信号のコードが必要になります
すべて事前にありがとう!

 
Matlabのpsdにspectrum.musicまたは関数を使用することができます。
Matlabのを入力する必要がある
spectrum.music
ドクまたは
ドクpsd

 
こんなに私はこのツールボックス
を発見 ...しかし僕は今...と思わMatlabので初心者には少し複雑なunderstant .. thatsなぜ..ガイドさせる必要があると語った
私は
、ツールボックス のにPSD信号... ...私は必要のような例を簡単に見ることができるたびには
、 N (サンプル
) と周波数を変更し
、 必要がある!
あなたの応答を
『 THX ! !

 
例えばこのようないくつかを簡単にすることができます:
レン= 1000 ;
fd = 10 ;
トン= 0:1 / fd : len/fd-1/fd ;
秒=罪(トン* 3 * 2 *円周率) ( 1 、レン)
randn ;
( psd 、
fd ) レン秒

 
ū salpツールボックスと同じ例を挙げることができる...私は
、 与えられたコードを用意して:
[イ、トン] = (ベータ版
は 、
F 、寸法、 <OSF> 、 <SNR> )
sg_cissoid%
%引数: -ベータ版:ベクター複雑な振幅
を% -金:ベクトル周波数
の 値
を ( -0.5 <金< 0.5 )
% -寸法:サンプルの数
% - OSF :オーバサンプリング係数(デフォルト: OSF = 1 )
% -のSNR :信号
の ノイズ比[ dBの] (デフォルト: SN比= infファイル)


の 出力: -イ:ベクトル信号のサンプル
で% -トン:時間軸

%例: sg_cissoid ( [ 1 0.3 0.5 0.7 ] 、 [ 0.1 0.2 0.25 0.48 ] 、 100 、 5 、 30 ) ;


、 適切なpsd
電話: [ PSD
は 、
F ] = psd_welch (イ、 <Win> 、 <N_Overlap> 、 <N_FFT> 、 <fs> )

%引数: -イ:ベクトル信号のサンプル
で% -ウィン: K_seg信号に分かれています重複
長さは
、 長さ
の セグメント
に 相当の割合
%は
、 ウィンドウとウィンドウを各セグメント
% (デフォルト:ハミングウィンドウで8セグメント)
% - N_Overlap :重複するサンプルの数のセグメントに
%セグメント(デフォルト: 50 %重複)
% - N_FFT : FFTのポイント数(デフォルト: N_FFT = 256 )
% -のFS : Hzでサンプリング周波数(デフォルト: fsの= 1 Hzの)


の 出力: - PSD :パワースペクトル密度
% -金:周波数軸

%例:イ= ( [ 1 0.5 0.25 ] 、 [ 0.2 0.3 0.4
sg_cissoid ] 、 100 、 1 、 10 ) ;
% psd_welch (イ、ハン( 50 ) 、 49 、 1024 ) ;

 

Welcome to EDABoard.com

Sponsor

Back
Top