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

nSamp = 1024;Fs = 1024e3;信噪比= 40;rng默认的t = (0:nSamp-1)'/Fs;x = chirp(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

生成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窗口。估计信号的中值频率，并在PSD图上标注它。

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

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

X2 = 2*sin(2*pi*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);流('Mean = %。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