c4_std5.c

mmo 信道估计算法 matlab 仿真 实现了无线信道环境下的信道估计

static const mxArray *c4_b_sf_marshall(void *c4_chartInstance, void *c4_u)
{
  const mxArray *c4_y = NULL;
  int32_T c4_i22;
  creal_T c4_b_u[52];
  const mxArray *c4_b_y = NULL;
  c4_y = NULL;
  for(c4_i22 = 0; c4_i22 < 52; c4_i22 = c4_i22 + 1) {
    c4_b_u[c4_i22].re = (*(creal_T (*)[52])c4_u)[c4_i22].re;
    c4_b_u[c4_i22].im = (*(creal_T (*)[52])c4_u)[c4_i22].im;
  }
  c4_b_y = NULL;
  sf_mex_assign(&c4_b_y, sf_mex_create(&c4_b_u, "y", 0, 1U, 1U, 1, 52));
  sf_mex_assign(&c4_y, c4_b_y);
  return c4_y;
}

static const mxArray *c4_c_sf_marshall(void *c4_chartInstance, void *c4_u)
{
  const mxArray *c4_y = NULL;
  int32_T c4_i23;
  real_T c4_b_u[52];
  const mxArray *c4_b_y = NULL;
  c4_y = NULL;
  for(c4_i23 = 0; c4_i23 < 52; c4_i23 = c4_i23 + 1) {
    c4_b_u[c4_i23] = (*(real_T (*)[52])c4_u)[c4_i23];
  }
  c4_b_y = NULL;
  sf_mex_assign(&c4_b_y, sf_mex_create(&c4_b_u, "y", 0, 0U, 1U, 2, 1, 52));
  sf_mex_assign(&c4_y, c4_b_y);
  return c4_y;
}

static const mxArray *c4_d_sf_marshall(void *c4_chartInstance, void *c4_u)
{
  const mxArray *c4_y = NULL;
  int32_T c4_i24;
  creal_T c4_b_u[52];
  const mxArray *c4_b_y = NULL;
  c4_y = NULL;
  for(c4_i24 = 0; c4_i24 < 52; c4_i24 = c4_i24 + 1) {
    c4_b_u[c4_i24].re = (*(creal_T (*)[52])c4_u)[c4_i24].re;
    c4_b_u[c4_i24].im = (*(creal_T (*)[52])c4_u)[c4_i24].im;
  }
  c4_b_y = NULL;
  sf_mex_assign(&c4_b_y, sf_mex_create(&c4_b_u, "y", 0, 1U, 1U, 2, 1, 52));
  sf_mex_assign(&c4_y, c4_b_y);
  return c4_y;
}

static const mxArray *c4_e_sf_marshall(void *c4_chartInstance, void *c4_u)
{
  const mxArray *c4_y = NULL;
  int32_T c4_i25;
  creal_T c4_b_u[52];
  const mxArray *c4_b_y = NULL;
  c4_y = NULL;
  for(c4_i25 = 0; c4_i25 < 52; c4_i25 = c4_i25 + 1) {
    c4_b_u[c4_i25].re = (*(creal_T (*)[52])c4_u)[c4_i25].re;
    c4_b_u[c4_i25].im = (*(creal_T (*)[52])c4_u)[c4_i25].im;
  }
  c4_b_y = NULL;
  sf_mex_assign(&c4_b_y, sf_mex_create(&c4_b_u, "y", 0, 1U, 1U, 2, 1, 52));
  sf_mex_assign(&c4_y, c4_b_y);
  return c4_y;
}

const mxArray *sf_c4_std5_get_eml_resolved_functions_info(void)
{
  const mxArray *c4_nameCaptureInfo = NULL;
  c4_nameCaptureInfo = NULL;
  sf_mex_assign(&c4_nameCaptureInfo, sf_mex_create(NULL, "nameCaptureInfo", 0,
    0U, 1U, 2, 0, 1));
  return c4_nameCaptureInfo;
}

static creal_T *c4_channel_estimate(void)
{
  return (creal_T *)ssGetOutputPortSignal(chartInstance.S, 1);
}

static creal_T *c4_freq_tr_syms(void)
{
  return (creal_T *)ssGetInputPortSignal(chartInstance.S, 0);
}

