高精度相位小存储量查表计算方法
【技术领域】
[0001] 本发明涉及雷达、通信和导航定位等包含大量相位计算的领域，特别是涉及利用 查找表进行相位计算的方法。
【背景技术】
[0002] 计算空间坐标点对应的相位值是常见的数字信号处理过程，如雷达探测中的角度 计算，通信中的调制星座点的相位计算和导航定位中的位置角度计算等。为了简化相位计 算过程，特别是简化在硬件（例如直接数字频率合成器）中实现相位计算的过程，通常采用 的相位计算方法是查表法。查表法将目标在直角坐标系中的坐标分别量化，并建立一个坐 标的量化数值到对应相位值的映射表，在进行相位计算时，直接用坐标的量化数值查相位 映射表就可以获得对应相位值。这种方法无需额外的计算过程，只需要一个存储器来存储 相位映射表，因此计算相位快捷。
[0003] 查表法的缺点在于，如果相位计算的精度要求较高，那么要求坐标的量化位宽较 宽，对应的相位映射表将会非常庞大，例如横坐标和纵坐标值量化位宽各为16比特，那么 完整的相位映射表的大小将达到2 32个存储单元，每个存储单元存储一个相位值，若相位值 采用16比特表述，那么存储这个相位映射表的存储器大小需要达到64M比特。显然，相位 映射表的大小与横坐标和纵坐标值的量化位宽的扩大呈4的幂次增长。对相位映射表缩减 的方法通常采用限制象限的方法，即将横坐标和纵坐标值都取绝对值，将相位映射表局限 在直角坐标系的一象限，计算相位时首先通过相位映射表获得在一象限的对应相位，而后 通过横坐标和纵坐标值的正负号将一象限的对应相位恢复到正确的象限。限制象限的相位 映射表缩减方法能够将相位映射表的大小缩减为完整相位映射表的四分之一，但对于高精 度的相位计算，仍然不能够解决相位映射表大小的膨胀问题。
【发明内容】
[0004] 本发明要解决的技术问题是，提供一种高精度相位小存储量查表计算方法，能够 在与现有查找表方法相位映射表存储量相同的情况下大幅提高相位计算精度，或者在保持 相同相位计算精度的情况下大幅度减少相位映射表规模，或者同时大幅提高相位计算精度 和大幅度减少相位映射表规模。
[0005] 本发明的技术方案是：一种高精度相位小存储量查表计算方法，其特征在于，包括 下述过程：
[0006] 第一步，预处理过程：
[0007] 已知相位精度△ Φ和相位映射表大小M，利用下式计算目标点坐标值的量化位宽 k ：
[0008] k = 0· 5 · Iog2M 式⑴
[0009] 构造一个包含M个存储单元的相位映射表，每个存储单元的比特位宽为B，地址为 IX2k+Q的存储单元存储的值为0，P即坐标为（I，Q)的目标点对应的相位值，其中：横坐 标I和纵坐标Q的值是对范围[0, 1]内的实数进行k比特均匀量化后的任意值，并且I多Q， 识=⑴7g/e ( /, ρ )，（函数angle (I，Q)表示横坐标为I和纵坐标为Q的目标点相位），可知 识G [0, ji/4]，识的量化位宽B的取值利用下式计算：
[0010] B = Iog2InKA · Δφ)\+\ 式（2)
[0011] 利用下式计算迭代次数值η :
[0012] n=cei/{ 5/4-1/2 · log2(2k · Δφ) } 式（3)
[0013] 上式中，函数ceil {χ}表示取不小于χ的最小整数。
[0014] 第二步，坐标迭代过程：
[0015] 设待计算的目标点坐标为（Itl, Qtl),利用下述步骤计算该目标点的查表坐标：
[0016] 令迭代次数 i = 1，Ii= I。，Qi= Q。；
[0017] 步骤①：记扒 mg/tM//，〇. }，IjP Q 满绝对值 abs (I i)和 abs (Qi)，IjP Q 满 正
[0018] 负号 Sw1和 Swq;
[0019] 步骤②：若 abs (Ii) = abs (Qi)，则 Ai= abs (I D，Bi= abs (Q D ;否贝丨」，令 Ai = max {abs (Ii)，
[0020] abs (Qi) }，Bi= min {abs (I J，abs (Qi) }，；
[0021] 步骤③：记沪 ',·= tmg/cK β,·}，
[0022] 若当前迭代次数i〈n，则利用下式计算I i+1和Q i+1:
[0023] /, I =cos(4^o ))=2 {2A,2-1}2-1 ； 式（4)
[0024] Qi I =^?η(4φ ))=4A ,5,· {2A ,2-1J ； 式（5)
[0025] 令Ii= I i+1，Qi= Q i+1，更新当前迭代次数i = i+1 ;并返回步骤②；
[0026] 若当前迭代次数 i = n，则 An= A i，Bn= B i，Ii= I n, Qi= Q n，Sw1= S ('，Swq = S(n)Q，并进入第三步查表过程；
[0027] 第三步，查表过程：
[0028] 利用在预处理过程构造的相位映射表中查坐标为（An，Bn)对应的相位值妒
[0029] 第四步，相位恢复过程：
[0030] 利用下述规则计算最终相位结果：
[0031 ]若 An= abs (I n)，Bn= abs (Q n)，则 V供 'η;
[0032] 若 Bn= abs (I n)，An= abs (Q η)，则
[0033] 其次：
[0034] 当Sw1为正且S (n)Q为正，则外
[0035] 当Sw1为正且S (n)Q为负，则
[0036] 当 Sw1 为负且 S (n)Q为正，则％=JI-V'》;
[0037] 当Sw1为负且S (n)Q为负，则厂Π ;
[0038] 令妒Vi=奶，按照上述规则，通过妒Vi获得仍-1，…，直至通过识'2获得奶，妁即最 终相位计算结果。
[0039] 本发明的有益效果是：本发明构造的相位映射表只包括第一象限的相位，并且相 位分辨率较低，因此有效地减少相位映射表的容量。通过迭代过程，能够不断恢复高精度相 位，实现了在小存储量相位映射表的情况下的高精度相位计算。
【附图说明】
[0040] 图1是本发明的高精度相位小存储量查表计算方法实现框图。
【具体实施方式】
[0041] 下面结合图1和一个实例来详细说明本发明的高精度相位小存储量查表计算方 法。
[0042] 图1中实现框图包含四个过程：预处理过程，坐标迭代过程，查表过程和相位恢复 过程。预处理过程利用相位精度参数和相位映射表大小M计算输迭代次数值η并构造 相位映射表；坐标迭代过程利用坐标I和Q，进行η次坐标迭代；查表过程接利用坐标迭代 的结果以及相位映射表获得迭代后相位；相位恢复过程利用迭代后相位恢复出最终相位计 算结果。
[0043] 实例如下：
[0044] 假设I、Q取值范围为[_1，+1]，16比特量化，相位计算精度要求为Λ Φ = 0. 1度， 相位映射表大小M = IK个存储单元。
[0045] 按照本发明的高精度相位小存储量查表计算方法，迭代的的次数为4次，实际得 到的相位计算精度为相位误差不大于〇. 0879度，满足相位计算精度要求为误差不大于0. 1 度的要求，同时满足相位映射表大小为IK个存储单元的要求。与现有的方法比较，按照普 通的直接查表法，达到0. 1度的相位精度，完整相位映射表的大小需要222个存储单元，即4Μ 个存储单元，无法满足实例中的相位映射表大小的要求。利用高精度相位小存储量查表计 算方法则能够同时满足相位精度和相位映射表大小的要求。
【主权项】
1. 一种高精度相位小存储量查表计算方法，其特征在于，包括下述过程： 第一步，预处理过程： 已知相位精度和相位映射表大小M，利用下式计算目标点坐标值的量化位宽k : k = 0· 5 · Iog2M (公式一) 构造一个包含M个存储单元的相位映射表，每个存储单元的比特位宽为B，地址为 IX2k+Q的存储单元存储的值为，供即坐标为（I，Q)的目标点对应的相位值，其中：横坐 标I和纵坐标Q的值是对范围[〇, 1]内的实数进行k比特均匀量化后的任意值，并且I多Q， 沪(/, 0)，函数angle(I，Q)表示横坐标为I和纵坐标为Q的目标点相位，识的量 化位宽B的取值利用下式计算： B = log2[s/(4 · zlf/9)]+l 利用下式计算迭代次数值η: n^ceil{ 5/4-1/2 · /og?(2/' · Δφ) } 上式中，函数ceil {χ}表示取不小于χ的最小整数； 第二步，坐标迭代过程： 设待计算的目标点坐标为（1〇, Qci)，利用下述步骤计算该目标点的查表坐标： 令迭代次数 i = 1，Ii= I cl，Qi = Qq; 步骤①：记扒=餅卜Ii和Q i的绝对值abs (I )和abs (Qi)，IjP Q i的正负 号 Sw1 和 S (i)Q; 步骤②：若 abs (Ii) = abs (Qi)，则 Ai= abs (I J，Bi= abs (Q J ;否贝丨J，令 Ai = max {abs (Ii)，abs (Qi) }，Bi= min {abs (I i)，abs (Qi) }，； 步骤③：记识'i= <37人?仏{為，及'}， 若当前迭代次数i〈n，则利用下式计算Ii+1和Qi+1: Ii+1= 2{2A /-Il2-I ； Qi+1= 4A AiZAi2-Il ； 令 Ii= 1 i+1，Qi= Q i+1，i = i+1，返回步骤②； 若当前迭代次数 i = n，则 An= A i，Bn= B i，Ii= I n，Qi= Q n，Sw1= S ('，Swq= S (n)Q， 并进入第三步查表过程； 第三步，查表过程： 利用在预处理过程构造的相位映射表中查坐标为（An，Bn)对应的相位值Φ ' n; 第四步，相位恢复过程： 利用下述规则计算最终相位结果： 若 An= abs (I n)，Bn= abs (Q n)，则 "'《二史 若 Bn= abs (I n)，An= abs (Q η)，则 V '"=Jl/2-供 其次： 当Sw1为正且S (n)Q为正，则％=W'"; 当Sw1为正且S (n)Q为负，则％ 当Sw1为负且S (n)Q为正，则外=J1-V',,； 当Sw1为负且S (n)Q为负，则外=<"-Jl; 令识Vl=奶，按照上述规则，通过识Vl获得奶-1，直至通过识'2获得供1，妁即最终相位 计算结果。
【专利摘要】本发明提供一种高精度相位小存储量查表计算方法。技术方案包括下述过程：第一步，预处理过程，计算目标点坐标值的量化位宽和迭代次数值，构造的相位映射表；第二步，坐标迭代过程，计算目标点的查表坐标；第三步，查表过程，利用相位映射表中查坐标对应的相位值；第四步，相位恢复过程，计算最终相位结果。本发明能够在与现有查找表方法相位映射表存储量相同的情况下大幅提高相位计算精度，或者在保持相同相位计算精度的情况下大幅度减少相位映射表规模，或者同时大幅提高相位计算精度和大幅度减少相位映射表规模。
【IPC分类】G06F19-00
【公开号】CN104636632
【申请号】CN201510104152
【发明人】杨军, 高凯, 朱江, 邹建彬
【申请人】中国人民解放军国防科学技术大学
【公开日】2015年5月20日
【申请日】2015年3月10日