第4章基于CCS的波束形成算法实现.docVIP

第4章 基于CCS的波束形成算法实现 CCS（Code Composer Studio）是TI公司推出的用于开发DSP芯片的集成开发环境，它采用Windows风格界面，集编辑、编译、链接、软件仿真、硬件调试以及实时跟踪等功能于一体，极大地方便了DSP芯片的开发与设计，是目前使用最为广泛的DSP开发软件之一。CCS有两种工作模式，即软件仿真器模式，可以脱离DSP芯片，在PC机上模拟DSP的指令集和工作机制，主要用于前期算法实现和调试；硬件在线编程模式，可以实时运行在DSP芯片上 [11]。 CCS的功能十分强大，支持C/C++和汇编的混合编程，本仿真程序全部使用C语言编写[12]。4.1 C语言实现复数矩阵基本变换 4.1.1 波束形成算法和矩阵运算分析 由上一章节中，利用Matlab实现几种波束形成算法仿真，分析可知，在算法求解过程中，对于矩阵的操作运算是整个波束形成算法的核心内容。分析几种波束形成算法可知，所涉及的矩阵运算主要有复矩阵共轭转置、复矩阵相乘、复矩阵求逆[13]。由于C语言无法直接的处理复数，没有相应的数据类型，所以只有间接的进行复数运算。把实虚部分开，分别做上标记放在不同的数组里面，处理的结果也分存在不同的数组里面。 本论文主要以MMSE波束形成算法为例，讨论该算法在CCS上的实现方法。分析波束形成中最优权值求解过程： 设最优权矢量为： （4-1） 其中阵元数。 由MMSE准则计算出的最优权向量为。 （4-2） 其中为输入矢量的自相关矩阵。 （4-3） 输入矢量与需要信号的互相关矩阵。 （4-4） 可见，在求解最优权矢量的过程中，复数矩阵的相乘运算、转置运算、共轭转置运算，以及矩阵求逆运算，是最优权值求解的核心。所以在这里，先进行波束形成算法中通用的矩阵基本操作的函数的设计[14]。 本论文中，CCS下待处理的数据由Matlab随机生成。考虑到存储空间问题，将阵列接收信号矢量阵X[M][N]，M为阵元数，N为快拍数，以及期望信号S[N]，N为快拍数作为本仿真程序的信号源数据直接在CCS下进行处理。 4.1.2 复矩阵共轭转置 求阶复数矩阵的共轭转置矩阵。其中 （4-5） （4-6） 将矩阵的行元素变成列元素，列元素变成行元素，虚部变成原来的相反数，存在矩阵中。 函数语句的形参说明见表4-1： 表4-1　形参说明 形参和函数类型 参数意义 double ar[m][n],ai[n][m] 存放矩阵的实部与虚部元素 double br[m][n],bi[n][m] 返回矩阵的实部与虚部元素 void contrans 过程 函数程序的主体部分如下： void contrans(ar,ai,br,bi) double ar[N][L],ai[N][L],br[L][N],bi[L][N]; { int i,j; double Tbi[L][N]; for(i=0;iN;i++) for(j=0;jL;j++) { br[j][i]=ar[i][j]; Tbi[j][i]=ai[i][j]; } for(i=0;iL;i++) for(j=0;jN;j++) bi[i][j]=-Tbi[i][j]; } 4.1.3 复矩阵相乘 求阶复数矩阵与阶矩阵的乘积矩阵。其中： （4-7） （4-8） （4-9） 乘积矩阵中的各元素为： （4-10） 其中两个复数相乘采用如下算法[15]： 假设 （4-11） 令 （4-12） 则有 （4-13） 函数语句的形参说