用噪声估计频响模型系统识别工具箱
打开Simulin金宝appk®,并指定要线性化模型的哪一部分:
load_system (“magball”);Io (1) = linio(“magball /期望的高度”1);Io (2) = linio(“磁球/磁球工厂”, 1“输出”);
计算稳态工作点,线性化模型:
Op = findop(“magball”operspec (“magball”),...findopOptions (“DisplayReport”,“关闭”));系统=线性化(“magball”io, op);
创建一个啁啾信号,并使用它来估计频率响应:
In =最先。唧唧声(“FreqRange”1000年[1],...“t”, 0.001,...“NumSamples”1 e4);[~,simout] = frestual (“magball”、io、op);
使用系统识别工具箱™软件来估计五阶状态空间模型。比较分析线性化和状态空间模型的结果:
input = generateTimeseries(in);output = simout{1}.Data;data = iddata(output,input.Data(:),in.Ts);Sys_id = n4sid(dettrend (data),5,“浸”,“没有”);bodemag (sys, ss (sys_id (“测量”)),“r”)传说('线性模型获得使用线性化',...“使用系统识别工具箱的状态空间模型”,...“位置”,“西南”)