24小时热门版块排行榜

返回列表

1229669946

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 120.2
帖子: 11
在线: 1.3小时
虫号: 2855571
注册: 2013-12-07
性别: GG
专业: 水工结构和材料及施工

[交流] 大型结构的随机减量法处理

现有一个重力坝的4个高程点的地震加速度记录，想用随机减量求加速度的自由衰减响应，用下列程序求，曲线总是不平滑，且只有一个测点有规律性，其他都是杂乱，具体情况见图片。有没有大神指导下，如何对实测记录进行处理呀？悬赏仅有的一个金币，谢谢啦！

for ii=1:mk; 测点个数
sf=100; %读入采样频率值
nt=length(x(ii,

);
t1=0:1/sf

nt-1)/sf; %原始数据的时间坐标
figure(1)
h1=plot(t1,x(ii,

);
xlabel('时间(s)');
ylabel('加速度幅值(cm/s2)');
grid on; %打印原始加速度时程曲线
filename1=['SE4.6' num2str(ii) '加速度时程曲线.jpg'];
saveas(h1,filename1)
close(1);

np=100; %输出数据长度
s=1.5*std(x(1,

);  %多次验算，取1.5倍的标准差自由曲线比较光滑
m=0;
y=zeros(1,np);
for k=2:nt-np
a=abs(x(1,k-1)-s);
b=abs(x(1,k)-s);
c=abs(x(1,k+1)-s); %触发条件为最接近S的点
if b<a&b<c
      y(1:np)=y(1:np)+x(ii,k:k+np-1);
      m=m+1;
end
end
y=y./m;
end