Сгенерируйте сигнал, который напоминает историю нагрузок, состоящую из синусоидальных полупериодов, соединяющих известные, равноудаленные реверсы. Дискретизация сигнала производится на частоте 512 Гц в течение 8 секунд. Постройте график экстремумов и сигнала.

fs = 512;

X = [-2 1 -3 5 -1 3 -4 4 -2];
lX = length(X)-1;

Y = -diff(X)/2.*cos(pi*(0:1/fs:1-1/fs)') + (X(1:lX)+X(2:lX+1))/2;
Y = [Y(:);X(end)];

plot(0:lX,X,'o',0:1/fs:lX,Y)

Вычислите счетчики цикла для данных. Отобразите матрицу отсчётов циклов.

[c,hist,edges,rmm,idx] = rainflow(Y,fs);

T = array2table(c,'VariableNames',{'Count','Range','Mean','Start','End'})

T=7×5 table
    Count    Range    Mean    Start    End
    _____    _____    ____    _____    ___

     0.5       3      -0.5      0       1 
     0.5       4        -1      1       2 
       1       4         1      4       5 
     0.5       8         1      2       3 
     0.5       9       0.5      3       6 
     0.5       8         0      6       7 
     0.5       6         1      7       8 

Отобразите гистограмму отсчета циклов как функцию от области значений напряжения.

histogram('BinEdges',edges','BinCounts',sum(hist,2))
xlabel('Stress Range')
ylabel('Cycle Counts')

Использование rainflow без выходных аргументов для отображения гистограммы циклов как функции от среднего цикла и области значений цикла.

rainflow(Y,fs)

Сгенерируйте сигнал, который напоминает историю нагрузок, состоящую из синусоидальных полупериодов, соединяющих известные, неравномерно расположенные реверсы. Дискретизация сигнала производится на частоте 10 Гц в течение 15 секунд. Постройте график экстремумов и сигнала.

fs = 10;

X = [0 1 3 4 5 6 8 10 13 15];
Y = [-2 1 -3 5 -1 3 -4 4 -2 6];

Z = [];
for k = 1:length(Y)-1
    x = X(k+1)-X(k);
    z = -(Y(k+1)-Y(k))*cos(pi*(0:1/fs:x-1/fs)/x)+Y(k+1)+Y(k);
    Z = [Z z/2];
end
Z = [Z Y(end)];

t = linspace(X(1),X(end),length(Z));
plot(X,Y,'o',t,Z)

Вычислите счетчики цикла для данных. Отобразите матрицу отсчётов циклов.

[c,hist,edges,rmm,idx] = rainflow(Z,t);

TT = array2table(c,'VariableNames',{'Count','Range','Mean','Start','End'})

TT=7×5 table
    Count    Range    Mean    Start    End
    _____    _____    ____    _____    ___

     0.5       3      -0.5      0       1 
     0.5       4        -1      1       3 
       1       4         1      5       6 
     0.5       8         1      3       4 
       1       6         1     10      13 
     0.5       9       0.5      4       8 
     0.5      10         1      8      15 

Использование rainflow без выходных аргументов для отображения гистограммы циклов как функции от среднего цикла и области значений цикла.

rainflow(Z,t)

Сгенерируйте случайный сигнал, дискретизированный на частоте 100 Гц в течение 100 секунд. Сохраните сигнал и его информацию о времени в расписании.

fs = 100;
t = seconds(0:1/fs:100-1/fs)';

x = randn(size(t));
TT = timetable(t,x);

Отобразите развороты и радужную матрицу сигнала.

rainflow(TT)

Сгенерируйте набор экстремумов, напоминающих развороты нагрузки. Постройте график данных.

X = [-2 1 -3 5 -1 3 -4 4 -2]';

plot(X)
xlabel('Sample Index')  
ylabel('Stress')

Вычислите счетчики цикла для данных. Укажите, что вход состоит из уже идентифицированных экстремумов.

[C,hist,edges] = rainflow(X,'ext');

Отобразите гистограмму отсчета циклов как функцию от области значений напряжения.

histogram('BinEdges',edges','BinCounts',sum(hist,2))
xlabel('Stress Range')  
ylabel('Cycle Counts')

Использование rainflow без выходных аргументов для отображения гистограммы циклов как функции от среднего цикла и области значений цикла.

rainflow(X,'ext')