用matlab2016a绘制风玫瑰图的代码

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

% 设置风速和风向数据
speed = [2 3 1 4 3 2 1 0 1 2 3 4];
direction = [135 135 135 135 135 135 135 135 135 135 135 135];

% 绘制风玫瑰图
polarhistogram(deg2rad(direction),12,'BinLimits',[0 pi],'BinWidth',pi/6,'FaceAlpha',0.5,'EdgeColor','none');
hold on;
pax = gca;
pax.ThetaZeroLocation = 'top';
pax.ThetaDir = 'clockwise';
pax.ThetaTick = [0:30:330];
pax.RTick = [0:5:max(speed)];
rticklabels(pax.RTick);
title('风玫瑰图');

解释：

• speed：风速数据，长度为12的一维向量，表示12个方向的风速。
• direction：风向数据，长度为12的一维向量，表示12个方向的风向（以正北为0度，正东为90度，正南为180度，正西为270度）。
• polarhistogram：绘制风玫瑰图的函数，第一个参数为弧度制的风向数据，第二个参数为方向的数量（12个方向），BinLimits参数设置极径的范围，BinWidth参数设置每个扇形的宽度，FaceAlpha参数设置扇形的透明度，EdgeColor参数设置扇形的边缘颜色。
• pax：获取当前的极坐标轴对象。
• pax.ThetaZeroLocation：设置0度方向的位置为顶部。
• pax.ThetaDir：设置极角的方向为顺时针。
• pax.ThetaTick：设置极角的刻度，每30度一个刻度。
• pax.RTick：设置极径的刻度，从0开始，每5个单位一个刻度。
• rticklabels：设置极径的标签，显示所有的刻度值。
• title：设置图像的标题。