ms_sfun_interleaver.m

CDMA2000 反向信道仿真 matlab

function [sys,x0,str,ts]=MS_SFUN_Interleaver(t,x,u,flag,interleaver_size)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% M文件S－函数主体部分
%函数名称: MS_SFUN_Interleaver
%主要功能: 根据输入参数flag的数值调用相应的函数
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
switch flag
    case 0
        %当flag等于0时调用mdlInitializeSize函数执行初始化
        [sys,x0,str,ts]=mdlInitializeSizes(interleaver_size);
    case 2
        %当flag等于2时调用mdlUpdate函数计算离散状态的数值
        sys=mdlUpdate(t,x,u,interleaver_size);
    case 3
        %当flag等于3时调用mdlOutputd函数计算输出信号的数值
        sys=mdlOutputs(t,x,u);
    case{1,4,9}
        %当flag等于1/4/9时没有相应的操作(没有连续状态)
        sys=[];
    otherwise
        %当flag等于其他数值时表示仿真过程出错
        error(['Unhandled flag=',num2str(flag)]);
end
%end MS_SFUN_Interleaver
%
%===========================================================
%mdlInitializeSizes
%S-函数的初始化
%向Simulink返回S-函数各种信号长度,初始设置和抽样时间设置
%===========================================================
function [sys,x0,str,ts]=mdlInitializeSizes(interleaver_size)
%调用simsize获得一个用于存放长度信息的结构
sizes=simsizes;
%设置连续状态的个数
sizes.NumContStates=0;
%设置离散状态的个数
sizes.NumDiscStates=2;
%设置输出信号的个数为动态确定
sizes.NumOutputs=-1;
%设置输入信号的个数为动态确定
sizes.NumInputs=-1;
%设置直接反馈的状态
%0表示没有直接反馈
%1表示存在直接反馈
sizes.DirFeedthrough=1;
%设置抽样时间的个数(>=1)
sizes.NumSampleTimes=1;
%通过simsizes把sizes结构返回给sys
sys=simsizes(sizes);
%
%设置S-函数的初始状态x0
[m J]=mdlInterleaverInit(interleaver_size);
x0=[m J]';
%
%设置S-函数的保留参数str(应该设置为空间向量[])
str=[];
%
%初始化抽样时间
ts=[-1 0];
%end mdlInitializeSizes
%
%==================================================================
%mdlUpdate
%更新S-函数的离散状态并且向Simulink返回这些状态的数值
%====================================================================
function sys=mdlUpdate(t,x,u,interleaver_size)
%计算S-函数的离散状态并且通过sys参数返回给Simulink
[m J]=mdlInterleaverInit(interleaver_size);
sys=[m J]';
%end mdlUpdate
%
%=================================================================
%mdlOutputs
%计算S-函数的输出信号并且返回给Simulink作为模块的输出
%===================================================================
function sys=mdlOutputs(t,x,u)
%计算S-函数的输出信号并且通过sys参数返回给Simulink
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;
%end mdlOutputs
%
%======================================================================
%mdlInterleaverInit
%初始化交织器长度
%=====================================================================
function [m,J]=mdlInterleaverInit(interleaver_size)
switch interleaver_size
    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