信号与系统课程实验3.doc

M仿真三 复杂温控系统调试仿真 1仿真目的: 学习系统仿真中参数调试与优化方法. 系统模型 如图3.1。被控对象含非最小相位环节和延迟环节.性能更复杂,控制难度更大. 控制器（c1*s+c2）/s 有两个参数，需要调试优化使系统动态响应最佳。 为此，建立动态响应性能指标 J=∫|偏差e|dt，J越小，动态性能越好。 用绝对值与积分模块实现并贴输出标签,便于调出数据分析。 图3.1 建立系统仿真模型文件 按图3.1选择模块并连接. 其中绝对值模块取自”Math Operation”包, 输出标签取自”Sinks”包. 基本参数设定: 模块/项目 设定参数 给定 阶跃值:3 限幅 输出限幅值: +1,-1 延迟环节 延迟时间τ=2s. 其它 如图3.1 仿真 时间?: 300s (3) 将仿真文件取名如mfzmx3保存. (4) 在命令窗口输入命令仿真 c1=0.2; c2=0.3;sim(‘mfzmx3’) . 并打开示波器,纪录仿真结果 改变参数, 再试. c1=1;c2=0.1;sim(mfzmx3) 表3.2 手工仿真记录 参数 c1=0.2 c2=0.3; 性能评价: 参数 c1=1 c2=0.1; 性能评价: 当M命令窗口信息过多时, 可用”Edit”清除. 4 计算机辅助参数优化 建立仿真参数优化程序： c1l=1.3504 ; c1h=1.3536; % 参数c1，c2的寻优范围 c2l=0.0112 ; c2h=0.0128; d1=(c1h-c1l)/10;c1x=c1l:d1:c1h; % 建立c1,c2 寻优向量 d2=(c2h-c2l)/10;c2x=c2l:d2:c2h; J=inf; % J为最小化指标,初值+∞ for i=1:length(c1x) % 仿真,求J最小的c1,c2 for k=1:length(c2x) % length(cx) 为cx的长度 c1=c1x(i);c2=c2x(k); [t,x,j]=sim(mfzmx3); % 对mfzmx3 模型仿真 % t为时间向量;x为仿真使用的状态向量;j为仿真模型标签的输出向量. Lj=length(j); if Jj(Lj) J=j(Lj);C1=c1;C2=c2;end end end J % 指标优化值 C1,[C1-d1,C1+d1] % 参数优化值与下次优化参考范围 C2,[C2-d2,C2+d2], c1=C1;c2=C2; % 用优化参数值再仿真 sim(mfzmx3); 取名“mfzcx3”保存。 程序将输入的c1和c2的选择范围,形成参数向量c1x和c2x, 并对全部参数组合进行仿真,求出优化值J和C1,C2及对应的仿真结果. 新寻优范围可用优化值前后2点,也可自行选择.继续优化,直到性能满意为止. 该程序简单, 有一定适应性. 专门的M函数优化方法请自选研究. 每次优化要作多次仿真, 当M命令窗口信息过多时, 可用”Edit”清除. 程序调试时,也应用”Edit”清数据库”Workspace”,以免受旧数据影响. 表3.2 计算机辅助参数优化过程 范围与优化值 C1:[ 0, 2](1.8 c2:[ 0,1 ](0.1 响应曲线 范围与优化值 c1:[1.6 2.0](1.6 c2:[ 0 0.2](0.02 响应曲线 范围与优化值 c1:[1.56 1.64]( 1.56 c2:[0 0.04 ]( 0.024 响应曲线 范围与优化值 c1:[1.552 1.568]( 1.552 c2:[ 0.02 0.028 ]( 0.0216 响应曲线 范围与优化值 c1:[1.55