Why i can't plot a spectrum of a rectangular window function?
조회 수: 15 (최근 30일)
이전 댓글 표시
Hi everyone! I have a problem, so could you explain to me? I knew that, to plot a spectrum of any windows function we can use "wvtool" command in Matlab or you can find a code to do this on google. But, why i can do this with my code? My plot (use my code) is different with that uses "wvtool". P/s:i understood that spectrum of any window (rect, hanning....) function is FFT of the windows. Is that right or wrong?. And my code:
t=0:100;
N=length(t)-1;
rect=rectwin(N+1); % create a rectangular windows function
taxis=-N/2:N/2; % time axis
plot(taxis,rect); % Plot rectangular windows in time domain.
fftrect=fft(rect); % calculate spectrum of the rectangular window
figure;
plot(taxis,20*log10(abs(fftrect))) % Plot rectangular windows in frequency domain.
wvtool(rect) % Plot rectangular windows in frequency domain use "wvtool" to compare
댓글 수: 0
채택된 답변
Greg Dionne
2018년 1월 29일
If you line up the frequency points, you'll see that your FFT is showing the zeros only. If you crank up the number of points in the FFT you are computing you'll see a similar picture:
plot(20*log10(abs(fft(rect,N*10))))
댓글 수: 3
Greg Dionne
2018년 1월 30일
I used 10 so you could see the lobes that are non-null. You could use any positive whole number.
추가 답변 (1개)
Honglei Chen
2018년 1월 29일
There is a mistake in your plotting command. Because the spectrum is in frequency domain, you need to have a frequency axis instead of time axis.
As to the curve it self, your concept is correct. However, the one in wvtool is computed with 1024 points, you can reproduce it using the following code snippet.
nfft = 1024;
fftrect=fft(rect,nfft); % calculate spectrum of the rectangular window
faxis = 0:2/nfft:2-2/nfft;
figure;
plot(faxis(1:nfft/2),20*log10(abs(fftrect(1:nfft/2)))) % Plot rectangular windows in frequency domain.
HTH
댓글 수: 3
Honglei Chen
2018년 1월 30일
That's how the frequency axis is defined in normalized frequency domain. Essentially it divides 0 to 2 to nfft equal bins.
HTH
참고 항목
카테고리
Help Center 및 File Exchange에서 Windows에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!