Вычислите функцию частотной характеристики для 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')