本文介绍如何使用xilinx的CORDIC 核生产一个sin和cos波形?
cordic的原理就不介绍了,百度一大堆,我们知道原理后,需要去使用它。
环境:xilinx ISE 14.5
cordic 4.0
modelsim
首先是IP核的选项设置,照例的有图才方便学习~
标注1:选择函数的类型,我们选择计算sin和cos值。
标注2:选择cordic的结构,是字串行还是并行,我们选择并行。
标注3:选择输出流水线类型,我们选择不要流水线。
关于cordic的计算结构,是并行还是串行,datasheet上面有介绍时序图:
看时序图可知,serial类型计算一个值,需要的clk比parallel需要的多,但是datasheet后面有说serial占用的资源少。
如果感兴趣可以看datasheet上介绍的二者的内部结构。
serial 占用的资源少,需要的clk多,也会使端口多出几个PIN,控制起来麻烦一点,如下图:
关于输出流水线类型,我们选择NO Pipelining ,如果选择其他的,那么输出的结果会延迟不定时的clk才产生,当然好处是节省资源。
选择NO Pipelining 时, 1个clk就出一个计算值,有点费资源。文档的最后有对二者的资源占用比较,可以参考。
如果对延时不在意,但是想节省资源的话,就选择优化。
datasheet相关介绍:
标注1:选择相位角的格式,我们选择Scaled Radians(意思是多少PI)。
标注2:选择数据宽度,这个和你需要多少位宽的数据有关,我们选择16位。
标注3:舍位模式,选择的是近似值。
对于相位角格式,datasheet的定义为:
数据的有效值的定义为:
相位角: 第一位是符号位,第二、三位是整数位,其他是小数位。
sin数据: 第一位符号位,第二位整数位,其他小数位。
cos数据:第一位符号位,第二位整数位,其他小数位。
数据的舍位模式:(具体可以看下表)
标注1:选择阶乘和精度,0表示根据数据来自动选择。
标注2:选择cordic算法的范围,将输出值控制在-pi/4 ~ +pi/4之间还是-pi ~ +pi之间。
标注3:选择PIN,本来sin和cos功能是可以不要RDY的,但是我需要RDY来查看什么时候输出数据才有效。
标注2的详解:
如果你要生成正弦余弦波,那么你的Coarse Rotation必须选上,让输出数据是整个圆平面。
这样设置好的cordic IP核只是第一步,下一步是来设置合适的相位phase,以产生完整的波形。
datasheet所示:
我设置的数据位是16位,那么+1.0 -> 16‘h2000; -1.0 -> 16'he000;
(暂时给定phase的步进为16’h00C0,至于怎么确定,下面会说到,cordic的时钟我设置的是100Mhz)
如果直接让phase从 16‘h2000 到 16'he000 ,那么波形就会变成这样:
原因是,16‘h2000 到 16'he000 之间有一长段的phase是大于+1和小于-1的,cordic IP核会当做+1和-1计算,所以才会一直是相同的值。
那么怎样才能产生完整的正弦波?
我们将phase从16'he000 到 16‘h2000 ,仿真结果为:
发现有的地方是完整的波形了,那么现在我们需要找到波形的衔接点:
我选取x_out的二个波谷做一个周期,那么找到二个波谷对应的phase即可,注意cordic计算是需要clk,当rdy为1时,数据才有效,所以我们要找到延时的clk数。
我用cnt来计算,大概11个clk,数据有效。(注意不同的cordic设置会产生不同的延时,所以cnt不是固定的)
那波谷对应的phase是向前数11个clk对应的phase:16‘hd980
同样找到另一个波谷的phase值为:16‘h1940
仿真产生完整的正弦余弦波:
测试的testbench:
那么步进的phase怎么计算呢?
我选择的是16’h00C0,十进制是192。
刚才提到了,+1 到 -1(16‘h2000 到 16'he000)之前有 :16'he000 - 16‘h2000 = 16‘hC000,也就是49152.
256 * 192 = 49152,那么当步进取192的时候,256个clk就可以从2000加到E000了。这样不会出现起始数据不定的状况。这里建议数值是49152的约数比较好,当然不是它的约数也行。
16‘h00C0 -> 0000 0000 1100 0000,按照datasheet的phase数据格式:000,0 0000 1100 0000,也就是十进制的0.0234375。而这个数的角度值为:0.0234375 * 180° = 4.21875°,也就是每个clk步进4.21875°。
一个圆有360°,那么从0°到360°一共步进360° / 4.21875° 约等于 85.3次。
我选择的clk为10ns,也就是85.3 * 10 ns输出一个完整的正弦余弦波。
所以正弦余弦波的周期为1 / 853 ns 约等于 1.17Mhz。
如果要产生自己想要的频率,那么需要更改步进,然后根据仿真找到衔接点的起始位置。
补充:怎样在modelsim里面看波形,而不是数据。
如果要从 波形形式 变回 数据形式,选择Literal格式即可。 |