基于MATLAB的氣象傳真信號的FM解調(diào)算法

李永鋒+寇學鋒+穆菲菲+楊啟超

摘 要 氣象傳真接收機是船舶在航行中獲取水文氣象資料的主要手段?，F(xiàn)在氣象傳真接收機大多采用模擬電路實現(xiàn)，導致接收機數(shù)字化程度不高，不便于氣象信息的有效利用。本文針對現(xiàn)在氣象傳真接收機問題，在基帶信號處理算法基礎上，結(jié)合了氣象傳真信號自身的特征，提出了FM的解調(diào)方案以及對應算法，并對解調(diào)算法在MATLAB平臺上做仿真和驗證，最后驗證了算法的科學性和可行性。

【關鍵詞】FM解調(diào) 氣象傳真信號 MATLAB

1 緒論

現(xiàn)階段通信技術出現(xiàn)了數(shù)字調(diào)制方式取代模擬調(diào)制方式的趨勢。但調(diào)頻調(diào)制方式因其具有抗噪聲性能好、易于實現(xiàn)的優(yōu)點而被廣泛應用于無線通信的各個領域。對于無線通信系統(tǒng)中，使用調(diào)頻調(diào)制方式時，調(diào)頻方式能克服電平波動的影響，因此無線氣象傳真機大都采用調(diào)頻方式進行調(diào)制。解調(diào)模塊同樣是接收機最核心的部分，解調(diào)方法及其實現(xiàn)過程直接影響解調(diào)效果，而解調(diào)效果又直接影響接收機的靈敏度和抗干擾等指標。

目前國內(nèi)主流氣象傳真接收機大多采用基于DSP的非相干解調(diào)算法，在實現(xiàn)上比較簡單，但其解調(diào)效果不理想，尤其是在信噪比較高時，其解調(diào)效果要比相干解調(diào)差。FM相干解調(diào)適用于窄帶調(diào)頻信號，抗干擾能力強。另外，F(xiàn)M解調(diào)器可采用模擬電路實現(xiàn)，也可采用數(shù)字信號處理方式實現(xiàn)。

本文采用了FM相干解調(diào)方法，結(jié)合數(shù)字正交解調(diào)原理對氣象傳真信號進行處理，并在Matlab平臺上實現(xiàn)解調(diào)算法，然后驗證算法的可行性以及觀察算法的解調(diào)效果。

2 氣象傳真信號特征及格式

氣象傳真通信由發(fā)送端和接收端構(gòu)成，首先發(fā)送端將一張圖像完整地送到接收端。在發(fā)送端，將要發(fā)送圖片通過光電轉(zhuǎn)換為相應強度的脈沖信號，然后將其調(diào)制成頻率相對中心頻率為1900Hz 400Hz的傳真信號，此信號調(diào)制后發(fā)送到接收端，接收端經(jīng)過相應的解調(diào)過程，最后恢復氣象傳真圖。傳真圖逐行掃描發(fā)送和接收，在合作指數(shù)576時，氣象傳真信號的發(fā)送和接收格式如圖1所示。

由圖1可知，整個氣象傳真圖由5部分構(gòu)成。當掃描速度采用120行/min時，即每行圖像信號發(fā)送時間為0.5s，接收端接收速度也為120行/min。表1為氣象傳真圖中各組成信號的說明。

3 解調(diào)方案及方法

本文采用直接調(diào)頻FM相干解調(diào)方法，結(jié)合數(shù)字正交分解技術，調(diào)頻信號的FM正交解調(diào)流程圖如圖2所示。

氣象傳真是短波信號經(jīng)過下變頻后得到的窄帶信號，其已調(diào)信號數(shù)學表達式：

要得到解調(diào)信號m（n）只需去除直流分量?ω即可。從整個解調(diào)算法的推導可知，該算法只有乘法和加減運算，便于實現(xiàn)和應用。

氣象傳真接收機通過天線接收氣象站臺發(fā)送的短波信號，短波信號頻率范圍為12～25MHz，通過下變頻處理變成基帶信號，下變頻用跟蹤環(huán)路來實現(xiàn)；基帶信號通過帶通濾波器（頻率范圍1.3～2.5kHz）處理，除掉頻段外不需要的信號，然后通過與本地載波信號相乘就得到解調(diào)后的信號；信號再通過低通濾波器濾除信號中的低頻，然后通過數(shù)據(jù)格式變換成圖像信號；再對圖像信號做打印前預處理，調(diào)整圖像顯示格式，最后對圖像中的噪聲信號做噪聲濾除處理。本設計在MATLAB中的調(diào)頻信號解調(diào)方案流程圖如圖3所示。

4 仿真與驗證

在Matlab平臺上對氣象傳真信號驗證FM正交解調(diào)算法的可行性和解調(diào)效果。仿真實驗輸入信號來自氣象傳真接收機接收到一段真實氣象信號，該信號經(jīng)過下變頻后得到基帶范圍內(nèi)的傳真信號。然后轉(zhuǎn)換為3min時長的音頻信號，通過MATLAB函數(shù)讀入這段音頻信號，實際氣象傳真信號數(shù)據(jù)波形如圖4所示。

由圖4可知實際傳真信號由有效傳真信號和噪聲信號混合而成。下變頻后的數(shù)據(jù)再經(jīng)過帶通濾波器（其上、下截止頻率分別為1.3kHz、2.5kHz）可得到去掉帶外信號的數(shù)據(jù)波形如圖5所示。

然后是對信號做正交解調(diào)去載波處理，解調(diào)出來的數(shù)據(jù)波形如圖6所示。

對解調(diào)出的傳真信號進行二值化圖像預處理。要將整個圖像呈現(xiàn)出明顯的黑白效果，需要采用不同的閾值對數(shù)據(jù)做門限判決。本文采用了30、65、95、120、150、175、200、227和255這9個數(shù)值對數(shù)據(jù)進行二值化調(diào)理后的數(shù)據(jù)如圖7所示。

對二值化圖像數(shù)據(jù)進行合并壓縮和結(jié)構(gòu)變形等處理后繪制出傳真圖。此處讀取音頻數(shù)據(jù)時長是3min，最終繪制出的氣象圖如圖8所示。

從最終解調(diào)的氣象圖可看出，圖中數(shù)字、字母和經(jīng)緯線等有用信息都非常清晰，但輸入信號本身含有噪聲，故解調(diào)出的氣象傳真圖也帶有噪聲。而去除氣象圖噪聲用圖片去噪方法即可，本文就不在贅述。

5 結(jié)論

本文對氣象傳真信號特征和格式做了詳細介紹，并且對FM正交解調(diào)算法進行了詳細推導。由MATLAB仿真結(jié)果知，F(xiàn)M解調(diào)算法能很好恢復出調(diào)制信號，解調(diào)出的氣象傳真圖能夠滿足接收機的性能指標要求。

參考文獻

[1]許開寧，曾連蓀.氣象傳真圖的實時接收[J].氣象與環(huán)境學報，1997，8（01）：30-35.

[2]顧兵.無線氣象傳真接收機數(shù)字解調(diào)系統(tǒng)設計[D].哈爾濱工程大學，2006.

[3]薛年喜.MATLAB在數(shù)字信號處理中的應用[M].北京：清華大學出版社，2003.

作者單位

陜西寶成航空儀表有限責任公司 陜西省寶雞市 721006