f4_5.m

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

%产生非平稳信号
%%暂态信号1
sig1=real(amgauss(128,45).*fmconst(128,0.25,45));
%%暂态信号2
sig2=real(amgauss(128,85).*fmconst(128,0.25,85));
sig=sig1+sig2;
%时域波形
figure(1)
plot(sig,'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');

%设置窗函数1
h1=window(65,'hamming');
%计算短时傅立叶变换
sig=hilbert(sig);
[tfr,t,f]=tfrstft(sig,1:128,128,h1);

%时频表示
figure(2)
contour(t,2*f(1:length(f)/2),abs(tfr));
xlabel('时间 t');
ylabel('频率 f');

%设置窗函数2
h2=window(17,'hamming');
%计算短时傅立叶变换
sig=hilbert(sig);
[tfr,t,f]=tfrstft(sig,1:128,128,h2);

%时频表示
figure(3)
contour(t,2*f(1:length(f)/2),abs(tfr));
xlabel('时间 t');
ylabel('频率 f');