stbc22.m

计算信道容量的matlab程序

function  Pb=stbc22(SNR_dB,a11,a12,a21,a22,a31,a32,a41,a42)
%function  Pb=stbc22(SNR_dB)
%------------------------------------------------------------------------
%本程序是对两发一收情况下采用空时分组码的性能分析
%星座映射采用的是BPSK映射
% SNR_dB 接收天线上的接收信噪比，在仿真的过程中把接收信号的功率归一化为1
% Pb是误比特率
%-------------------------------------------------------------------------
SNR=10^(SNR_dB/10);
Cons=sqrt(1/2)*[1 -1]; % BPSK星座
E=1; %接收端的信号功率
No=E/(SNR);%高斯白噪声的功率谱密度
Frame=10000;
ErrorNum=0;
for i=1:1:Frame
    x1=Cons((rand(1)>=1/2)+1); % 随机产生发送数据
    x2=Cons((rand(1)>=1/2)+1);
    h1=sqrt(1/2)*(a11+j*a12);  % 产生信道 两个发送天线到接收天线的信道
    h2=sqrt(1/2)*(a21+j*a22);
    h3=sqrt(1/2)*(a31+j*a32);
    h4=sqrt(1/2)*(a41+j*a42);
%    h1=sqrt(1/2)*(randn(1)+j*randn(1));  % 产生信道 两个发送天线到接收天线的信道
%    h2=sqrt(1/2)*(randn(1)+j*randn(1));
%    h3=sqrt(1/2)*(randn(1)+j*randn(1));
%    h4=sqrt(1/2)*(randn(1)+j*randn(1));
   
    n=sqrt(No/2)*(randn(1,1)+j*randn(1,1)); %接收端的噪声
    R1=h1*x1+h2*x2+n; % 第一根天线接收的信号，在连续两个符号周期中
    R2=-h1*conj(x2)+h2*conj(x1)+n;
    R3=h3*x1+h4*x2+n; % 第二根天线接收的信号，在连续两个符号周期中
    R4=-h3*conj(x2)+h4*conj(x1)+n;
    
    %-----------------------------------------------进行空时解码
    A=h1*conj(h1)+h2*conj(h2)+h3*conj(h3)+h4*conj(h4);
    X1=(conj(h1)*R1+h2*conj(R2)+conj(h3)*R3+h4*conj(R4))/A; %统计结果
    X2=(conj(h2)*R1-h1*conj(R2)+conj(h4)*R3-h3*conj(R4))/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=ErrorNum+(x1~=X1)+(x2~=X2); %误码数
end
Pb=ErrorNum/(Frame*2); %误码率