my.asv

来自「经典数字信号处理滤波器的源代码 重点是利用巴特沃斯模拟滤波器转而设计其它数字滤」· ASV 代码 · 共 27 行

 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)');