打开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”)传奇(“使用线性化获得的线性模型”,...“使用系统识别工具箱的状态空间模型”,...“位置”,“西南”)