【参数识别】频域MIMO模态参数辨识（Matlab实现）

“在代码的海洋里，有无尽的知识等待你去发现。我就是那艘领航的船，带你乘风破浪，驶向代码的彼岸。

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

💥1 概述

 频域 MIMO 模态参数辨识是一种先进的系统动态特性分析技术。MIMO（多输入多输出）意味着系统具有多个输入和多个输出，这增加了系统的复杂性和分析的难度。 在频域中进行 MIMO 模态参数辨识时，通过对多个输入信号和相应的多个输出信号在频率域上的测量和分析，来确定系统的模态参数。这种方法充分考虑了输入输出之间的相互关系，能够更全面、准确地描述系统的动态行为。 它通常需要复杂的数学算法和强大的计算能力来处理大量的频域数据。常见的技术包括频响函数矩阵估计、奇异值分解、子空间辨识等，以提取系统的固有频率、阻尼比和模态振型等关键参数。 频域 MIMO 模态参数辨识在众多领域具有重要应用价值。例如，在航空航天领域，可用于飞机结构的振动特性分析，提高飞行安全性和舒适性；在汽车工程中，有助于优化车辆的悬挂系统和底盘设计；在机械制造领域，能够帮助检测和诊断复杂机械装备的故障。 然而，该方法也面临一些挑战，如测量误差对结果的影响、复杂系统建模的准确性以及对高性能计算资源的需求等。但随着技术的不断进步，频域 MIMO 模态参数辨识将在系统动态特性研究中发挥越来越重要的作用。 

📚2 运行结果

主函数部分代码：

clear all
close all
clc

% file test using 4 dofs system with proportional damping
% for MPIFD_MIMO algorithm
%
% modal parameters of the analytical 4 dofs system
%  -- f --   - xi -
%   0.0494   0.0177
%   0.1203   0.0104
%   0.2037   0.0103
%   0.3059   0.0122
%
% load frf in matrix form: one column = one FRF [3001x4]
load('FRF_tot.mat')
% load angular frequency in rad/s and in vector form [3001x1]
load('omg_tot.mat')

[nbl, nbc] = size(frf) ;

% selection of the FRF range using frequency betwen 0.02 and 0.38 Hz
[w1, frf1] = select_frf(w, frf, 0.02,  0.38) ;

% FRFs plot with frequency range selection, all the next plot will be plotted in
% figure 1 (stabilization chart, modal FRF, poles).
figure(1)
hold on
frf1log = 20*log10(abs(frf1)) ;
plot(w1/2/pi, frf1log,'LineWidth', 2)


% identification using stabilization chart between order 2 and 20
order = 2:20 ;
[fn, xin, frfnumtot, FST, FF, XIXI, MATHP] = lscf(w1, frf1, order) ;
% plot stabilization chart
plotstabchart(w1, frf1, FST, FF, XIXI, MATHP, order) ;

% selection of stable poles using stabilization chart, returns poles repeated 
% nbit = 6 times in the tolerance interval tolf = 0.01.
[fns, xins] = select_stabchart(fn, xin, FST, 0.01, 6) ;

% calcul of residues and modal frf using selected poles and FRF between 
% force and displacement ('disp') for the modes specified in idx
idx = 1:4 ;
[rntot, LRtot, URtot, frfmtot] = lsfd(w1, frf1, fns(idx), xins(idx), 'disp', 1) ;

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]姚艳春. 复杂激励下玉米收获机框架结构振动特性与模态参数辨识方法研究[D].中国农业大学,2018.

[2]周思达,曹博远,周小陈.结构应变模态参数辨识的最小二乘复频域方法[J].噪声与振动控制,2017,37(06):17-22.

🌈4 Matlab代码实现