生成1024个采样的1024khz的啁啾。啁啾的初始频率为50 kHz，在采样结束时达到100 kHz。添加高斯白噪声，使信噪比为40 dB。

nSamp = 1024;Fs = 1024e3;信噪比= 40;t = (0:nSamp-1)'/Fs;x = chirp(t,50e3,nSamp/Fs,100e3);x = x + randn(size(x))*std(x)/db2mag(SNR);

估计信号的3-dB带宽，并在功率谱密度(PSD)图上标注它。

powerbw (x, Fs)

Ans = 4.4386e+04

产生另一个啁啾。指定初始频率为200khz，最终频率为300khz，振幅为第一个信号的两倍。添加高斯白噪声。

x2 = 2*chirp(t,200e3,nSamp/Fs,300e3);x2 = x2 + randn(size(x2))*std(x2)/db2mag(SNR);

将这些啁啾串联起来，产生一个双通道信号。估计每个通道的3db带宽。

y = powerbw([x x2]，Fs)

y =1×2104× 4.4386 9.2208

在psd图上标注两个通道的3db带宽。

powerbw (x2 [x]、Fs);

把这两个通道加起来形成一个新信号。绘制PSD，标注3db带宽。

powerbw (x + x2, Fs)

Ans = 9.2243e+04

生成1024khz采样的100.123 kHz正弦波的1024个样本。添加高斯白噪声，使信噪比为40 dB。重置随机数发生器可重现的结果。

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

使用周期图函数计算信号的功率谱密度(PSD)。指定与信号长度相同且形状因子为38的Kaiser窗口。估计信号的3-dB带宽，并在PSD图上标注它。

[Pxx,f] = periodogram(x,kaiser(nSamp,38)，[]，Fs);powerbw (Pxx f);

生成另一个正弦波，频率为257.321 kHz，振幅是第一个正弦波的两倍。添加高斯白噪声。

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

串联正弦波以产生一个双通道信号。估计每个通道的PSD，并使用结果来确定3-dB带宽。

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

y =1×2103.× 3.1753 3.3015

在psd图上标注两个通道的3db带宽。

powerbw (Pyy组,f);

把这两个通道加起来形成一个新信号。预估PSD，标注3db带宽。

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

生成一个信号，其PSD类似于88阶带通FIR滤波器的频率响应，具有归一化的截止频率 $$0．25\pi$$rad /样本和 $$0．45\pi$$rad /样品。

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

计算信号占用的3db带宽。指定中间波段的平均功率作为参考电平 $$0．2\pi$$rad /样本和 $$0．6\pi$$rad /样品。绘制PSD并标注带宽。

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

输出带宽，其下限和上限，以及带宽功率。的采样率 $$2\pi$$就相当于不设定速率。

[bw,弗洛,fhi,权力]= powerbw (d, 2 *π,(0.2 - 0.6)*π);流('bw = %。3.f*pi, flo = %.3f*pi, fhi = %.3f*pi \n'，.．.[bw Flo fhi]/pi

Bw = 0.200*pi, Flo = 0.250*pi, fhi = 0.450*pi

流('power = %。总额的1f%% '、电力/ bandpower (d) * 100)

电力=总数的96.9%

添加第二个具有标准化截止频率的通道 $$0．5\pi$$rad /样本和 $$0．8\pi$$Rad /sample和振幅是第一个通道的十分之一。

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

计算双通道信号的6db带宽。指定频谱的最大功率电平作为参考电平。

powerbw (d, [] [], 6);

输出每个通道的6db带宽和上下限。

[bw,flo,fhi] = powerbw(d，[]，[]，6);BDS = [bw;flo;fhi];流(“1:bw = %。3.f*pi, flo = %.3f*pi, fhi = %.3f*pi \n'bds(: 1) /π)

一:bw = 0.198*pi, flo = 0.252*pi, fhi = 0.450*pi

流(二:bw = %。3.f*pi, flo = %.3f*pi, fhi = %.3f*pi \n'bds(: 2) /π)

二:bw = 0.294*pi, flo = 0.503*pi, fhi = 0.797*pi