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

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

fs = 600;
t1 = 5;
t = 0:1/fs:t1;

f0 = 10;
f1 = 40;
rpm = 60*linspace(f0,f1,length(t));

Сигнал состоит из четырех гармонично связанных щебетов с порядками 1, 1/2, √2, и 2. Амплитуды щебетов равняются 1, 1/2, √2, и 2, соответственно. Чтобы сгенерировать щебеты, используйте метод трапеций, чтобы выразить фазу как интеграл скорости вращения.

ord = [1 0.5 sqrt(2) 2];
amp = [1 0.5 sqrt(2) 2];

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

x(1,:) = amp(1)*cos(ord(1)*ph);
x(2,:) = amp(2)*cos(ord(2)*ph);
x(3,:) = amp(3)*cos(ord(3)*ph);
x(4,:) = amp(4)*cos(ord(4)*ph);

xsum = sum(x);

Восстановите формы волны временного интервала, которые составляют сигнал.

xrec = orderwaveform(xsum,fs,rpm,ord);

Визуализируйте результаты. Увеличьте масштаб временного интервала, происходящего после того, как переходные процессы затухнут.

for kj = 1:4
    subplot(2,2,kj)
    plot(t,x(kj,:),t,xrec(:,kj))
    title(['Order = ' num2str(ord(kj))])
    xlim([2 3])
end

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

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

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

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

x = sum(y,2);

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

rpmfreqmap(x,fs,rpm1)

Восстановите формы волны временного интервала, которые составляют сигнал. Используйте алгоритм Волд-Кальмана, чтобы разъединить пересекающиеся порядки.

xrec = orderwaveform(x,fs,[rpm1 rpm2],[1.2 0.8],[1 2],'Decouple',true);

Постройте исходные и восстановленные формы волны.

for kj = 1:2
    figure(kj)
    subplot(2,1,1)
    plot((0:nsamp-1)/fs,y(:,kj))
    legend('Original')
    title(['Motor ' int2str(kj)])
    subplot(2,1,2)
    plot((0:nsamp-1)/fs,xrec(:,kj))
    legend('Reconstructed')
end