基于角度計(jì)算轉(zhuǎn)換的瞬時(shí)相位實(shí)時(shí)生成方法

張 強(qiáng) ，齊 明

（1. 航天科工防御技術(shù)研究試驗(yàn)中心 北京 100854；2. 中國人民解放軍駐二三二廠軍事代表室 北京 100086）

基于角度計(jì)算轉(zhuǎn)換的瞬時(shí)相位實(shí)時(shí)生成方法

張 強(qiáng)1，齊 明2

（1. 航天科工防御技術(shù)研究試驗(yàn)中心 北京 100854；2. 中國人民解放軍駐二三二廠軍事代表室 北京 100086）

針對高數(shù)據(jù)率復(fù)信號的瞬時(shí)相位難以實(shí)時(shí)求解問題，提出了基于角度計(jì)算轉(zhuǎn)換的求解方法。該方法先把瞬時(shí)相位的反正切運(yùn)算轉(zhuǎn)換成反正/余弦運(yùn)算，再利用查找表代替具體的運(yùn)算。證明了在角度計(jì)算轉(zhuǎn)換前后，復(fù)信號I/Q兩路通道噪聲對瞬時(shí)相位求解結(jié)果影響的等效性。微波著陸系統(tǒng)(Microwave Landing System)實(shí)測數(shù)據(jù)處理結(jié)果表明，采用該方法可以在不明顯降低瞬時(shí)相位求解精度的前提下，有效地降低計(jì)算量，并可完成實(shí)時(shí)處理。

復(fù)信號；瞬時(shí)相位；角度計(jì)算轉(zhuǎn)換；查找表

當(dāng)今通信、導(dǎo)航和雷達(dá)系統(tǒng)的發(fā)展有兩個(gè)明顯的特點(diǎn)。一方面信號的時(shí)寬帶寬積逐漸增大；另一方面系統(tǒng)的ADC(Analog to Digital Conversion) 越來越靠近射頻前端。其主要原因是，增大時(shí)寬帶寬積可以有效提高通信系統(tǒng)的信道容量、提高導(dǎo)航系統(tǒng)的導(dǎo)航定位精度和提高雷達(dá)系統(tǒng)的作用距離及分辨率；而ADC靠近射頻前端可以有效提高系統(tǒng)的信號處理靈活性（例如軟件無線電技術(shù)[1]）。上述兩個(gè)特點(diǎn)都會造成系統(tǒng)信號數(shù)據(jù)率的明顯增高。以FPGA (Field Programmable Gate Array)和DSP (Digital Signal Processor)為代表的微電子技術(shù)的迅速發(fā)展在很大程度上減輕了高速數(shù)字信號實(shí)時(shí)處理的壓力，但其仍難以完成一些復(fù)雜的信號處理[2]。為此，在不明顯降低運(yùn)算精度的前提下研究實(shí)時(shí)高效的信號處理方法在當(dāng)前仍然具有很重要的意義?，F(xiàn)有的一些通信、導(dǎo)航和雷達(dá)系統(tǒng)在對接收信號做中頻數(shù)字化處理之后，往往還需要進(jìn)一步求解復(fù)信號的瞬時(shí)幅度和瞬時(shí)相位[3-5]。其中瞬時(shí)相位的求解涉及到反正切運(yùn)算，因而在高數(shù)據(jù)率條件下難以實(shí)時(shí)實(shí)現(xiàn)。這也給后續(xù)的信號處理帶來很大的限制。為了克服該問題，本文提出了一種基于角度計(jì)算轉(zhuǎn)換的瞬時(shí)相位獲取方法。該方法先通過角度計(jì)算轉(zhuǎn)換將求解瞬時(shí)相位的反正切運(yùn)算轉(zhuǎn)換成反正弦運(yùn)算和反余弦運(yùn)算，接著根據(jù)瞬時(shí)相位的求解精度需求設(shè)計(jì)相應(yīng)的反正弦表和反余弦表，最后通過查找表的方式實(shí)現(xiàn)反正弦和反余弦運(yùn)算，即通過查找表的方式實(shí)現(xiàn)瞬時(shí)相位的獲取。具體內(nèi)容請見下文。

1 瞬時(shí)相位的實(shí)時(shí)獲取

若記復(fù)信號的同相分量(In-Phase)和正交分量(Quadrature-Phase)分別為xI和xQ，則其瞬時(shí)相位可定義為：

顯然上式需要求解反正切函數(shù)，因此計(jì)算量較大。對于數(shù)據(jù)率很高的場合，其實(shí)時(shí)處理的壓力很大。為此，下面采用簡化的求解方法，具體實(shí)現(xiàn)步驟如下：

1 ) 角度計(jì)算的轉(zhuǎn)換

為了能夠利用查找表來求解瞬時(shí)相位，首先需要將式(1)的反正切函數(shù)求解變?yōu)榉从嘞液瘮?shù)求解，即：

