自适应数字波束形成(adbf)是一种目前常采用的抗干扰处理方法，能在恶劣的敌方干扰和电磁兼容环境中大大提高系统的抗干扰能力，广泛应用于通信、雷达等领域。其中，并行性能较好的qr分解(正交分解)算法已经在实际系统中得到了应用。文献[3]提出了一种比传统qr分解adbf算法性能更优越的mqr分解(混合qr分解)smi(采样矩阵求逆)算法，他能避免传统qr分解算法需要通过前、后向回代才能得到自适应权向量的缺点，能做到真正意义上的实时权向量提取。在这些处理方式中一般采用标准givens旋转来实现qr分解或mqr分解，标准givens旋转包含开方和除法运算，保证足够精度及稳定性的开方和除法运算的运算量相当大。我们提出了一种基于cordic(坐标旋转数字计算机)技术的无开方无除法的mqr阵分解方法，并应用于自适应阵抗干扰处理，取得了良好的效果。

　　cordic算法最初由volder提出，最早用于三角函数的计算。其基本思想是：若要求平面矢量进行θ角度的旋转，则将此角度值θ分解，用一组预先规定好的基本角度的线性组合去逼近。该方法的巧妙之处在于基本角度的选取恰好使每次矢量以基本角度值旋转后，新矢量坐标值的计算只需要简单的移位和进行加法操作就能完成。我们讨论的就是利用cordic方法以避免开方和除法运算来实现givens旋转，从而进行mqr分解，完成自适应波束形成。

2 mqr-smi阵结构

　　图1给出了mqr-smi(混合qr分解采样矩阵求逆)阵结构(以三阵元为例)。阵中包含两种单元，边界单元(引导单元)aii和内部单元aij，bij和vi。随着自适应阵列接收数据从顶端输入，阵中各个单元存储数据依次被更新，其中边界单元aij产生旋转因子，后面的内部单元按照相同的旋转方式进行更新，最后在阵的底部实时地得到权向量w。

3 cordic技术

　　coric有两种计算模式：旋转(rotation)模式和求模(vectoring)模式。这里以求模模式来实现mqr阵分解。如图2所示，假设坐标为(x，y)的点旋转到(x′，y′)，两向量之间的夹角为θ，则有：

　　cordic技术的核心是采用一系列特殊的基本角度θi去逼近θ，其中θi满足：tan(θi)=±2-i。亦即把任意角度θ表示成一系列微小角度θi的前向或后向旋转的叠加(i=0，1，2，…)。这样就有：。其中ρi∈{1，-1，0)来代表第i次旋转的方向(向上、向下或不旋转)。于是图2所示的旋转可通过下面的过程来实现：

　　cordic旋转并非完美的旋转，他改变了矢量的模长。因而循环完成后要进行一个校正才能保证原矢量模长不变。

　　在求模模式中，目的是把坐标(x，y)旋转为(x′，0)。我们并不需要知道旋转的角度θ是多少，只要能通过该角进行旋转即可，亦即真正需要的是决定每次微旋转的方向的控制量ρi，然后通过相同的ρi去旋转其他坐标。ρi的确定：若x=0，则ρi=0；否则，ρi=sgn(x)sgn(y)。即旋转的目的是把矢量旋转到x轴上。

4 cordic更新mqr阵

　　cordic电路可以用在一个并行的systolic阵流水系统中来进行上述mqr阵的分解。在自适应波束形成中需要处理的数据为复数，这可以采用一个mqr处理单元执行3次cordic操作来解决。mqr阵每一行有一个引导单元，在引导单元中执行“求模”变换，内部单元则执行相同角度的“旋转”变换。具体说来，把他分为两种变换：θ变换和φ 变换。θ变换是一个相位变换，使引导单元变为实数，后面的内部单元进行相同的旋转变换；φ变换即为一个旋转，把复数坐标通过一