Вычислить и отобразить карту порядка - RPM

Создайте сигнал, который состоит из двух линейных чирп и квадратичной чирп, все дискретизированные при 600 Гц в течение 5 секунд. Система, которая выдает сигнал, увеличивает скорость его вращения с 10 до 40 оборотов в секунду в течение периода испытаний.

Создайте показания тахометра.

fs = 600;
t1 = 5;
t = 0:1/fs:t1;
f0 = 10;
f1 = 40;
rpm = 60*linspace(f0,f1,length(t));

Линейные чирки имеют порядки 1 и 2,5. Компонент с порядком 1 имеет вдвое большую амплитуду, чем другой. Квадратичная чирп начинается с порядка 6 и возвращается к этому порядку в конце измерения. Его амплитуда составляет 0,8. Создайте сигнал, используя эту информацию.

o1 = 1;
o2 = 2.5;
o6 = 6;
x = 2*chirp(t,o1*f0,t1,o1*f1)+chirp(t,o2*f0,t1,o2*f1) + ...
    0.8*chirp(t,o6*f0,t1,o6*f1,'quadratic');

Вычислите карту порядка-оборотов сигнала. Используйте амплитуду пика в каждой измерительной ячейке. Укажите разрешение 0,25 порядка. Окно данных с чебышёвским окном, затухание боковых зон которого составляет 50 дБ.

[map,or,rp] = rpmordermap(x,fs,rpm,0.25, ...
    'Amplitude','peak','Window',{'chebwin',50});

Нарисуйте карту order-RPM как график водопада.

[OR,RP] = meshgrid(or,rp);
waterfall(OR,RP,map')
view(-15,45)
xlabel('Order')
ylabel('RPM')
zlabel('Amplitude')

Использовать rpmordermap без выходных аргументов для отображения карты. Задайте разрешение 0,2 порядка и 80% перекрытия между соседними сегментами. Установите затухание боковины окна Чебышева на 80 дБ.

rpmordermap(x,fs,rpm,0.2, ...
    'Amplitude','peak','OverlapPercent',80,'Window',{'chebwin',80})

См. также

orderspectrum | ordertrack | orderwaveform | rpmfreqmap | rpmordermap | tachorpm