命題：在角度計(jì)算轉(zhuǎn)換前后，由通道噪聲和量化誤差造成的瞬時(shí)相位求解誤差是等效的。

證明：由式(1)得到的瞬時(shí)相位計(jì)算誤差為：

其中，δ xI和δ xQ表示由通道噪聲和量化引起的 xI和xQ的測量誤差。并且

同樣，由式（4）得到的瞬時(shí)相位計(jì)算誤差為：

因此，上述兩種角度求解方法具有相同的瞬時(shí)相位估計(jì)誤差 δφ ，證畢。

同樣地，也可以將反正切函數(shù)的求解變?yōu)榍蠓凑液瘮?shù)的求解。

同理可以證明，上式給出的求解方法具有相同的角度求解誤差，因此可以認(rèn)為該角度求解方法的轉(zhuǎn)換是等效的。

2 ) 縮小角度查找范圍

在利用查表的方法求解瞬時(shí)相位時(shí)，縮小角度的查找范圍是很有必要的，縮小查找范圍意味著在相同的求解精度下，可以減小查找表的存儲空間，提高查找效率。已知由IQ數(shù)據(jù)得到的瞬時(shí)相位的取值范圍是00～3600，通過象限變換，其它三個(gè)象限的角度值可通過求解第一象限的角度值間接得到，所以為了縮小查表的范圍，只需對求解第一個(gè)象限的瞬時(shí)相位。現(xiàn)設(shè)定中間變量p，使得

顯然， 2-0.5≤p≤1。為了進(jìn)一步縮小查表的范圍，可采取如下的瞬時(shí)相位求解方式：

不難看出上述3種方法的估計(jì)誤差是等效的。因此，第一象限的瞬時(shí)相位求解可以轉(zhuǎn)換成00～450的反余弦函數(shù)求解。這樣，所欲求解的角度范圍便可以縮小一半。

當(dāng)然，也可以采用求反正弦函數(shù)的方式實(shí)現(xiàn)瞬時(shí)相位的求解，先設(shè)定中間變量：

顯然，0≤q≤2-0.5。由此可得到瞬時(shí)相位的求解：

3 )查找表的實(shí)現(xiàn)

和正弦表和余弦表的查找方式不同，為了保證在瞬時(shí)相位的求解過程中，查找反正弦表和反余弦表所得到的角度值在整個(gè)查找范圍內(nèi)是均勻取值的。用于查找表的中間變量的取值就不能是均勻的，由此造成了查找反正弦表和反余弦表的困難。

首先考慮查找反余弦表的情況，用于查找表的中間變量 p的取值范圍為： 2-0.5≤p≤1，讓 p以0.001的步進(jìn)遍歷其取值范圍，則總共有294個(gè)數(shù)據(jù)。因?yàn)閜的均勻取值會造成 θ取值的不均勻，所以為了保證θ的精度，需要在p接近1的時(shí)候減小p取值的步進(jìn)，或者是在整個(gè)遍歷范圍內(nèi)減小p取值的步進(jìn)。顯然，前一種方法不利于查找，而后一種方法需要增大取值的數(shù)據(jù)點(diǎn)個(gè)數(shù)。為了克服查找反余弦表過程中存在的缺點(diǎn)，可以采用“二次查表”法作彌補(bǔ)，主要思想是設(shè)置反余弦表和余弦表兩張查找表，首先保持p的遍歷步進(jìn)不變，獲得了一張反余弦表，將反余弦表的內(nèi)容由原來的角度值替換成下一次需要在第二張查找表中查找的起始位置和查找范圍，根據(jù)該信息，中間變量p可以很方便地在接下來的余弦表中找到與自己最為接近的余弦值，該余弦值所對應(yīng)的角度值即為瞬時(shí)相位的估算結(jié)果。

據(jù)此可見，由p的量化誤差引起的φ值估計(jì)誤差為：

設(shè)q以步進(jìn)為0.001遍歷其取值范圍，則反正弦查找表中共有708個(gè)數(shù)據(jù)。當(dāng)q的取值接近2-0.5時(shí)，此時(shí)由量化引起的φ值估計(jì)誤差最大：

4 ) 象限恢復(fù)

求得第一象限的φ 后，將其取值擴(kuò)展到4個(gè)象限：

至此，有關(guān)瞬時(shí)相位的簡化求解算法已經(jīng)介紹完畢?？偟目磥恚鋵?shí)現(xiàn)步驟可以概括如下：

Step2.將所求得的瞬時(shí)幅度代入式，求得中間變量q；

Step3.根據(jù)q的取值，在00～450的反正弦表中查得φ值。因?yàn)閝的查找范圍是均勻量化的，所以可以很方便地實(shí)現(xiàn)查找；

Step4.將φ值的取值擴(kuò)展到4個(gè)象限，得到所要求解的瞬時(shí)相位。

