1024千赫でサンプリングされたチャープのサンプルを1024個生成します。チャープの初期周波数は50千赫でサンプリングの最後には100 kHzに到達します。S / N比が40 dBとなるホワイトガウスノイズを付加します。再現可能な結果が必要な場合は,乱数発生器をリセットします。

nSamp = 1024;Fs = 1024年e3;信噪比= 40;rng默认的t = (0: nSamp-1) / Fs;x =唧唧声(t, 50 e3, nSamp / Fs 100 e3);x = x + randn(大小(x)) *性病(x) / db2mag(信噪比);

信号の占有帯域幅を推定し,パワースペクトル密度(PSD)のプロット上で推定値に注釈を付けます。

obw (x, Fs)

ans = 5.5377 e + 04

別のチャープを生成します。初期周波数200 kHz,最終周波数300 kHzおよび最初の信号の振幅の2倍の振幅を指定します。ホワイトガウスノイズを付加します。

x2 = 2 *唧唧声(t, 200年e3, nSamp / Fs 300 e3);x2 = x2 + randn(大小(x2)) *性病(x2) / db2mag(信噪比);

チャープを連結して2チャネル信号を生成します。各チャネルの占有帯域幅を推定します。

y = obw([x x]， f)

y =1×2105×0.5538 - 1.0546

PSDのプロット上で2つのチャネルの占有帯域幅に注釈を付けます。

obw (x2 [x]、Fs);

2つのチャネルを追加して新しい信号を作成します。PSDをプロットして占有帯域幅に注釈を付けます。

obw (x + x2, Fs);

1024千赫でサンプリングされた100.123 kHzの正弦波のサンプルを1024個生成します。S / N比が40 dBとなるホワイトガウスノイズを付加します。再現可能な結果が必要な場合は,乱数発生器をリセットします。

nSamp = 1024;Fs = 1024年e3;信噪比= 40;rng默认的t = (0: nSamp-1) / Fs;x =罪(2 *π* t * 100.123 e3);x = x + randn(size(x))*std(x)/db2mag(SNR);

周期图を使用して,信号のパワースペクトル密度(PSD)を計算します。信号と同じ長さおよび形状係数38をもつカイザーウィンドウを指定します。信号の占有帯域幅を推定し,PSDのプロット上で推定値に注釈を付けます。

[Pxx f] =周期图(x, kaiser (nSamp 38), [], Fs);obw (Pxx f);

別の正弦波を生成します。この正弦波の周波数は257.321千赫で振幅は最初の正弦波の2倍です。ホワイトガウスノイズを付加します。

x2 = 2 * sin(2 *π* t * 257.321 e3);x2 = x2 + randn(size(x2))*std(x2)/db2mag(SNR);

正弦波を連結して2チャネル信号を生成します。各チャネルのPSDを推定し,結果から占有帯域幅を特定します。

(Pyy组,f) =周期图(x2 [x], kaiser (nSamp 38), [], Fs);y = obw (Pyy组,f)

y =1×2103.×7.2001 - 7.3777

PSDのプロット上で2つのチャネルの占有帯域幅に注釈を付けます。

obw (Pyy组,f);

2つのチャネルを追加して新しい信号を作成します。PSDを推定して占有帯域幅に注釈を付けます。

[Pzz f] =周期图(x + x2, kaiser (nSamp 38), [], Fs);obw (Pzz f);

正規化されたカットオフ周波数が $$0．25\pi$$ラジアン/サンプルおよび $$0．45\pi$$ラジアン/サンプルである88次のバンドパス冷杉フィルターの周波数応答に類似したPSDをもつ信号を生成します。

D = fir1(88，[0.25 0.45]);

$$0．2\pi$$ラジアン/サンプルと $$0．6\pi$$ラジアン/サンプル間の信号の占99%有帯域幅を計算します。PSDをプロットして占有帯域幅と測定間隔に注釈を付けます。

obw (d,[],[0.2 - 0.6] *π);

占有帯域幅,その上限と下限および占有帯域パワーを出力します。サンプルレート $$2\pi$$の指定は,レートを設定しない場合と同じです。

[bw,flo,fhi,power] = obw(d,2*pi，[0.2 0.6]*pi);流(' bw = %。3.f*pi, flo = %.3f*pi, fhi = %.3f*pi \n', (bw弗洛fhi) / pi)

Bw = 0.217*pi, Flo = 0.240*pi, fhi = 0.458*pi

流(“权力= %。1 f % %的总、电力/ bandpower (d) * 100)

功率=总功率的99.0%

正規化されたカットオフ周波数が $$0．5\pi$$ラジアン/サンプルと $$0．8\pi$$ラジアン/サンプルであり,最初のチャネルの1/10の振幅をもつ2番目のチャネルを追加します。

D = [D;fir1(88，[0.5 0.8])/10]';

$$0．3.\pi$$ラジアン/サンプルと $$0．9\pi$$ラジアン/サンプル間の信号の占50%有帯域幅を計算します。PSDをプロットして占有帯域幅と測定間隔に注釈を付けます。

obw (d,[],[0.3 - 0.9] *π,50);

各チャネルの占有帯域幅を出力します。 $$\pi$$で除算します。

Bw = obw(d，[]，[0.3 0.9]*pi,50)/pi

bw =1×20.0705 - 0.1412