短时傅里叶变换中的trfstft的归一化处理完如何用程序将坐标在...

大天强

clear;close all;clc;
load cengliu.txt;
y=load('cengliu.txt');
fs=2500;
ts=1/fs;
nfft=10240;
N=150000;
t=(0:N-1)/fs;
figure(1);
plot(t,y)
xlabel('时间 t');
ylabel('幅值 y');
title('时域波形');
%计算功率谱
nfft=10240;
dsp=abs(fft(y,nfft));
%绘制功率谱图
figure(2);
plot((0:nfft/2-1)/nfft*fs,dsp(1:nfft/2),'LineWidth',2);
xlabel('频率 f');
ylabel('功率谱 PSD');%计算短时傅里叶变化
%%窗函数
N = 65;
h=window(@hamming,N);
sig=hilbert(y);
[tfr,t,f]=tfrstft(sig',1:60,1024,h);
%%绘制等高线图
figure(3);
%[tfr,t,f]=tfrstft(A);
[X,Y] = meshgrid(t,f);
contour(t,f(1:1024),abs(tfr(1:1024,:)));
%contour(t/2500,f(1:1024),abs(tfr(1:1024),:));
% contour(t(1:128)/256*fs,abs(tfr(1:128,:)).^2)
%contour(X,Y,abs(tfr(1:1024,:)).^2);
xlabel('时间 t');
ylabel('频率 f');
%axis([0 60 -0.1 0.1]);
%grid on