用非平稳Kanai-Tajimi模型模拟了地面加速度记录。
非平稳地面加速度记录的模拟是基于Guo et al.[3]的论文中提出的一个例子,我发现这个例子解释得很好。然而,该方法本身较旧,参见例[1,2]。此外,本文还包含了一个拟合非平稳Kanai-Tajimi模型到地面加速度记录的Matlab函数。优化工具箱用于拟合程序(lsqcurves fit)。但是,也可以使用其他函数。信号处理工具箱在这里通过函数filtfilt使用。欢迎任何对改进代码的评论、问题或建议。
参考文献
林永康,杨永勇(1987)。演化的Kanai-Tajimi地震模型。工程力学学报,113(8),1119-1137。
[2] Rofooei, f.r., Mobarake, A., & Ahmadi, G.(2001)。用非平稳Kanai-Tajimi模型生成人工地震记录。工程结构学报,23(7),827-837。
[3] Guo, Y., & Kareem, A.(2016)。基于时频盲源分离的非平稳数据系统辨识。声与振动学报,371,110 -131。
嗨,克里斯多夫,
您使用的Matlab版本可能与我使用的版本不兼容。目前的提交是用Matlab R2013b编写的,应该兼容最新版本和版本R2012
嗨Cheynet:
首先,谢谢你回答我的问题,我真的很感激jejeje cd,但再次我试图使用的例子。m文件你加入提交,它仍然给我过滤后的NaN值。我没有对你的程序做任何改动,我只是下载了你的文件并运行了它。
我甚至把自由的检查您所使用的过滤器和尝试与他人低通滤波器(椭圆和chebysheb i型)但问题仍然是相同的,我甚至narginchk遇到了麻烦,因为它说,函数是一个未定义的函数或方法参数类型的两倍,和“y”和“t”值矩阵。
我在很类似,您正在使用您的程序,但我使用它来计算爆破振动产生的但我一直有一些问题来计算加速度记录的参数,因为我无法适应信封和计算精确的参数。我真的认为你的程序解决了这个问题,但我不能使它正常工作,我不知道你是否会这么甜来帮助我与这个…
嗨Cristopher,
感谢您的反馈!
如果我使用已加入提交的示例,在过滤加速记录的包膜后,我看不到任何NaN。你从真实的地面加速度记录中得到NaN值了吗?如果是:
检查你的数据是否包含任何NaN值
-这可能是因为你的记录的特征与我所附的例子明显不同。拟合程序可能对拟合参数的初始值敏感。这意味着对于给定的加速度记录,变量“猜信封”没有理由与我给出的示例完全相同。
最好的问候,
艾蒂安
你有没有注意到,当你试图计算“Y = filtfilt(d1. sosmatrix,d1。ScaleValues, abs(hilbert(y)));
%%获得包络参数”你有一些NaN值??我认为这可能是你在程序开始时使用的参数,或者你认为这可能是问题,我试图使用它与一些地面加速度记录,我有,它不工作....
可能好。让我们检查。
好的代码!