How to calculate mean wind direction
29 次查看(过去 30 天)
显示 更早的评论
Hello!
I need help figuring out how to calculate mean wind direction when my data is in degrees (0-360). I just realized my current program does not take into account that the data is circular, and the mean of 355 and 5 will be 180, instead of 0. Any help is greatly appreciated! I am a beginner when it comes to MATLAB programming
采纳的回答
José-Luis
2014-5-5
average = mod(sum(data),360)./numel(data)
Please accept an answer if it helped you.
12 个评论
Naseef Muhammed
2022-9-12
@Walter Roberson I'm refering to the function 'windir_avg' you provided above. In that function, if we given 2 anlges from third or fourth, it will give a wrong answer. for example, if we put 345 and 355 the answer will be 190. similarly, 240 and 250 will result in 295. I hope it can be resolved if we edit your script as below:
function [windir_avged] = windir_avg(windir)
windir_avged=180/pi*angle(sum(exp(1i.*windir.*pi/180))/numel(windir));
if windir_avged < 0
windir_avged=windir_avged+360;
end
Walter Roberson
2022-9-12
I assumed that @giacomo labbri provided the correct formula in https://www.mathworks.com/matlabcentral/answers/128381-how-to-calculate-mean-wind-direction#comment_896547
更多回答(3 个)
Jenna Marie
2014-5-5
1 个评论
José-Luis
2014-5-5
No worries.
doc numel
Counts the number of elements in the matrix.
Robert Daly
2021-6-16
I needed a solution that would ignore NAN values in the data.
Converts the direction data into X & Y vector components, averages those, then converts back to direction.
function [windir_avged] = windir_avg(windir)
[x,y] = pol2cart(deg2rad(windir),ones(size(windir)));
x=mean(x,'omitnan');
y=mean(y,'omitnan');
[windir_avged,~]=cart2pol(x,y);
windir_avged = rad2deg(windir_avged);
end
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!