脉动风时程matlab程序4页.doc

根据风的记录，脉动风可作为高斯平稳过程来考虑。观察个具有零均值的平稳高斯过程，其谱密度函数矩阵为： （9）将进行Cholesky分解，得有效方法。 （10）其中， （11）为的共轭转置。根据文献8，对于功率谱密度函数矩阵为的多维随机过程向量，模拟风速具有如下形式： （12）其中，风谱在频率范围内划分成个相同部分，为频率增量，为上述下三角矩阵的模，为两个不同作用点之间的相位角，为介于和之间均匀分布的随机数，是频域的递增变量。文中模拟开孔处的来流风，因而只作单点模拟。即式（4）可简化为： （13）本文采用Davenport水平脉动风速谱： （14） 式中，脉动风速功率谱；脉动风频率(Hz)；地面粗糙度系数；标准高度为10m处的风速(m/s)。Matlab程序:N=10;d=0.001;n=d:d:N;%频率区间（0.0110）v10=16;k=0.005;x=1200*n/v10;s1=