國 立 交 通 大 學
電信工程學系碩士班
碩士論文
用於數位影像廣播系統之以多重路徑干擾消除為基礎之通
道等化方法的設計與模擬
Design and Simulation of A Multipath Interference Cancellation
Based Channel Equalization Method in DVB-T Systems
研 究 生 : 林家輝
指導教授 : 黃家齊 博士
用於數位影像廣播系統之以多重路徑干擾消除為基礎之通
道等化方法的設計與模擬
Design and Simulation of A Multipath Interference Cancellation
Based Channel Equalization Method in DVB-T Systems
研究生:林家輝 Student: Chia-Hui Lin
指導教授:黃家齊 博士 Advisor: Dr. Chia-chi Huang
國 立 交 通 大 學
電 信 工 程 學 系 碩 士 班
碩士論文
A Thesis
Submitted to Department of Communication Engineering
College of Electrical Engineering and Computer Science
National Chiao Tung University
in Partial Fulfillment of the Requirements
for the Degree of
Master of Science
In
Communication Engineering
July 2004
Hsinchu, Taiwan, Republic of China
用於數位影像廣播系統之以多重路徑干擾消除為基礎之通
道等化方法的設計與模擬
研究生:林家輝 指導教授:黃家齊 博士
國立交通大學
電信工程學系碩士班
摘要
在現今的社會中,多媒體以及電腦通訊扮演著一個重要的角色。垂直 正交多工(OFDM)調變具有高速資訊傳輸的能力。目前已應用在新一代 的數位廣播傳輸技術,如數位音訊廣播(DAB)及數位視訊廣播之地面廣 播(DVB-T)等系統中。DVB-T 的系統擁有比傳統類比電視更高畫質的收 視品質(HDTV、SDTV)。 本論文將介紹用於 DVB-T 系統的 OFDM 通道等化方法。在接收端利 用估計的通道計算路徑的延遲,重建多重路徑的干擾。利用反覆演算的機 制來改善通道等化器的效能。另外,使用通道編碼(華士碼)產生傳送資 料的關連性進而增加編碼增益。以電腦模擬的方式來驗證並分析所設計之 通道等化方法對於 DVB-T 系統性能之提升。Design and Simulation of A Multipath Interference Cancellation
Based Channel Equalization Method in DVB-T Systems
Student: Chia-Hui Lin Advisor: Dr. Chia-chi Huang
Institute of Communication Engineering
National Chiao Tung University
Abstract
Nowadays, multimedia and computer communication are playing an important role. Orthogonal Frequency Division Multiplexing(OFDM), which makes information transmitted at a high speed, has applied to the digit of new generation broadcasts transmission technology, such as the DAB system and the DVB-T system. The DVB-T system has higher quality (HDTV, SDTV) than the traditional analogue TV system.
This thesis will introduce an OFDM channel equalization method which is used in the DVB-T system. Utilizing the estimated channel response in receiver to compute the path delay, the multi path interference can be reconstructed. By the implementation of iterative algorithm, it will improve the performance of the channel equalization. In addition, we can use the channel coding to produce the correlation between the transmitting data, and then possess the coding gain. The computer simulation method will verify and analyze the effect of the channel equalization method to the DVB-T
誌 謝
首先我要特別感謝我的指導老師-黃家齊博士,謝謝他幾年來在
學業上及生活上的耐心和細心指導,使我得以順利完成碩士學位。
感謝我的父親、母親,謝謝他們的養育之恩並且給予我無限的支
持,使我得以專心致力於我的研究。姐姐慧娟,謝謝她陪我渡過成長
的歲月。同時也謝謝所有關心我的家人。
謝謝實驗室所有的學長姐及同學們,感謝他們在研究學涯裡旳照
顧與提攜以及在學業上的教學相長使我獲益良多。
最後,僅以此論文獻給所有關心我的人,願大家順心如意。
目錄
第一章 正交分頻多工系統簡介(OFDM) ... 1 1.1 系統簡介... 1 1.2 基於快速傅立葉轉換之正交分頻多工系統... 4 1.3 護衛間隔的放置... 6 1.4 正交分頻多工系統之優缺點... 7 第二章 正交分頻多工系統之應用-數位影像廣播之地面廣播系統(DVB-T System) ... 9 2.1 系統簡介... 9 2.2 傳輸訊號及碼框結構... 13 2.3 參考訊號... 14 2.3.1 連續領航訊號(Continual Pilots) ... 152.3.2 傳輸參數訊號(Transmission Parameter Signaling, TPS) ... 16
2.3.3 散射領航訊號(Scattered Pilots)... 18 第三章 傳輸通道之模型及估測... 20 3.1 衰變通道... 21 3.2 通道模型... 22 3.3 通道的估計... 23 第四章 以多重路徑干擾消除為基礎之通道等化器... 25 4.1 通道等化器架構... 25 4.2 通道估測... 27 4.2.1 一維線性內插法... 27 4.2.2 基於快速傅利葉轉換之一維內插法(頻域上的 Sinc 內插)... 28
4.3 逼零等化器... 35 4.4 路徑干擾重建... 37 4.5 時間位移與最大比例合併... 39 第五章 電腦模擬與討論... 40 5.1 模擬環境及參數... 40 5.2 模擬的結果及討論... 41 5.2.1 通道等化器 without MPIC ... 41 5.2.1.1 AWGN 通道... 41 5.2.1.2 行動接收通道... 42 5.2.2 MPIC... 44 第六章 結論與未來方向... 50 參考文獻... 51
表目錄
表 2-1 8K 及 2K 兩種傳輸模式在 8MHz 通道的相關參數表 ... 14 表 2-2 傳輸參數訊號所攜帶之 TPS 位元的功能表 ... 17 表 5-1 模擬的環境與參數 ... 40
圖目錄
圖 1.1 (a)傳統分頻多工系統所使用之次載波 (b)正交分頻多工系統所使用之次 載波 ... 2 圖 1.2 正交分頻多工次載波示意圖... 2 圖 1.3 (a)正交分頻多工系統之發射機 (b)正交分頻多工系統之接收機 ... 4 圖 1.4 以快速傅立葉轉換為基礎的正交分頻多工系統... 5 圖 1.5 護衛間隔的放置... 6 圖 1.6 護衛間隔對抗多重路境延遲示意圖... 7 圖 2.1 數位影像廣播之地面廣播系統方塊圖... 10 圖 2.2 階層式 64 點二維振幅調變訊號 (a)α =1和 (b)α =4之星象圖 ... 11 圖 2.3 三種參考訊號次載波和資料訊號次載波位置分佈示意圖... 15 圖 2.4 連續領航訊號次載波和其他訊號次載波位置分佈示意圖... 16 圖 2.5 傳輸參數訊號次載波和其他訊號次載波位置分佈示意圖... 17 圖 2.6 散射領航訊號次載波和其他訊號次載波位置分佈示意圖... 19 圖 3.1 (a)頻率選擇性衰變通道 (b)平緩衰變通道 ... 20 圖 3.2 多重路徑時變通道脈衝響應模型... 21 圖 3.3 在延遲時間上無法區分之窄頻複數訊號... 22 圖 3.4 多路徑衰變通道之模型... 23 圖 4.1 正交分頻多工系統架構下傳送端及通道的影響... 25 圖 4.2 正交分頻多工系統架構下接收端資料重建... 26 圖 4.3 數位影像廣播之地面廣播系統之散設領航訊號... 27 圖 4.4 一維線性內插法示意圖... 28 圖 4.5 (a)完整的 Least Square 通道估計 (b) Least Square 通道估計之取樣 (c)取圖 4.6 應用於 DVB-T 系統之基於快速傅利葉轉換內插通道估計法 ... 30 圖 4.7 通道之最大延遲小於等於護衛間隔... 31 圖 4.8 時域上通道脈衝響應重複的出現... 33 圖 4.9 (a)DVB-T 的符元架構 (b)額外的干擾 ... 33 圖 4.10 (a)完整的通道估計值|hˆLS|(延遲為 50) (b)僅使用 1705 點估計出的通道 ' ˆ |hLS |... 34 圖 4.11 (a)完整的通道估計值|hˆLS |(延遲為 500) (b)僅使用 1705 點估計出的通道 35 圖 4.12 傳送訊號經通道至接收端簡圖... 36 圖 4.13 通道估計所估計出的時域通道特性... 37 圖 4.14 傳送資料經過編碼的機制... 38 圖 5.1 逼零等化器在 AWGN 通道下的 BER... 41 圖 5.2 逼零等化器在行動接收通道下的 BER(車速 30km/hr) ... 42 圖 5.3 逼零等化器在行動接收通道下的 BER(車速 90km/hr) ... 42 圖 5.4 逼零等化器在兩條路徑下行動接收通道下的 BER... 43 圖 5.5 逼零等化器在三條路徑下行動接收通道下的 BER... 43 圖 5.6 AWGN 通道使用相同能量比例下的 BER(a)完美通道估計 (b)通道 估計 ... 44 圖 5.7 AWGN 通道使用相同路徑延遲下的 BER(a)完美通道估計 (b)通道 估計 ... 45 圖 5.8 行動接收通道相同能量比例使用完美通道估計下的 BER(a)車速 30 公 里 (b)車速 120 公里 ... 46 圖 5.9 行動接收通道相同路徑延遲使用完美通道估計下的 BER(a)車速 30 公 里 (b)車速 120 公里 ... 47
圖 5.10 行動接收通道相同能量比例使用通道估計下的 BER(a)車速 30 公里 (b)車速 120 公里 ... 48 圖 5.11 行動接收通道相同路徑延遲使用通道估計下的 BER(a)車速 30 公里
第一章 正交分頻多工系統簡介(OFDM)
正交分頻多工技術是由多載波的調變方式發展而來,大約在 1960 年代被提出。 傳統的單一載波系統是將一個頻寬內的資料,全部載在一個載波上傳輸,因此,當 信號的傳輸速率越大時,所需的頻寬越大,但其傳輸的符元長度也越短,當其通過 一多重路徑延遲的通道時,碼際干擾 (Inter Symbol Interference) 會因符元的長度變 短而顯得更為嚴重。而多載波傳輸的原理則是將原本傳送速率快、寬頻的信號切分 成許多個傳送速率慢、窄頻的訊號載在不同的載波頻率上,形成多個次載波訊號, 然後再合併一起傳送[1],因此它的符元長度變長,能夠減輕多重路徑所造成的碼際 干擾。也由於其次載波間彼此正交的特性,因而提高頻譜的使用效率。
1.1 系統簡介
一個分頻多工系統在頻譜的配置上把使用的寬頻頻道分成許多窄頻頻道,為了 使每個次載波頻道間所載的資料彼此不干擾,直覺上所有的頻道必須分開來,甚至 在次頻道與次頻道之間要預留一些頻寬以確保不會有次頻道間訊號重疊。在正交分 頻多工系統中為了能更進一步節省頻寬,因此選擇彼此正交的載波頻率,如圖 1.1 (a)、(b)所示。利用載波正交性質允許頻道重疊,以節省頻寬使用,即 2 2 * 0 0 1 ( )( ) , 0 T j fnt j fnt k n k e e f f k f n k π π =⎧ = = + ∆ ⎨ ≠ ⎩∫
其中 f 為第k個次載波頻率;k f 為中心頻率;0 ( )• 為共軛複數 (Complex conjugate)。 * (a) (1.1)(b) 圖1.1 (a)傳統分頻多工系統所使用之次載波 (b)正交分頻多工系統所使用之次載波 k = 0 1 2 3 4 圖 1.2 正交分頻多工次載波示意圖 一個正交分頻多工系統之發射機與接收機架構如圖 1.3(a)、(b)。圖 1.3 中經過 資料編碼(Data Encoder)後為複數資料,可表示為 ( )d k =a k( )+ jb k( ),假設有N 個次 載波,則每個次載波上的訊號可表示成d tk( )=a tk( )+ jb tk( ),則發射機所發射的訊 號可表示為
{
} {
}
-1 0 -1 0 1 ( ) ( ) 1( ) cos( ) ( )sin( ) ( )sin( ) ( ) cos( )
k N jw t k k N k k k k k k k k k D t d t e N a t w t b t w t j a t w t b t w t N = = = = − + +
∑
∑
若每隔 t∆ 取樣一次,wk =2π fk,則(1-2)式可寫成 Frequency Saving of bandwidth (1.2)N-1 2 n 0 1 D(n )= (n ) N k j f t k k t d t e π ∆ = ∆
∑
∆ 若 f =0,c ∆f=N t1∆ 則 k k N t f = ∆ ,所以(1.3)式可化簡為與反離散傅立葉轉換(Inverse Discrete Fourier Transform,IDFT)相同的等式,即[ ]
-1[ ]
2 0 1 nk N N j k k D n d n e N π = =∑
所以我們可將反離散傅立葉轉換應用於OFDM 上,使其容易實現。 Data Encoder S/P Multiplex Input d(k) {a(k)+jb(k)} a(0) b(0) t w0 cos t wN 1 cos − t wN 1 sin − t w0 sin a(N-1) b(N-1) D(t) Channel k c 1 f f f , f N t k = + ∆ ∆ = ∆ (a) (1.3) (1.4)Data Decoder P/S Output t w0 cos t wN 1 cos − t wN 1 sin − t w0 sin Channel Integrator Integrator Integrator Integrator ˆ(0) a ˆ(0) b ˆ( 1) a N− ˆ( 1) b N− ˆ( ) d k (b) 圖1.3 (a)正交分頻多工系統之發射機 (b)正交分頻多工系統之接收機
1.2 基於快速傅立葉轉換之正交分頻多工系統
前一節我們提到可以使用快速傅利葉轉換來產生正交分頻多工系統之傳輸信 號,在本節我們將說明一個以快速傅利葉轉換為基礎而設計的正交分頻多工系統 [2]。 在圖 1.4 中,輸入的序列位元資料首先經由一個串列轉並列的轉換器(S/P Converter)分成 N 個複數資料,作為 N 個正交分頻多工系統次載波的輸入,每個複 數資料都包含數個位元,以便符合複數資料在訊號空間(Signal Space)上的位置, 而這些位置和正交分頻多工系統所使用的調變技術有關,例如可使用四相位移鍵調 變(QPSK)或 16 點二維振幅調變(16-QAM)。此時的複數資料可以視為在頻域的訊 號。接著這N 個複數資料經由反快速傅立葉轉換轉換到時域,而反快速傅立葉轉換 的輸出經由一個並列轉序列的轉換器(P/S Converter)變成序列形式。並在每個有 效符元的前端我們插入護衛間隔(Guard Interval)用以避免由多路徑通道所產生的碼際干擾。然後這些離散的資料經由數位轉類比轉換器(D/A Converter)轉成類比形 式且經過低通濾波以限定其傳輸頻寬,再將訊號傳送出去。接收機的工作大體上為 發射機的反運作,唯有其中單軌等化器(One-Tap Equalizer)的作用為修正每個次 載波通道上所產生的失真,因此每個等化器上的相乘係數(Tap Coefficient)也就和 每個次載波通道的特性有關。使用快速傅立葉轉換來實現離散傅立葉轉換的優點為 運算量可由N 個相乘計算次數減少為2 N log N⋅ 2 個相乘計算次數,大量的減少了在 實作上的複雜度。
S/P MapperSignal IFFT P/S
Guard Interval Insertion D/A LPF Up Converter Channel Down Converter P/S MapperSignal FFT S/P Guard Interval Removal LPF A/D One-tap Equalizer Serial Data Input Serial Data Output x bits d(0) d(1) d(N-1) x bits 圖1.4 以快速傅立葉轉換為基礎的正交分頻多工系統
1.3 護衛間隔的放置
在多重路徑的通道下,為了避免碼際干擾的發生,在正交分頻多工系統中我們 可以用增長符元時間的方法來避免這樣的問題。其作法就是在每個訊號的前面加個 護衛間隔,也就是TSYM =TU+ 。Tg T 是有效符元(Useful Symbol)時間,U T 是護衛間g 隔。護衛間隔中所放置的是有效符元循環延伸(Cyclic Extension)的訊號。在每個
正交分頻多工的訊號之前就是它自己的週期延伸,圖 1.5 表示出一個插有護衛間隔
的正交分頻多工符元。因為在時域中的循環平移(Cyclic Shift)會在頻域上造成線 性相位偏移,而線性相位偏移可用差分檢測技術(Differential Detection Technique) 消除,或者此線性相位偏移效應將合併於通道響應中再利用領航次載波一起估計並 加以補償之,所以在不受碼際干擾的區間內開始取樣,其對應之有效符元時間內的 資料都可以用來正確地解調。
Guard interval Useful symbol Copy
圖1.5 護衛間隔的放置
前置的護衛間隔使得符元與通道由原本的線性旋積(Linear Convolution)變為環旋積 (Circular Convolution),只要接收訊號取樣在非碼際干擾區間(ISI Free Region),則取 樣結果在時域上相當於符元和通道的環旋積,經過離散傅利葉轉換在頻域上等效於 符元和通道的乘積,在兩個相鄰符元通道不變或變化不大的假設下,利用差分偵測 消除通道效應,或利用通道估計,經過等化器消除通道對訊號的影響。但是若多重
路徑延遲大於護衛間隔,則在解調的時候就會有碼際干擾的問題出現,如圖 1.6(b)
Undesired Data Path1 Path 2 Path 1 Path 2 τ g (a) : τ <T g (b) : τ >T g T U T One OFDM Symbol
U
T One OFDM Symbol τ
Sampling Window
Sampling Window ISI Free Region
g
T
圖1.6 護衛間隔對抗多重路境延遲示意圖
1.4 正交分頻多工系統之優缺點
正交分頻多工系統將一個寬頻信號分成許多個窄頻信號傳送,因此每個次載波 上的窄頻訊號,可視為經過平緩衰變通道(Flat Fading Channel),只需要經過簡單的 等化器就可被還原,同時,符元的長度增加,並加上護衛間隔的放置,使得符元對 於碼際干擾的效應有較佳的抵抗能力,因此能在低複雜度接收機架構下達到高速傳 輸的實現。除此之外,利用彼此正交的次載波,更能有效的節省頻寬。
因為次載波間彼此的正交特性,使得系統對於頻率的同步要求非常嚴格,只要 有些許的頻率漂移(Frequency offset),頻率之間的正交特性即被破壞而產生次載波間 干擾(Inter-carrier interference, ICI)。系統的另一個缺點為容易造成非線性失真。這是
因為多個載波訊號輸出時相加在一起,使得其輸出功率的尖峰 對平均之比值 (Peak-to-average power ratio, PAPR)變得很大,系統容易出現非線性失真,亦或使得 功率放大器的功率效率(Power efficiency)降低及數位/類比轉換器(D/A converter)與 類比/數位轉換器(A/D converter)所需的位元數提高,這些缺點對系統都有極大的影
第二章
正交分頻多工系統之應用-數位影像廣播之地面廣播
系統(DVB-T System)
歐規的數位影像廣播之地面廣播系統(DVB-T System)是採用正交分頻多工 傳輸技術來傳送訊號。如同在第一章所述,由於正交分頻多工傳輸技術的特性, 特別是可以抵抗碼際干擾,全世界許多國家在數位影像廣播之地面廣播系統的規 格上都採用了此歐規的系統,其中台灣也將在2004 年六月一日正式開播,相較 於傳統類比電視系統,數位電視提供了更好的收視品質(MPEG-2),同時也適用 於行動接收的各種應用,如車用視訊等。本章我們將對此系統做簡單的敘述:首 先,我們將介紹這個系統的系統方塊圖,包含說明每個方塊的功能;接著,我們 將展示這個系統的傳輸訊號以及碼框結構,包含各種傳輸模式的參數表以及各種 參考訊號的說明。2.1 系統簡介
圖2.1 是數位影像廣播之地面廣播系統方塊圖[3]。原始影像、聲音及一般資料在經過編碼壓縮(MPEG-2 coded TV signals)及傳輸層多工(Source Coding and Transport Multiplexing)產生一個個長 188 個位元組的封包(Packet),其中每個封包
含有1 個同步位元組(Sync-Word Byte=47HEX)和 187 個資料位元組(Data Bytes)。
每八個封包為一組,其中資料位元組的部份被乘上一個假隨機二進位序列 (Pseudo Random Binary Sequence),此假隨機序列的暫存器初始值亦每八個封包
重設一次。每八個封包的第一個封包其同步位元組由47HEX轉換成B8HEX。這些
動作稱作傳輸多工調整暨能量分散(MUX Adaptation,Energy Dispersal),其目的
Splitter Video Coder Audio Coder Data Coder Programme MUX Transport MUX 2 n 1 MPEG-2 Source coding and Multiplexing
MUX Adaptation, Energy Dispersal Outer Coder MUX Adaptation, Energy Dispersal Outer Coder Outer Interleaver Outer Interleaver Inner Coder Inner Coder Inner Interleaver Mapper Frame Adaptation Pilot & TPS Signals OFDM Guard Interval Insertion
D/A Front End
Terrestrial Channel Adapter To Aerial
圖2.1 數位影像廣播之地面廣播系統方塊圖
此系統使用階層式傳輸(Hierarchical Transmission),可以依照資料的優先順 序來傳送。首先利用分離器(Splitter)將傳輸資料流(Transport Stream)分為高優先 權(High Priority)資料及低優先權(Low Priority)資料;在解調時先判斷訊號座落在 訊號空間上的象限即可先解調出高優先權資料,再經由判斷訊號在某一象限中所 對應的位置來解調出低優先權資料。這兩筆資料各自有不同的通道編碼及調變方
式。圖2.2 是階層式 64 點二維振幅調變訊號α =1和α =4的星象圖。由此圖可知
當α 值愈大,這四個象限內的訊號點離原點愈遠,愈不容易將訊號點所在的象限 判斷錯誤,所以高優先權資料愈不會有錯誤。
100000100010 101010 101000 100001100011 101011 101001 100101100111 101111 101101 100100100110 101110 101100 110100110110 111110 111100 110101110111 111111 111101 110001110011 111011 111001 110000110010 111010 111000 001000001010000010 000000 001001001011000011 000001 001101001111000111 000101 001100001110000110 000100 011100011110010110 010100 011101011111010111 010101 011001011011010011 010001 011000011010010010 010000 Re{Z} Im{Z} 1 3 5 7 9 1 3 5 7 9 -1 -3 -5 -7 -9 -1 -3 -5 -7 -9 1 ) a ( α= 4 10 8 6 -4 -10 -8 -6 4 6 8 10 -4 -10 -8 -6 100000 100010101010 101000 100001 100011101011 101001 100101 100111101111 101101 100100 100110101110 101100 110100 110110111110 111100 110101 110111111111 111101 110001 110011111011 111001 110000 110010111010 111000 001000 001010000010 000000 001001 001011000011 000001 001101 001111000111 000101 001100 001110000110 000100 011100 011110010110 010100 011101 011111010111 010101 011001 011011010011 010001 011000 011010010010 010000 Re{Z} Im{Z} 4 ) b ( α= 圖2.2 階層式 64 點二維振幅調變訊號 (a)α =1和 (b)α =4之星象圖 系統共有兩層通道編碼器及交錯器。外層通道編碼是使用里德-所羅門短碼 (Reed-Solomon Shorten Code,RS(204,188,t=8)),使每個封包由 188 個位元組變
為204 個位元組,並且可以校正最長為 8 個位元組的連續錯誤﹔外層交錯器是一
個以位元組為單位的迴旋交錯器(Convolutional Interleaver),這個交錯器有 12 個
的時間延遲量(深度)是 1122 個位元組(即是 5.5 個封包長)為一般區塊交錯器 (Block Interleaver)延遲量的一半。內層通道碼是編碼率為 1/2 的迴旋碼,藉由打 孔法可以將編碼率提高至2/3、3/4、5/6 或 7/8。內層交錯器是由一個深度為 126 位元的位元式區塊交錯器以及一個深度為一個正交分頻多工符元所擁有的資料 次載波數目的區塊交錯器。以64 點的二維振幅調變訊號為例,其位元訊號分成 六個分支,每個分支經過一種深度為126 位元的位元式區塊交錯器(共六種),然 後這六個分支上的位元訊號合併對映成一個個64 點二維振幅調變的複數資料訊 號,這些複數資料訊號再經過一個以複數資料為單位的區塊交錯器(其深度以 2K 模式而言是1512,以 8K 模式而言是 6048)。交錯器的功能是用來降低傳輸通道 以及加成性白高斯雜訊對訊號造成的錯誤率,特別是當訊號干擾是隨機產生時它 的效果最好。 經過內層交錯器的資料根據不同的調變方式將資料對映成複數資料,可使用 的調變方式共有四相位移鍵調變、16 點二維振幅調變、64 點二維振幅調變 (64-QAM)、不均勻的 16 點二維振幅調變(Non-uniform 16-QAM)以及不均勻的 64 點二維振幅調變(Non-uniform 64-QAM)。所謂的”不均勻”是指階層式調變技術中 α 值不為 1 的情形,如圖 2.2(b)所示。 這個系統提供兩種傳輸模式:2K 模式以及 8K 模式。根據不同的模式使用 不同長度(2048 or 8192)的反快速傅立葉轉換以得到正交分頻多工的有效符元,然 後在每個有效符元之前會加上一段護衛間隔,形成一個完整的符元。接著這個數 位訊號經過數位/類比轉換器(Digital-to-Analog Converter,DAC)轉成類比訊號, 然後通過前級(Front End),最後將傳輸訊號傳送出去。除了一些載有資料的次 載波之外還有一些特定的次載波上載著已知的參考訊號(Reference Signal)。參考 訊號共分為三種:連續領航(Continual Pilots)訊號、散射領航(Scattered Pilots)訊號 以及傳輸參數訊號(Transmission Parameter Signaling,TPS)。整個系統屬於哪一 種傳輸模式、護衛間隔的比例、是否使用階層式調變技術、內層通道編碼的編碼
位元區塊當中。這三種參考訊號的放置的方式、量值以及其用途將會在本章的第 三小節說明。
2.2 傳輸訊號及碼框結構
每個傳輸的碼框包含68 個正交分頻多工符元,每四個碼框再構成一個大碼 框(Superframe)。在時間 t 所傳的訊號可由下列數學式表示: max min K 67 2 , , F SYM 0 K S(t)=Re{ j f tc ( T T )} m l k k m l k e π Z g t m l ∞ =−∞ = = ⋅ − ⋅ − ⋅∑ ∑ ∑
(2-1) 其中 2 K K , 68T T , ) T Rect( ) ( max min SYM F SYM ) T ( T 2 g U ⋅ = ′= − − = − ′ k k t e t g t k j k π (2-2) 每個符元由[K , Kmin max] 之間次載波上所載的複數資料Zm l k, , (為第m個碼框的第l 個符元上第k個次載波所載的資料)經反快速傅立葉轉換後得到。 f 是載波的中c 心頻率;TSYM是一個完整符元的長度,包含一個有效符元的長度(T )及一個護U 衛間隔的長度(T ),g TF是一個碼框的長度﹔ SYM Rect( ) T t 是一個寬度為TSYM的方 波。 在頻寬為8MHz 的系統中,其載波的中心頻率位置為︰ 470MHz+4MHz+i 8MHz, i=0,1,2,3,...× (2-3) 兩種傳輸模式相關的參數如表 2-1 所示,其基本取樣週期(Elementary Sampling Period : T)皆為 7 64µs。由此表可知,因為8K 模式的快速傅立葉轉換長度(N) 是 2K 模式的四倍,所以其有效符元長度為 2K 模式的四倍而次載波間隔為 2K 模式的四分之一。又因為8K 模式的主動次載波個數(K)約為 2K 模式的四倍, 所以兩個傳輸模式的有效頻寬皆約為7.61MHz。每一種傳輸模式都可根據需要選 擇適當的護衛間隔長度(為有效符元長度的 1/4、1/8、1/16 或 1/32)。另外還有頻寬為7MHz 及 6MHz 的系統其基本取樣週期(T)分別為1 8µs及 7 48µs,頻寬愈小其相對的基本取樣週期愈長。這兩種系統也同樣各有其 2K 模 式和8K 模式,而他們的相關參數中,全部次載波的個數(N)、主動次載波個數 (K)、Kmin和Kmax的量值和頻寬為8Mz 的系統一樣,如表 2-1 中所示。由此可 以算出頻寬為7MHz 及 6MHz 的系統其有效頻寬約為 6.66MHz 及 5.71MHz。同 樣的,每一種傳輸模式都可根據需要選擇適當的護衛間隔長度(為有效符元長度 的1/4、1/8、1/16 或 1/32)。 sec 896µ U 1/T U T = ⋅N T K max K min K max min and K K carriers between Spacing (K-1)/TU
Parameter 8K mode 2K mode
Number of active carriers Vaule of carrier number Vaule of carrier number Duration Carrier spacing 6817 1705 0 0 6816 1704 1116 Hz 4464 Hz 7.61 MHz 7.61 MHz sec 224µ N
Number of total carriers 8192 2048
g U T / T The ratio of toTg TU 1 1 1 1 , , or 4 8 16 32 1 1 1 1 , , or 4 8 16 32 表2-1 8K 及 2K 兩種傳輸模式在 8MHz 通道的相關參數表
2.3 參考訊號
在數位影像廣播-地面廣播系統中,訊號同步、傳輸方式的相關參數以及通 道估計等相關資訊均可利用已知的參考訊號以及護衛間隔來求得。而其他的TPS 位元將提供此系統其他相關參數的訊息。通道頻率響應的估計則是先利用散射領 航訊號(Scattered Pilots)來得到通道頻率響應的取樣訊號,再將此取樣訊號做內插 來求得。接下來,我們將分別敘述這三種參考訊號的分佈位置、量值。 圖2.3 是三種參考次載波訊號和資料次載波訊號的示意圖,實心點是三種參分別為連續領航訊號(Continual Pilots)、傳輸參數訊號(Transmission Parameter Signaling, TPS)以及散射領航訊號(Scattered Pilots)。此圖的橫軸標示次載波的順
序,而縱軸則標示時間順序。次載波的順序由第Kmin個次載波到第Kmax個次載 波。時間順序是指這個次載波所在的時間位置是在哪一個碼框中的哪一個符元, 一個碼框有68 個符元,所以符元順序每 68 個循環一次。每一列的實心點及空心 點即是某一碼框中的某一符元在頻率軸上其次載波分佈的情形。 min K Kmax 圖2.3 三種參考訊號次載波和資料訊號次載波位置分佈示意圖
虛擬隨機二位元序列(Pseudo Random Binary Sequence,PRBS)w 是由隨機的k
0 或 1 組成,這個序列被拿來當作三種參考訊號的量值。它的值只隨次載波的位 置而變,不會因為碼框時間或符元時間的不同而改變。連續領航次載波訊號及散 射領航次載波訊號的量值由此序列產生,此序列也決定傳輸參數訊號在第一個符 元的量值。 2.3.1 連續領航訊號(Continual Pilots) 連續領航次載波訊號上的值Zm,l,k由虛擬隨機二位元序列決定,其值為 4 3 ± , 即是 ) 5 . 0 ( 2 3 4 , ,lk k m w Z = ⋅ ⋅ − 圖2.4 特別標示出連續領航次載波訊號。如圖 2.4 所示,實心點為連續領航次載 (2.4)
波訊號,而空心點則為其它次載波訊號。連續領航次載波訊存在於每個符元的固 定次載波位置上,同時在所有符元的相同連續領航次載波上的值相同。就8K 模 式而言,一個符元中含有177 個連續領航次載波訊號;而就 2K 模式而言,一個 符元中含有45 個連續領航次載波訊號。這些連續領航次載波訊號的位置有一個 特定的組合,而整數載波頻率偏移則會造成這個位置組合整體的循環平移。因 此,偵測這個連續領航次載波位置組合的循環平移量即可得到整數載波頻率偏移 的估計值。 symbol 67 symbol 3 symbol 2 symbol 1 symbol 0 min
K continual pilot others Kmax
圖2.4 連續領航訊號次載波和其他訊號次載波位置分佈示意圖
2.3.2 傳輸參數訊號(Transmission Parameter Signaling, TPS)
傳輸參數訊號次載波上的值為± 。第一個符元之傳輸參數訊號次載波上的1 值由虛擬隨機二位元序列決定,即 ) 5 . 0 ( 2 , 0 , k k m w Z = ⋅ − 而其它符元上之傳輸參數訊號次載波的值則為差分編碼前一符元上的值所得來 且帶有一個值為0 或 1 的 TPS 位元S ,即 l 67 ,...., 2 , 1 ) 1 ( , 1, , , = − ⋅Z − l= Z ml k S k l m l ) 圖2.5 特別標示出傳輸參數訊號次載波。如圖 2.5 所示,實心點為傳輸參數訊號 次載波,而空心點則為其它次載波訊號。就 8K 模式而言,一個符元中含有 68 個傳輸參數訊號次載波;而就 2K 模式而言,一個符元中含有 17 個傳輸參數訊 (2.5) (2.6)
號次載波。 symbol 67 symbol 3 symbol 2 symbol 1 symbol 0 min K TPS others Kmax 圖2.5 傳輸參數訊號次載波和其他訊號次載波位置分佈示意圖 每一個符元上的所有傳輸參數訊號次載波所攜帶的TPS 位元皆相同,如等式(2.6) 所示,所以每個碼框共有67 個 TPS 位元。表 2-2 為這 67 個 TPS 位元的功能表, 我們接下來將會簡述這67 個 TPS 位元的使用目的。
Bit number
Purpose/Content
24 23
, s
s
29 28 27,
s
,
s
s
32 31 30,
s
,
s
s
35 34 33,
s
,
s
s
26 25, s
s
37 36, s
s
39 38, s
s
16 1s
s
22 17s
s
4047
s
s
67 54s
s
Synchronization word Length indicator Frame number Constellation Hierarchy informationCode rate, HP stream Code rate, LP stream
Guard interval transmission mode
Cell identifier Error protection
48
53
s
s
Reserved for future use16 synchronization bits
37 information bits
14 redundancy bits
表2-2 傳輸參數訊號所攜帶之 TPS 位元的功能表
(S1,S2,...,S16)有固定的組合,匹配這個組合可以找出符元的順序,即是碼框時 間。37 個資訊位元(S17,S18,...,S53)是用來控制系統的相關傳輸參數,例如調變 的方式、階層式調變的資訊、高優先權資料的編碼率、低優先權資料的編碼率、 護衛間隔和有效符元長度的比例,以及使用哪一種傳輸模式等等。最後14 個多 餘位元(S54,S55,...,S67)是用來保護之前的 53 個 TPS 位元,經由 BCH(67,53,t=2) 編碼將53 個 TPS 位元延伸為 67 個 TPS 位元,並增加保護的能力。 2.3.3 散射領航訊號(Scattered Pilots) 散射領航次載波訊號上的值Zm,l,k由虛擬隨機二位元序列決定,其值為 4 3 ± , 即 ) 5 . 0 ( 2 3 4 , ,lk k m w Z = ⋅ ⋅ − 圖2.6 特別標示出散射領航次載波訊號。如圖 2.6 所示,實心點為散射領航次載 波訊號,而空心點則為其它次載波訊號。由此圖可知,在一個符元中,每12 個 次載波便會放置一個散射領航次載波訊號。而在相鄰的符元中,散射領航次載波 訊號的位置組合會相對的間隔3 個次載波間隔。一個碼框中第l個符元的散射領 航次載波訊號其次載波的位置k為:
min min max
K 3 ( mod 4) 12 0, k [K ,K ]
k= + × l + ×p p≥ ∈
就2K 模式而言,Kmin=0, Kmax =1704;就8K 模式而言,Kmin =0, Kmax =6816。
在接收端取出存在於散射領航次載波位置上的值並且將其它次載波位置上的值 設成0,我們就可以估計出傳輸通道頻率響應(Frequency Response)的取樣訊號。 接著,我們可以依次針對時域及頻域使用內插法(Interpolation),如此我們就可以 得到傳輸通道頻率響應的估計值。 (2.7) (2.8)
symbol 67 symbol 3 symbol 2 symbol 1 symbol 0 min
K scattered pilot others Kmax
圖2.6 散射領航訊號次載波和其他訊號次載波位置分佈示意圖
第三章 傳輸通道之模型及估測
在無線通訊中,通道的效應是影響效能的重要關鍵,多重路徑的通道 (Multi-path),是造成接收訊號失真的主要因素,它是由於電波在空中傳播時碰到 物體產生反射所造成。如果接收端和傳送端有相對的運動,則會造成通道隨著時 間而改變,這樣的情形稱為通道的衰變(Fading)。衰變通道一般可分為頻率選擇 性衰變通道,與平緩衰變通道。當傳輸信號頻寬大於其通道的同調頻寬(Coherent Bandwidth)時,造成在頻譜上每個頻率的衰變都不同,如圖3.1(a)示,因此稱之為 頻率選擇性衰變(Frequency Selective Fading)通道,反之,我們稱之為平緩衰變通 道(Flat Fading)。圖3.1 (a) 頻率選擇性衰變通道 (b) 平緩衰變通道
通道的同調頻寬和通道延遲方均根值(Root Mean Square Value of Delay Spread)成 反比,也就是當通道延遲方均根值越大時,其同調頻寬則越小,反之則相反。同 調 頻 寬 經 常 被 定 義 為 兩 個 頻 域 上 不 同 頻 率 的 訊 號 彼 此 能 夠 維 持 相 關 係 數 (Correlation Coefficient)大於0.9 的最大頻率差。
3.1 衰變通道
對一個寬頻傳輸訊號的多重路徑衰變通道,可以表示成: 1 ( ) 0 1 0 ( , ) ( ) ( ) ( ) ( ) l L j t l l l L l l l h t A t e r t φ τ δ τ τ δ τ τ − = − = = − = −∑
∑
t 為時間,τl為第l個路徑的延遲時間,A t 為第l( ) l個路徑的增益,φl( )t 為第l個 路徑的隨機相位(Random Phase), ( )r t 所表示的是第l l個路徑的增益及其隨機相 位,如圖3.2 所示。 ( , ) h tτ
0 ( )t τ 1 ( )tτ
2 ( )t τ 3 ( )t τ 0 τ τ1 τ2 τl τL 0 t 1 t 2 t 3 t 圖3.2 多重路徑時變通道脈衝響應模型 接收機在一小區域內所收到的每個路徑,可視為來自各方向在的窄頻複數訊 號之和。每個路徑以複數訊號表示如下: , ( ) ( ) , 1 ( ) ( ) ( ) l l n j t l l N j t l n n r t A t e A t e φ φ = = =∑
N 表示在第l個路徑中,來自各方向在延遲時間上無法區分之窄頻複數訊號的數 (3.1) (3.2)目(如圖 3.3 所示)。當 N 值很大時,根據中央極限定理(Central Limit Theory), ( )
l
r t 將近似於一個複數高斯隨機程序(Complex Gaussian Random Process),因此
複數訊號之振幅會形成雷利(Rayleigh Distribution)分佈,其相位會形成均勻分佈 (Uniform Distribution)。因為 Jake’s Model[4]可以產生一個近似雷利分佈的序列,
適合用來描述一個窄頻訊號,所以可用來表示式中的r t 。 l( )
圖3.3 在延遲時間上無法區分之窄頻複數訊號
在 Jake’s model 中,一個等效基頻 (Equivalent Baseband)的窄頻衰變模型
(Narrow Band Fading Profile)可以表示成來自 N 個方向之複數訊號和:
( ) 0 ( ) 1 cos 0 ( ) ( ) , 2 l m n n j t l N j w t E m m N n T t A t e e w f φ α θ π − + = = =
∑
= Eo 是 T(t)的平均功率,fm是最大都卜勒(Doppler)頻率偏移量,αn是其入射角 度,θn為此方向入射訊號之隨機相位,假設αn和θn在[0 , 2π]內均勻分佈(Uniform Distribution)。3.2 通道模型
利用Jake’s model 我們可以來模擬所要的多重路徑衰變通道。假設共有 N 個 (3.3)的增益及相位都可由Jake’s Model 產生,如式(3-3)所示,且不同路徑的增益及相 位彼此互不相關。圖3.4 為多路徑衰變通道的模型圖。 0( ) 0( ) j t A t eφ 1( ) 1( ) j t A t eφ 1( ) 1( ) N j t N A t eφ − − 2( ) 2( ) j t A t eφ ∑ 1 τ τ2 τN−1 AWGN INPUT OUTPUT 圖3.4 多路徑衰變通道之模型
3.3 通道的估計
對一個同調調變的系統而言,在接收機完成同步的動作後接著就必須經過等 化器以消除通道所造成的振幅變動以及隋機相位偏移。因此我們必需作通道的估 計以提供等化器所須的傳輸通道資訊。 在正交分頻多工系統中,我們通常使用領航訊號(Pilots)來估計通道,所謂領 航符元(Pilot Symbols)即是在一個 OFDM 符元中,在頻域上間隔特定次載波數目 插入已知訊號。在接收端利用已知的 Pilot,我們就可以估計出通道。許多種的估計方法已 被提出,在此我們先介紹一種最基本的方法:最小平方誤差法(Least Square estimation )[5]。
首先,我們假設在正交分頻多工系統下,滿足下列的線性模型: Y = XH + Z 其中 0 1 0 1 0 1 0 1 [ , ,..., ] 0 0 0 0 0 0 0 0 0 0 0 0 [ , ,..., ] [ , ,..., ] T N N T N T N Y Y Y Y X X X X H H H H Z Z Z Z = ⎡ ⎤ ⎢ ⎥ ⎢ ⎥ = ⎢ ⎥ ⎢ ⎥ ⎣ ⎦ = = % 在此線性模型中,Y 為接收機收到的信號、X 為已知的 Pilot,H 為待估計的通道 頻率響應而Z 則為雜訊。估計的方法是找出一組解,HˆLS,使得其平方誤差最小 即:
將(3.6)式對HˆLS作偏微分,並令其為零即可解出HˆLS(Least Square solutions):
式(3.7)的概念是將接收到的信號轉到頻域,直接除以已知的 Pilot 就可以估計出 通道的頻率響應,這個方法的好處是它很容易實現,但由式(3.7)可知,收到的信 號 Y 是包含有雜訊的,所以這個方法並不能有效的抑制雜訊的干擾。所以,我 們通常會利用Least Square 的方法初步的估計出HˆLS,接著再利用HˆLS做進一步 的處理,細估出更精確的通道。 (3.4) (3.5) 2 min{ } min{( - ) ( -T )} e = Y XH Y XH (3.6) 2 1 1 ˆ 2 ( ) ˆ ˆ ( ) T LS LS T T LS e X Y XH H H X X X Y X Y − − ∂ = − ∂ ⇒ = = (3.7)
第四章
以多重路徑干擾消除為基礎之通道等化器
我們考慮在正交分頻多工系統的架構下,利用通道估計可以估出通道的特 性,加上等化器及決策器所估出的傳送資料,可以估計出傳送資料經過通道所受 到的效應,分別取出每一個路徑的訊號以排除多路徑的干擾,再以最大比例匹配 (Maximum ratio combining)重建出較為準確的接收訊號,繼而可以估計出較為 準確的傳送訊號。如此又可以反覆的進行,再得到更新之後的每一個路徑的訊 號,而又重建出更準確的接收訊號,最後得到更準確的傳送訊號。
4.1 通道等化器架構
IFFT insertionGI h[n]+
(L-path) [ ] y n [ ] w n [ ] x n [ ] d n [ ] D n 圖4.1 正交分頻多工系統架構下傳送端及通道的影響 圖 4.1 說明正交分頻多工系統架構下輸入的資料(D n )經由反快速[ ] 傅立葉轉換轉換到時域內(d n )[ ] ,而反快速傅立葉轉換的輸出加上護衛區間 (Guard Interval),得到我們的發送訊號(x n )[ ] ,經過通道之後加上雜訊的影響, 即為接收訊號(y n )[ ] 。 在接收端收到的訊號 y n 除去護衛區間後經過快速傅立葉轉換得到頻域的[ ] 訊號,通道估計所得到的通道特性 ˆ[ ]h n ,首先用在通道等化器,使用逼零等化器 消除通道的效應而得到估計的傳送資料 ˆD n ,其經過反快速傅立葉轉換可獲得[ ]時域訊號 ˆ[ ]d n ,而估計的時域傳送資料加上護衛區間後又可以得到估計的發送資 料 ˆ[ ]x n 。估計的通道特性 ˆ[ ]h n 加上估計的發送之料 ˆ[ ]x n 可重建出每一路徑的資 料。 ˆ[ ] h n [ ] y n d nˆ[ ] 1[ ] i n 2[ ] i n [ ] L i n 1[ ] i n 2[ ] i n [ ] L i n [ ] y n + + + − − − [ ] d n 1 [ 1] ˆ1[ ] a x t−τ +w n 2 [ 2] ˆ2[ ] a x t−τ +w n ˆ [ ] [ ] L L L a x t−τ +w n ˆ[ ] y n ˆ[ ] x n [ ] x n 圖4.2 正交分頻多工系統架構下接收端資料重建
4.2 通道估測
前一章我們提到利用已知的 Pilot 可以求出其 Lease square solutions,但這只
限於在放置 Pilot 的次載波位置,實際的正交分頻多工系統不可能所有的次載波 都放置 Pilot,圖 4.3 即為數位影像廣播之地面廣播系統的 Pilot 放置方法,其中 黑點的部分為已知的Pilot,在橫軸也就是頻域上,Pilot 的放置是每隔 12 個次載 波間隔,這種間隔特定次載波數目而放置 Pilot 的方法就稱為領航符元(Pilot Symbol)。 symbol 67 symbol 3 symbol 2 symbol 1 symbol 0 min
K scattered pilot others Kmax
圖4.3 數位影像廣播之地面廣播系統之散設領航訊號
利用這些已知的 Pilot,我們可以求相對應次載波位置上的 Least Square 解,
ˆ LS H ,其他沒有放置 Pilot 的次載波位置可以用內插法來求得,詳細方法如下面 的章節所述。
4.2.1 一維線性內插法[6]
如同式(4.1)所示,利用兩已知的值,可以內插出中間未知的值,這個方法的 好處是它很簡單,但是,當通道在頻域變化得很快的時候,Pilot 放置的密度必 須增加,否則這個方法將不適用。 ( ) [ ] [ ] { [( 1) ] [ ]} k iT , ( 1) H k H iT H i T H iT iT k i T T − = + + − < < + (4.1)iT k (i+1)T H[iT] H[k] H[(i+1)T] 圖4.4 一維線性內插法示意圖
4.2.2 基 於 快 速 傅 利 葉 轉 換 之 一 維 內 插 法 ( 頻 域 上 的 Sinc 內
插)[7][8][9]
如同圖 4.3 所示,如果我們利用 Least square 的方法求出放置 Pilot 次載波位
置的通道頻率響應值,剩下沒有放置Pilot 的次載波位置其值設為零,那麼就可 視為是對HˆLS做取樣: 其中T 為 Pilot 放置的週期,N 為所有的次載波數目。 將HˆLS在頻域上做取樣,經反傅利葉轉換到時域時,hˆLS將在時域上重複出現, 重複的週期為T/N,即: ˆS {ˆ [0],0,..., ˆ [ ],0..., ˆ [2 ],0,..., ˆ [( 1) ],0...,0} LS LS LS LS LS N H H H T H T H T T = − (4.2) 2 1 0 2 1 0 ^ ( ) [ ] 1 [ ] ( ) ˆ ( ) ( ) , , 0,1,... 1 N j kn N n N j kn N k S LS LS H k h n e h n H k e N N H k H k if k qT q T π π − − = − = = = = = = −
∑
∑
(4.3)由取樣定理得知,只要通道的最大延遲小於N/T,那麼,在時域上重複出現的hˆLS 就不會發生交疊(Aliasing)的情況,只要取出 ˆ [ ]S LS h n 中0 n N T ≤ ≤ 的部份,並乘上T 補償因取樣造成的振幅衰減,如此以來及可還原出完整的hˆLS,如圖4.5 所示。 ˆ |hLS[ ] |n | ˆ ( ) | LS H k (a) Time ˆ |hLS[ ] |n
⊗
Time N/T ˆ | S ( ) | LS H k (b) 2 1 0 1 2 ( ) 0 1 2 2 1 ( ) ( ) 0 0 2 1 ( ) 1 0 0 1 ˆ [ ] ˆ ( ) 1 ˆ = ( ) 1 ˆ = [ ] 1 ˆ = [ ] 1 ˆ = [ ] [( N j kn S S N LS LS k N T j qT n N LS q N N T j qT x j qT n N N LS q x N j q x n n N T N T LS x q LS h n H k e N H qT e N h x e e N h x e N N N h x x n N T T π π π π π δ − = − = − − − = = − − − − = = = − −∑
∑
∑∑
∑
∑
1 1 0 0 1 0 ) ] 1 ˆ = [( ) ] N T N x z T LS N z z N h n z T T − − = = − = +∑
∑
∑
(4.4)ˆ | S [ ] |
LS
h n
(c)
圖4.5 (a)完整的 Least Square 通道估計 (b) Least Square 通道估計之取樣 (c)取
樣後經IFFT,hˆLS在時域重複出現,振幅縮為1/T 圖 4.6 描述了一個應用於 DVB-T 系統的基於快速傅利葉轉換內插通道估計 法, LS 的通道估計是無法有效濾除雜訊的,因此我們必需使用符元平均(Symbol by Symbol Average)與路徑選擇的方法來抑制雜訊,在做符元平均時要選擇幾個 符元來平均有下列兩個考量,第一,既然是使用平均的方法來抑制雜訊,那麼, 理論上當然是平均越多個符元效果越好,但第二個考量是:如果平均太多個符 元,在這些符元時間內,通道可能已經改變了,反而造成估計誤差。因此,必需 有一個協調,在模擬中我們將使用2 個符元的長度來平均。 圖4.6 應用於 DVB-T 系統之基於快速傅利葉轉換內插通道估計法
4.2.3 通道路徑的選擇(Path Selection)
在這一節中我們將介紹”路徑選擇”的方法,路徑選擇的目地有兩個,第一個 是能有效的濾除雜訊,另一個則是降低運算的複雜度。它的方法描述如下: 在一 般的正交分頻多工系統中,我們往往假設通道的最大延遲會小於時間護衛的長 度,如此一來方可避免碼際干擾(Inter Symbol Interference)的效應,因此我們可以
進一步將通道中大於護衛間隔的值都視為是雜訊並令之為零,如圖4.7 所示。 圖4.7 通道之最大延遲小於等於護衛間隔 但在4.2 節中所提到的基於快速傅立葉轉換內插法中,估出的通道在時域是重複 出現的,如同圖4.8 所示: ˆ | S [ ] | LS h n 圖4.8 時域上通道脈衝響應重複的出現 因此,在做路徑選擇時,第一步我們必需先將t>N T 的點其值設為0,儘管如此, 通道中t<N T 的值還是有可能包含雜訊,因此,我們可以用一個簡單的方法來判 別,首先,在t<N T 的點內找出一最大能量值的路徑,利用這最大能量值我們設
訊,反之則視為通道,因此,式(4.20)可以改寫成
其中L 為選中為 Path 的總數,l則為選中Path 的 index。因此將通道路徑的選擇
應用在式(4.19)中,我們只要計算選中為 Path 次載波位置的加權值,如此即可簡 化計算: 4. 2.4 通道延遲對於估計效能的影響 在 DVB-T 的規範中,2048 點的次載波只使用了 1705 點來傳送資料(包含 Pilot),如圖 4.9(a),但是,當我們在做基於快速傅利葉轉換的內插法時,用的卻 是2048 點的 IFFT,這樣的效應等效是 ˆS LS H 乘上一個方波,經過2048 點的 IFFT 後等效是 ˆS LS h 摺積一個sinc 波,就如同圖 4.9(b)所示,這樣的效應造成兩個影響, 第一是振幅衰減成1705 2048,第二是波瓣(side-lobe)的干擾,這可以透過路徑的選擇 來消除。 (a) 0 1 1 1 1 0 2 2 2 2 2, 2 2,..., 2 2 ˆ ˆ ˆ L L T LS LS LS l l l l LS LS LS l z l z l z h h h k h h h
σ
σ
σ
− − ⎡ ⎤ ⎢ ⎥ = ⎢ + + + ⎥ ⎣ ⎦ 2 2 1 ˆ ˆ ( ) / z LS n l h n N T L σ ∉ = −∑
(4.5) (4.6)ˆ | S ( ) | LS H k |HˆLSS '( ) |k
×
ˆ | [ ] | LS S h n⊗
(b) 圖4.9 (a)DVB-T 的符元架構 (b)額外的干擾 圖 4.10 顯示了雙重路徑下,不同延遲所受的影響,當兩條路徑的延遲較小時, 波瓣(side-lobe)的干擾較為嚴重,當兩條路徑的延遲較大時(圖 4.11),波瓣 (side-lobe)的干擾可以忽略。 (a) 'ˆ
|
S|
LSh
=
(b)
圖4.10 (a)完整的通道估計值|hˆLS |(延遲為 50) (b)僅使用 1705 點估計出的通道
'
ˆ |hLS |
(b) 圖4.11 (a)完整的通道估計值|hˆLS |(延遲為 500) (b)僅使用 1705 點估計出的通道
4.3 逼零等化器[10]
由於 DVB-T 為 OFDM 系統,傳輸通道被切割成許多細小子通道,因此每 個次載波的頻寬將相對變窄,使得每個次載波的傳輸通道為平坦衰變(flat fading),也就是訊號只有振幅衰減,而不會有波形失真的情形。此時,在接收 端就只需要用一個簡單的頻域等化器即能將訊號的振幅失真調整回來,大大的降 低了系統運算的複雜度。而頻域等化器的主要功能就是補償每個子通道增益上的 損失。就每一個次載波而言,即為單閥等化器(one-tap equalizer)。 由於,DVB-T 接收機經傅立葉轉換之後得頻域訊號,所以我們以頻域觀點 來設計通道等化器設計,圖4.12 顯示一系統傳送訊號經通道至接收端簡圖。
Channel H X N Y=XH+N 圖4.12 傳送訊號經通道至接收端簡圖 X 表示為傳送端的傳送訊號,通道的頻域響應為 H,N 為可加性白色高斯雜訊 (AWGN )。則接收端所收到的訊號 Y 為傳送訊號乘上通道的頻域響應再加上 可加性白色高斯雜訊即。 Y=XH+N 通道等化器的主要目的是為了等化傳送訊號所受到通道影響的振幅失真。 所以,我們可經由反乘通道參數,而求得傳送訊號,即如(4.8)。 H Y X = 而通道參數可以藉由第二章所提的幾個方法估得,所以實際上的通道等化器工作 公式如(4.9)。 H Y X ˆ = 由(4.9)可看出,逼零等化器的缺點是忽略可加性雜訊的影響,因此,它會造 成重大的雜訊增強效果,當通道的頻域響應 H 或 Hˆ 很小,則逼零等化器的補償 帶來了大的增益,同時,雜訊也被增強了許多。 (4.7) (4.8) (4.9)
4.4 路徑干擾重建[11]
1 τ τ2 τ3 τL−1 τL 0 1 a 2 a L a 1 L a − 3 a 圖4.13 通道估計所估計出的時域通道特性 接收訊號y n 為發送訊號乘上通道的特性加上雜訊,如(4.10)[ ] : 1 1 2 2 1 1 [ ] [ ] [ ] ... L [ L ] L [ L] [ ] y n =a x t−τ +a x t−τ + +a −x t−τ − +a x t−τ +w n (4.10) 重建每個路徑的干擾: ] [ ˆ ˆ ] [ ˆ ˆ ... ] [ ˆ ˆ ] [ ˆ ˆ ] [ ] [ ˆ ˆ ] [ ˆ ˆ ... ] [ ˆ ˆ ] [ ˆ ˆ ] [ ] [ ˆ ˆ ] [ ˆ ˆ ... ] [ ˆ ˆ ] [ ˆ ˆ ] [ ] [ ˆ ˆ ] [ ˆ ˆ ... ] [ ˆ ˆ ] [ ˆ ˆ ] [ 1 1 2 2 2 2 1 1 2 2 2 2 1 1 1 1 1 3 3 2 2 2 1 1 3 3 2 2 1 − − − − − − − − − − − − + − + + − + − = − + − + + − + − = − + − + + − + − = − + − + + − + − = L L L L L L L L L L L L L L L L L L t x a t x a t x a t x a n i t x a t x a t x a t x a n i t x a t x a t x a t x a n i t x a t x a t x a t x a n i τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ # (4.11) 將接收訊號分別扣掉每個路徑的干擾,則可得到每條路徑的傳送資料,且只受到 雜訊的影響而沒有其他路徑的影響: ] [ ˆ ] [ ] [ ] [ ] [ ˆ ] [ ] [ ] [ ] [ ˆ ] [ ] [ ] [ 2 2 2 2 1 1 1 1 n w t x a n i n y n w t x a n i n y + − = − + − = − + − = − τ τ τ # (4.12)其中 ( ) ( ) 1 2 2 2 2 ˆ[ ] [ ] ˆ ˆ[ ] ... L [ L] ˆ ˆL [ L] [ ] w n = a x t−τ −a x t−τ + + a x t−τ −a x t−τ +w n
(
) (
)
(
)
2 1 1 1 1 3 3 3 3 ˆ [ ] [ ] ˆ ˆ[ ] [ ] ˆ ˆ[ ] ... ˆ ˆ [ ] [ ] [ ] L L L L w n a x t a x t a x t a x t a x t a x t w n τ τ τ τ τ τ = − − − + − − − + + − − − + (4.13) #(
1 1 1 1)
(
1 1 1 1)
ˆL[ ] [ ] ˆ ˆ[ ] ... L [ L ] ˆ ˆL [ L ] [ ] w n = a x t−τ −a x t−τ + + a −x t−τ − −a −x t−τ − +w n 在得到每條路徑的傳送資料後,經過時間延遲調整以及最大比例匹配的增 益,可以得到估計的接收訊號 ˆ[ ]y n ,此估計的接收訊號已經消除了多路徑的相互 干擾,同樣的通過去除護衛區間、快速傅立葉轉換及等化器後得到消除了多路徑 相互干擾的估計傳送訊號D n[ ],經過反快速傅立葉轉換再加上護衛區間得到發 送訊號x n[ ]。 估計出發送訊號x n[ ]後,可代回接收端架構中取代原本的估出的發送訊號 ˆ[ ] x n ,藉此重建出更為準確的每個路徑的干擾,繼而接收訊號再一次扣掉每個路 徑的干擾而得到每條路徑的傳送資料,最後又估出更為準確的發送訊號,不斷的 反覆後估計得以完全的消除多路徑之間彼此的干擾。 想要得到更好的效能表現,可以先將估計得到的傳送資料( ˆD n or [ ][ ] D n ) 先經過解碼及編碼的動作,可以有編碼的增益(Coding gain),得到修正的傳送 資料( ˆD n′[ ]orD n′[ ]),如下圖所示,再代回上面的接收端架構重建資料的機制。ˆ[ ]
D n
D n
ˆ [ ]
′
D n[ ] D n′[ ] 圖4.14 傳送資料經過編碼的機制4.5 時間位移與最大比例合併[12]
當接收端收到 L 個分集訊號,而每個分集訊號的強度為r ,則根據 MRC,l 只要將各個分集訊號的相位調成相同,而各個分集再乘上與訊號強度成正比的增 益,再把這 L 個分集相加在一起來做資料決策,則可得到最大的 SNR 比。我們可 以把MRC 的觀念應用在接收端的架構上,因為傳送使用到展頻的方式,我們可 以把 L 個切片看成是 L 個分集訊號,而這 L 個分集訊號為複數的形式。如果我們 使用MRC 的通道匹配方式來做決策,則須一方面把每個切片訊號的相位調成相 同,另一方面再乘上與每個切片訊號的振幅成正比的增益,最直觀的做法就是在 每個切片乘上 *( ) k i H 而後再將這 L 個切片相加,即 ( ) ( ) *( ) 1 ( ) ( ) *( ) 1 2 ( ) ( ) 1 (1) (1) (2) (2) [ ] [ ] [ ] [ ] [ ] {
L k k k i c m m i m L k k k M M c i i i m m m k k M M c i i m m M m c i m i m L m D m E c s H m E e H m H m I N E e H m I N where I E e c e c e = = = ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ + + ⋅ ⋅ + + = ⋅ ⋅ + ⋅ + +
=
=
∑
∑
∑
" ( 1) ( 1) ( 1) ( 1) ( ) ( ) ( ) ( ) *( ) (1) (1) (2) (2) ( 1) ( 1) ( 1) ( 1) 2 ( ) ( ) ( ) ( ) } [ ] [ ] { } [ ] k k k k i m i m K K k k k i m i m i k k k k c i m i m i m i m K K k k i m i m M m c c e c e c H m c H m E e c e c e c e c e c H m c N E n − − + + − − + + ⋅ + ⋅ + + ⋅ ⋅ ⋅ ⋅ = ⋅ ⋅ + ⋅ + + ⋅ + ⋅ + + ⋅ ⋅ ⋅ = ⋅ " " " ( ) *( ) [ ] k k [ ] m m i m c⋅ ⋅ ⋅s H m (4.14)第五章 電腦模擬與討論
5.1 模擬環境及參數
Number of Sub-carriers 2048 Number of Active Sub-carriers 1705
Sampling Interval 7/64µ sec Useful Symbol Duration TU = 224µsec
Sub-carrier Frequency Spacing 4464Hz Effective Bandwidth 7.61Mhz
Guard Interval 1/4 TU = 56µsec
Maximum Delay 500TS =500*7/64µsec
Coherent Bandwidth 4.816 carrier spacing Carrier Frequency
f
C 474MHzMaximum Mobile Velocity 300Km/hr Maximum Doppler Frequency fd ~=131.67Hz
Total Pilot Power 176×(16/9)
Modulation QPSK
Simulation Time 650 DVB-T Frames
表5-1 模擬的環境與參數
我們模擬的環境如同表 5-1 所示,用的是 DVB-T 標準中的 2K 模式,護衛
5.2 模擬的結果及討論
5.2.1 通道等化器 without MPIC
5.2.1.1 AWGN 通道
5.2.1.1 行動接收通道
圖5.2 逼零等化器在行動接收通道下的 BER(車速 30km/hr)
圖5.4 逼零等化器在兩條路徑下行動接收通道下的 BER
(車速30km/hr and 90km/hr)
圖5.5 逼零等化器在三條路徑下行動接收通道下的 BER
5.2.2 MPIC
(a)
(b)
(a)
(b)
圖5.7 AWGN 通道使用相同路徑延遲下的 BER(a)完美通道估計(b)通道估
(a)
(b)
圖5.8 行動接收通道相同能量比例使用完美通道估計下的 BER(a)車速 30 公
(a)
(b)
圖5.9 行動接收通道相同路徑延遲使用完美通道估計下的 BER(a)車速 30 公
(a)
(b)
圖5.10 行動接收通道相同能量比例使用通道估計下的 BER(a)車速 30 公里
(a)
(b)
圖5.11 行動接收通道相同路徑延遲使用通道估計下的 BER(a)車速 30 公里(b)
第六章 結論與未來方向
近年來,電子科技的快速發展為人類的文明與生活帶來進步,科技日新月 異,電子系統的設計也越多樣化,如何把多種系統整合到同一系統架構,已成為 未來電子業發展的趨勢,當相似標準及架構的系統可以整合起來,代表電子生產 的成本可以降低,產品的價值與效能可以提高,在市場上將可以更有競爭力。此 份報告介紹接收端等化器設計的方法,設計出其相對應之方塊圖及數學演算法, 並且以電腦模擬的方式作系統架構的驗證。利用通道估計估出通道的特性,其配 合等化器及決策器的運作可估計出傳送資料,此傳送資料可加上估計出的通道效 應而構成估計出的接收訊號,其可用來區分不同路徑的訊號及排除多路徑干擾的 效應,系統再以最大比例匹配(Maximum ratio combining)重建出較為準確的接 收訊號,如此反覆進行得到更準確的傳送訊號。在未來繼續加入編碼(華士碼)使得資料載波具有相關性,再利用多重路徑 干擾消除扣掉多重路徑的影響恢復其正交性,亦藉此提升多重路徑干擾消除的成 效,可得到更佳的成效。
參考文獻
[1] Richard van Nee, “OFDM for Wireless Multimedia Communications”, Artech House, 2000
[2] S. Weinstein and P. Ebert, “Data Transmission by Frequency-Division
Multiplexing Using the Discrete Fourier Transform”, Communications, IEEE
Transactions on [legacy, pre - 1988] , Volume: 19 Issue: 5, pp.628 –634, October 1971
[3] Draft EN 300 744 v1.4.1(2001-01). “Digital Video Broadcasting (DVB);
Framing structure, channel coding and modulation for digital terrestrial television”.
[4] Theodore S. Rappaport, “Wireless communications: principles and practice”, Prentice Hall PTR, Inc. 2002
[5] Louis L. Scharf, Statistical Signal Processing, Addison-Wesley,1991
[6] Rinne, J. and Renfors, M., “Pilot spacing in orthogonal frequency division
multiplexing systems on practical channels
”,
IEEE Transactions on , Volume:42 , Issue: 4 , Nov. 1996 Pages:959 – 962
[7] Ye Li, “Pilot-symbol-aided channel estimation for OFDM in wireless systems.”
IEEE Transactions on Vehicular Technology, Volume: 49 Issue: 4 , July 2000
[8] Che-Shen Yeh and Yinyi Lin, “Channel Estimation Using Pilot Tones in OFDM
Systems” IEEE Transactions on Broadcasting, Vol. 45, No. 4, December 1999.
[9] Hiaing Minn, Vijay K. Bhargava, “An Investigation into Time-Domain Approach
for OFDM Channel Estimation” IEEE Transactions on Broadcasting, Vol. 46,
No 4, December 2000.
[11] Kenichi Higuchi, Akihiro Fujiwara, and Mamoru Sawahashi, “Multipath
Interference Canceller for High-Speed Packet Transmission With Adaptive
Modulation and Coding Scheme in W-CDMA Forward Link”,IEEE JOURNAL
ON SELECTED AREAS IN COMMUNICATIONS, VOL. 20, NO. 2, FEBRUARY 2002
[12] L. Hanzo, T.H. Liew, and B.L. Yeap, “Turbo Coding, Turbo Equalisation and
Space-Time Coding for Transmission over Fading Channels” IEEE Press, pp.