clc;clear;close all
b=[21,23,11,14,26,126,79,25]/328;%按照北、东北、东、东南、南、西南、西、西北顺时针顺序；
a=2/328;%静风频率
n=length(b);
po=b;
m=5;%圈数
w=ceil(max(po)/m)/100
dir_ang=(0:360/n:360)*pi/180;
wn=1/m;
po=po/w/100*wn;
po(n+1)=po(1);
po=po*wn/max(po)*m;
px=(po).*sin(dir_ang);
py=(po).*cos(dir_ang);
ang_cyc=0:0.01:2*pi;
for i=1:n+1
plot([-1-0.08]*cos(dir_ang(i)),[-1 -0.08]*sin(dir_ang(i)),'k');
hold on
plot([0.08 1]*cos(dir_ang(i)),[0.08 1]*sin(dir_ang(i)),'k');
hold on
end
for i=1:m
plot((i*wn)*cos(ang_cyc),(i*wn)*sin(ang_cyc),'k:')
% text((i-0.4)*wn,0.04,[num2str(i*w*100),'%']);
text((i)*wn,0.04,num2str(i*w*100));
hold on
end
plot(0.08* cos(ang_cyc),0.08*sin(ang_cyc) ,'k')
hold on
plot(px,py,'-','LineWidth',2);
plot(px,py,'.','LineWidth',4);
%text(-0.04,0,[num2str(a),'%']);
text((m*wn+0.1)*sin(2/180*pi),(m*wn+0.1)*cos(0/180*pi),'N','rotation',180)
%text((m*wn+0.1)*sin(17.5/180*pi),(m*wn+0.1)*cos(17.5/180*pi),'NNE','rotation',337)
text((m*wn+0.1)*sin(42/180*pi),(m*wn+0.1)*cos(42/180*pi),'NE','rotation',315)
%text((m*wn+0.1)*sin(62.5/180*pi),(m*wn+0.1)*cos(62.5/180*pi),'ENE','rotation',292.5)
text((m*wn+0.1)*sin(88/180*pi),(m*wn-0.8)*cos(88/180*pi),'E','rotation',0)
%text((m*wn+0.05)*sin(108.5/180*pi),(m*wn+0.6)*cos(108.5/180*pi),'ESE','rotation',67.5)
text((m*wn)*sin(132/180*pi),(m*wn+0.25)*cos(132/180*pi),'SE','rotation',45)
%text((m*wn-0.25)*sin(152.5/180*pi),(m*wn+0.2)*cos(152.5/180*pi),'SSE','rotation',22.5)
text((m*wn+0.1)*sin(178/180*pi),(m*wn+0.1)*cos(178/180*pi),'S','rotation',180)
%text((m*wn+0.9)*sin(195.5/180*pi),(m*wn+0.03)*cos(195.5/180*pi),'SSW','rotation',-22.5)
text((m*wn+0.3)*sin(220/180*pi),(m*wn)*cos(220/180*pi),'SW','rotation',-45)
%text((m*wn+0.2)*sin(241.5/180*pi),(m*wn-0.3)*cos(241.5/180*pi),'WSW','rotation',-67.5)
text((m*wn+0.15)*sin(268/180*pi),(m*wn-0.8)*cos(268/180*pi),'W','rotation',0)
%text((m*wn+0.1)*sin(286.5/180*pi),(m*wn+0.1)*cos(286.5/180*pi),'WNW','rotation',67.5)
text((m*wn+0.1)*sin(312/180*pi),(m*wn+0.1)*cos(312/180*pi),'NW','rotation',45)
%text((m*wn+0.1)*sin(331.5/180*pi),(m*wn+0.1)*cos(331.5/180*pi),'NNW','rotation',22.5)
axis off
axis equal
