% Script : plot the LF model waveform and the spectrum
%

 figure(Gwdis_f);
 
 lens=100;
 noise=0;

 gp=lfpuls(Tp,Te,Ta,Tc,Ee,lens,noise);

 subplot(121);
   plot(gp,'k');
   xlabel('samples');
   title('Quasi-periodic Waveform');
 subplot(122);
   fgp=fft(gp);
   fgp=20*log(abs(fgp(1:50)));
   plot(0:5000/49:5000,fgp,'k');
   xlabel('frequency (Hz)');
   ylabel('dB');
   title('Quasi-periodic Spectrum')
 clear fgp;