生成1024个样本的唧唧声在1024千赫采样。指定的唧唧声,有一个初始频率50千赫,达到100千赫采样结束的时候。添加高斯白噪声,信噪比是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)和注释中值频率。

medfreq (x, Fs)

ans = 7.4998 e + 04

生成另一个唧唧喳喳。指定一个初始频率的200 kHz,最后300 kHz的频率,振幅,第一信号的两倍。添加高斯白噪声。

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

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

y = medfreq (x2 [x], Fs)

y =1×2105×0.7500 - 2.4999

情节PSDs的两个渠道和注释中值频率。

medfreq (x2 [x]、Fs);

添加两个渠道,形成一个新的信号。情节PSD和注释中值频率。

medfreq (x + x2, Fs)

ans = 2.3756 e + 05

产生1024样品100.123 kHz的正弦信号在1024千赫采样。添加高斯白噪声,信噪比是40 dB。重置的随机数字生成器可重复的结果。

nSamp = 1024;Fs = 1024年e3;信噪比= 40;rng默认的t = (0: nSamp-1) / Fs;x =罪(2 *π* t * 100.123 e3);x = x + randn(大小(x)) *性病(x) / db2mag(信噪比);

使用周期图函数计算功率谱密度(PSD)的信号。指定一个Kaiser窗长度相同的信号和38的形状系数。估计信号的中值频率和注释PSD的情节。

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

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

x2 = 2 * sin(2 *π* t * 257.321 e3);x2 = x2 + randn(大小(x2)) *性病(x2) / db2mag(信噪比);

将正弦信号产生一个双通道的信号。估计每个通道的PSD和使用结果来确定平均频率。

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

y =1×2105×1.0012 - 2.5731

注释中值频率的两个渠道PSDs的阴谋活动。

medfreq (Pyy组,f);

添加两个渠道,形成一个新的信号。估计PSD和注释中值频率。

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

生成一个信号的PSD像88阶带通滤波器的频率响应与归一化截止频率 $$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、电力)= medfreq (d, 2 *π,(0.3 - 0.6)*π);流('意味着= %。3f*pi, power = %.1f%% of total \n',…mdf /π,电源/ bandpower (d) * 100)

意味着权力= = 0.371 *π,总额的77.4%

添加第二个通道与归一化截止频率 $$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和注释中值频率和测量时间间隔。

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

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

mdf = medfreq (d,[],[0.3 - 0.9] *π)/π

mdf =1×20.3706 - 0.6500