Сгенерируйте сигнал, который напоминает историю загрузки, состоя из соединения полупериодов синусоиды известные, равномерно расположенные реверсирования. Сигнал производится на уровне 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 матрицу сигнала.

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')