mqam_mrc_twotr.m

MIMO系统

%MQAM_MRC---Sfunction
function [sys,x0,str,ts]=MQAM_MRC_twoTR(t,x,u,flag,BitRate)
j=sqrt(-1);
M=16;
d=1;
% 对应的星座图
mapping=[-3*d+j*3*d;
         -d+j*3*d;
         d+j*3*d;
         3*d+j*3*d;
         -3*d+j*d;
         -d+j*d;
         d+j*d;
         3*d+j*d;
         -3*d-j*d;
         -d-j*d;
         d-j*d;
         3*d-j*d;
         -3*d-j*3*d;
         -d-j*3*d;
         d-j*3*d;
         3*d-j*3*d];

switch flag,
    case 0
        [sys,x0,str,ts]=mdlInitializeSizes(BitRate);%初始化模块
    case 3
        sys=mdlOutputs(t,x,u,BitRate,M,mapping);%计算输出
    case {1,2,4,9}
        sys=[];%无用的标志
    otherwise
        error(['Unhandled flag =',num2str(flag)]);%错误处理
end
%模块结束
%%%%%%%初始化状态、采样周期、Size结构%%%%%%%%%
function [sys,x0,str,ts]=mdlInitializeSizes(BitRate)
%调用simsizes创建sizes结构
sizes=simsizes;
%填写sizes结构
sizes.NumContStates=0;
sizes.NumDiscStates=0;
sizes.NumOutputs=2*BitRate;
sizes.NumInputs=3*BitRate;
sizes.DirFeedthrough=1;
sizes.NumSampleTimes=1;
sys=simsizes(sizes);
x0=[];%没有连续状态
str=[];%没有状态次序
ts=[0 0];%继承性的采样周期

%%%%%模块输出%%%%%%%%%%
function sys=mdlOutputs(t,x,u,BitRate,M,mapping)

%形成信道
delta=u(1:BitRate);%[(BitRate)*1]

%还原输出矩阵信号
r1=u(BitRate+1:2*BitRate)+j*u(2*BitRate+1:3*BitRate);%[(BitRate)*1]

%%%%%%%%%%%%%MRC rule%%%%%%%%%%%%%%
for i=1:M
    metric(i,:)=((abs(r1-delta*mapping(i))).^2)';%every branch's distance[1*BitRate]
end
[min_metrics num]=min(metric);%%select the smallest

decei_symbol=mapping(num);%(BitRate*1)
real_decei_symbol=real(decei_symbol);%实部
imag_decei_symbol=imag(decei_symbol);%虚部

sys=[real_decei_symbol;imag_decei_symbol];%((2*BitRate)*1)