Документация

В mdlInitializeSizes, укажите количество векторов DWork с помощью ssSetNumDWork макро. Например, чтобы указать, что S-функция содержит два вектора DWork, используйте команду

ssSetNumDWork(S, 2);

Хотя mdlInitializeSizes метод сообщает движку Simulink ®, сколько векторов DWork будет использовать S-функция, в данный момент движок не выделяет память для векторов DWork.

S-функция может отложить задание количества векторов DWork до тех пор, пока не будет доступна вся информация о входах S-функции путем передачи значения. DYNAMICALLY_SIZED в ssSetNumDWork макро. Если S-функция откладывает задание количества векторов DWork в mdlInitializeSizes, он должен обеспечить mdlSetWorkWidths для настройки векторов DWork.

Если S-функция не предоставляет mdlSetWorkWidths способ, mdlInitializeSizes задает любые применимые атрибуты для каждого вектора DWork. Например, следующие строки инициализируют ширину и типы данных векторов DWork, инициализированных на предыдущем шаге.

ssSetDWorkWidth(S, 0, 2);
ssSetDWorkWidth(S, 1, 1);

ssSetDWorkDataType(S, 0, SS_DOUBLE);
ssSetDWorkDataType(S, 1, SS_BOOLEAN);

В следующей таблице перечислены атрибуты, которые можно задать для вектора DWork, и показан пример макроса, который его задает. Посмотрите ssSetDWorkRTWStorageClass список поддерживаемых классов хранения.

В mdlStartинициализируйте значения всех векторов DWork, которые должны быть установлены только в начале моделирования. Используйте ssGetDWork для извлечения указателя на каждый вектор DWork и инициализации значений. Например, следующее: mdlStart инициализирует первый вектор DWork.

static void mdlStart(SimStruct *S)
{
    real_T *x = (real_T*) ssGetDWork(S,0);

    /*  Initialize the first DWork vector */
    x[0] = 0;
    x[1] = 2;
}

Модуль Simulink выделяет память для вектора DWork перед вызовом mdlStart способ. Потому что mdlStart метод вызывается только один раз в начале моделирования, не используйте его для данных или состояний, которые необходимо повторно инициализировать, например, при повторном включении отключенной подсистемы, содержащей S-функцию.

В mdlInitializeConditionsинициализируйте значения всех векторов DWork, которые необходимо повторно инициализировать в определенных точках моделирования. Модуль выполняет mdlInitializeConditions в начале моделирования и в любое время снова включается включенная подсистема, содержащая S-функцию. См. раздел mdlStart пример на предыдущем шаге для команд, используемых для инициализации значений вектора DWork.

В mdlOutputs, mdlUpdateи т.д., используйте ssGetDWork для извлечения указателя на вектор DWork и использования или обновления значений вектора DWork. Например, для вектора DWork, хранящего два дискретных состояния, следующее mdlOutputs и mdlUpdate методы вычисляют выходные данные и обновляют значения дискретного состояния.

Ранее определенная S-функция U(element) как (*uPtrs[element]), A, B, C, и D в качестве матриц состояния-пространства для дискретной системы состояния-пространства.

/* Function: mdlOutputs ==============================================
 * Abstract:
 *      y = Cx + Du 
 */
static void mdlOutputs(SimStruct *S, int_T tid)
{
    if( ssGetDWorkUsageType(S, 0) == SS_DWORK_USED_AS_DSTATE) {
        real_T            *y      = ssGetOutputPortRealSignal(S,0);
        real_T            *x      = (real_T*) ssGetDWork(S, 0);
        InputRealPtrsType uPtrs   = ssGetInputPortRealSignalPtrs(S,0);
        
        UNUSED_ARG(tid); /* not used in single tasking mode */
        
    /* y=Cx+Du */
        y[0]=C[0][0]*x[0]+C[0][1]*x[1]+D[0][0]*U(0)+D[0][1]*U(1);
        y[1]=C[1][0]*x[0]+C[1][1]*x[1]+D[1][0]*U(0)+D[1][1]*U(1);
    }
}



#define MDL_UPDATE
/* Function: mdlUpdate ===============================================
 * Abstract:
 *      xdot = Ax + Bu
 */
static void mdlUpdate(SimStruct *S, int_T tid)
{
    real_T            tempX[2] = {0.0, 0.0};
    real_T            *x       = (real_T*) ssGetDWork(S, 0);
    InputRealPtrsType uPtrs    = ssGetInputPortRealSignalPtrs(S,0);

    UNUSED_ARG(tid); /* not used in single tasking mode */

    /* xdot=Ax+Bu */
    tempX[0]=A[0][0]*x[0]+A[0][1]*x[1]+B[0][0]*U(0)+B[0][1]*U(1);
    tempX[1]=A[1][0]*x[0]+A[1][1]*x[1]+B[1][0]*U(0)+B[1][1]*U(1);
 
    x[0]=tempX[0];
    x[1]=tempX[1];
}