%樣本點數
N=256;
%線性調頻信號
x1=fmlin(N,0,0.2);
x2=fmlin(N,0.3,0.5);
%疊加
x=x1+x2;
%瞬時頻率
ifr=instfreq(x);
%群延遲
f=0:0.01:0.5;
gd=sgrpdlay(x,f);
%顯示
figure(1)
plot(ifr,'LineWidth',2);
xlabel('時間 t');
ylabel('歸一化的頻率 f');
figure(2)
plot(gd,f,'LineWidth',2);
xlabel('時間 t');
ylabel('歸一化的頻率 f');