请教，地震波处理的命令
目的把地震波（加速度时程）转换成速度时程有没有现成的命令可用。
求地震波的反应谱又没有现成的命令？ http://forum.vibunion.com/forum/viewthread.php?tid=9253
http://forum.vibunion.com/forum/viewthread.php?tid=13598 地震波反应谱计算程序
function =rspsp(damping,eqw,T)
%
% function =rspsp(damping,eqw,T)
% Purpose: To calculate response spectrum of EQW
% damping阻尼比
% eqw 地震波，地面加速度时程. 应具有两列
% 当 eqw 是向量时，则认为eqw(1)是时间间隔
% T 向量，周期
% x 计算结果，共三列，分别为位移，速度，加速度反应谱
% Note: What we calulate are called Pseudo-velocity response
% spectrum and Pseudo-acceleration response spectrum.
% Author:hjdwg@0451.com98-11-25 15:57
%
if nargin ~= 3 ,error('Nargin ~=3');end
=size(eqw);
if s2>2,error('Illegal eqw format');end
if s2==2,
u=eqw(:,2);
t=eqw(:,1);
end
if s2==1,
t=eqw(1)*(0:s1-2);
u=eqw(2:s1,1);
end
x=zeros();
w=2*pi./T;
for i=1:length(w)
SYS=tf(1,);
x(i,1)=max(abs(lsim(SYS,u,t)));
end
x(:,2)=x(:,1)*2*pi./T(:);
x(:,3)=x(:,2)*2*pi./T(:);
return
页:
[1]