非平稳的地震加速度波视为平稳随机过程和考虑非平稳特性的外包线函数的乘积[1]:
a(t)=f(t)a(t)
其中,t1、t2 分别为控制强震平稳段的首末时刻; c 是衰减因子,控制下降段衰减的快慢; td 表示地震动总持时。一般地,用强震动平稳段的持时 ts 来反应参数 t2 的特性,需满足关系式: ts = t2 - t1。[2] 按该表取值:
本文使用随机数生成平稳相位波与外包线函数相乘得到非平稳相位波,只需将a(t)函数进行替换即可:
Python代码为:
import numpy as np
import math
import matplotlib.pyplot as plt
plt.style.use("classic")
plt.rcParams['font.sans-serif'] = ['SimSu