声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 5063|回复: 12

[综合讨论] 如何把地震波时程转化为反应谱曲线啊?

[复制链接]
发表于 2007-4-24 16:45 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 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=[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)

评分

1

查看全部评分

发表于 2007-4-24 21:17 | 显示全部楼层
请问el是什么?没看明白……
发表于 2007-4-24 21:30 | 显示全部楼层
发表于 2007-4-24 21:40 | 显示全部楼层
load d:\EL.txt
EL;
这个EL我没看出有什么用来。
发表于 2007-4-24 21:44 | 显示全部楼层
原帖由 后知后觉 于 2007-4-24 21:40 发表
load d:EL.txt
EL;
这个EL我没看出有什么用来。



的确无特别作用,可能是去除分号后调试用的。
 楼主| 发表于 2007-12-16 19:05 | 显示全部楼层
EL指的是ELCENTRO地震波
发表于 2008-3-14 17:06 | 显示全部楼层
请教,我运行了一下,怎么在长周期段(3-6秒),纵坐标值反而很大呢,而且不是下降趋势的?
程序是不是有问题啊
发表于 2008-4-8 10:02 | 显示全部楼层
貌似不是很正确哦,模拟的曲线的走向都不正确
发表于 2010-4-3 21:21 | 显示全部楼层
为什么我的提示程序有错误呢?我对MATLAB掌握的不怎么好 可能是没运行好!
发表于 2011-4-24 13:09 | 显示全部楼层
这个程序要求输入的 以行的形式输入,大家可以再试试看
发表于 2011-11-17 16:56 | 显示全部楼层
老说没定义是怎么回事呀
发表于 2012-6-19 17:15 | 显示全部楼层
运行中,出现下面两行,是怎么回事?

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

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-25 01:51 , Processed in 0.157899 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表