时频分析程序哪里错了？求助：

鸭鸭112

200110120145MHS-UD.txt (1.41 KB, 下载次数: 7)

2008-7-30 11:28 上传

点击文件名下载附件

x=load('200110120145MHS-UD.txt');%读入信号
fs=100; %采样率
n=length(x);
t=0:1/fs:(n-1)/fs;
nfft=2^nextpow2(n);%傅氏谱点数
x1=detrend(x,'linear');%去倾斜
x2=fft(x1,nfft);
a=load('Gs(MHSUD).txt');%读入仪器系统响应
n1=length(a);
nfft1=2^nextpow2(n1);
a=fft(a,nfft);
x2=x2./a';%扣除仪器响应
x2=ifft(x2,nfft);
x2=real(x2(1:n));
tfrsp(x2);
能否用X2的实部来进行时频分析？望指教，谢谢 ！