MATLAB_Simulink建模与仿真[张德丰][程序源代码]

/* 文件名为xiu_exp11.c * 利用C MEX模板文件编写S函数，实现单位延迟（1/z）积分（1/s）的连接离散混合系统 */ #define S_FUNCTION_NAME mixedm #define S_FUNCTION_LEVEL 2 #include "simstruc.h" #include U(element)(*uPtrs[element]) /*====================* * S-函数方法 * *====================*/ /* 函数 mdlInitializeSizes ====================== * Simulink利用该信息决定S函数模块的特性（如输入、输出和状态等变量数目等） */ static void mdlInitializeSizes(SimStruct *S) { ssSetNumSFcnParams(S, 0); /*需要参数的数目*/ if (ssGetNumSFcnParams(S) != ssGetSFcnParamsCount(S)) { /*当需要参数与实际参数数目相等时，则返回*/ return; } ssSetNumContStates(S, 1); ssSetNumDiscStates(S, 1); ssSetNumRWork(S,1); /*对积分结果进行延迟操作*/ if (!ssSetNumInputPorts(S, 1)) return; ssSetInputPortWidth(S, 0, 1); ssSetInputPortDirectFeedThrough(S, 0, 1); ssSetInputPortSampleTime(S,0,CONTINUOUS_SAMPLE_TIME); ssSetInputPortOffsetTime(S,0,0.0); if (!ssSetNumOutputPorts(S, 1)) return; ssSetOutputPortWidth(S, 0, 1); ssSetOutputPortSampleTimes(S,0,1); ssSetOutputPortOffsetTime(S,0,0); ssSetNumNonsampledZCs(S, 2); ssSetOptions(S, (SS_OPTION_EXCEPTION_FREE_CODE| SS_OPTION_PORT_SAMPLE_TIMES_ASSIGNED)); } /* 函数mdlInitializeSampleTimes ========================= * 指定S函数的连接和离散系统的采样时间都为1.0 */ static void mdlInitializeSampleTimes(SimStruct *S) { ssSetSampleTime(S, 0, CONTINUOUS_SAMPLE_TIME); ssSetOffsetTime(S, 0, 0.0); ssSetSampleTime(S,1,1.0); ssSetOffsetTime(S, 1, 0.0); ssSetModelReferenceSampleTimeDefaultInheritance(S); } #define MDL_INITIALIZE_CONDITIONS /*改变到#undef以移除函数*/ #if defined(MDL_INITIALIZE_CONDITIONS) /* 函数 mdlInitializeConditions ============================ *初始化连续和离散系统的状态为1 */ static void mdlInitializeConditions(SimStruct *S) { real_T*xC0=ssGetContStates(S); real_T*xD0=ssGetRealDiscStates(S); xCO[0]=1.0; xD0[0]=1.0; } /*函数mdlOutputs ============================================= * 计算输出y=xD,同时更新积分器的内部输出结果 */ static void mdlOutputs(SimStruct *S, int_T tid) { /*更新积分器内部输出结果"zoh"*/ if(ssIsContinuousTask(S,tid)){ real_T*zoh=ssGetRWork(S); real_T*xC=ssGetContStates(S) } } /*计算y=xD*/ if(ssIsSampleHit(S,1,tid)){ real_T*y=ssGetOutputPortRealSignal(S,0); real_T*xD=ssGetRealDiscStates(S); y[0]=xD[0]; } } #define MDL_UPDATE /*改变#undef以移除函数*/ #if defined(MDL_UPDATE) /*函数mdlUpdate =========================================== *计算xD=xC */ static void mdlUpdate(SimStruct *S, int_T tid) { UNUSED_ARG(tid); %计算xD=xC if(ssIsSamleHit(S,1,tid)){ real_T*xD=ssGetRealDiscStates(S); real_T*zoh=ssGetRWorks(S); xD[0]=*zoh; } } #define MDL_DERIVATIVES /*改变#undef以移除函数*/ #if defined(MDL_DERIVATIVES) /* 函数mdlDerivatives ========================================= *计算连续模型系统状态的导数xdot=U */ static void mdlDerivatives(SimStruct *S) { real_T*dx=ssGetdX(S); InputRealPtrsType uPtrs=ssGetInputPortRealSignalPtrs(S,0); /*计算xdot=U*/ dx[0]=U[0]; } /*函数mdlTerminate ========================================== *处理仿真结 */ static void mdlTerminate(SimStruct *S) { UNUSED_ARG(S); /*没有使用的输入参数*/ } #ifdef MATLAB_MEX_FILE /* 该文件被编译为MEX文件? */ #include "simulink.c" /*MEX文件接口机制*/ #else #include "cg_sfun.h" /*代码产生注册函数*/ #endif