Создайте симулированный сигнал, произведенный на уровне 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, 0.5, 4, и 6. Щебет порядка 4 имеет дважды амплитуду других. Чтобы сгенерировать щебеты, используйте метод трапеций, чтобы описать фазу как интеграл скорости вращения.

o1 = 1;
o2 = 0.5;
o3 = 4;
o4 = 6;

ph = 2*pi*cumtrapz(rpm/60)/fs;

x = [1 1 2 1]*cos([o1 o2 o3 o4]'*ph);

Визуализируйте карту об/мин порядка сигнала.

rpmordermap(x,fs,rpm)

Анализируйте симулированные данные из акселерометра, помещенного в кабину вертолета.

Загрузите вертолетные данные. Вибрационные измерения, vib, производятся на уровне 500 Гц в течение 10 секунд. Контроль данных показывает, что это имеет линейный тренд. Удалите тренд, чтобы препятствовать тому, чтобы он ухудшил качество оценки порядка.

load('helidata.mat')

vib = detrend(vib);

Постройте нелинейный профиль об/мин. Запуски ротора вплоть до него достигают максимальной скорости вращения приблизительно 27 600 оборотов в минуту и затем двигаются вперед без усилий.

plot(t,rpm)
xlabel('Time (s)')
ylabel('RPM')

Вычислите карту об/мин порядка. Задайте разрешение порядка 0,015.

[map,order,rpmOut,time] = rpmordermap(vib,fs,rpm,0.015);

Визуализируйте карту.

imagesc(time,order,map)
ax = gca;
ax.YDir = 'normal';
xlabel('Time (s)')
ylabel('Order')

Повторите расчет с помощью более прекрасного разрешения порядка. Постройте карту с помощью встроенной функциональности rpmordermap. Низшие порядки разрешены более ясно.

rpmordermap(vib,fs,rpm,0.005)

Сгенерируйте сигнал, который состоит из двух линейных щебетов и квадратичного щебета, все произведенные на уровне 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 порядков. Окно данные с Окном Чебышева, затухание бокового лепестка которого составляет 80 дБ.

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

Чертите карту об/мин порядка как график водопада.

[OR,RP] = meshgrid(or,rp);
waterfall(OR,RP,map')

view(-15,45)
xlabel('Order')
ylabel('RPM')
zlabel('Amplitude')

Постройте интерактивную карту об/мин порядка путем вызова rpmordermap без выходных аргументов.

load helidata.mat
rpmordermap(detrend(vib),fs,rpm,0.005)