%時間
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('瞬時頻率','群延遲');