zxwflyer 发表于 2007-4-24 16:45

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

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

laohuren 发表于 2007-4-24 20:08

我用过的

这里有源程序，调试一下
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=;
b=1;
=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)

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

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

eight 发表于 2007-4-24 21:30

原帖由 后知后觉 于 2007-4-24 21:17 发表 http://forum.vibunion.com/forum/images/common/back.gif
请问el是什么？没看明白……


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

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

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

eight 发表于 2007-4-24 21:44

原帖由 后知后觉 于 2007-4-24 21:40 发表 http://forum.vibunion.com/forum/images/common/back.gif
load d:EL.txt
EL;
这个EL我没看出有什么用来。


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

zxwflyer 发表于 2007-12-16 19:05

EL指的是ELCENTRO地震波

banywong 发表于 2008-3-14 17:06

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

happyren526 发表于 2008-4-8 10:02

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

liujiashun 发表于 2010-4-3 21:21

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

pansidong311 发表于 2011-4-24 13:09

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

beauty33333 发表于 2011-11-17 16:56

老说没定义是怎么回事呀

yzsldj 发表于 2012-6-19 17:15

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

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

页: [1]
查看完整版本: 如何把地震波时程转化为反应谱曲线啊？