cordic ip核 vivado_Xilinx Vivado CORDIC IP求解atan 反正切

赛灵思官方提供了cordic(coordinate rotational digital computer) ip核实现直角坐标极坐标变化，三角函数的操作。我介绍下它进行反正切求解的使用：

新建个简单工程：bd如下

进行ip设置，选择运算位反正切后，ip端口回自动变为上图，再引出2个总线和时钟，xilinx的ip核不少是基于AXI4-Stream总线，这里使用并不复杂，默认只有2个信号，一个数据线tdata，一个握手信号tvalid，tvalid拉高时数据信号有效。

确定了输入输出位宽后，系统会自动求解出需要的延时latency。这里说一句，AXI-Stream总线收发基于byte的，无论输入什么位宽最后都是8的整数倍。所以需要参照UG，放置好有效输入位和分解输出位。

如果输入输出位宽恰好不是8的倍数，那么tdata位宽是大于设定的值的。所以需要上图提到进行填充(PAD)，比如我输入20位，x和y加起来40位，那么输入S_AXIS的位宽是48，我要在0-19放x，20-23放填充数据(随意，我用0)，24-43放y，44-47放填充数据。必须保证0-19和24-25分别是我的x和y。

atan=y/x，要求输入必须是[-1,1]，所以如果数据不是这个区间还需要进行归一化处理(可以用除法器div_gen ip核，这里不做介绍了)。上图是一个简单的示例，介绍输入输出的定点小数格式。整数部分的第一位是符号位，小数部分没有符号位是正的。※Q格式：小数点位于第 n 位元之右侧，称为Qn 格式。

这里关于二进制小数表示，可以看下这个https://blog.csdn.net/AaricYang/article/details/87882868。整数的每一位权重是2^(n-1) 111就是2^2+2^1+2^0=7，小数部分的权重是2^(-n) 0.111就是2^(-1)+2^(-2)+2^(-3)=0.875

这里我觉得还有一个细节问题。如果输入的数值x或y是负数小数，比如1100，常规默认符号位第一位，第二位开始就是小数部分了，所以1010就是1.010= -1+0.25= -0.75；这种符号位后一位就是小数部分的第一位的格式就是vivado 除法器ip输出小数的格式。也可以说这种是0Qn format，即符号位后0位就是小数。但是cordic ip输入要求 1Qn format，符号位结束后1位才是小数格式(为了兼容恰好等于±1的情况)，那么刚刚1010就不可以直接作为cordic ip核的输入，1010按照1Qn format，第3位才是小数部分，1010就是10.10= -1.5 因为整数部分是10，反码01，补码10为2，因为是负数所以-2，小数部分是0.5，最终结果是-2+0.5=1.5。这样就不符合cordic ip核输入要求了。从0Qn改为1Qn为了满足，如果是绝对值小于1的数只要在符号位后重复符号位就可以了，1 010改为11 010，这样就可以正常输入到cordic中了。

`timescale 1 ns / 1psmodulecordic_tb_top;regaclk;wire [31:0]M_AXIS_DOUT_0_tdata;wireM_AXIS_DOUT_0_tvalid;reg [47:0]S_AXIS_CARTESIAN_0_tdata;regS_AXIS_CARTESIAN_0_tvalid;integer handle;

initial beginaclk=0;

S_AXIS_CARTESIAN_0_tvalid=1'b0;

S_AXIS_CARTESIAN_0_tdata='b0;

end

initial handle = $fopen("D:/FPGAcode/_file/cordic.txt");//打开文件

always #10 aclk=~aclk;localparam PAD=(48-20*2)/2;always@(posedge aclk) beginS_AXIS_CARTESIAN_0_tvalid<=1'b1;

S_AXIS_CARTESIAN_0_tdata<={{PAD{1'b0}},{3'b000,17'b0},{PAD{1'b0}},{3'b111,17'b0}};end

always@(posedge aclk) begin

if(M_AXIS_DOUT_0_tvalid)

$fdisplay(handle,"%b",M_AXIS_DOUT_0_tdata);//写数据

enddesign_1_wrapper design_1_i

(.M_AXIS_DOUT_0_tdata(M_AXIS_DOUT_0_tdata),