创建由三个嘈杂的正弦曲线和啁啾组成的信号,以200kHz采样0.1秒。正弦波的频率为1 kHz,10 kHz和20 kHz。正弦曲线具有不同的幅度和噪声水平。无噪声啁啾具有以20 kHz开头的频率,并在采样期间将线性地增加至30 kHz。
fs = 200e3;FC = [110 20]'* 1E3;ns = 0.1 * fs;t =(0:ns-1)/ fs;X = [1 1/10 10] * SIN(2 * PI * FC * T)+ [1/200 1/2000 1/20] * RANDN(3,NS);x = x + chirp(t,20e3,t(端),30e3);
计算信号的Welch PSD估计和最大保持和最小保持光谱。绘制结果。
[pxx,f] = pwelch(x,[],[],[],fs);pmax = pwelch(x,[],[],[],fs,'maxhold');pmin = pwelch(x,[],[],[],fs,'minhold');绘图(F / 1000,POW2DB(PXX))保持在绘图(F / 1000,POW2DB([PMAX PMIN]),':') 抓住离开Xlabel('频率(kHz)')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]),':') 抓住离开Xlabel('频率(kHz)')ylabel('Power(DB)') 传奇('pwelch'那'maxhold'那'minhold') 标题(“中心的功率谱估计”) 网格