- 积分
- 12511
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-7-4
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
%% 计算风向
for i = 1:length(U)
if U(i,1)>0&V(i,1)>0
dir(i,1) = 270-atan(V(i,1)./U(i,1))*180/pi;
elseif U(i,1)<0&V(i,1)>0
dir(i,1) = 90-atan(V(i,1)./U(i,1))*180/pi;
elseif U(i,1)<0&V(i,1)<0
dir(i,1) = 90-atan(V(i,1)./U(i,1))*180/pi;
elseif U(i,1)>0&V(i,1)<0
dir(i,1) = 270-atan(V(i,1)./U(i,1))*180/pi;
elseif U(i,1)==0&V(i,1)<0
dir(i,1) = 180;
elseif U(i,1)==0&V(i,1)<0
dir(i,1) = 0;
elseif U(i,1)>0&V(i,1)==0
dir(i,1) = 270;
elseif U(i,1)<0&V(i,1)==0
dir(i,1) = 90;
else U(i,1)==0&V(i,1)==0
dir(i,1) = nan;
end
end
%% 风向反向处理
for i = 1:length(dir)
if dir(i,1)>=180
dir(i,1) = dir(i,1)-180;
else
dir(i,1) = dir(i,1)+180;
end
end
|
|