生成1024个样本的啁啾采样在1024千赫。指定啁啾，使其初始频率为50khz，并在采样结束时达到100khz。添加高斯白噪声，使信噪比为40 dB。重置随机数生成器以获得可重复的结果。

nSamp = 1024;Fs = 1024e3;信噪比= 40;rng默认的t = (0: nsample -1)'/Fs;x =啁啾(t,50e3,nSamp/Fs,100e3);x = x+randn(size(x))*std(x)/db2mag(SNR);

估计啁啾的中值频率。绘制功率谱密度(PSD)并标注中位数频率。

medfreq (x, Fs)

Ans = 7.4998e+04

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

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

把这些啁啾连在一起就能产生双通道信号。估计每个通道的中位数频率。

y = medfreq([x x2]，Fs)

y =1×2105× 0.7500 2.4999

绘制两个通道的psd，并标注它们的中位数频率。

medfreq (x2 [x]、Fs);

把这两个通道加起来形成一个新信号。绘制PSD并标注中位数频率。

medfreq (x + x2, Fs)

Ans = 2.3756e+05

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

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

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

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

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

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

把正弦信号连接起来，就会产生双通道信号。估计每个通道的PSD，并利用结果确定中位数频率。

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

y =1×2105× 1.0012 2.5731

在psd图上标注两个通道的中位数频率。

medfreq (Pyy组,f);

把这两个通道加起来形成一个新信号。估计PSD并标注中位数频率。

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

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

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

计算信号之间的中值频率 $$0．3.\pi$$rad /样本和 $$0．6\pi$$rad /样品。绘制PSD图，标注中位数频率和测量间隔。

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

输出测量间隔的中位数频率和带功率。的抽样速率 $$2\pi$$相当于不设定利率。

[mdf,power] = medfreq(d,2*pi，[0.3 0.6]*pi);流(“平均值= %。”3.f*pi, power = %.1f%% of total \n'，.．.mdf /π,电源/ bandpower (d) * 100)

平均值= 0.371*pi，功率=总功率的77.4%

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

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

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

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

输出每个通道的中位数频率。除以 $$\pi$$．

MDF = medfreq(d，[]，[0.3 0.9]*pi)/pi

mdf =1×20.3706 - 0.6500