利用噪声估计频率响应模型系统识别工具箱

打开Simulin万博1manbetxk®模型,并指定要线性化的模型部分:

加载系统(“magball”);io(1)=linio(“磁球/所需高度”,1);io(2)=linio(“磁球/磁球装置”1.“输出”);

计算稳态工作点,并将模型线性化:

op=findop(“magball”,operspec(“magball”),...findopOptions(“显示报告”,“关”)); sys=线性化(“magball”,io,op);

创建一个啁啾信号,并使用它来估计频率响应:

in=菲斯特啁啾(“频率范围”,[1 1000],...“Ts”,0.001,...“NumSamples”,1e4);[~,simout]=frestimate(“magball”,io,op,in);

使用系统识别工具箱™ 估算五阶状态空间模型的软件。比较分析线性化和状态空间模型的结果:

input=generateTimeseries(in);output=simout{1}.Data;Data=iddata(output,input.Data(:),in.Ts);sys_id=n4sid(detrend(Data),5,“cov”,“没有”); bodemag(系统,ss(系统id(“测量”)),“r”)传奇(“使用线性化获得的线性模型”,...“使用系统识别工具箱的状态空间模型”,...“位置”,“西南”)