8-Channel Models Estimation
8-Channel Models Estimation
The impulse
response changes
with time
1. Collect Data and partition in blocks of length
:
x
y
1/
MAX
s D
N T F <<
max D
s
F
F
N <<
X=reshape(x,N,length(x)/N);
Y=reshape(x,N,length(y)/N);
[ ]
N
X,Y =
0 = n N n = N n 2 =
N n =
B
N
N N n
B
=
Within each block the channel is
almost time invariant
2. Estimate impulse response in each block :
h(:,i)=xcorr(Y(:,i),X(:,i))/N;
h =
[
B
N
]
2 1 N
h =
[
]
B
N
2 1 N
Each row is an impulse response
taken at different times
t
s
F N / /
s
N F +
Take the transpose:
plot((-N+1:N-1)/Fs, abs(h(:,i)));
B
N
2 1 N
[
B
N
]
2 1 N
time
time
s
T t =
( )
s
t n NT =
S =
time
s
T t =
Freq.
s
B
F
F k
N N
=
4. Take the sum over rows for Doppler Spread and sum over columns for
Time Spread (fftshift each vector to have zero term (sec or Hz) in the
middle
s
F m t / =
B
s
N
N F
k f
) / (
=
St
Sf
Frequency Resolution:
Hz
N N
F
F
B
S
) length(sec data total
1
=
= A
Therefore if we want to a resolution in the doppler spread of (say) 1Hz,
we need to collect at least 1 sec of data.
Time Resolution:
S
F t / 1 = A
Example: % channel
Fs=10^6; % sampling freq. In Hz
P=[0,-2,-3]; % attenuations in dB
T=[0, 10, 15]*10^(-6); % time delays in sec
fd=70; %doppler shift in Hz
y
To Workspace1
x
To Workspace
Rectangul ar
QAM
Rectangul ar QAM
Modul ator
Baseband
Rayl ei gh
Fadi ng
Mul ti path Rayl ei gh
Fadi ng Channel
Bernoul l i
Bi nary
Bernoul l i Bi nary
Generator
test_scattering.mdl
Channel Output (Magnitude) with a QPSK Transmitted Signal:
t (sec)
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
0
0.5
1
1.5
2
2.5
3
-1.5 -1 -0.5 0 0.5 1
x 10
-4
0
1
2
3
4
5
6
7
8
9
x 10
-3
Time Spread
time (sec)
-1000 -800 -600 -400 -200 0 200 400 600 800 1000
0
0.2
0.4
0.6
0.8
1
1.2
x 10
-3
Frequency Spread
frequency (Hz)
Time Spread
Hz 70 + Hz 70
Frequency Spread
sec 15
sum(S)/NB; sum of each column
sum(S)/(2N-1); ave. of each row