How to calculate mean wind direction

20 次查看（过去 30 天）
Jenna Marie
Jenna Marie 2014-5-5
评论： Soeren Bilges 2023-2-10
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
  2 个评论
Jenna Marie
Jenna Marie 2014-5-5
Ah yes thank you for catching that! Accidentally mixed up 360 degrees in a circle with 365 days in a year. It's been a long day....

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

José-Luis
José-Luis 2014-5-5
average = mod(sum(data),360)./numel(data)
Please accept an answer if it helped you.
  12 个评论
显示 10更早的评论
Naseef Muhammed
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

请先登录，再进行评论。

更多回答（3 个）

Jenna Marie
Jenna Marie 2014-5-5
Thank you so much for your response!
What does numel mean?
Walter Roberson
Walter Roberson 2017-7-9
See unwrap() but you will need to convert to radians
Robert Daly
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 中查找有关 Shifting and Sorting Matrices 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by