f3_4.m

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

%时间
t=2:255;
%信号1产生
sig1=amgauss(256,128,90).*fmlin(256,0,0.5);
%计算时宽
[tm,T1]=loctime(sig1);
%计算频宽
[fm,B1]=locfreq(sig1);
TB1=T1*B1
%计算瞬时频率
ifr1=instfreq(sig1,t);

%频率
f1=linspace(0,0.5-1/256,256);
%计算群延迟
gd1=sgrpdlay(sig1,f1);
%显示
plot(t,ifr1,'b-','LineWidth',2);
hold on;
plot(gd1,f1,'o--','LineWidth',2);
xlabel('时间 t');
ylabel('归一化的频率 f');
legend('瞬时频率','群延迟');