tmp617.m

信号分析与处理的一些源代码

%MATLAB PROGRAM 6-17
%Create a Rectangular window
clf
Nfft=512;
for i=1:4
  switch i
  case 1  
     Nwin=10;
  case 2
     Nwin=20;
  case 3
     Nwin=50;
  case 4
     Nwin=100;
  end   
   
w=boxcar(Nwin);
[y,f]=freqz(w,1,Nfft);
mag=abs(y);
posplot=['22'int2str(i)];
subplot(posplot)
plot(f/pi,20*log10(mag/max(mag)));
xlabel('Normalized Frequency');
ylabel(' Magnitude dB');
stext=['N='int2str(Nwin)];
title(stext)
grid
end