s_spreading.c

这是一个c＋＋编写的WCDMA链路采用RAKE接收的方针源代码

/*
 |
 | Copyright disclaimer:
 |   This software was developed at the National Institute of Standards
 |   and Technology by employees of the Federal Government in the course
 |   of their official duties. Pursuant to title 17 Section 105 of the
 |   United States Code this software is not subject to copyright
 |   protection and is in the public domain.
 |
 |   We would appreciate acknowledgement if the software is used.
 |
*/

/*
 | Project:     WCDMA simulation environment
 | Module:      
 | Author:      Maarit Melvasalo
 | Date:        February 23, 1999
 |
 | History:
 |             June 10, 1999 Maarit Melvasalo
 |                  Function rewritten
 *
 *  File    : s_spreading.c
 *  Abstract:
 *
 * USER GIVEN PARAMETERS
 *
 * 1: B  = the input block size  
 * 2: C  = the spreading code 
 * 3: SF = spreading factor
 * 4: nSlot = Number of slots in a frame (=used for sample timing)
 *
 * 
 * DEPENDS ON FILES 
 *            spreading.c 
 *            mealy.c  
 *            bitroutines.c
 *            conversions.c  
 *            utility.c      
 */

#define S_FUNCTION_NAME s_spreading
#define S_FUNCTION_LEVEL 2

#include <math.h>
#include "simstruc.h"
#include "tmwtypes.h"
#include "config_wcdma.h"
#include "wcdma_simulink.h"


/* USER GIVEN PARAMETERS AND DEFINITIONS /**/

/* General Parameters /**/
#define NPARAMS 4                                 /* Number of parameters /**/


/* Spreading parameters /**/
#define C_PARAM(S) ssGetSFcnParam(S,1)            /* Code */
#define nC   mxGetN(C_PARAM(S))                   /* Code length */ 

#define SF_PARAM(S) ssGetSFcnParam(S,2)            /* Block size*/
#define SF (int_T)(mxGetPr(SF_PARAM(S))[0])   /* Number of inputs*/

#define SLOTS(S)  ssGetSFcnParam(S,3)
#define nSlots (int_T)(mxGetPr(SLOTS(S))[0]) 

/* Input and Output Size parameters /**/
#define SIZE_P(S) ssGetSFcnParam(S,0)            /* Block size*/
#define nIn (int_T)(mxGetPr(SIZE_P(S))[0])   /* Number of inputs*/
#define nOut (int_T)(nIn * SF )        /* Number of outputs /**/
#define U(element) (*uPtrs[element])              /* Pointer to Input Port0 */

#define td  TD_FRAME / (real_T)nSlots  

/*====================*
 * S-function methods *
 *====================*/

/* Function: mdlInitializeSizes ===============================================
 * Abstract:
 *    The sizes information is used by Simulink to determine the S-function
 *    block's characteristics (number of inputs, outputs, states, etc.).
 */
static void mdlInitializeSizes(SimStruct *S)
{
  ssSetNumSFcnParams(S, NPARAMS);  /* Number of expected parameters */

  if (ssGetNumSFcnParams(S) != ssGetSFcnParamsCount(S)) 
 {
      return;
 }
  
  if (!ssSetNumInputPorts(S, 1)) return;  
    ssSetInputPortWidth(S, 0, nIn);
    ssSetInputPortDirectFeedThrough(S, 0, 1);

    if (!ssSetNumOutputPorts(S, 1)) return; 
    ssSetOutputPortWidth(S, 0,nOut);

    ssSetNumSampleTimes(S, 1);
    ssSetNumIWork(S, nC);

    ssSetSFcnParamNotTunable(S,0); 
    ssSetSFcnParamNotTunable(S,1);  
    ssSetSFcnParamNotTunable(S,2);  
    ssSetSFcnParamNotTunable(S,3);  

    /* Take care when specifying exception free code - see sfuntmpl.doc */
    ssSetOptions(S, SS_OPTION_EXCEPTION_FREE_CODE);
}   /* do not use mxCalloc or  mexErrMsgTxt functions */


/* Function: mdlInitializeSampleTimes =========================================
 * Abstract:
 *    Specifiy that sampletime is given by the user
 */
static void mdlInitializeSampleTimes(SimStruct *S)
{   
     ssSetSampleTime(S, 0, td); 
     ssSetOffsetTime(S, 0, 0.0);
}

#define MDL_INITIALIZE_CONDITIONS
/* Function: mdlInitializeConditions ========================================
 * Abstract:
 *           Stores the spreading code to IWork vector 
 */
static void mdlInitializeConditions(SimStruct *S)
{
    int_T             *iwork = ssGetIWork(S);
    real_T            *cpr  = mxGetPr(C_PARAM(S));
    int_T  i;

    for (i = 0; i < nC; i++){
      iwork[i] = (int_T)cpr[i];
    }
}

/* Function: mdlOutputs =======================================================
 * Abstract:
 * 
 */
static void mdlOutputs(SimStruct *S, int_T tid)
{
    real_T            *y    = ssGetOutputPortRealSignal(S,0);
    InputRealPtrsType uPtrs = ssGetInputPortRealSignalPtrs(S,0);
    int_T             antipodal[nOut]; /**/
    int_T             inputs[nIn]; /**/
    int_T             out[nOut]; /**/
    int_T             *iwork = ssGetIWork(S); /*spreading code/**/
    int_T             lp,tmp;
    
    for (lp=0; lp<nIn; lp++){
      inputs[lp]=(int_T)U(lp);
    }  

    bin2antipodal(inputs, nIn, antipodal);
    wcdma_spread(antipodal, nIn, iwork, nC, SF, out);
  

    for (lp=0; lp < nOut; lp++){ 
      y[lp] = (real_T)out[lp];/**/
    }/**/
}


/* Function: mdlTerminate =====================================================
 * Abstract:
 */
static void mdlTerminate(SimStruct *S)
{ 

}

#ifdef  MATLAB_MEX_FILE    /* Is this file being compiled as a MEX-file? */
#include "simulink.c"      /* MEX-file interface mechanism */
#else
#include "cg_sfun.h"       /* Code generation registration function */
#endif