%?BER?for?QPSK?modulation?in?a?Rayleigh?fading?channel?using?multiple?taps?Taps
%with?and?without?transmit?beamforming

clear
N?=?10^5%?number?of?bits?

%?Transmitter

ip?=?（2*（rand（1N）>0.5）-1）?+?j*（2*（rand（1N）>0.5）-1）;?%symbols?or?constellation?points

s?=?（1/sqrt（2））*ip;?%?normalization?of?energy?to?1

Th=input（‘Enter?the?number?of?Channel?Taps?:‘）;%take?input?from?user?1234....

?
Eb_N0_dB?=?[0:35];?%?multiple?Eb/N0?values

for?ii?=?1:length（Eb_N0_dB）
???
???n?=?1/sqrt（2）*[randn（1N）?+?j*randn（1N）];?%?white?gaussian?noise?

???h?=1/sqrt（2）*[randn（ThN）?+?j*randn（ThN）];?%?Rayleigh?fading?channel
???

???signal?=?（1/sqrt（Th））*kron（ones（Th1）s）;?
???
???%?Channel?and?noise?Noise?addition
???phase=?exp（j*angle（h））;%calculate?the?phase?of?channel
???phase_conjugate=conj（phase）;%conjugate?of?phase?
???bf_ch=?h.*phase_conjugate;
????y1?=?sum（h.*signal1）?+?10^（-Eb_N0_dB（ii）/20）*n;?%signal?received?without?Beamforming
????y2?=?sum（bf_ch.*signal1）?+?10^（-Eb_N0_dB（ii）/20）*n;?%receivd?with?Beamforming
???


???%?Receiver?side-Equalization
??
???y1Hat?=?y1./sum（h1）;?
???y2Hat?=?y2./sum（bf_ch1）;?

??y_re1?=?real（y1Hat）;?%?real
y_im1?=?imag（y1Hat）;?%?imaginary
ipHat1（find（y_re1?ipHat1（find（y_re1?>