Сгенерируйте сигнал, который состоит из двух линейных щебета и квадратичного щебета, все дискретизированные на частоте 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