Создайте сигнал, состоящий из трех шумных синусоидов и щебета, дискретизированный при 200 кГц в течение 0,1 секунды. Частоты синусоидов составляют 1 кГц, 10 кГц и 20 кГц. Синусоиды имеют различные амплитуды и уровни шума. Бесшумный щебет имеет частоту, которая начинается с 20 кГц и увеличивается линейно до 30 кГц во время дискретизации.
Fs = 200e3; Fc = [1 10 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(end),30e3);
Вычислите оценку PSD Welch и спектры максимального удержания и минимального удержания сигнала. Постройте график результатов.
[pxx,f] = pwelch(x,[],[],[],Fs); pmax = pwelch(x,[],[],[],Fs,'maxhold'); pmin = pwelch(x,[],[],[],Fs,'minhold'); plot(f/1000,pow2db(pxx)) hold on plot(f/1000,pow2db([pmax pmin]),':') hold off xlabel('Frequency (kHz)') ylabel('PSD (dB/Hz)') legend('pwelch','maxhold','minhold') grid
Повторите процедуру, на этот раз вычисляя оценки центрированного спектра степени.
[pxx,f] = pwelch(x,[],[],[],Fs,'centered','power'); pmax = pwelch(x,[],[],[],Fs,'maxhold','centered','power'); pmin = pwelch(x,[],[],[],Fs,'minhold','centered','power'); plot(f/1000,pow2db(pxx)) hold on plot(f/1000,pow2db([pmax pmin]),':') hold off xlabel('Frequency (kHz)') ylabel('Power (dB)') legend('pwelch','maxhold','minhold') title('Centered Power Spectrum Estimates') grid
