1、内容简介
matlab137-模拟风场的随机脉动风
可以交流、咨询、答疑
2、内容说明
略
模拟风场的随机脉动风，并进行相关的统计分析和计算，包括风速谱、空间相关性、自谱、互谱、以及POD（Proper Orthogonal Decomposition）分解等。下面是对各部分的中文解释：
1. 数据图形变量清除
clear all;
close all;
clc;
tic;
1.clear all：清除工作空间中的所有变量。
2.close all：关闭所有已打开的图形窗口。
3.clc：清除命令窗口中的所有内容。
4.tic：启动计时器，用于测量代码执行时间。
2. 生成样本的参数
这里设置了空间模拟和频率相关的参数。
Npoint = 17; % 空间模拟点数
fu = 2.0; % 频率范围上限
fl = 0.0001; % 频率范围下限
Wu = fu * 2 * pi; % 圆频率上限
Wl = fl * 2 * pi; % 圆频率下限
Nfrq = 1024; % 频段数
Dw = (Wu - Wl) / Nfrq; % 频段宽度
Nxx = 2; % 频段倍数
Mfrq = Nfrq * Nxx; % 总频段数
5.Npoint：设置空间点数，这里为17个。
6.fu 和 fl：设定频率范围，上限为2 Hz，下限为0.0001 Hz。
7.Mfrq：计算后的总频段数。
通过这些参数，计算出每个频段的中心频率 (Wll) 和其他频率相关的参数。
3. 随机脉动风的参数
这部分设置了风速、湍流和风谱的参数。
8.用于计算平均风速：
matlab
U10 = 24; % 10米处的平均风速
AA = 0.16; % 风速指数
9.用于计算湍流：
matlab
Zk = 0.4; % 摩擦速度
Z0 = 0.02; % 地表粗糙度
AAIu = 0.15; % 湍流动量通量
Zg = 350.0; % 测量高度
Uaa = 1.544747; % 测量高度的风速
10.用于计算湍流的功率谱：设置了湍流各个方向的功率谱参数，如Sdu，Sdv，Sdw，Yu，Yv等。
这些参数用于模拟和计算湍流和风速变化，以及风场的统计特性。
4. 生成空间点的坐标
生成一个空间点阵，表示风场中的不同点：
Point0 = [0, 0, 100; 5, 0, 100; ...]; % 每个点的坐标 (x, y, z)
这些点用于模拟风场在不同位置的风速特性。
5. 计算沿着高度方向的风参数
使用TurbPara函数计算风场的湍流参数：
Turb0 = TurbPara(Npoint, Point1, Turbpara1, Turbpara2, NumPSDchoice);
这一步根据先前设置的参数计算每个空间点的风速和湍流特性。
6. 计算自谱
自谱表示在每个空间点的频率谱分布：
PSDNum = 3;
PSD00 = PSDxx(PSDNum, Point1, Turb0, Turbpara1, Turbpara2, Sampara);
通过指定的参数计算出每个空间点的自谱。
7. 计算空间相关函数
空间相关函数用于描述不同位置点之间的风速相互影响：
COH00 = cohxy(PSDNum, Point1, Sampara, Corrpara, Turb0);
该函数计算了空间点之间的风速相关性。
8. 计算互谱函数
互谱描述了不同空间点之间的频率关系：
PSD12 = PSDxy(Sampara, PSD00, COH00);
该函数基于自谱和空间相关函数计算了互谱。
9. POD分解
这部分通过POD方法对互谱进行分解，提取出主要的模态：
for LL = 1:Sampara(4)
AS1 = tril(PSD12(:,:,LL), -1) + triu(PSD12(:,:,LL)', 0);
[EigVec, EigVal] = eig(AS1, 'vector');
end
POD分解用于将风场中的风速场表示为有限个主要模态，用于分析风场的主要特征。
。
3、仿真分析
略
4、参考论文
略