f3_5.m

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

%时间
t=2:255;
%信号1产生
sig2=amgauss(256,128,30).*fmlin(256,0.2,0.4);
%计算时宽
[tm,T2]=loctime(sig2);
%计算频宽
[fm,B2]=locfreq(sig2);
TB2=T2*B2
%计算瞬时频率
ifr2=instfreq(sig2,t);

%频率
f2=linspace(0.2,0.4,256);
%计算群延迟
gd2=sgrpdlay(sig2,f2);
%显示
plot(t,ifr2,'b-','LineWidth',2);
hold on;
plot(gd2,f2,'o--','LineWidth',2);
xlabel('时间 t');
ylabel('归一化的频率 f');
legend('瞬时频率','群延迟');