时频分析中怎样使坐标反映真实时间和频率

s18990

我想用wigner-ville变换处理信号，信号是包括时间和幅值的10000*2数据，编了一段程序，坐标轴不能反映实际中的时间和频率，哪位高人给解决一下，万分感谢。程序如下：
clear;
[fname,pname]=uigetfile('*.csv','选择数据文件');
fp=fopen(strcat(pname,fname),'r');
M=dlmread(strcat(pname,fname),',',1,0); fclose(fp);
t_original=M(:,3); y_original=M(:,4);
% figure(1); plot(t_original,y_original); grid on; axis([-inf,inf,-inf,inf]);

y_original(1:3000)=0; y_original(7001:10000)=0;
for n=1:2000, y1(n)=y_original(n.*5); t1(n)=t_original(n.*5); end
for n=2001:2048, y1(n)=0; t1(n)=0; end

figure; y2=hilbert(y1); x=y2'; [tfr,t,f]=tfrwv(x); contour(tfr,10)
axis([600,1200,1960,2000]); box on; grid on; set(gcf,'color','w')

[ 本帖最后由 ChaChing 于 2010-3-10 09:43 编辑 ]