running mean of circular data (angles)

15 次查看（过去 30 天）
Jakob Sievers
Jakob Sievers 2011-10-5
评论： Walter Roberson 2020-5-7
Hi there
Im trying to do a running mean of a quite large wind-direction data-set but I'm having trouble finding solutions that work. Does anyone know of a simple solution?
Thanks in advance

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

采纳的回答

Mark Brandon
Mark Brandon 2011-10-5
You could look at the excellent "Circular Statistics Toolbox"
http://www.mathworks.com/matlabcentral/fileexchange/10676
  5 个评论
显示 3更早的评论
Ana Soph
Ana Soph 2020-5-6
  1. how can i use mean angle?
  2. M = xlsread('Dic2019','Dic2019mat')
  3. ten_minute_avarege = meanangle(reshape(M, 10, []), 'omitnan')
  4. xlswrite('datosexcel1.xlsx', ten_minute_avarege, 'Hoja1', 'A1');
Walter Roberson
Walter Roberson 2020-5-7
I am not aware of any meanangle() function that has 'omitnan' as an option.

请先登录，再进行评论。

更多回答（2 个）

David Young
David Young 2011-10-5
Have a look at
doc unwrap
Jakob Sievers
Jakob Sievers 2011-10-5
The Circular Statistics Toolbox doesnt appear to have a running mean function though. I have used its circ_mean function to make my own running mean but it takes hours to complete its calculation due to the size of the dataset (basically a long vector containing 10 million points).
Does anyone have any other suggestions?
  4 个评论
显示 2更早的评论
David Young
David Young 2011-10-5
Note that unwrap assumes your angles are in radians. If you're using degrees (and from the 360 in your code I think you are), multiply your data by pi/180 before using unwrap. You can convert back afterwards if necessary by multiplying by 180/pi.
Also, perhaps there's a problem with runmean. Maybe you could explain a little more about what this function does, or if it's very short even copy it into your question. (Not into a comment or answer though.)
Jakob Sievers
Jakob Sievers 2011-10-5
Hi David
I just talked with a colleague who said that a 30min mean would be sufficient, meaning that I can make a simple running mean that runs in a 20 seconds or so myself using the circ_mean as described above. Thanks for the help though!

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by