my.asv

MATLAB的部分源代码,与MATLAB数字信号处理与应用一书相配套

 N=512; 
 [fm,am,iflaw]=doppler(N,200,65,10,50); 
 figure;

 plot(real(am.*fm));  
 title( 'Signal in time')
 xlabel('Time');
 ylabel('real part');
 figure;
 
 plot(iflaw); 
 title( 'Instantaneous frequency estiwmation')
 xlabel('Time');
 ylabel('Nomized Frequency (Hz)');
 
 sig=sigmerge(am.*fm,noisecg(N),20);
 figure;
 title( 'Instantaneous frequency estiwmation')
 xlabel('Time');
 ylabel('Nomized Frequency (Hz)');
 
 [ifhat,t]=instfreq(sig,11:502,10);
 figure;
 hold on; plot(t,ifhat,'g')
 title( 'Instantaneous frequency estiwmation')
 xlabel('Time');
 ylabel('Nomized Frequency (Hz)');