时频分析程序哪里错了?求助:
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的实部来进行时频分析?望指教,谢谢 !
页:
[1]