static void init_dsm_address_info(void)
{
}

/* SFunction Glue Code */
void sf_c4_std5_get_check_sum(mxArray *plhs[])
{
  ((real_T *)mxGetPr((plhs[0])))[0] = (real_T)(1710200764U);
  ((real_T *)mxGetPr((plhs[0])))[1] = (real_T)(607161690U);
  ((real_T *)mxGetPr((plhs[0])))[2] = (real_T)(745703221U);
  ((real_T *)mxGetPr((plhs[0])))[3] = (real_T)(1379764459U);
}

mxArray *sf_c4_std5_get_autoinheritance_info(void)
{
  const char *autoinheritanceFields[] =
  {"checksum","inputs","parameters","outputs"};
  mxArray *mxAutoinheritanceInfo =
  mxCreateStructMatrix(1,1,4,autoinheritanceFields);
  {
    mxArray *mxChecksum = mxCreateDoubleMatrix(4,1,mxREAL);
    double *pr = mxGetPr(mxChecksum);
    pr[0] = (double)(202484585U);
    pr[1] = (double)(1218703978U);
    pr[2] = (double)(4233831050U);
    pr[3] = (double)(1992993190U);
    mxSetField(mxAutoinheritanceInfo,0,"checksum",mxChecksum);
  }
  {
    const char *dataFields[] = {"size","type","complexity"};
    mxArray *mxData = mxCreateStructMatrix(1,1,3,dataFields);
    {
      mxArray *mxSize = mxCreateDoubleMatrix(1,2,mxREAL);
      double *pr = mxGetPr(mxSize);
      pr[0] = (double)(52);
      pr[1] = (double)(2);
      mxSetField(mxData,0,"size",mxSize);
    }
    {
      const char *typeFields[] = {"base","aliasId","fixpt"};
      mxArray *mxType = mxCreateStructMatrix(1,1,3,typeFields);
      mxSetField(mxType,0,"base",mxCreateDoubleScalar(10));
      mxSetField(mxType,0,"aliasId",mxCreateDoubleScalar(0));
      mxSetField(mxType,0,"fixpt",mxCreateDoubleMatrix(0,0,mxREAL));
      mxSetField(mxData,0,"type",mxType);
    }
    mxSetField(mxData,0,"complexity",mxCreateDoubleScalar(1));
    mxSetField(mxAutoinheritanceInfo,0,"inputs",mxData);
  }
  {
    mxSetField(mxAutoinheritanceInfo,0,"parameters",mxCreateDoubleMatrix(0,0,mxREAL));
  }
  {
    const char *dataFields[] = {"size","type","complexity"};
    mxArray *mxData = mxCreateStructMatrix(1,1,3,dataFields);
    {
      mxArray *mxSize = mxCreateDoubleMatrix(1,2,mxREAL);
      double *pr = mxGetPr(mxSize);
      pr[0] = (double)(52);
      pr[1] = (double)(1);
      mxSetField(mxData,0,"size",mxSize);
    }
    {
      const char *typeFields[] = {"base","aliasId","fixpt"};
      mxArray *mxType = mxCreateStructMatrix(1,1,3,typeFields);
      mxSetField(mxType,0,"base",mxCreateDoubleScalar(10));
      mxSetField(mxType,0,"aliasId",mxCreateDoubleScalar(0));
      mxSetField(mxType,0,"fixpt",mxCreateDoubleMatrix(0,0,mxREAL));
      mxSetField(mxData,0,"type",mxType);
    }
    mxSetField(mxData,0,"complexity",mxCreateDoubleScalar(1));
    mxSetField(mxAutoinheritanceInfo,0,"outputs",mxData);
  }
  return(mxAutoinheritanceInfo);
}

