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

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

我这有个初程序不知对你有没有用 ^1\[hyZ!  
clear s$3WJ'yr  
dt=;   %时间间隔 yhsbso,5 a  
s=   %信号数组 <)]j;Tl  
len=128; (?vKe5  
F=2*pi*(0:len/2)/len/dt; 1nZ7xCDK98  
fs=fft(s,len); O<\h_   
p=fs.*conj(fs)  _(_U=  
p=p(1:len/2+1); By;{Y[@rS  
plot(F,p) cx
_$`H  
ylabel('power spectrum density ') $5r1Si)  
xlabel('sequence ') V9\y*6#Y,  
title('03051171') @xkI?vK6  
参数自已改吧

这个有点强悍