关于人工合成地震波 [复制链接]

人工合成地震波方面，目前我正在研究，但是精度不是很理想。迭代后0.03s-4s之间的平均误差基本在5~10%。看咱论坛上有不少讨论，可能关于人造波的拟合问题是比较复杂，并不是一下子能达到一个比较好结果吧。我把我的几个小程序贴出来看看，希望高手给予一定的指导。 Mo`7YS-Y  
我这里一共有三个程序 X]d["  
第一个是关于计算反应谱的，使用方法是纽马克线性加速度法，程序是从别人的matlab程序稍微改编后进行移植成Fortran程序。 >(+g:p  
第二个是关于线性插值的，生成造波的时候所需要的周期所对应的谱值。 Q.:SIBP  
第三个是造波用的，我是把别的论坛里“老先生”所提供的程序进行了一定的改编，输入文件为ASEW.DAT。主要是频率输入更加精细（以前提供的是分2个圆频率变化值，1-163的增量为0.2，164-257增量为1，而我的程序里所有输入圆频率的增量都不相同），反应谱的计算方法用的也是纽马克线性加速度法。结果文件一共有WAVEACompare.resp（反应谱值对比），WAVEAFINAL.OUT（迭代完成后生成的地震波程序），WAVEAFINAL.resp（迭代完成后生成的地震波谱值），WAVEA.OUT（所有中间生成的地震波），WAVEA.resp（所有中间生成的地震波反应谱）。 `ZL~k  
PS。不知道站内高手什么时候上线，我把这几个fortran的源代码贴出来，向高手请教。对于精度问题，我想肯定有比较好的解决办法，只是似乎文献上写的都不是很详细。不过地震局做的安评报告拟合的就相当不错，而且北京波谱的Vib'SYS程序做到精度也比较不错，不过软件收费，而且输入比较麻烦。

不错，谢谢！

是地震分析的好资料的，我刚参与地震分析，有很多还需要学习

先看看再说了

ricozhou挺能钻研，支持一下，呵呵！ IB~`Ht8 b  
我08年的时候也用MATLAB写过一个人工合成地震动的程序，方法用的最经典的三角级数叠加法，不过好像当时对初始相位的选取做过一点尝试，不是用的0-2pi均匀分布的假设。计算的精度也有些问题。不过我当时好像用的目标谱是实际计算出来的频率谱，且精度控制点非常多，不是安评时候用的50-60个控制点的对数周期谱，所以也可能是这个原因导致不是所有控制点的精度都能满足5%的要求。 ci]IH]x  
有时间我把那个程序找出来改改，然后贴出来供大家参考！

shantuyue用随机振动理论合成地震动的程序上传后，我也尝试着用了一下，但是总出错误，不知道是什么原因.不知道是不是少文件。

谢谢楼主分享

引用第5楼ricozhou1231于2010-05-11 11:39发表的  : /< OoZf+[  
shantuyue用随机振动理论合成地震动的程序上传后，我也尝试着用了一下，但是总出错误，不知道是什么原因.不知道是不是少文件。 **oaR  

P^wDt14>  
什么错误？我最近太忙，等有时间给你看一下。 V30Om3C  
.u+ZrA#  
那个程序是因为看见有人上传那本书但没上程序，就顺手补了一下。我们现在都不用那方法的，过时了。

不知道shantuyue现在都用什么方法合成地震波呢？请教中……

不错，谢谢了 u0s25JY.%  

谢谢楼主分享。。。。

我当时编出来的程序误差也比较大，后来在迭代中改进了下，是有所好转，但是没有达到很多论文中那样