f3_3.m

matlab时频分析配套源码第3章的程序

%线性调频信号
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');