stbc.m

空时码中的STBC技术编码方案

function  BER=stbc(SNR_dB)
%------------------------------------------------------------------------
%本程序是对两发一收情况下采用空时分组码的性能分析
%星座映射采用的是BPSK映射
% SNR_dB 接收天线上的接收信噪比，在仿真的过程中把接收信号的功率归一化为1
% BER是误比特率
%-------------------------------------------------------------------------
SNR_dB=[0:5]
SNR=10.^(SNR_dB/10);
Cons=sqrt(1/2)*[1 -1]; % BPSK星座
E=1; %接收端的信号功率
No=E./(SNR);%高斯白噪声的功率谱密度
Frame=10000;
ErrorNum=zeros(1,6);
BER=zeros(1,6);
for i1=1:length(SNR_dB)   
    tic
    fprintf('@@@@@@@@@@@@@@@@@@@@@@@@@Now the SNR is %fdb.@@@@@@@@@@@@@@@@@@@@@@@@@@\n',SNR_dB(i1));  
for i2=1:1:Frame
    x1=Cons((rand(1)>=1/2)+1); % 随机产生发送数据
    x2=Cons((rand(1)>=1/2)+1);
    X=[x1 x2 ;-conj(x2) conj(x1)]; % 进行空时编码
    h1=sqrt(1/2)*(randn(1)+j*randn(1));  % 产生信道 两个发送天线到接收天线的信道
    h2=sqrt(1/2)*(randn(1)+j*randn(1));
    H=[h1;h2];
    N=sqrt(No(i1)/2)*(randn(2,1)+j*randn(2,1)); %接收端的噪声
    R=X*H+N;
    %-----------------------------------------------进行空时解码
    A=h1*conj(h1)+h2*conj(h2);
    X1=(conj(h1)*R(1)+h2*conj(R(2)))/A;
    X2=(conj(h2)*R(1)-h1*conj(R(2)))/A;
    if real(X1)>=0
        X1=Cons(1);
    else X1=Cons(2);
    end
    if real(X2)>=0
        X2=Cons(1);
    else X2=Cons(2);
    end
    ErrorNum(i1)=ErrorNum(i1)+(x1~=X1)+(x2~=X2);
end
BER(i1)=ErrorNum(i1)/(Frame*2);
save BER BER;
save ErrorNum ErrorNum;
end
toc
figure;
plot(SNR_dB,BER,'b-');
end