主题:[讨论]抑制频谱泄露是这样的吗?
t = 0:0.001:0.6;
x = 3*sin(2*pi*50*t)+4*sin(2*pi*120*t);
y = x + 2*randn(size(t));
subplot(311)
plot(1000*t(1:50),y(1:50)) % 初始信号图
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
Y = fft(y,512);
Pyy=abs(Y)*2/512;
f = 1000*(0:256)/512;
subplot(312)
plot(f,Pyy(1:257)) %未经处理的信号频谱图
title('Frequency content of y')
xlabel('frequency (Hz)')
Yh=fft(y(1:512).*hanning(512)',512); % 对信号进行加窗,防止频谱泄露?
Pyy=abs(Yh)*2/512;
f = 1000*(0:256)/512;
subplot(313)
plot(f,Pyy(1:257)) %加窗后的频谱图,是否也是抑制了频谱泄露的频谱图?
title('Frequency content of y')
xlabel('frequency (Hz)')
是这样抑制频谱泄露的么?
x = 3*sin(2*pi*50*t)+4*sin(2*pi*120*t);
y = x + 2*randn(size(t));
subplot(311)
plot(1000*t(1:50),y(1:50)) % 初始信号图
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
Y = fft(y,512);
Pyy=abs(Y)*2/512;
f = 1000*(0:256)/512;
subplot(312)
plot(f,Pyy(1:257)) %未经处理的信号频谱图
title('Frequency content of y')
xlabel('frequency (Hz)')
Yh=fft(y(1:512).*hanning(512)',512); % 对信号进行加窗,防止频谱泄露?
Pyy=abs(Yh)*2/512;
f = 1000*(0:256)/512;
subplot(313)
plot(f,Pyy(1:257)) %加窗后的频谱图,是否也是抑制了频谱泄露的频谱图?
title('Frequency content of y')
xlabel('frequency (Hz)')
是这样抑制频谱泄露的么?