Сгенерируйте диаграмму стабилизации для модального анализа
modalsd( генерирует диаграмму стабилизации на текущей фигуре. frf,f,fs)modalsd оценивает естественные частоты и коэффициенты затухания от 1 до 50 режимов и генерирует схему с помощью алгоритма наименьших квадратов комплексной экспоненциальной (LSCE). fs - частота дискретизации. Частота, f, - вектор с количеством элементов, равным количеству строк функции частотной характеристики, frf. Можно использовать эту схему, чтобы дифференцировать вычислительный и физический режимы.
modalsd( задает опции, используя аргументы пары "имя-значение".frf,f,fs,Name,Value)
fn = modalsd(___)fn, идентифицированный как стабильный между последовательными порядками моделей. i-й элемент содержит длинно- i вектор собственных частот устойчивых полюсов. Поляки, которые не являются стабильными, возвращаются как NaNs. Этот синтаксис принимает любую комбинацию входов из предыдущих синтаксисов.
Вычислите функции частотной характеристики для системы с двумя входами/двумя выходами, возбужденными случайным шумом.
Загрузите файл данных. Вычислите функции частотной характеристики с помощью 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] Брандт, Андерс. Анализ шума и вибрации: анализ сигналов и экспериментальные процедуры. Chichester, UK: John Wiley & Sons, 2011.
[2] Оздемир, Ахмет Арда и Суат Гумуссой. Оценка передаточной функции в System Identification Toolbox™ через Вектор Подбора кривой. Материалы XX Всемирного конгресса Международной федерации автоматического управления, Тулуза, Франция, июль 2017 года.
[3] Vold, Hovard, John Crowley, and G. Thomas Rocklin. Новые способы оценки функций частотной характеристики. Звук и вибрация. Том 18, ноябрь 1984, стр. 34-38.