创建一个信号组成的三个嘈杂的正弦曲线和一个唧唧喳喳,取样0.1秒的200千赫。正弦信号的频率是1 kHz, 10 kHz, 20 kHz。正弦信号有不同的振幅和噪音水平。无声的唧唧声频率,始于20 kHz和线性增加到30千赫采样期间。
Fs = 200年e3;Fc = 10 20 [1]“* 1 e3;Ns = 0.1 * Fs;t = (0: Ns-1) / Fs;x = [1 1/10 10] * sin(2 *π* Fc * t) + (1/200 1/2000 1/20) * randn (Ns);x = x +唧唧喳喳(t, 20 e3, t(结束),30 e3);
计算韦尔奇功率谱估计和maximum-hold minimum-hold光谱的信号。策划的结果。
[pxx f] = pwelch (x, [] [], [], Fs);pmax = pwelch (x, [] [], [], Fs,“maxhold”);pmin = pwelch (x, [] [], [], Fs,“minhold”);情节(f / 1000, pow2db (pxx))在情节(f / 1000, pow2db ([pmax pmin]),“:”)举行从包含(的频率(赫兹))ylabel (“PSD (dB / Hz)”)传说(“pwelch”,“maxhold”,“minhold”网格)
重复这个过程,这段时间计算为中心的功率谱估计。
[pxx f] = pwelch (x, [] [], [], Fs,“中心”,“权力”);pmax = pwelch (x, [] [], [], Fs,“maxhold”,“中心”,“权力”);pmin = pwelch (x, [] [], [], Fs,“minhold”,“中心”,“权力”);情节(f / 1000, pow2db (pxx))在情节(f / 1000, pow2db ([pmax pmin]),“:”)举行从包含(的频率(赫兹))ylabel (“权力(dB)”)传说(“pwelch”,“maxhold”,“minhold”)标题(集中功率谱估计的网格)