第7章S函数论述.ppt

　　运行仿真，在示波器中观察结果波形，图7-15(b)中示波器上面的窗口显示的是输入信号波形，下面的窗口显示的是S-函数的输出波形。 (a) 　　 (b) 图7-15 　　用户也可以封装S-Function模块，这样可以在模块图标上标明该函数的功能，如在封装编辑器的Icon选项页内输入绘制命令：disp(Unit Delay)，那么所得S-Function模块如图7-16所示。 图7-16 　　sys = simsizes(sizes); 　　 　　x0 = 0; 　　str = [ ]; 　　ts = [0.1 0]; % Sample period of 0.1 seconds (10Hz) 　　 　　% end mdlInitializeSizes 　　%============================================================ 　　% mdlUpdate 　　% Handle discrete state updates, sample time hits, and major time step requirements. 　　%============================================================ 　　function sys = mdlUpdate(t, x, u) 　　sys = u; 　　%end mdlUpdate　　%=============================================　　% mdlOutputs　　% Return the output vector for the S-function　　%=============================================　　function sys = mdlOutputs(t, x, u)　　sys = x;　　%end mdlOutputs 　　例7-6 人口系统S-函数。　　编写一个S-函数描述人口的动态变化。设繁殖率为r，资源为k，初始人口数量为init，则人口变化规律为 系统模型如图7-17所示。 图7-17 　　修改S-函数的第一行：　　function [sys, x0, str, ts] = sfun_population(t, x, u, flag, r, k, init)　　对子函数的调用为：　　case 0,　　 [sys, x0, str, ts] = mdlInitializeSizes(init);　　case 2,　　 sys = mdlUpdate(t, x, u, r, k);　　在这个例子中，p为状态，输出等于状态，于是，更新函数和输出函数为：　　function sys = mdlUpdate(t, x, u, r, k)　　sys = [r*x*(1-x/k)];　　function sys = mdlOutputs(t, x, u)　　sys = [x]; 　　打开S-Function对话框，设置S-function name参数为sfun_population，设置S-function parameter参数值为1.05、1e6、100000，如图7-18(a)所示。选择变步长ode45求解器，仿真100个时间单位，运行仿真后得到的结果曲线如图7-18(b)所示。 (a) 图7-18 (b) 图7-18 　　下面是sfun_population函数文件完整的M代码：　　function [sys,x0,str,ts] = sfun_population(t, x, u, flag, r, K, init)　　%population dynamics. Rate(r), resources(K) and initial population　　%(init) are user inputs.　　switch flag,　　 case 0,　　 [sys,x0,str,ts] = mdlInitializeSizes(init);　　 case 1,　　 sys = mdlDerivatives(t, x, u);　　 case 2,　　 sys = mdlUpdate(t, x, u, r, K);　　 case 3,　　 sys = mdlOutputs(t, x, u);　　 case 4, 　　 sys = mdlGetTimeOfNext