Calculating the directional and magnitude frequency of wind at spec... (2024)

Hey Matt,

I appreciate the help. I am 'relatively new' to coding despite taking one class during my undergrad. Therefore, I spend time writing and figuring out what works and doesn't work while developing something.

The attempt with frequency is to find how often specific wind directions occur (work demands). Any suggestions?

I also need to extract time from the file (datetime, I think would be good), any suggestions?

I am trying to do one thing. How can we extract data from a specific set of rows and plot it using wind rose? For example, use data from row 2751:2798). I tried it just now and got an error ("index must not exceed 18 or index in position 2 exceeds array bounds")

(Don't ask about the sum code, it was an experimental attempt).

I also appreciate the help with the code you provided as I found the wind rose function and couldn't quite get it to work.

hello again

ok, I understand you wanted to create your own code and not rely on someone else code (including the wind rose function)

now this is a work that is maybe a bit tough for someone just starting coding so I wonder what kind of result you need to provide

IMO what you have to do is to create a histogram , but here it would be nice to have it in polar coordinates and that's exactly what polarhistogram can provide

so I wonder why someone needs to reinvent this function or if you are allowed to use it

also I still don't have the answer about where does your meanfreq function comes from ? seems not in the standard matlab package (the only one I am aware comes from the Signal Processing Toolbox and has only one input argument : Mean frequency - MATLAB meanfreq - MathWorks France

Yeah, it's a challenge for sure. I'm actually being forced to learn Python as I go (no excuse is allowed here at my particular position).

Before I found the wind rose code, is why I tried to develop my own. However, if I can find someone else's code that does what I need, I will gladly use it.

Ultimately, what I need is the number of times and count of (frequency) specific wind directions occurred (e.g., 330 degrees occurred x amount of times at time x on the 1st of January) at a given sector. Do that for the entire 24-hr period of January 1st. The follow-up will then do that for every single day and plot each one.

Can you verify my results from this code? I rewrote the main body to the following:

x = readmatrix('WindRose.txt')

WNDIR = x(2:2798,6);

WSPD = x(2:2798,7);

[C,ia] = unique(WNDIR)

W = WNDIR(ia,:)

[S,ia1] = unique(WSPD)

M = WSPD(ia1,:)

wind_rose(W,M)

I'm hoping that the unique function allows me to plot the windrose given the number of times a specific DIR and magnitude occurs beginning with the first row (row 2). On my end, the code works. I just would like a professional overview. Thanks!

hello again

a few suggestions to improve your code (don't forget also the ; at the end of each line , otherwise you fill your command window with the content of each variable that appears in the line)

I wonder if applying unique is a good thing , why only plot unique events ? if the wind blows often at the same speed and / or direction, why hide this is the plot ?

Thank you for the tip! I was hoping that the unique function would do the following:

• count the number of times a specific wind direction occurred (count) at a specific zone (e.g. 0 or 360 is zone 1)
• calculate the frequency those winds occurred
• plot the wind rose
• disp or print the chart of the results

well that's exactly the contrary, unique will remove all duplicates so in terms of statistics your are not taking into account the most frequent events and you will miss that in your plot

what you need to do is an histogram (see histcount) and that's what is doing the wind_rose function with polarhistogram

Hi Sakshi,

That's awesome, I appreciate the help. I gotta figure out how to frequencies now.