Introduction To Fading Channel: Ece 4202 Communication Engineering Lab Iv
Introduction To Fading Channel: Ece 4202 Communication Engineering Lab Iv
Week Topics
1 Introduction to Computer Simulation and MATLAB
2 Monte Carlo Simulations of Communication Systems
3 –Simulation
AWGN Channel
of MPSK
4 Simulation of MQAM
5 Simulation of Multipath Fading Channel
6 Performance Evaluation of Physical /Link Layer
7-11 Mini Project
12-13 Evaluation and Assessment
14 Lab Test
ECE4202 Communication Engineering Lab IV
Agenda
5
Large-scale propagation
Behavior of radio waves over ‘long’ distance (several
tens to hundreds of meters)
Small-scale propagation
Behavior of radio waves over ‘short’ distance (5 to 40
wavelengths)
THIS LAB IS ABOUT SMALL-SCALE FADING
sc
at
te r
in
g
SMALL SCALE
Non-resolvable;
Seen as a single path by receiver
>>sumsig=wave1+wave2;
>>subplot(311)
>>plot(10*log10(abs(wave1)))
>>subplot(312)
>> plot(10*log10(abs(wave2)))
>>subplot(313)
>>plot(10*log10(abs(sumsig)))
Rayleigh fading
No Line of Sight (LoS) component presents
Rician fading
Line of Sight (LoS) component presents
Deterministic component, contributes significantly to the total
received power
-5
-15
-20
-25
Deep fading
(weak signal)
-30
-35
0 200 400 600 800 1000
#sample
ECE 4202 COMM ENG LAB IV
Fading Envelope
13
-2
10
BER
-3
10
-4
10
-5
10
-6
10
0 5 10 15 20 25 30 35 40
Eb/N0
ECE 4202 COMM ENG LAB IV
What happen to the tx signal?
17
Scatter plot S c atter plot
2 2
1.5 1.5
1 1
0.5 0.5
Quadrature
Quadrature
0 0
-1
n(t) -1
-1.5
TXSIG -1.5 AWGN
-2 -2
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2
h(t)
FADING
1.5 1.5
1 1
0.5 0.5
Quadrature
Quadrature
AWGN
0 0
n(t)
-0.5 -0.5
-1 -1
FADING
ECE 4202 COMM ENG LAB IV
-1.5 -1.5
FADING+AWGN
-2 -2
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2
In-Phase In-Phase
Cont’d
18
hMod = modem.pskmod(M);
hDemod = modem.pskdemod(hMod);
txsig=modulate(hMod,data);
scatterplot(txsig);
scatterplot(fadesig);
%add noise
rxsig=awgn(fadesig,15);
%demodulate
recsig=demodulate(hDemod,rxsig);
1.5
0.5
Quadrature
0
-0.5
-1
-1.5
-2
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
ECE 4202 COMM ENG LAB IV In-Phase
Cont’d
22
y (t ) = s (t ) * h(t )
Note that, when h(t)=d(t), the convolution reduces to
multiplication
Frequency domain representation (after Fourier
transform)
Y ( f ) = S( f ) • H ( f )
Y( f )
H( f ) =
S( f )
ECE 4202 COMM ENG LAB IV
Channel Estimation
23
N p −1
~ 1
h (t ) =
Np
∑ y ( n) • p ( n)
n =0
~
h (t ) : Channel estimate
Np : number of pilot symbols
y(n) : received signal samples
p(n) : n-th pilot symbol
%add noise
rxsig=awgn(fadesig,15);
scatterplot(rxsig);
%demodulate
recsig=pskdemod(modemin,2);
%calculate the number of errors
sum(recsig~=data)