2 基于實(shí)測數(shù)據(jù)的求解精度驗(yàn)證

為了進(jìn)一步驗(yàn)證本文所提方法的性能，下面結(jié)合微波著陸系統(tǒng)(Microwave Landing System, MLS)實(shí)測數(shù)據(jù)開展實(shí)驗(yàn)。

在實(shí)驗(yàn)中，利用MLS內(nèi)場模擬器產(chǎn)生所需的信號（如圖1所示），采用NI (National Instrument)公司的PXI (PCI eXtensions for Instrumentation)總線[6]對MLS信號進(jìn)行中頻數(shù)字化處理，得到的I/Q通道的信號分別如圖2和圖3所示。

由I/Q通道的兩路信號也可以得到未經(jīng)簡化的瞬時(shí)相位，其結(jié)果如圖4所示。本文所采用的近似求解算法的求解結(jié)果如圖5所示，由近似帶來的誤差如圖6所示。由此可見，由近似帶來的誤差最大不超過0.08deg，該精度是能夠滿足MLS機(jī)載接收機(jī)信號解調(diào)的要求。這也主要得益于本文的角度計(jì)算轉(zhuǎn)換方法和查找表設(shè)計(jì)方法。

圖1 MLS的發(fā)射信號Fig.1 The transmit signal of MLS

圖2 MLS信號中頻數(shù)字化之后的I路信號Fig.2 The in-phase component of MLS signal after intermediate frequency digitalization

圖3 MLS信號中頻數(shù)字化之后的Q路信號Fig.3 The quadrature-phase component of MLS signal after intermediate frequency digitalization

圖4 MLS信號中頻數(shù)字化之后的瞬時(shí)相位（未作近似）Fig.4 The instant-phase of MLS signal after intermediate frequency digitalization without any approximation

圖5 瞬時(shí)相位近似求解算法的求解結(jié)果Fig.5 The instant-phase derived by the approximated calculation algorithm

圖6 簡化瞬時(shí)相位算法的求解誤差Fig.6 The calculation error for the simplified instant-phase calculation algorithm

3 結(jié) 論

基于角度計(jì)算轉(zhuǎn)換和查找表方法的瞬時(shí)相位簡化求解算法可以在近似誤差的允許范圍內(nèi)，有效地減小計(jì)算量，因此其適合于高數(shù)據(jù)率條件下信號的實(shí)時(shí)處理。該簡化求解算法也可以推廣到其他的中頻數(shù)字化瞬時(shí)相位求解。

[1]楊小牛，樓才義，徐建良. 軟件無線電原理與應(yīng)用[M]. 北京：電子工業(yè)出版社, 2002.

[2]Kumardeb B, Bivas D, Kalyan M. Novel FPGA-based LVDT signal conditioner[C] //2013 IEEE International Symposium on Industrial Electronics (ISIE),2013:1-6.

[3]Landis G P, Galysh I, Petsopoulos T. New Digital Phase Measurement System[C] //Naval Research Lab., Washington,DC. 2009.

[4]Han D,Yuanjin Zheng.A GFSK demodulator based on instant phase computation and adaptive multi-threshold quantization[C]//IEEE Asian Solid-State Circuits Conference,2009: 249-252.

[5]Vyroubal D.Optical method for instant estimate of vibration signature based on spectrum analysis of phase-modulated light pulses[J].IEEE Transactions on Instrumentation and Measurement,2004,53(1):181-185.

[6]Gutterman, Loofie. PXI for military test applications[C] // IEEE Systems Readines Technology Conference,2002:215-226.

A real-time instant phase generating method based on angle calculation transforming

ZHANG Qiang1, QI Ming2
（ 1. Aerospace Science and Industry Defense Technology Research and Test Center, Beijing 100854, China;2. PLA Military Representative Office in No.232 Factory, Beijing 100086, China）

Because of the big computational burden of arc tangent with high data rate, the real-time calculation of instant phase can be a problem for complex signal processing. To solve this problem, a new method based on angle calculation transforming is proposed in this paper. The arc tangent computation is firstly transformed to the arc sine and/or arc cosine computation, and then the look-up tables can be utilized to replace the computation of angles. The influences of I/Q channel noises on the calculation results before and after transforming are proved to be equivalent to each other. The experiment results of Microwave Landing System (MLS) are given that, with the calculation precision not degraded obviously, the method proposed can effectively reduced the computational burden of instant phase, and thus, is very suitable to real-time signal processing.

complex signals; instant phase; angle calculation transforming; look-up table

TN967.1

A

1674-6236(2014)14-0043-03

2013-09-23 稿件編號：201309165

張 強(qiáng)（1983—），男，安徽合肥人，碩士，工程師，技術(shù)研究員。研究方向：無線電導(dǎo)航檢測與校驗(yàn)、系統(tǒng)可靠性設(shè)計(jì)。

book=0,ebook=429