亲爱的师傅们，帮忙看看我的程序时频三维图该怎么画

huangwanbo

第一次发帖，实在是弄不懂了。任务是做振动数据的短时傅里叶变换，可是我的三维图怎么也画不出来，还是黑色的。帮帮忙。程序如下：
clc;
clear;
close all;
t=0:0.0001:0.3499;
N=10000;
fs=10000;
n=1:3500;
sampledata=load('E:\资料\matlab学习\data(wangfeng)\7.20\2\4500.TXT');
y=sampledata(:,2);
y=y-mean(y);
p=y(1:3500);
% figure(1)
% plot(n,p);
% title('原始信号');
h1=window(625,'hamming');
[tfr,T,F]=tfrstft(y,1:3500,3500,h1);
figure(1)
subplot(2,1,1);
imagesc(t,F(1:100)*20000,abs(tfr(1:100,:)));
title('窗的长度取655');
h2=window(1655,'hamming');
[tfr1,T1,F1]=tfrstft(y,1:3500,3500,h2);
subplot(2,1,2);
imagesc(t,F1(1:100)*20000,abs(tfr1(1:100,:)));
title('窗的长度取1655');
figure(2)   %画时频三维图
b=F1(1:100)*20000;
c=abs(tfr1(1:100,:));
surfl(t,b,c);

kiefer0107

只这一行tfrstft(y,1:3500,3500,h1);用它自带的作图程序先试试。