Вычислите функции частотной характеристики для системы с двумя входами/двумя выходами, возбужденными случайным шумом.

Загрузите файл данных. Вычислите функции частотной характеристики с помощью 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')