How can I specify the width of Bin Edges when calculating the rainflow cycle counts matrix 'c = rainflow(x)'

Dear Madam/Sir,

How can I specify the width of Bin Edges when calculating the rainflow cycle counts matrix 'c = rainflow(x)' ?

After computing the cycle counts for the data, using the command: [c,hist,edges,rmm,idx] = rainflow(Y) , I plot a histogram of cycle counts as a function of stress range, using the command: histogram('BinEdges',edges','BinCounts',sum(hist,2)) , as indicated in the example of the Matlab Documentation: https://ch.mathworks.com/help/signal/ref/rainflow.html#d120e139233

It seems that the 'Bin Edges' are assigned an arbitrary width value. Instead, I would like to assign a specific width value to the Bin Edges, while ploting the histogram of the cycle counts, resulting from the cycle counts matrix 'rainflow(x)'.

Could you please help me, how can I specify the width of Bin Edges when calculating the rainflow cycle counts matrix 'c = rainflow(x)' ?

Thank you very much in advance for your help

Best regards,

Geri

[Giuseppe Naselli]

The answer was already given to a similar question

https://uk.mathworks.com/matlabcentral/answers/446260-assign-width-of-bin-edges-within-rainflow-matrix-c-rainflow-x#answer_380775

[Cyril Besson]

No the question has not been answered!! (u can see in the link)

[Sandro Günter]

The rainflow function outputs the exact cycle range. You can count them by yourself using the histcounts function.

Here's an example:

[count,hist,edges,rmm,idx] = rainflow(T_J);  % Rainflow counting
my_edges = [1 5 10 20 30 40 50 60 80 100];  % define own bin edges
N = histcounts(count(:,2),my_edges);  % count according to own bin width
histogram('BinEdges',transpose(my_edges)','BinCounts',N)  % plot

[Sandro Günter]

However, it counts the half cycles as full ones.

Another way is to change the code of the rainflow function itself as suggested here: https://stackoverflow.com/questions/52513809/reduce-bin-size-for-matlab-rainflow-counting

[Sandro Günter]

In order to count the half cycles correctly, you can add the following code:

index = find(count(:,1)==0.5);  % find half cycles
N = N - 0.5 * histcounts(count(index,2),my_edges);  % subtract half cycles