| 查看: 1054 | 回复: 1 | |||
1229669946新虫 (初入文坛)
|
[交流]
大型结构的随机减量法处理
|
|
现有一个重力坝的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  SE4.61加速度时程曲线.jpg  SE4.62加速度时程曲线.jpg  SE4.63加速度时程曲线.jpg  SE4.64加速度时程曲线.jpg  SE4.61随机减量后的自由响应曲线.jpg  SE4.62随机减量后的自由响应曲线.jpg  SE4.63随机减量后的自由响应曲线.jpg  SE4.64随机减量后的自由响应曲线.jpg |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有14人回复
-
这样的filecode谁见过
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】求孔隙率和比表面积的测试方法
已经有7人回复
1229669946
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 120.2
- 帖子: 11
- 在线: 1.3小时
- 虫号: 2855571
- 注册: 2013-12-07
- 性别: GG
- 专业: 水工结构和材料及施工
2楼2015-06-14 20:16:25
100