mqam_cornaka_twotrmain.m

MIMO系统

%%%%%%%%%%%%%%%%%%MIMO系统MQAM通过相关Nakagami衰落信道进行各种分集技术的仿真%%%%%%%%%%%%%%%%%%%%%
clear;clc
% hold on
tic;
%Rou_gamma_Tx: correlation coefficient matrix of transmitter
%Rou_gamma_Rx: correlation coefficient matrix of receiver
%m: fading parameter
%omegap:  mean fading power omegap=E(x^2)
%N: number of transmit branches
%L: number of receive branches
%R: symbol transmission rate
%sel_num: select sel_num out of L diversity branchs
%parameters:
snr=11:19;%%平均比特信噪比
m=2;
Rou_gamma_Tx=[1 0;0 1];
N=2;%确定
Rou_gamma_Rx=[1 0.795 0.605 0.375; 0.795 1 0.795 0.605;0.605 0.795 1 0.795;0.375 0.605 0.795 1];
L=4;
sel_num=2;

x=10.^(snr/10);
% 信噪比(db转化)
R=1;
epsilon=2;
d=1;
BitRate=100000;
% 比特率(一帧中采样点的数目)
omegap=1;
for i=1:length(snr)
%    SNR=snr(i);
%     sgma=sqrt(20*d^2*omegap/(2*x(i)));
sgma=sqrt(d^2*omegap*5*L*epsilon/(4*R*x(i)));
    sim('MQAM_CORNAKA_TwoTR');
%     仿真模块
    y1(i)=mean(BitErrorRate);
%     输出平均错误概率
end
semilogy(snr,y1,'*');
toc;