bpsk_cooper_rayleigh_ber1.m

协作通信中

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
%%%%%BPSK瑞利信道下的仿真%%%%%%%%%%%
clear all
clc
omega=1;%%%%%瑞利信道的参数模值平方值
N=1000;%仿真点数
N2=1000;
SNRindB=10;%%%%%高斯噪声设置
SNR=10.^(SNRindB/10);
noise_var=1./10.^(SNRindB/10);
m1=length(noise_var);
bit_error=zeros(1,m1);
for k=1:m1
    for k2=1:N2
B=[1 -1];%BPSK的可能取值
A=randsrc(1,N,B);
%%%%%%%%%%%%%产生信道衰落H值%%%%%%%%%
hI=normrnd(0,sqrt(omega/2),3,N);
hQ=normrnd(0,sqrt(omega/2),3,N);
h0=hI(1,:)+j*hQ(1,:);
h1=hI(2,:)+j*hQ(2,:);
h2=hI(3,:)+j*hQ(3,:);
%alpha=abs(h);
%%%%%%%%%%%%%%%产生高斯噪声%%%%%%%%%%%%%%%%
NoiseI=normrnd(0,sqrt(noise_var(k)),3,N);
NoiseQ=normrnd(0,sqrt(noise_var(k)),3,N);
z0=NoiseI(1,:)+j*NoiseQ(1,:);
z1=NoiseI(2,:)+j*NoiseQ(2,:);
z2=NoiseI(3,:)+j*NoiseQ(3,:);
%%%%%%%%%%%%%%%%接收信号%%%%%%%%%%%%%%%%%%%%%
y0=h0.*A+z0;
y1=h1.*A+z1;
a=sqrt(1./(abs(h1).^2+noise_var(k)));
y2=h2.*a.*y1+z2;
y=conj(h0).*y0+conj(h1).*conj(h2).*a.*y2;
%y=conj(h0).*h0.*A+conj(h0).*z0;
%%%%%%%%%%%%%%判决变量%%%%%%%%%%%%%%%%%%%%%
%z=conj(h)*y;
%%%%%%%%%%%%%判决%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:N
    if y(i)>0
        ybit(i)=1;
    else
        ybit(i)=-1;
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
bit_error(k)=N-sum(ybit==A)+bit_error(k);
%the(k)=0.5*(1-sqrt(SNR(k)/(SNR(k)+1)));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    end
    bit_error(k)=bit_error(k)/N/N2;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%the
bit_error
%semilogy(SNRindB,the,SNRindB,bit_error,'r')