SISO.zip_siso_zip

clc; close all; clear all; N=1000; BPSK=randsrc(1,N); for Eb_No_dB=0:1:12 sig_info=sqrt(2*10.^(Eb_No_dB/10)).*BPSK; for m=1:1:1 %----channel impulse response--------- X=randn(1,N); Y=randn(1,N); unit_X=X/sqrt(var(X)); unit_Y=Y/sqrt(var(Y)); alpha=sqrt((unit_X.^2)+(unit_Y.^2)); alpha_norm=alpha/sqrt((mean(alpha.^2))); phase=rand(1,N); uni_phase=2*pi*phase; uni_phase_norm=uni_phase/sqrt(var(uni_phase)); e_ph=complex(cos(uni_phase_norm),-sin(uni_phase_norm)); h=alpha_norm.*e_ph; %----AWGN-------------------------- noise1=randn(1,N); noise2=randn(1,N); uni_var_noise1=noise1/sqrt(var(noise1)); uni_var_noise2=noise2/sqrt(var(noise2)); cmplex_noise=complex(uni_var_noise1,uni_var_noise2); %-----Recieved signal--------------- rxsig=h.*sig_info+cmplex_noise;%equalising for channel charaecteristics Rx_eq=rxsig.*conj(h); for i=1:1:N if(real(Rx_eq(i))>0) decision(i)=1; else decision(i)=-1; end end error=0; for i=1:1:N if decision(i)~=BPSK(i) error=error+1; end end ber(m)=error/N; end Pb_SISO(Eb_No_dB+1)=mean(ber); end Eb_No_dB=0:1:12; semilogy(Eb_No_dB,Pb_SISO,'r--'); grid on; hold on; axis([0 12 10^-4 0.5]); legend('SISO'); xlabel('SNR(dB)'); ylabel('BER'); title('SISO (MT=1,MR=1) in Rayleigh fading channel');