Генерировать сигнал, который напоминает историю нагрузки, состоящую из синусоидальных полупериодов, соединяющих известные равноудаленные реверсии. Сигнал дискретизируется при частоте 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')