[热点讨论]求教抗震设计功率谱密度 [复制链接]

请问各位大侠如何能由地震波数据得到相应的功率谱密度(PSD)值？传说有相关的MATLAB程序可以实现，不知能否附上。。万分感谢！！

我这有个初程序不知对你有没有用 |K.I%B  
clear OsOfo({I_  
dt=;   %时间间隔 +wj}x?ZeV  
s=   %信号数组 OTYkJEC8\N  
len=128; UK6x]tE  
F=2*pi*(0:len/2)/len/dt; _E9[4%f  
fs=fft(s,len); @Ov}X]ELi  
p=fs.*conj(fs) 7b~uU@L`  
p=p(1:len/2+1); s58dHnj5+  
plot(F,p) >JNdtP8s/1  
ylabel('power spectrum density ') CL7_3^2qI  
xlabel('sequence ') \6AM?}v
 
title('03051171') !}}
)f/  
参数自已改吧

这个有点强悍