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