北邮信息工程信号与信息处理综合实验FPGA实验二报告（CORDIC算法）.docxVIP

信号与信息处理综合实验报告 第 PAGE 8页 信息与通信工程学院 信号与信息处理综合实验报告 （FPGA部分） 班 级： 姓 名： 学 号： 序 号： 日 期： 实验二 CORDIC算法 一、实验目的 掌握FPGA设计中的流水线技术； 掌握Cordic算法的基本原理及其实现方法； 了解通过在片内生成ROM的方式进行在板模块测试的方法。 实验内容 按实验指导书所给出的步骤，在FPGA上实现Cordic算法用于计算sin(x)。 修改程序使其能够用于计算x2+ 三、程序基本信息 （一）、程序模块描述： 本次实验里涉及了两个程序，一个是原始程序，另一个是对原始程序作简单修改后的功能拓展，代码本身基本是类似的。以下只对原始程序作较为详细的说明，而功能拓展只简单说明修改的情况。 1、CORDIC算法计算sin(x)（6次迭代）： 模块名：cordic。 以下介绍程序中的输入输出参数和临时变量，在这里首先要说明一点，本程序主要进行的是二进制位运算，因此所有角度和函数值也都是采用二进制表示，具体描述如下： 角度值：均采用8位二进制数表示，也就11111111（十进制0-255），对应角度值中的0-2π，因应π（180°）应π2（90°），依此类推。每一最低位表示360 函数值：程序中输出的三角函数值同样采用8位二进制数表示，但三角函数有正有负，因此实际采用的是8位补码，第一位是符号位（0为正，1为负），后面7位为数据位，01111111（0-127）对应0-+1，11111111（128-255）对应-1到-1256，每一最低位表示 输入及输出参数： clk：时钟信号； rst：重置信号； ena： 启动信号，该信号置1时系统开始运算过程； phase_in（8位二进制数）：输入的角度值； sin_out（8位二进制数）：输出的sin函数值 函数变量（仅在运行中调用）： cos_out（8位二进制数）：计算得到的cos函数值，与sin_out同时产生但不输出； phase_in_reg（8位二进制数）：程序中实际参与运算的角度值，下面会进行说明； x0-x6、y0-y6、z0-z6（均为8位二进制数）：在迭代算法中使用的临时参数，x和y参数存储三角函数值，z存储角度值。x0、y0、z0存储初值，其它参数存储6次迭代的临时值。 quadrant（8个1位存储器）：存储输入角度值的相位信息，具体用处见下面说明。 PIPELINE：流水线，必须和数据位长相同。本程序使用8级流水线。 主要工作流程： 以下先用流程图简单的说明工作过程。 可见程序本身的流程并不复杂，但是某些流程涉及到算法本身，下面作详细说明。 x0 这里置x0 x 在本程序中，x0取值用8位数表示为1An×27≈ 迭代运算： 此段代码是程序的主要部分，简单的说就是实现CORDIC算法的旋转模式： 它采用二进制位运算，以下对一些复杂地方作简要说明： eq \o\ac(○,1)zi（di）正负的判断：通过判断zi符号位（首位）； eq \o\ac(○,2)2-i的实现：即每次迭代时对xi或yi作除2操作，具体实现方法为把原二进制值向右移动一位，移出的空位补0 eq \o\ac(○,3)每次迭代时，z的值要减去的角度值满足公式tan-12-i。 根据原输入角度值所在象限对迭代后的结果进行处理： 利用三角函数的性质，由于结果是在角度位于第一象限的情况下获得的，因此可作如下处理。 （输入角度在）第二象限：sinx值不变，cos 第三象限：sinx和cos 第四象限：cosx值不变，sin 取负方法：所有二进制位取反加一（补码）。 2、CORDIC算法计算sin(x)（10次迭代）： 鉴于8次和10次迭代都是原程序（6次迭代）的拓展，且10次迭代的修改更多，因此本次实验中我只做了10次迭代的功能拓展。 这个程序是对6次迭代的程序进行修改后获得的，基本代码与工作流程完全相同，以下只说明修改的地方。 （1）数据位长度：由上述公式可得，每次xi、yi和zi的变化范围随迭代次数增加而减小，若仍采用原先的8位二进制数存储数据，到了较高次迭代时会出现数据过小而无法完全表示的情况，因此数据位长需要增加。本程序中，xi、yi和zi位长度都增加到10（zi按理论值应增加到 （2）增加迭代次数及对应临