my.m

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

 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),10);
 figure;
 plot(real(sig));
 title( 'Signal in time')
 xlabel('Time');
 ylabel('real part');
 
 [ifhat,t]=instfreq(sig,11:502,10);
 figure;
 plot(t,ifhat)
 title( 'Instantaneous frequency estiwmation')
 xlabel('Time');
 ylabel('Nomized Frequency (Hz)');