First, we introduce the format notations. In complex divider, the equation can be expressed by:
currence step based algorithm [21]. In recurrence step algorithm, it only requires an adder (substracter)
Motivation
-T/H systems for 8Mhz channels is about 109ns. Due to the long cycle time and the high hardware cost of long digits dividers, we propose a low cost architecture to implement the equivalent divider.
s and n is at structure is show
are (m1, n1 ll produce some
intermediate 1, 2n1).
α. When
the output da α. The complex
division includes two real ore, in order
to get n2 bits in decimal ac t n2 bits left. The dividend beco
igh.
Fig. 3.1 The format (m, n) structure
First, we define the format (m, n) of the signed number, where m is the total bit
the bits of decimal. The form n in Fig.3.1. The formats of inputs a, b, c, d, ), and the formats of outputs e, f, are (m2, n2). In this process, it wi
values like (ac + bd), (bc - ad), and (c2 + d2), which formats are (2m Since the output can only present in dynamic range P, we define saturation point
ta is out of the range P, it will be saturated at saturation point number divisions with 2m1 bits word-length. Furtherm
curacy, so the dividend should shif
mes (2m1 + n2) bits and the divisor is (2m1) bits. And we can find in Fig.3.2, we can find the output of single cycle division is 2m1+n2 bits which is bigger than m2 bits. These bits present the saturation cases and over design here. It needs (2m1+n2) bits subtractor. So if we implement the divider by single cycle division model directly, the hardware will be cost h
Fig. 3.2 The bits presentation of division
In proposed separate one cycle
into many cy lify the divider to
subtraction by hardware reus
it will detect if state 3, else goes to
state 2. In state 2, it will do m ration and getting the
quotient by iterations. At last, it ate 3.The state diagram is shown in Fig.3.3.
design, the multi-cycle division,
t
he basic concept is to cles and get quotient by iterative subtraction. So we can simpe and raising clock rate.
First, the procedure of proposed division scheme consists of three states [21]. In state 1, the result is saturated or not. If saturation occurs, it goes to
ain function including subtraction ope will output the result of division in st
detect saturation
output the result do main function
Next operation Not saturated
Saturation
Fig. 3.3 The state diagram
As mention before, because of the format (m , n ) of output so we can detect that if the
so it will not be affected. In this method, the relative position of dividend and divisor is
2 2
turated in the beginning. If the result is saturated, we can get the diately. In that case, other logics will be idle for power saving.
First, we define A is the minuend, B is the subtrahend, and sub is the result of subtraction.
ine the quotient by subtracting B from A. In this algorithm, the quotient can only or ‘0’, so if A is larger than two times of B, the quotient can’t represent.
reason, we must make sure that A can not be larger than two times of B in state 2.
t normalize to saturation poin
ting (m2-n2-1) bits right, and let B be the divisor. The quotient will shift back at last,
exactly to get the MSB of the quotient, and we can detect the saturation case easily. We only take one cycle for saturation detection, and we can save (2m1+n2-m2) cycles.
Here, A is the dividend >> (m2-n2-1), and B is the divisor. If A is larger than B, it will represent that the quotient is saturation, and the quotient will be the saturation point α. If A is smaller than B, it will make sure that A can not be larger than two times of B, and it will go to state 2 to get the quotient. The flow of state 1 is shown in Fig.4.
Fig. 3.4 The flow of the state 1
In state 2, we can use subtraction to determine the quotient q[k] is ‘1’
akes sure A can not be larger than two times of B e can determine it by the sign bit of the subtraction and q[k] is ‘1’, A < B, the sign bit of sub is ‘1’
sub. Then, we update A depends on the sign bit of sub to get the next bit of quotient. We can get
2-1) cycles in state 2. The flow
between parameters is listed in
or ‘0’ by binary
property. Because it m . When A ≥ B, the q[k]
will be ‘1’ else not, w result. If A ≥ B, the
sign bit of sub is ‘0’ and q[k] is ‘0’. We can
find q[k] is the inverse of the sign bit of by (sub<<1) or (A<<1) the quotient of format (m2, n2) one by one bit through (m of the main function in state 2 is shown in Fig.3.5, and the relationship table 3-1.
Fig. 3.5 The flow of the main function in state 2 Table 3-1 the parameter relationship in state2
Sub>=0 Sub<0
Sign bit of sub “0” “1”
q[k] “1” “0”
Remainder (A) Sub<<1 A<<1
In state 3, it is output stage. e last cycle of state 2, and goes to next new division operation.
Furthermore, we should determine the clock rate of this architecture which depends on
the cycles of one operati e is for state 1 to detect
if it is saturation. The following (m2-1) cycles are for compu (m2-1) bits of quotient excluding sign bit. The sign bit can be dete before state 2 will output the complete quotient stably in the last . So there are cles needed in operation, in other word the ratio between mbol : m2}. Th g diagram is shown in Fig.3
It can output the whole quotient at th
on. It will take m2 cycles totally. The first cycl ting the
rmine . It
cycle m2 cy one
clock rate and sy rate is {1 e timin .6. The hardware implementation will mention in section 5.3.2.
Fig. 3.6 Timing diagram
Chapter 4 .
imulation and Performance Analysis
ter, the overall simulation platform built for DVB-T/H system will be illustrated odel and some other distortion source such as Doppler delay spread and iscussed later. Finally, the performance analysis of the proposed channel ualizer scheme and comparison with state of the art will be performed.
In order to verify the performance of the proposed channel equalizer scheme, a complete VB-T/H baseband simulation platform is constructed in Matlab. The block diagram of the
S
In this chap
first. The channel m SCO model will be d eq
4.1 Simulation Platform
D
overall simulation platform is shown as Fig. 4.1.
Fig. 4.1 Overall DVB-T/H platform
As shown in Fig. 4.1, the blocks with dotted line is the specific function blocks for DVB-H
system. By adding support of 4k IFFT/FFT, in-depth interleaving, and additional TPS information, the developed DVB-T system platform can share most of the function blocks with DVB-H system at the same time. The platform is composed of transmitter, channel, and receiver. A typical transmitter that receives data from MPEG2 encoder or IP datagram is completely established. The transmitter consist the full function of FEC blocks and OFDM modulation blocks. The coding rate, interleaving mode, constellation mapping mode, IFFT length, and guard-interval length are all parameterized and able to be selected while simulation. An oversampling and pulse shaping filter is added before data entering channel to simulate discrete signal as far as continuously. The oversampling rate is also parameterized and can be chosen according to the simulation accuracy. The roll-off factor of the pulse-shaping filter is chosen as a normal value α =0.15 because it is not defined in the DVB-T/H standard.
Various distortion models are adopted in the channel model to simulate real mobile environment such as multipath fading, Doppler spread, AWGN, CFO, and SCO. In practically, there are still some imperfect effects which contain co-channel interference, adjacent-channel interference, phase noise, and common phase error caused by imperfect front-end receiving.
However, the distortion of these imperfect effects is relatively smaller compared with effective time-varying channel response caused by Doppler spread, CFO, and SCO. Therefore these effects are neglected in our simulation platform.
The baseband receiver in our system platform can be divided into inner receiver and outer receiver as Fig. 4.2 shows. The inner receiver includes all of the timing and frequency synchronization function, FFT demodulation, ch ation, equalization, and pilot remove blocks. The outer receiver consists of
annel estim
other functional blocks that following the de-mapping. The transmission parameters extracted by TPS decoder such as constellation mapping mode and Viterbi code rate will be sent the relative blocks as control parameters.
Besides, the extracted TPS parameter such as guard interval length and IFFT/FFT mode
should be checked all the time during online receiving to prevent synchronization error. Once TPS check fail occurs, the acquisition and tracking of inner receiver must be shut down and then restart all the synchronization schemes. As for bit-error-rate (BER) measurement, the DVB-T standard defines quasi error-free condition, which means less than one uncorrelated error event per hour, while the BER of the output of the Viterbi decoder is equal to 2 10× −4. Therefore, in order to verify the overall system performance, the BER after Viterbi decoder should be measured.
Fig. 4.2 The baseband receiver design
ζ^
^
εF
^
εI ε^R δ^
Fig. 4.3 Functional blocks of inner receiver
Fig. 4.3 shows the detail functional blocks of the inner receiver. The main functional blocks consists of symbol timing offset synchronization, carrier frequency offset synchronization, SCO synchronization, channel estimation, and equalizer, respectively. Th acquisition parts (gray color) only operate in the beginning of the receiving and then are turned off when the tracking parts work, and the tracking parts works all the time until th receiver is turned off or TPS check error occurs. In this thesis, we only focus on the performance analysis of the channel estimation and equalization scheme. The detailed discussion of other functional blocks such as timing synchronization and CFO synchronization will be neglected in this work and can be found in [22].
e
e
4.2 Channel Model
The typical baseband equivalent channel model for DVB-T/H system platform is shown as in Fig. 4.4. The transmitted data will pass through multipath fading, Doppler delay spread, CFO, SCO, and AWGN in turn. The effects of co-channel interference, adjacent-channel interference, phase noise, and common phase error are neglected in our simulation. In the following sections, the detailed effect of each channel distortion will be illustrated.
⊕
⊗
2 j ft
e πΔ (1+ζ)fs
.4 Channel model of DVB-T/H system Fig. 4
4.2.1 Multipath Fading Channel Model
In wireless communication transmission, the multipath fading is caused by the reception through different paths with different time delay and power decay. In DVB-T standard, two
types of multipath fading channel model are specified [3]. The fixed reception condition is modeled
by Rayleigh floating
the
following equations whe pectively
modeled by Ricean channel (Ricean factor = 10dB) while the portable reception is channel. The full 20-tap Ricean and Rayleigh channel was used with
point tap magnitude and phase values with tap delay accuracies rounded to within 1/2 of duration for practical discrete simulation. The channel models can be generated from
re x(t) and y(t) are input and output signals res
Rayleigh: path, and τi is the relatively delay of the i-th path, respectively. The detailed value of these parameters is listed in table B.1 of [3]. The rms delay of Rayleigh and Ricean channel is 1.4426 sμ (about 13 samples) and 0.4491 sμ (about 4 samples). From the above two equations, we can find that the major difference between Ricean and Rayleigh channel is the ma
(the ratio of the power of the direct path to the reflected path) and can be expressed as
in path (the sight way). In Ricean channel, a main path is defined with the Ricean factor K
2
ain path in Rayleigh ch
(4-3)
However, there is no m annel. Hence the received signals consist of several reflected signals with similar power d bring serious synchronization error. The impulse response and frequency response of the two types of channel when K=10dB are shown in Fig. 4.5. As we can see there is a significant direct path in the impulse response of the Ricean channel. In the impulse response of the Rayleigh channel, there is no any direct
an
path and all the paths have similar magnitude. Therefore, the frequency selective fading effect in the frequency response of the Rayleigh channel is more serious than that of the Ricean channel.
0 200 400 600 800 1000 1200 1400 1600 1800 0
(a) Impulse response of Rayleigh channel (b) Frequency response of Rayleigh channel
ud
200 400 600 800 1000 1200 1400 1600 1800 Subcarrier index
Amplitude
(c) Impulse response of Ricean channel (d) Frequency response of Ricean channel Fig. 4.5 Channel response of Rayleigh and Ricean (K=10dB) channel
4.2.2 Mobile channel model
In DVB-T standard, it only provides two static channel models which is described in section 4.2.1. However, the applications in DVB-T/H systems are not only fir fix reception, but also for mobile reception. Therefore, we refer to the channel models Typical Urban 6 (TU6) and Rural Area 6 (RA6) in GSM COST 207 project [23]. The two single-transmitter
profiles come from the set defined by the COST 207 project (GSM transmission). The technical specification of COST207 describes the equipment and techniques used to measure the channel characteristics over typical bandwidths of 10~20 MHz at near 900MHz. Therefore, the COST207 profiles are applicable to the DVB-T transmission situations. The detailed value of these parameters is listed in table 4-1 and table 4-2. The Fig4.6 shows the TU6 model, the Doppler spectrum filter will introduce in next section.
Fig. 4.6 TU6 model
Table 4-1 Typical Urban Reception (TU6) channel model
Tap number Delay(us) Power(dB) Doppler spectrum
1 0 -3 Rayleigh 2 0.2 0 Rayleigh 3 0.5 -2 Rayleigh 4 1.6 -6 Rayleigh 5 2.3 -8 Rayleigh 6 5.0 -10 Rayleigh
Table 4-2 Rural Area Reception (RA6) channel model
Tap number Delay(us) Power(dB) Doppler spectrum
1 0 0 Rice
oppler spec types
Doppler
In DVB-T/H system, the reception ability in mobile environment is necessary. Hence a mobile radio channel including Doppler spread must be constructed. A simplified Doppler spread model is shown in Fig. 4.7 [24]. In the beginning, we assume a channel with a known and fixed number of paths P such as Rayleigh or Ricean with a Doppler frequency d( )k
Fig. 4.7 Doppler spread model
f ,
attenuation ρ(k)ejθ( )k , and time delay τ( )k . All the parameters are fixed as described in section 4.2.1 except the Doppler frequency. Since each path has its own Doppler frequency, the decision of the statistic distribution of f is very important. There are two commonly d sed Doppler frequency PDFs, uniform and classical, where the former exploits uniform d, and the later uses Jake’s Doppler spectrum [25],
respectively e Jak s essed as
u
distribution to model Doppler sprea
. The PDF of th e’s Doppler pread can be expr
2
After transformation of random variable, each f can be obtained by the following equation d cos(2 (1)) max
d d
f = π⋅rand ⋅f (4-5) The type of Doppler spread (uniform or Jake’s) affects the system performance enormously.
Because each path gets different f in each simulation case with different d fdmax, the value of fdmax should be fixed for each simulation and comparison.
B. Rayleigh fading
In wireless communication, the multi path effect will cause the frequency selective fading problems. If the transmitting environment is without the main direct path, according to central limit theory the amplitude will be Rayleigh distribution and the phase will be uniform distribution. So here we based on the Jake’s model to simulate the Doppler effects, in concept,
in-phase) and imaginary components of each fader are zero-mean independent Gaussian the Jake’s model will produce the complex signal with the same statistic character s [25-27].
Simulating a multipath fading channel often requires the generation of multiple independent Rayleigh faders which can be odeled as complex-valued random processes.
Ideally the Rayleigh faders should conform to the following criteria. First, the real (or istic
m
processes with identical power spectra or auto-correlation functions. As a result, the envelope is Rayleigh distributed and the phase is uniformly distributed. Second, the cross-correlation
yleigh fading waveform
deterministically or statistically. The Rayleigh Doppler spectrum generator is shown in Fig.
4.8. The Rayleigh fading Doppler spectrum is generating by Jake’s model which Doppler between any pair of faders should be zero. Ra s can be generated
frequency is 100Hz is shown in Fig. 4.9.
Fig. 4.8 Rayleigh Doppler spectrum generator
-400 -300 -200 -100 0 100 200 300 400
0 20
-20 -40 -60 -80 -100 -120
Frequency (Hz)
eensity (dPowr Spectrum D
Fig. 4.9 Rayleigh fading Doppler spectrum by jake’s model
B)
Rayleigh Fading Spectrum Doppler Frequency : 100Hz
C. Rice fading
The Rice fading spectrum is similar to Rayleigh fading spectrum, but there are one direct fects is like path between the transmitter and the receiver. This direct path with Doppler ef
pure Doppler spectrum, but the other path is similar to Rayleigh fading spectrum.
The Rice fading Doppler spectrum is generating by Jake’s model which Doppler frequency is 100Hz is shown in Fig. 4.10.
-400 -300 -200 -100 0 100 200 300 400 Doppler Frequency : 100Hz
4.2.4
The additive requency
domain is defined by :
Fig. 4.10 Rice fading Doppler spectrum by jake’s model
Additive White Gaussian Noise (AWGN)
white noise is added after multipath fading channel. The SNR in f
2
the variance of AWGN. Therefore, the signal in real part and in imaginary part should be
added AWGN which variance is : 2 2 1
discussed t of sinc
ediate value that the samp
Carrier Frequency Offset and Sampling Clock Offset model
The detailed signal model of CFO is already described in section 2.1.1 and will not be repeatedly in this section. The model of SCO is built based on the concep
interpolation. The input digital signals can be exploited to interpolate the interm between two successive samples by using the shifted value of sic function. Assume
ling period is T and SCO is s ζ . Then the sampling phase can be represented as
The performa the proposed channel estimation scheme is illustrated in this
ation algorithm in chapter 2. The 3. First, we do the channel estimator for pilot signal, cond, we do the transform domain processing finally do the interpolation process. First we ke the overall respec
nalysis
nce analysis of
section. Except the verification of the proposed algorithms, some performance or computational complexity comparison between different methods and proposed scheme will be made too. We introduced the proposed channel estim
proposed scheme is depicted in Fig. 2.
se
discuss the interpolation process result, then discuss each block finally ma
syste
4.3.1 Interpolati
After interpolation in tim ts, we can get estimated CFR every three subcarriers. lation results to verify the algorithm
ic channel (Gaussian, Ricean , Rayleigh) which cified in DVB-T standard [3] to compare the MSE results. Here we use the MSE
he frequency interpolators.
m performance. The simulation environment is 2k mode, GI=1/4, 64-Qam, code rate=2/3.
A good compromise between bandwidth efficiency and robustness is 64-QAM with code rate 2/3 (mode used in UK). Furthermore, the performance of the overall system is based on the quasi-error free. The quasi-error free (QEF) condition which is corresponds to 2x10-4 BER after Viterbi decoder will be the system erformance target. p
on in frequency domain
e domain between scattered pilo Here we will make simu
mentioned in section 2.5.2.
Fig.4.11~Fig.4.13 are the simulation results of the different interpolation methods under different channel condition. First we use stat
are spe
criterion to discuss the performance of t
0 10 20 30 40 50
100 Gaussian channel awgn only
SNR(dB)
Fig. 4.11 MSE of different frequency interpolator in Gaussian channel
0 10 20 30 40 50 10-6
10-5 10-4 10-3 10-2 10-1
100 Ricean channel awgn only
SNR(dB)
MSE
Cubic interpolator Linear interpolator Parabolic interpolator Average interpolator Second-order interpolator
44 45 46 47 48 49 50
10-5
Ricean channel awgn only
SNR(dB)
MSE
Cubic interpolator Linear interpolator Parabolic interpolator Average interpolator Second-order interpolator
Crossover point
Fig. 4.12 MSE of different frequency interpolator in Ricean channel
0 10 20 30 40 50 10-6
10-5 10-4 10-3 10-2 10-1 100
101 Rayleigh channel awgn only
SNR(dB)
MSE
Cubic interpolator Linear interpolator Parabolic interpolator Average interpolator Second-order interpolator
37 38 39 40 41 42 43 44 45
10-4
Rayleigh channel awgn only
SNR(dB)
MSE
Cubic interpolator Linear interpolator Parabolic interpolator Average interpolator Second-order interpolator
Crossover point
Fig. 4.13 MSE of different frequency interpolator in Rayleigh channel
Here, we define the MSE is the sum of MSE expectation at all subcarriers which is expressed
Here, we define the MSE is the sum of MSE expectation at all subcarriers which is expressed