matlab 模拟风,matlab谐波合成法的脉动风模拟

clc

clear

m=10;                   %模拟风速个数

N=2^8;                   %频率采样个数

dt=0.5 ;                   %时间间隔

omegaup=2*pi;                %上限频率

%风速谱参数设定

L=1000;                         %跨度

z=50;                           %模拟点离地面高度

z0=0.03;                          %地面粗糙度

Uz=40;                         %50米处的平均风速

delta=100;                      %模拟点间距(m)

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

K=0.4;                          %Kaman常熟0.4

M=2*N;

V=zeros(m,M*m);

t=0.5*(0:1:(M*m-1));

domega=omegaup/N;

D=zeros(m,m,N);

U=K*Uz/log(z/z0);                    %U为摩阻速速

%形成目标谱

omegal=omegaup/N:omegaup/N: omegaup;

Sw1=200*U^2.*z/Uz./(1+50.*omegal.*z/(2*pi*Uz)).^(5/3);

for j=1:m