如何把地震波时程转化为反应谱曲线啊？

zxwflyer

各位高手，如何用MATLAB把地震波时程转化为反应谱曲线啊？:handshake

laohuren

这里有源程序，调试一下
load d:\EL.txt
EL;
amax=max(abs(EL));
dt=0.02
ff=EL;
n=0;
zeta=0.02;
for w=100:-0.2:1;
x=0;
v=0;
alhp=35;
acc=0;
a=[1,2*zeta*w,w^2];
b=1;
[r,p]=residue(b,a);
t=0:dt:length(EL)*dt;
h=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)
x=conv(h,ff)*dt;
v=diff(x)/dt;
acc=diff(v)/dt;
n=n+1;
beta(n)=max(abs(acc))*alhp/amax;
T(n)=2*pi/w;
end;
plot(T,beta)

后知后觉

请问el是什么？没看明白……

eight

原帖由 后知后觉 于 2007-4-24 21:17 发表

                               
登录/注册后可看大图

请问el是什么？没看明白……

应该是原始数据的结构，即 txt 文件存储的是结构体变量EL

后知后觉

load d:\EL.txt
EL;
这个EL我没看出有什么用来。

eight

原帖由 后知后觉 于 2007-4-24 21:40 发表

                               
登录/注册后可看大图

load d:EL.txt
EL;
这个EL我没看出有什么用来。

的确无特别作用，可能是去除分号后调试用的。

zxwflyer

EL指的是ELCENTRO地震波

banywong

请教，我运行了一下，怎么在长周期段（3－6秒），纵坐标值反而很大呢，而且不是下降趋势的？
程序是不是有问题啊

happyren526

貌似不是很正确哦,模拟的曲线的走向都不正确

liujiashun

为什么我的提示程序有错误呢？我对MATLAB掌握的不怎么好 可能是没运行好！

pansidong311

这个程序要求输入的 以行的形式输入，大家可以再试试看

beauty33333

老说没定义是怎么回事呀

yzsldj

运行中，出现下面两行，是怎么回事？

??? Error using ==> conv
A and B must be vectors.