static void chart_debug_initialization(SimStruct *S, unsigned int
 fullDebuggerInitialization)
{
  if(ssIsFirstInitCond(S) && fullDebuggerInitialization==1) {
    /* do this only if simulation is starting */
    if(!sim_mode_is_rtw_gen(S)) {
      {
        unsigned int chartAlreadyPresent;
        chartAlreadyPresent = sf_debug_initialize_chart(_std5MachineNumber_,
         4,
         1,
         1,
         2,
         0,
         0,
         0,
         0,
         &(chartInstance.chartNumber),
         &(chartInstance.instanceNumber),
         ssGetPath(S),
         (void *)S);
        if(chartAlreadyPresent==0) {
          /* this is the first instance */
          sf_debug_set_chart_disable_implicit_casting(_std5MachineNumber_,chartInstance.chartNumber,1);
          sf_debug_set_chart_event_thresholds(_std5MachineNumber_,
           chartInstance.chartNumber,
           0,
           0,
           0);

          {
            unsigned int dimVector[1];
            dimVector[0]= 52;
            _SFD_SET_DATA_PROPS(1,2,0,1,SF_DOUBLE,1,&(dimVector[0]),0,0,0,0.0,1.0,0,"channel_estimate",1,c4_b_sf_marshall);
          }
          {
            unsigned int dimVector[2];
            dimVector[0]= 52;
            dimVector[1]= 2;
            _SFD_SET_DATA_PROPS(0,1,1,0,SF_DOUBLE,2,&(dimVector[0]),0,0,0,0.0,1.0,0,"freq_tr_syms",1,c4_sf_marshall);
          }
          _SFD_STATE_INFO(0,0,2);
          _SFD_CH_SUBSTATE_COUNT(0);
          _SFD_CH_SUBSTATE_DECOMP(0);
        }
        _SFD_CV_INIT_CHART(0,0,0,0);
        {
          _SFD_CV_INIT_STATE(0,0,0,0,0,0,NULL,NULL);
        }

        _SFD_CV_INIT_TRANS(0,0,NULL,NULL,0,NULL);

        /* Initialization of EML Model Coverage */
        _SFD_CV_INIT_EML(0,1,0,0,0,0,0,0);
        _SFD_CV_INIT_EML_FCN(0,0,"eML_blk_kernel",2,-1,376);
        _SFD_TRANS_COV_WTS(0,0,0,1,0);
        if(chartAlreadyPresent==0)
        {
          _SFD_TRANS_COV_MAPS(0,
           0,NULL,NULL,
           0,NULL,NULL,
           1,NULL,NULL,
           0,NULL,NULL);
        }
        _SFD_SET_DATA_VALUE_PTR(1U, c4_channel_estimate());
        _SFD_SET_DATA_VALUE_PTR(0U, c4_freq_tr_syms());
      }
    }
  } else {
    sf_debug_reset_current_state_configuration(_std5MachineNumber_,chartInstance.chartNumber,chartInstance.instanceNumber);
  }
}

static void sf_opaque_initialize_c4_std5(void *chartInstanceVar)
{
  chart_debug_initialization(chartInstance.S,0);
  initialize_params_c4_std5();
  initialize_c4_std5();
}

static void sf_opaque_enable_c4_std5(void *chartInstanceVar)
{
  enable_c4_std5();
}

static void sf_opaque_disable_c4_std5(void *chartInstanceVar)
{
  disable_c4_std5();
}

static void sf_opaque_gateway_c4_std5(void *chartInstanceVar)
{
  sf_c4_std5();
}

static void sf_opaque_terminate_c4_std5(void *chartInstanceVar)
{
  if (sim_mode_is_rtw_gen(chartInstance.S)) {
    sf_clear_rtw_identifier(chartInstance.S);
  }
  finalize_c4_std5();
}

static void mdlProcessParameters_c4_std5(SimStruct *S)
{
  int i;
  for(i=0;i<ssGetNumRunTimeParams(S);i++) {
    if(ssGetSFcnParamTunable(S,i)) {
      ssUpdateDlgParamAsRunTimeParam(S,i);
    }
  }
  initialize_params_c4_std5();
}

