nodiversity.m

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

% 无分集，一发一收
function  BER=stbc(SNR_dB)
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(1)+j*randn(1)); %接收端的噪声
    R=x1*h1+N;
    %-----------------------------------------------进行空时解码
%     A=h1*conj(h1)+h2*conj(h2);
    X1=conj(h1)*R;
%     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);
end
BER(i1)=ErrorNum(i1)/Frame;
save BER BER;
save ErrorNum ErrorNum;
end
toc
figure;
plot(SNR_dB,BER,'b-');
end