DVB_T系统的matlab调制解调仿真代码

/*%**************************************************************************** %* All rights reserved (C). 2008 AICIA (Spain) %* This file may not be used, copied, distributed, corrected, modified, %* translated, transmitted or assigned without AICIA's prior authorization. %**************************************************************************** %**************************************************************************** %* Project : DVBT2 %* Date : $Date: 2008-02-01 15:31:38 $ %* Version : $Revision: 1.1 $ %* Author(s) : SIDSA %* Modifier(s) : %* Description : Fast Convolution. MEX file %* %* This function outperformes the matlab conv function speed. %* %* INPUTS: %* prhs[0] : Complex Baseband Channel Impulse Response %* prhs[1] : Complex (IQ) Transmitted Data %* %* OUTPUTS: %* prls[0] : Complex (IQ) Received Data %***************************************************************************/ #include "mex.h" void mexFunction( int Nlhs , /* Number of left-hand and side args */ mxArray *Plhs[] , /* Pointer to left-hand and side args */ int Nrhs , /* Number of right-hand and side args */ const mxArray *Prhs[] /* Pointer to right-hand and side args */ ) { int len_X ; /* Lenght of conv. input X */ int len_Y ; /* Lenght of conv. output Y */ int len_H ; /* Lenght of CIR output H */ int n ; /* Sample number */ double *X_re ; /* X input samples (real) */ double *X_im ; /* X input samples (imag.) */ double *H_re ; /* H CIR samples (real) */ double *H_im ; /* H CIR samples (imag.) */ double acc_re ; /* Convolution accumulator (real) */ double acc_im ; /* Convolution accumulator (imag.) */ int k ; /* Convolution iteration */ int lowlimit ; /* Limit to control first samples */ int uplimit ; /* Limit to control last samples */ double *y_re ; /* Convolution result (real) */ double *y_im ; /* Convolution result (imag.) */ double *k_out ; /* Convolution result (num samples) */ H_re = mxGetPr(Prhs[0]); H_im = mxGetPi(Prhs[0]); len_H = mxGetN(Prhs[0]); X_re = mxGetPr(Prhs[1]); X_im = mxGetPi(Prhs[1]); len_X = mxGetN(Prhs[1]); len_Y = len_X + len_H - 1; Plhs[0] = mxCreateDoubleMatrix(1,len_Y,mxCOMPLEX); y_re = mxGetPr(Plhs[0]); y_im = mxGetPi(Plhs[0]); /* Convolution loop */ for(k = 0;k < len_Y;k++) { acc_re = 0; acc_im = 0; /* To control first samples */ lowlimit = k-len_X+1; if (lowlimit < 0) lowlimit = 0; /* To control last samples */ uplimit = k; if (uplimit > len_H-1) uplimit = len_H-1; /* Complex convolution */ for (n=lowlimit;n <= uplimit;n++) { acc_re = acc_re + H_re[n]*X_re[k-n] - H_im[n]*X_im[k-n]; acc_im = acc_im + H_re[n]*X_im[k-n] + H_im[n]*X_re[k-n]; } y_re[k] = acc_re; y_im[k] = acc_im; } }