f7_1.m

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

%产生暂态信号
transig=amexpo1s(64).*fmconst(64);
sig=[zeros(100,1);transig;zeros(92,1)];
%加噪声
sign=sigmerge(sig,noisecg(256),-5);

%绘制时域波形
figure(1);
plot(real(sign),'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');

%计算功率谱
nfft=1024;
dsp=abs(fft(sign,nfft)).^2;
%绘制频谱图
figure(2);
plot((0:nfft/2-1)/nfft,dsp(1:nfft/2));
xlabel('频率 f');
ylabel('功率谱 PSD');

%计算信号的谱图
[tfr,t,f]=tfrsp(sign);
%绘图
%%绘制等高线图
figure(3);
contour(t,f(1:128),abs(tfr(1:128,:)).^2);
xlabel('时间 t');
ylabel('频率 f');
%%三维图

figure(4);
mesh(t,f(1:128),abs(tfr(1:128,:)).^2);
xlabel('时间 t');
ylabel('频率 f');
zlabel('幅值 A');