function  Pb=stbc1(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=10e5;
ErrorNum=0;
for i=1:1:Frame
    x1=Cons((rand(1)>=1/2)+1); % 隨機產生發送數據
%     x2=Cons((rand(1)>=1/2)+1);
    X=[x1]; % 進行空時編碼
    h1=sqrt(1/2)*(randn(1)+j*randn(1));  % 產生信道兩個發送天線到接收天線的信道
%     h2=sqrt(1/2)*(randn(1)+j*randn(1));
    H=[h1];
    N=sqrt(No/2)*(randn(1,1)+j*randn(1,1)); %接收端的噪聲
    R=X*H+N;
    %-----------------------------------------------進行空時解碼
    A=h1*conj(h1);
    X1=(conj(h1)*R(1))/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);
end
Pb=ErrorNum/(Frame);