谐波合成法matlab,基于Kaimal谱采用谐波合成法生成脉动风场

%********************* 谐波叠加法模拟风速时程 修改1.0(采用Kaimal谱)**************

clc

clear

%************************* 风速时程参数设定 *************************

m=32;         %模拟风速点数

N=2^11;       %频率采样点数

dt=0.25;      %时间间隔

omegaup=3*pi; %上限频率

%************************  设定风速谱参数 *************************

L=514.5;      %斜拉桥跨度

z=45;         %风速测点离地面高度

z0=0.03;      %地面粗糙度

Uz=20;        %平均风速

lambda=10;    %空间相关函数中的衰减系数

K=0.4;        %Kaimal常数=0.4

M=2*N;        %nfft傅里叶变换长度，取采样个数的2倍

%形成风速时程矩阵

v=zeros(m,M*m);    %创建m行，M*m列的时程矩阵

u=zeros(m,M*m);

v1=zeros(M*m,m);   %创建m行，M*m列的时程矩阵

u2=zeros(M*m,m);

t=dt*(0:1:(M*m-1));%创建时程横坐标时间点

domega=(omegaup-0.001)/N;%频率间距

D=zeros(m,m,N);