产生1024采样啁啾在1024 kHz。啁啾的初始频率为50 kHz，采样结束时达到100 kHz。加入高斯白噪声，使信噪比为40 dB。重置随机数生成器以获得可重复的结果。

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

估计啁啾的平均频率。绘制功率谱密度（PSD）并注释平均频率。

meanfreq (x, Fs)

ans = 7.5032 e + 04

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

x2=2*chirp（t，200e3，nSamp/Fs，300e3）；x2=x2+randn（大小（x2））*std（x2）/db2mag（信噪比）；

将啁啾连接起来产生双通道信号。估计每个信道的平均频率。

y=平均频率（[x x2]，Fs）

y=1×2105.× 0.7503 2.4999

绘制两个通道的PSD并注释其平均频率。

平均频率（[x2]，Fs）；

添加两个通道以形成新信号。绘制PSD并注释平均频率。

平均频率（x+x2，Fs）

ans=2.1496e+05

生成1024个采样100.123 kHz的正弦信号。加入高斯白噪声，使信噪比为40 dB。重置随机数生成器以获得可重复的结果。

nSamp=1024；Fs=1024e3；信噪比=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);meanfreq (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 = meanfreq (Pyy组,f)

y=1×2105.× 1.0013 2.5732

在PSD图上标注两个通道的平均频率。

meanfreq (Pyy组,f);

添加两个通道以形成新信号。估计PSD并注释平均频率。

[Pzz，f]=周期图（x+x2，kaiser（nSamp，38），[]，Fs）；平均频率（Pzz，f）；

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

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

计算之间信号的平均频率 $$0．3.\pi$$rad/样本和 $$0．6.\pi$$rad/样本。绘制PSD并注释平均频率和测量间隔。

meanfreq (d,[],[0.3 - 0.6] *π);

输出测量间隔的平均频率和频带功率。指定的采样率 $$2.\pi$$相当于不设定速率。

[mnf，功率]=平均频率（d，2*pi，[0.30.6]*pi）；fprintf('意味着= %。3.F*pi, power = %.1f%% of total \n',...mnf/pi，功率/频带功率（d）*100）

平均值= 0.373*pi，功率= 75.6%

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

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

计算之间信号的平均频率 $$0．3.\pi$$rad/样本和 $$0．9\pi$$rad /样品。绘制PSD，标注每个通道的平均频率和测量间隔。

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

输出每个通道的平均频率。除以 $$\pi$$．

mpf = meanfreq(d，[]，[0.3 0.9]*pi)/pi

延长=1×20.3730 0.6500