Сгенерируйте диаграмму стабилизации для модального анализа
modalsd( генерирует диаграмму стабилизации в текущей фигуре. frf,f,fs)modalsd оценивает собственные частоты и коэффициенты затухания от 1 до 50 режимов и генерирует схему с помощью алгоритма комплексной экпоненты наименьших квадратов (LSCE). fs частота дискретизации. Частота, f, вектор со многими элементами, равными количеству строк функции частотной характеристики, frf. Можно использовать эту схему, чтобы дифференцироваться между вычислительными и физическими режимами.
modalsd( задает аргументы пары "имя-значение" использования опций.frf,f,fs,Name,Value)
fn = modalsd(___)fn, идентифицированный как являющийся устойчивым между последовательными порядками модели. i th элемент содержит вектор длины-i из собственных частот устойчивых полюсов. Полюса, которые не устойчивы, возвращены как NaNs. Этот синтаксис принимает любую комбинацию входных параметров от предыдущих синтаксисов.
Вычислите функцию частотной характеристики для two-input/two-output системы, взволнованной случайным шумом.
Загрузите файл данных. Вычислите функцию частотной характеристики с помощью окна Hann с 5000 выборками и 50%-го перекрытия между смежными сегментами данных. Укажите, что выходные измерения являются смещениями.
load modaldata winlen = 5000; [frf,f] = modalfrf(Xrand,Yrand,fs,hann(winlen),0.5*winlen,'Sensor','dis');
Сгенерируйте диаграмму стабилизации, чтобы идентифицировать до 20 физических режимов.
modalsd(frf,f,fs,'MaxModes',20)
Повторите расчет, но теперь сожмите критерии устойчивости. Классифицируйте данный полюс, столь же устойчивый в частоте, если ее собственная частота изменяется меньше чем на 0,01% как увеличения порядка модели. Классифицируйте данный полюс, столь же устойчивый в затухании, если оценка коэффициента затухания изменяется меньше чем на 0,2% как увеличения порядка модели.
modalsd(frf,f,fs,'MaxModes',20,'SCriteria',[1e-4 0.002])
Ограничьте частотный диапазон между 0 и 500 Гц. Ослабьте критерии стабильности к 0,5% для частоты и 10% для затухания.
modalsd(frf,f,fs,'MaxModes',20,'SCriteria',[5e-3 0.1],'FreqRange',[0 500])
Повторите расчет с помощью алгоритма рациональной функции наименьших квадратов. Ограничьте частотный диапазон от 100 Гц до 350 Гц и идентифицируйте до 10 физических режимов.
modalsd(frf,f,fs,'MaxModes',10,'FreqRange',[100 350],'FitMethod','lsrf')
[1] Брандт, Андерс. Шум и анализ вибрации: анализ сигнала и экспериментальные процедуры. Чичестер, Великобритания: John Wiley & Sons, 2011.
[2] Ozdemir, Ахмет Арда и Суэт Гумассой. "Оценка Передаточной функции System Identification Toolbox™ через Подбор кривой Вектора". Продолжения 20-го Мирового Конгресса Международной федерации Автоматического управления, Тулузы, Франция, июль 2017.
[3] Vold, Håvard, Джон Кроули и Г. Томас Роклин. “Новые Способы Оценить Функции Частотной характеристики”. Звук и Вибрация. Издание 18, ноябрь 1984, стр 34–38.