static void mdlSetWorkWidths_c4_std5(SimStruct *S)
{
  if(sim_mode_is_rtw_gen(S)) {
    int_T chartIsInlinable =
      (int_T)sf_is_chart_inlinable("std5",4);
    ssSetStateflowIsInlinable(S,chartIsInlinable);
    ssSetEnableFcnIsTrivial(S,1);
    ssSetDisableFcnIsTrivial(S,1);
    ssSetNotMultipleInlinable(S,sf_rtw_info_uint_prop("std5",4,"gatewayCannotBeInlinedMultipleTimes"));
    if(chartIsInlinable) {
      ssSetInputPortOptimOpts(S, 0, SS_REUSABLE_AND_LOCAL);
      sf_mark_chart_expressionable_inputs(S,"std5",4,1);
      sf_mark_chart_reusable_outputs(S,"std5",4,1);
    }
    if (!sf_is_chart_instance_optimized_out("std5",4)) {
      int dtId;
      char *chartInstanceTypedefName =
        sf_chart_instance_typedef_name("std5",4);
      dtId = ssRegisterDataType(S, chartInstanceTypedefName);
      if (dtId == INVALID_DTYPE_ID ) return;
      /* Register the size of the udt */
      if (!ssSetDataTypeSize(S, dtId, 8)) return;
      if(!ssSetNumDWork(S,1)) return;
      ssSetDWorkDataType(S, 0, dtId);
      ssSetDWorkWidth(S, 0, 1);
      ssSetDWorkName(S, 0, "ChartInstance"); /*optional name, less than 16 chars*/
      sf_set_rtw_identifier(S);
    }
    ssSetHasSubFunctions(S,!(chartIsInlinable));
    ssSetOptions(S,ssGetOptions(S)|SS_OPTION_WORKS_WITH_CODE_REUSE);
  }

  ssSetChecksum0(S,(3006270258U));
  ssSetChecksum1(S,(2175242051U));
  ssSetChecksum2(S,(2685313568U));
  ssSetChecksum3(S,(322155487U));

  ssSetExplicitFCSSCtrl(S,1);
}

static void mdlRTW_c4_std5(SimStruct *S)
{
  if(sim_mode_is_rtw_gen(S)) {
    sf_write_symbol_mapping(S, "std5", 4);
    ssWriteRTWStrParam(S, "StateflowChartType", "Embedded MATLAB");
  }
}

static void mdlStart_c4_std5(SimStruct *S)
{
  chartInstance.chartInfo.chartInstance = NULL;
  chartInstance.chartInfo.isEMLChart = 1;
  chartInstance.chartInfo.chartInitialized = 0;
  chartInstance.chartInfo.sFunctionGateway = sf_opaque_gateway_c4_std5;
  chartInstance.chartInfo.initializeChart = sf_opaque_initialize_c4_std5;
  chartInstance.chartInfo.terminateChart = sf_opaque_terminate_c4_std5;
  chartInstance.chartInfo.enableChart = sf_opaque_enable_c4_std5;
  chartInstance.chartInfo.disableChart = sf_opaque_disable_c4_std5;
  chartInstance.chartInfo.mdlRTW = mdlRTW_c4_std5;
  chartInstance.chartInfo.mdlStart = mdlStart_c4_std5;
  chartInstance.chartInfo.mdlSetWorkWidths = mdlSetWorkWidths_c4_std5;
  chartInstance.chartInfo.restoreLastMajorStepConfiguration = NULL;
  chartInstance.chartInfo.restoreBeforeLastMajorStepConfiguration = NULL;
  chartInstance.chartInfo.storeCurrentConfiguration = NULL;
  chartInstance.S = S;
  ssSetUserData(S,(void *)(&(chartInstance.chartInfo))); /* register the chart instance with simstruct */

  if(!sim_mode_is_rtw_gen(S)) {
    init_dsm_address_info();
  }
  chart_debug_initialization(S,1);
}

void c4_std5_method_dispatcher(SimStruct *S, int_T method, void *data)
{
  switch (method) {
   case SS_CALL_MDL_START:
    mdlStart_c4_std5(S);
    break;
   case SS_CALL_MDL_SET_WORK_WIDTHS:
    mdlSetWorkWidths_c4_std5(S);
    break;
   case SS_CALL_MDL_PROCESS_PARAMETERS:
    mdlProcessParameters_c4_std5(S);
    break;
   default:
    /* Unhandled method */
    sf_mex_error_message("Stateflow Internal Error:\n"
     "Error calling c4_std5_method_dispatcher.\n"
     "Can't handle method %d.\n", method);
    break;
  }
}