请教，地震波处理的命令

dongzhijun

目的
把地震波（加速度时程）转换成速度时程有没有现成的命令可用。
求地震波的反应谱又没有现成的命令？

kkkttt

http://forum.vibunion.com/forum/viewthread.php?tid=9253

http://forum.vibunion.com/forum/viewthread.php?tid=13598

happy

地震波反应谱计算程序

1. function [x]=rspsp(damping,eqw,T)
   
2. %
   
3. % function [x]=rspsp(damping,eqw,T)
   
4. % Purpose: To calculate response spectrum of EQW
   
5. %    damping  阻尼比
   
6. %    eqw      地震波，地面加速度时程. 应具有两列[t,u]
   
7. %             当 eqw 是向量时，则认为eqw(1)是时间间隔
   
8. %    T        向量，周期
   
9. %    x        计算结果，共三列，分别为位移，速度，加速度反应谱
   
10. % Note:    What we calulate are called Pseudo-velocity response
    
11. %          spectrum and Pseudo-acceleration response spectrum.
    
12. % Author:  hjdwg@0451.com  98-11-25 15:57   
    
13. %
    
14. 
    
15. if nargin ~= 3 ,error('Nargin ~=3');end
    
16. [s1,s2]=size(eqw);
    
17. if s2>2,error('Illegal eqw format');end
    
18. if s2==2,
    
19.     u=eqw(:,2);
    
20.     t=eqw(:,1);
    
21. end
    
22. if s2==1,
    
23.     t=eqw(1)*(0:s1-2);
    
24.     u=eqw(2:s1,1);
    
25. end
    
26. x=zeros([length(T),3]);
    
27. w=2*pi./T;
    
28. for i=1:length(w)
    
29. SYS=tf(1,[1,2*damping*w(i),w(i)*w(i)]);
    
30. x(i,1)=max(abs(lsim(SYS,u,t)));
    
31. end
    
32. x(:,2)=x(:,1)*2*pi./T(:);
    
33. x(:,3)=x(:,2)*2*pi./T(:);
    
34. return

复制代码