产生1024采样啁啾在1024 kHz。啁啾的初始频率为50 kHz，采样结束时达到100 kHz。加入高斯白噪声，使信噪比为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

生成另一个唧唧喳喳。指定初始频率为200khz，最终频率为300khz，振幅为第一个信号的两倍。加入白色高斯噪声。

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

将啁啾连接起来产生双通道信号。估计每个信道的占用带宽。

y = obw([x x]， f)

y =1×2105×0.5538 - 1.0546

在psd的图上标注两个信道的占用带宽。

obw (x2 [x]、Fs);

添加两个通道，形成一个新的信号。绘制PSD图并标注占用的带宽。

obw (x + x2, Fs);

生成1024个采样100.123 kHz的正弦信号。加入高斯白噪声，使信噪比为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的Kaiser窗口。估计信号占用的带宽，并将其标注在PSD的图上。

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

生成另一个正弦信号，这个频率为257.321 kHz，振幅是第一个正弦信号的两倍。加入白色高斯噪声。

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

连接正弦波以产生双通道信号。估计每个信道的PSD，并使用结果来确定占用带宽。

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

y =1×2103.×7.2001 - 7.3777

在psd的图上标注两个信道的占用带宽。

obw (Pyy组,f);

添加两个通道，形成一个新的信号。估计PSD并标注占用带宽。

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

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

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

计算99%之间占用的信号带宽 $$0．2\pi$$rad /样本和 $$0．6\pi$$rad /样品。绘制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$$rad /样本和 $$0．8\pi$$Rad /采样，振幅为第一个通道的十分之一。

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

计算50%占用的带宽之间的信号 $$0．3.\pi$$rad /样本和 $$0．9\pi$$rad /样品。绘制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