ms_sfun_interleaver.m

利用matlab/simulink仿真了CDMA2

function[sys,x0,str,ts]=MS_SFUN_Interleaver(t,x,u,flag,xInerleaverSize)
switch flag
    case 0
        [sys,x0,str,ts]=mdlInitializeSize(xInerleaverSize);
    case 2
        sys=mdlUpdate(t,x,u,xInerleaverSize);
    case 3
        sys=mdlOutputs(t,x,u);
    case{1,4,9}
        sys=[];
    otherwise
        error(['Unhandled flag=',num2str(flag)']);
end


function [sys,x0,str,ts]=mdlInitializeSize(xInerleaverSize)
sizes=simsizes;
%设置连续状态的个数
sizes.NumContStates=0;
%设置离散状态的个数
sizes.NumDiscStates=2;
%设置输出信号的个数为动态确定
sizes.NumOutputs=-1;
%设置输入信号的个数为动态确定
sizes.NumInputs=-1;
%设置直接反馈的状态
sizes.DirFeedthrough=1;
%设置抽样时间的个数（大于等于1）
sizes.NumSampleTimes=1;
sys=simsizes(sizes);
%设置s函数的初始状态x0
[m J]=mdlInterleaverInit(xInerleaverSize);
x0=[m J]';
str=[];
ts=[-1 0];


function sys=mdlUpdate(t,x,u,xInerleaverSize)
[m J]=mdlInterleaverInit(xInerleaverSize);
sys=[m J]';


function sys=mdlOutputs(t,x,u)
m=x(1);
J=x(2);
exp_m=2^m;
A=zeros(exp_m*J,1);
for i=0:length(u)-1
    quotient=fix(i/J);
    BRO=0;
    for k=1:m
        if bitget(quotient,k)==1
            BRO=bitset(BRO,m-k+1);
        end
    end
    A(exp_m*mod(i,J)+BRO+1)=u(i+1);
end
sys=A;


function [m,J]=mdlInterleaverInit(xInerleaverSize)
switch xInerleaverSize
    case 384
        m=6;
        J=6;
    case 768
        m=6;
        J=12;
    case 1536
        m=6;
        J=24;
    case 3072
        m=6;
        J=48;
    case 6144
        m=7;
        J=48;
    case 12288
        m=7;
        J=96;
    case 576
        m=5;
        J=18;
    case 2304
        m=6;
        J=36;
    case 4608
        m=7;
        J=36;
    case 9216
        m=7;
        J=72;
    case 18432
        m=8;
        J=72;
    case 36864
        m=8;
        J=144;
    otherwise
        error('Error:Invalid Interleaver Size in MS_SFUN_Interleaver');
end