clear;close;
load('ysc_197510.mat');
windspeed=d(:,7);winddirection=d(:,8);
q=1;b(1:17)=0;
for i=0:22.5:360
[e,j]=find(winddirection==i);
if ~isempty(e)%分析其是否为空
i=d(e,7);
b(q)=max(i);
else
b(q)=0;
if i==360
b(q)=b(1)
end
end
q=q+1;
end
a=0:22.5*pi/180:360*pi/180;
h=polar(a,b);
view(90,-90);
hold on
patch(get(h,'xdata'),get(h,'ydata'),'r','EraseMode','xor');%利用补丁
填充颜色
title('max windspeed at yinshuichuan station
197510','fontsize',14,'fontname','times new roman');
评论0