那位大侠帮忙看看我编制的地震信号处理的程序和结果是否正确？

zzglyp

我选择el－centro地震波，分析其功率谱，编制的程序是这样的，请问正确吗？恳请各位大侠指点
function suru()
load ele1.dat
fs=50;T=1/fs;
ug=(ele1)';M=length(ug);
t=[0:M-1]*T;
figure(1);plot(t,ug);
ylabel('ug(t)');title('ug of input signal ');
xlabel('time(s)');
%compute Power spectrum and Rxx
N=M*2-1;
X=T.*fft(ug,N);
Xpower=X.*conj(X);
Rxx=1/T.*ifft(Xpower,N)
t1=[-M+1:M-1]*T;
figure(2);
plot(t1,real(fftshift(Rxx)));
ylabel('Rxx(t)');title('auto-correlation function');
xlabel('time(s)');
f=[-M+1:M-1]*fs/N;
figure(3);
plot(f,abs(fftshift(Xpower)))
ylabel('P(f)');title('power spectrum');
xlabel('Frequency(HZ)');

zzglyp

得到的谱图是这样的。

郁闷，文件大小有限制发不上来

Radiant

尽管功率谱的学术定义来自自相关函数的FFT，但是很少有人使用这个定义去计算功率谱，一般都是使用周期图法计算功率谱。