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

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

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. Амплитуды щебёнок составляют 1, 1/2, √ 2 и 2 соответственно. Чтобы сгенерировать щебет, используйте метод трапеций, чтобы выразить фазу как интеграл скорости вращения.

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

a1 = 1;
a2 = 0.5;
a3 = sqrt(2);
a4 = 2;

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

x = [a1 a2 a3 a4]*cos([o1 o2 o3 o4]'*ph);

Извлеките и визуализируйте величины порядков.

ordertrack(x,fs,rpm,[o1 o2 o3 o4])

Создайте имитированный сигнал вибрации, состоящий из двух порядков пересечения, соответствующих двум разным двигателям. Дискретизация сигнала производится на частоте 300 Гц в течение 3 секунд. Первый двигатель увеличивает скорость вращения с 10 до 100 оборотов в секунду (или, что эквивалентно, с 600 до 6000 оборотов в минуту) во время измерения. Второй двигатель увеличивает скорость вращения с 50 до 70 оборотов в секунду (или с 3000 до 4200 оборотов в минуту) за тот же период.

fs = 300;
nsamp = 3*fs;

rpm1 = linspace(10,100,nsamp)'*60;
rpm2 = linspace(50,70,nsamp)'*60;

Измеренный сигнал имеет порядок 1.2 и амплитуду 2√2 относительно первого двигателя. Что касается второго двигателя, то сигнал имеет порядок 0.8 и амплитуду 4√2.

x = [2 4]*sqrt(2).*cos(2*pi*cumtrapz([1.2*rpm1 0.8*rpm2]/60)/fs);

Заставьте первый двигатель возбуждать резонанс в середине частотной области значений.

rs = [1+1./(1+linspace(-10,10,nsamp).^4)'/2 ones(nsamp,1)];

x = sum(rs.*x,2);

Визуализируйте порядки с помощью rpmfreqmap.

rpmfreqmap(x,fs,rpm1)

Вычислите амплитуды порядка для обоих двигателей как функцию от об/мин. Используйте алгоритм Vold-Kalman, чтобы развязать порядки пересечения.

ordertrack(x,fs,[rpm1 rpm2],[1.2 0.8],[1 2],'Decouple',true)

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

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

load('helidata.mat')

vib = detrend(vib);

Вычислите карту order-RPM. Задайте разрешение порядка 0,005.

[map,order,rpm,time,res] = rpmordermap(vib,fs,rpm,0.005);

Вычислите и постройте график среднего спектра порядка сигнала. Найдите три самых высоких peaks спектра.

[spectrum,specorder] = orderspectrum(map,order);

[~,pkords] = findpeaks(spectrum,specorder,'SortStr','descend','Npeaks',3);

findpeaks(spectrum,specorder,'SortStr','descend','Npeaks',3)

Отслеживайте амплитуды трех самых высоких peaks.

ordertrack(map,order,rpm,time,pkords)