Matlab风玫瑰图绘制代码示例

以下是使用Matlab绘制风玫瑰图的代码示例：

% 生成风速数据
wind_speeds = [3 5 6 7 8 10 11 11 12 13 14 14 15 15 16 17 17 18 18 20 22 24 25 25 26 26 27 27 28 29 30 32 34 35 35 36 36 37 37 38 39 40];

% 计算风向角度
wind_directions = [0:10:350];
angles = deg2rad(wind_directions);

% 统计每个角度的风速频率
freq = histcounts(mod(wind_speeds, 360), [0:10:360]);
freq(end) = freq(end) + freq(1);
freq(1) = [];

% 绘制风玫瑰图
figure
ax = polarplot(angles, freq);
ax.ThetaZeroLocation = 'top';
ax.ThetaDir = 'clockwise';
ax.ThetaTick = [0:30:330];
ax.RLim = [0 ceil(max(freq)/10)*10];
ax.RTick = [0:10:ax.RLim(2)];
ax.RTickLabel = ax.RTick;
thetaticklabels({'N','NE','E','SE','S','SW','W','NW'});
title('Wind Rose Diagram');

解释：

1. 首先生成了一个包含风速数据的向量 wind_speeds。

2. 然后定义了一个包含风向角度的向量 wind_directions，并将其转换为弧度制的角度 angles。

3. 使用 histcounts 函数统计每个角度的风速频率，并将其存储在向量 freq 中。在计算频率时，使用 mod(wind_speeds, 360) 将风速值转换为 0 到 360 度之间的角度值。

4. 最后，使用 polarplot 函数绘制风玫瑰图。在绘制图形之前，需要对极坐标轴进行一些设置，使其适合风玫瑰图的显示。例如， ax.ThetaZeroLocation = 'top' 将极坐标轴的起始点设置为顶部， ax.ThetaDir = 'clockwise' 将极坐标轴的方向设置为顺时针， ax.ThetaTick = [0:30:330] 设置角度刻度线的位置， ax.RLim = [0 ceil(max(freq)/10)*10] 设置极径的范围， ax.RTick = [0:10:ax.RLim(2)] 设置极径刻度线的位置， ax.RTickLabel = ax.RTick 将极径刻度线的值设置为刻度线位置， thetaticklabels({'N','NE','E','SE','S','SW','W','NW'}) 设置角度刻度线的标签， title('Wind Rose Diagram') 设置图形标题。