f3_3.m

来自「《Matlab时频分析技术及其应用》(人民邮电出版社)中第三章的随书程序。」· M 代码 · 共 18 行

%线性调频信号
sig=fmlin(256);
%频率
f=0:0.05:0.5;
%时间
t=3:256;
%瞬时频率估计
ifr=instfreq(sig);
%群延迟估计
gd=sgrpdlay(sig,f);
%显示结果
plot(t,ifr','LineWidth',2);
hold on
plot(gd,f,'o--','LineWidth',2);
legend('瞬时频率','群延迟');
xlabel('时间 t');
ylabel('归一化的频率 f');