北航MATLAB第六次作业

1、1、 当电源电压，求解表达式的波形，且，其中。用simulink动态仿真和S函数两种方法，建立系统仿真模型，求解的波形。（1）（2）function sys,x0,str,ts,simStateCompliance = sfun1(t,x,u,flag) switch flag, case 0, sys,x0,str,ts,simStateCompliance=mdlInitializeSizes; case 1, sys=mdlDerivatives(t,x,u); case 2, sys=mdlUpdate(t,x,u); case 3, sys=mdlOutputs(t,x,u); ca

2、se 4, sys=mdlGetTimeOfNextVarHit(t,x,u); case 9, sys=mdlTerminate(t,x,u); otherwise DAStudio.error('Simulink:blocks:unhandledFlag', num2str(flag); end function sys,x0,str,ts,simStateCompliance=mdlInitializeSizes sizes = simsizes; sizes.NumContStates = 1;sizes.NumDiscStates = 0;sizes.NumOutpu

3、ts = 1;sizes.NumInputs = 1;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1; sys = simsizes(sizes); x0 = 1; str = ; ts = 0 0; simStateCompliance = 'UnknownSimState' function sys=mdlDerivatives(t,x,u) sys =4*100*sin(100*pi*t+30)+1-3*x; function sys=mdlUpdate(t,x,u) sys = ; function sys=mdlOu

4、tputs(t,x,u) sys = x; function sys=mdlGetTimeOfNextVarHit(t,x,u) sampleTime = 1; sys = t + sampleTime; function sys=mdlTerminate(t,x,u) sys = ;2、 已知末系统的状态方程如下：试用simulink建立仿真模型，且其输入为单位阶跃响应曲线。（1）（2）functionsys,x0,str,ts=sfun3(t,x,u,flag)switch flag,case 0,sys,x0,str,ts=mdlInitializeSizes;case 1,sys= m

5、dlDerivatives(t,x,u);case 2,sys= mdlUpdate(t,x,u);case 3,sys= mdlOutputs(t,x,u);case 4,sys=mdlGetTimeOfNextVarHit(t,x,u);case 9,sys=mdlTerminate(t,x,u);otherwiseDAStudio.error('Simulink:blocks:unhandledFlag',num2str(flag);endfunctionsys,x0,str,ts=mdlInitializeSizessizes=simsizes;sizes.NumCon

6、tStates = 3;sizes.NumDiscStates = 0;sizes.NumOutputs = 1;sizes.NumInputs = 1;sizes.DirFeedthrough = 0;sizes.NumSampleTimes = 1;sys=simsizes(sizes);x0=0;0;0;str=;ts=0 0;function sys=mdlDerivatives(t,x,u)x(1)=1*x(1)-1*x(2)-2*x(3)+u;x(2)=2*x(2)-1*x(3);x(3)=3*x(3)+u;sys=x;function sys=mdlUpdate(t,x,u)sy

7、s=;function sys=mdlOutputs(t,x,u)sys=1*x(1)-2*x(2)+3*x(3);function sys=mdlGetTimeOfNextVarHit(t,x,u)sampleTime=1sys=t+sampleTime;function sys=mdlTerminate(t,x,u)sys=;3、 食饵-捕食者模型：设食饵（如鱼，兔等）数量为，捕食者数量为，有设r=1, d=0.5, a=0.1, b=0.02, .用simulink动态仿真和S函数两种方法，建立系统仿真模型，绘制，的图形。function sys,x0,str,ts,simStateCo

8、mpliance = sfuntmpl2(t,x,u,flag) switch flag, case 0, sys,x0,str,ts,simStateCompliance=mdlInitializeSizes; case 1, sys=mdlDerivatives(t,x,u); case 2, sys=mdlUpdate(t,x,u); case 3, sys=mdlOutputs(t,x,u); case 4, sys=mdlGetTimeOfNextVarHit(t,x,u); case 9, sys=mdlTerminate(t,x,u); otherwise DAStudio.er

9、ror('Simulink:blocks:unhandledFlag', num2str(flag); end function sys,x0,str,ts,simStateCompliance=mdlInitializeSizes sizes = simsizes; sizes.NumContStates = 2;sizes.NumDiscStates = 0;sizes.NumOutputs = 2;sizes.NumInputs = 0;sizes.DirFeedthrough = 0;sizes.NumSampleTimes = 1; sys = simsizes(si

10、zes); x0 = 25,2; str = ; ts = 0 0; simStateCompliance = 'UnknownSimState' function sys=mdlDerivatives(t,x,u) dx=zeros(2,1);dx(1)=x*(r-a*y);dx(2)=y*(-b+d*x);sys =dx; function sys=mdlOutputs(t,x,u) sys = x;4、 已知，求该方程的解。用simulink动态仿真和S函数两种方法，建立系统仿真模型，绘制，的图形。function sys,x0,str,ts,simStateCompli

11、ance = sfuntmpl2(t,x,u,flag) switch flag, case 0, sys,x0,str,ts,simStateCompliance=mdlInitializeSizes; case 1, sys=mdlDerivatives(t,x,u); case 2, sys=mdlUpdate(t,x,u); case 3, sys=mdlOutputs(t,x,u); case 4, sys=mdlGetTimeOfNextVarHit(t,x,u); case 9, sys=mdlTerminate(t,x,u); otherwise DAStudio.error(

12、'Simulink:blocks:unhandledFlag', num2str(flag); end function sys,x0,str,ts,simStateCompliance=mdlInitializeSizes sizes = simsizes; sizes.NumContStates = 1;sizes.NumDiscStates = 0;sizes.NumOutputs = 2;sizes.NumInputs = 1;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1; sys = simsizes(sizes); x0 = 25; str = ; ts = 0 0; simStateCompliance = 'UnknownSimState' function sys=mdlDerivatives(t,x,u) sys =4*100*sin(100*pi*t+30)+1-3*x; function sys=mdlUpdate(t,x,u) sys = ; function sys=mdlOutputs(t,x,u) sys = x; function sys=mdlGetT