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

В 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];
}