1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
|
/* file: CreateUmtsInterleaver.c
Description: Produce an interleaver according to the UMTS spec.
The calling syntax is:
[alpha] = CreateUmtsInterleaver( K )
alpha = the interleaver in a length K vector
K = the size of the interleaver
Copyright (C) 2005-2006, Matthew C. Valenti
Last updated on June 10, 2006
Function CreateUmtsInterleaver is part of the Iterative Solutions
Coded Modulation Library. The Iterative Solutions Coded Modulation
Library is free software; you can redistribute it and/or modify it
under the terms of the GNU Lesser General Public License as published
by the Free Software Foundation; either version 2.1 of the License,
or (at your option) any later version.
This library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with this library; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
*/
#include <math.h>
#include <mex.h>
#include <Matrix.h>
#include <stdlib.h>
/* library of functions */
#include "./include/interleaver.h"
/* Input Arguments */
#define INPUT prhs[0]
/* Output Arguments */
#define OUTPUT plhs[0]
/* main function that interfaces with MATLAB */
void mexFunction(
int nlhs,
mxArray *plhs[],
int nrhs,
const mxArray *prhs[] )
{
int DataLength;
double *output_p;
int *alpha_code; /* interleaver */
int *interleaver_input; /* Temporary array used to initialize the interleaver. */
int i;
/* Check for proper number of arguments */
if (nrhs < 1) {
mexErrMsgTxt("[alpha] = CreateUmtsInterleaver( K )");
} else if (nlhs > 2) {
mexErrMsgTxt("[alpha] = CreateUmtsInterleaver( K )");
}
/* initialize the input data */
DataLength = (int) *mxGetPr(INPUT);
if ( (DataLength < 40)|( DataLength > 5114) )
mexErrMsgTxt("CreateUmtsInterleaver: Input must be between 40 and 5114");
/* Create the interleaver */
alpha_code = (int*)calloc( DataLength, sizeof(int) );
interleaver_input = (int*)calloc( DataLength, sizeof(int) );
for (i=0;i<DataLength;i++)
interleaver_input[i] = i;
CreateUmtsInterleaver( DataLength, interleaver_input, alpha_code );
/* Output encoded data */
OUTPUT = mxCreateDoubleMatrix(1, DataLength, mxREAL);
output_p = mxGetPr(OUTPUT);
for (i=0;i<DataLength;i++) {
output_p[i] = alpha_code[i];
}
free( alpha_code );
free( interleaver_input );
return;
}
|
