最新のリリースでは、このページがまだ翻訳されていません。 このページの最新版は英語でご覧になれます。
medfreq
周波数中央値
構文
freq = medfreq(x)freq = medfreq(x,fs)freq = medfreq(pxx,f)freq = medfreq(sxx,f,rbw)freq = medfreq(___,freqrange)[freq,power]
= medfreq(___)medfreq(___)説明
出力引数なしで medfreq(___) を使用すると、PSD またはパワー スペクトルがプロットされて周波数中央値に注釈が付けられます。
例
チャープの周波数中央値
1024 kHz でサンプリングされたチャープのサンプルを 1024 個生成します。チャープの初期周波数が 50 kHz で、サンプリングの最後には 100 kHz に到達するようにチャープを指定します。S/N 比が 40 dB となるホワイト ガウス ノイズを付加します。再現可能な結果が必要な場合は、乱数発生器をリセットします。
nSamp = 1024;
Fs = 1024e3;
SNR = 40;
rng default
t = (0:nSamp-1)'/Fs;
x = chirp(t,50e3,nSamp/Fs,100e3);
x = x+randn(size(x))*std(x)/db2mag(SNR);チャープの周波数中央値を推定します。パワー スペクトル密度 (PSD) をプロットして周波数中央値に注釈を付けます。
medfreq(x,Fs)
ans = 7.4998e+04
別のチャープを生成します。初期周波数 200 kHz、最終周波数 300 kHz および最初の信号の振幅の 2 倍の振幅を指定します。ホワイト ガウス ノイズを付加します。
x2 = 2*chirp(t,200e3,nSamp/Fs,300e3); x2 = x2+randn(size(x2))*std(x2)/db2mag(SNR);
チャープを連結して 2 チャネル信号を生成します。各チャネルの周波数中央値を推定します。
y = medfreq([x x2],Fs)
y = 1×2
105 ×
0.7500 2.4999
2 つのチャネルの PSD をプロットして周波数中央値に注釈を付けます。
medfreq([x x2],Fs);
2 つのチャネルを追加して新しい信号を作成します。PSD をプロットして周波数中央値に注釈を付けます。
medfreq(x+x2,Fs)
ans = 2.3756e+05
正弦波の周波数中央値
1024 kHz でサンプリングされた 100.123 kHz の正弦波のサンプルを 1024 個生成します。S/N 比が 40 dB となるホワイト ガウス ノイズを付加します。再現可能な結果が必要な場合は、乱数発生器をリセットします。
nSamp = 1024;
Fs = 1024e3;
SNR = 40;
rng default
t = (0:nSamp-1)'/Fs;
x = sin(2*pi*t*100.123e3);
x = x + randn(size(x))*std(x)/db2mag(SNR);periodogram を使用して、信号のパワー スペクトル密度 (PSD) を計算します。信号と同じ長さおよび形状係数 38 をもつカイザー ウィンドウを指定します。信号の周波数中央値を推定し、PSD のプロット上でそれに注釈を付けます。
[Pxx,f] = periodogram(x,kaiser(nSamp,38),[],Fs); medfreq(Pxx,f);
別の正弦波を生成します。この正弦波の周波数は 257.321 kHz で、振幅は最初の正弦波の 2 倍です。ホワイト ノイズを付加します。
x2 = 2*sin(2*pi*t*257.321e3); x2 = x2 + randn(size(x2))*std(x2)/db2mag(SNR);
正弦波を連結して 2 チャネル信号を生成します。各チャネルの PSD を推定し、結果から周波数中央値を特定します。
[Pyy,f] = periodogram([x x2],kaiser(nSamp,38),[],Fs); y = medfreq(Pyy,f)
y = 1×2
105 ×
1.0012 2.5731
PSD のプロット上で 2 つのチャネルの周波数中央値に注釈を付けます。
medfreq(Pyy,f);
2 つのチャネルを追加して新しい信号を作成します。PSD を推定して周波数中央値に注釈を付けます。
[Pzz,f] = periodogram(x+x2,kaiser(nSamp,38),[],Fs); medfreq(Pzz,f);
帯域制限された信号の周波数中央値
正規化されたカットオフ周波数が 
ラジアン/サンプルおよび 
ラジアン/サンプルの 88 次のバンドパス FIR フィルターの周波数応答に類似した PSD をもつ信号を生成します。
d = fir1(88,[0.25 0.45]);
ラジアン/サンプルと 
ラジアン/サンプル間の信号の周波数中央値を計算します。PSD をプロットして周波数中央値と測定間隔に注釈を付けます。
medfreq(d,[],[0.3 0.6]*pi);
測定間隔の周波数中央値と帯域パワーを出力します。サンプルレート 
の指定は、レートを設定しない場合と同じです。
[mdf,power] = medfreq(d,2*pi,[0.3 0.6]*pi); fprintf('Mean = %.3f*pi, power = %.1f%% of total \n', ... mdf/pi,power/bandpower(d)*100)
Mean = 0.371*pi, power = 77.4% of total
正規化されたカットオフ周波数が 
ラジアン/サンプルと 
ラジアン/サンプルであり、最初のチャネルの 1/10 の振幅をもつ 2 番目のチャネルを追加します。
d = [d;fir1(88,[0.5 0.8])/10]';
ラジアン/サンプルと 
ラジアン/サンプル間の信号の周波数中央値を計算します。PSD をプロットして各チャネルの周波数中央値と測定間隔に注釈を付けます。
medfreq(d,[],[0.3 0.9]*pi);
各チャネルの周波数中央値を出力します。
で除算します。
mdf = medfreq(d,[],[0.3 0.9]*pi)/pi
mdf = 1×2
0.3706 0.6500
入力引数
出力引数
参考
findpeaks | meanfreq | periodogram | plomb | pwelch
