f3_26.m

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

%样本点数
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');