ms_sfun_turbopuncture.m

用matlab的Simulink对cdma2000中信道编码进行仿真 主要是turbo码部分。可能要用matlab7。0打开。

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


function [sys,x0,str,ts]=mdlInitializeSizes(interleaver_size,code_rate)
sizes=simsizes;
sizes.NumContStates=0;
len=(interleaver_size+6)*6;
switch code_rate
    case 1
        out_len=len/3;
    case 2
        out_len=len/2;
    case 3
        out_len=len*2/3;
    otherwise
        error('Error:Invalid Code Rate in MS_SFUN_TurboPuncture');
end
sizes.NumOutputs=out_len;
sizes.NumInputs=len;
sizes.DirFeedthrough=1;
sizes.NumSampleTimes=1;
sys=simsizes(sizes);
x0=[];
str=[];
ts=[-1 0];

function sys=mdlOutputs(t,x,u,interleaver_size,code_rate)
len=(interleaver_size+6)*6;
index1=1;
index2=len/2+1;
A=zeros(1,len);
for i=1:len
    if mod(i-1,6)<3
        A(i)=u(index1);
        index1=index1+1;
    else
        A(i)=u(index2);
        index2=index2+1;
    end
end

switch code_rate
    case 1
    pattern =[1 1 0 0 0 0 1 0 0 0 1 0];
    pattern1=[1 1 0 0 0 0];
    pattern2=[0 0 0 1 1 0];
    out_len=len/3;
    case 2
        pattern=[1 1 0 0 1 0 1 1 0 0 1 0];
        pattern1=[2 1 0 0 0 0];
        pattern2=[0 0 0 2 1 0];
        out_len=len/2;
    case 3
        pattern=[1 1 1 0 0 1 1 1 0 0 1 1];
        pattern1=[2 1 1 0 0 0];
        pattern2=[0 0 0 2 1 1];
        out_len=len*2/3;
    otherwise
        error('Error:Invalid Code Rate in MS_SFUN_TurboPuncture');
end

B=zeros(1,out_len);
index=1;
for i=1:(len-36)
    if pattern(mod(i-1,12)+1)==1
        B(index)=A(i);
        index=index+1;
    end
end

for i=(len-36)+1:len-18
    for k=1:pattern1(mod(i-1,6)+1)
        B(index)=A(i);
        index=index+1;
    end
end

for i=(len-18)+1:len
    for k=1:pattern2(mod(i-1,6)+1)
        B(index)=A(i);
        index=index+1;
    end
end

sys=B;