%
clear all;
clc;
N=500; %采样点数
wu=4*2*pi; %截断频率
dm=wu/N; %频率步长
dt=3.2*pi/(2*wu); %时间步长0.2
k=0.00464; %地面粗糙度系数地面粗糙度等级A B C D:K= 0.00129 0.00215 0.00464
0.01291
d=0.001;
f=d:d:10; %时间从0.001到10s,步进值为0.001
v10=26.8; %设计风速 26.8m/s――50年一遇十分钟平均风速最大值
x=1200*f/v10;
s=4*k*v10*v10.*x.^2./f./(1+x.^2).^(4/3);
z1=10; %取第一点为10米高度
z2=52.08; %取第二点为52.08米高度
r=0.2; %考虑地表粗糙度影响的无量纲幂指数,按中国规范取0.22-c类
v5=8; %计算n米高处的平均风速――52.8m处平均风速期望值
C=10; %指数衰减系数(取经验值)
v1=zeros(2*N,1); %产生一个全零矩阵
v2=zeros(2*N,1);
thta1=rand(N,1);
thta2=rand(N,1); %随机矩阵
node=1;
for K=1:node
for j=1:2*N
sum1=0;
sum2=0;
for l=1:N
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余2页未读,立即下载