基于cordic ip核实现arctan

最新推荐文章于 2024-04-19 21:28:53 发布
最新推荐文章于 2024-04-19 21:28:53 发布
阅读量5.3k 收藏 74
点赞数 12
分类专栏： 小欧的FPGA进阶日记 文章标签： verilog fpga
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_42014414/article/details/108685941
版权

基于cordic ip核实现arctan@TOC

基于cordic ip核实现arctan

最近在研究相位补偿电路，需要用到反正切函数求解相位差，因此学习了该ip核的使用，写一篇小记录，希望能帮到有需要的盆友。

IP核介绍

简单粗暴地介绍一下IP核的使用。
在这里插入图片描述

  1. 首先将functional selection选为arc tan，随后设置输入位宽和输出位宽。
    需要注意的是，输入数据默认格式为“有2bit整数位的有符号定点小数”，输出数据默认格式为“有3bit整数位的有符号定点小数”，依据需要制定位宽即可。
    在这里插入图片描述

  2. 图中设置的input width为16bit，即cos和sin值都用16bit来表示，因此可以看到数据端口的输入数据总位宽为32bit。在实际写入数据时，将sin值写在[31:16]处，将cos值写在[15:0]处。
    需要注意的是，输入数据的总位宽一定为8的倍数。例如如果input width都为10bit，则有效数据为20bit，实际数据端口位宽为32bit。在实际写入数据时，需要将sin值写在[25:16]处，将cos值写在[9:0]处，[15:10]和[31:26]处补0即可。
    在这里插入图片描述

  3. 一通设置完以后，可以在左边栏看到延迟的时钟周期，如图所示延迟了20个时钟。
    在这里插入图片描述

仿真实现

为了验证该ip核的功能，写了一个.v文件，依次写入30°、45°、60°、60°、45°、30°、-30°、-45°、-60°、-60°、-45°、-30°… 的cos值和sin值，求取反正切的结果。
代码如下，较为基础，不做解释。

// A code block
module test_arctan(
input                   clk,
input                   rst_n,
output reg [15:0]      arctan_out_tdata,
output wire            arctan_out_tvalid
    );
 
 reg [2:0]                   state;
 reg [2:0]                   next_state;
 reg [4:0]                   cnt;
 reg [15:0]                  cos,sin;
 reg [31:0]                  arctan_in_tdata;
 reg                         arctan_in_tvalid; 
 wire[15:0]                  arctan_dout;
 
 localparam                  S_IDLE      =3'b001;
 localparam                  S_START     =3'b101;
 localparam                  S_CPT       =3'b100;
 localparam                  S_STOP      =3'b110;
 
 always@(posedge clk or negedge rst_n)
 begin
     if(rst_n==1'b0)
         state<=S_IDLE;
     else
         state<=next_state;
 end
 
 always@(*)
 begin
     case(state)
         S_IDLE:
             next_state<=S_START;
         S_START:
             next_state<=S_CPT;
         S_CPT:
             if(arctan_out_tvalid)
                 next_state<=S_STOP;
             else
                 next_state<=S_CPT;
         S_STOP:
             if(cnt<5'd23)
                 next_state<=S_START;
             else
                 next_state<=S_STOP;
        default:
            next_state<=state;
     endcase
 end
 
 always@(posedge clk or negedge rst_n)
 begin
     if(rst_n==1'b0)
         cnt<=0;
     else if(state==S_STOP && next_state!=state)
               cnt<=cnt+1;
           else
               cnt<=cnt;
 end
 
 always@(posedge clk or negedge rst_n)
 begin
     if(rst_n==1'b0)
         arctan_in_tvalid<=0;
     else if(state==S_START)
               arctan_in_tvalid<=1;
           else
               arctan_in_tvalid<=0;
 end
 
 always@(posedge clk or negedge rst_n)
 begin
     if(rst_n==1'b0)
         arctan_in_tdata<=0;
     else if(state==S_START)
               arctan_in_tdata<={sin,cos};
 end

 always@(posedge clk or negedge rst_n)
 begin
     if(rst_n==1'b0)
         arctan_out_tdata<=0;
     else if(next_state==S_STOP)
          begin
              if(cos[15]==1'b0)
                  arctan_out_tdata<=arctan_dout;
              else
                  arctan_out_tdata<=arctan_dout-16'b0110_0100_1000_0111;
          end
     else
         arctan_out_tdata<=arctan_out_tdata;
 end

 always@(*)
 begin
     case(cnt)
             5'd0: begin //30°
                       cos<=16'b0011_0111_0110_1100;//根号3/2
                       sin<=16'b0010_0000_0000_0000;//1/2                    
                   end    
             5'd1: begin //45°
                       cos<=16'b0010_1101_0100_0001;
                       sin<=16'b0010_1101_0100_0001;                   
                   end  
             5'd2: begin //60°
                       cos<=16'b0010_0000_0000_0000;
                       sin<=16'b0011_0111_0110_1100;                  
                   end 
             5'd3: begin //60°
                       cos<=16'b0010_0000_0000_0000;
                       sin<=16'b0011_0111_0110_1100;                  
                   end  
             5'd4: begin //45°
                       cos<=16'b0010_1101_0100_0001;
                       sin<=16'b0010_1101_0100_0001;                   
                   end                                                         
             5'd5: begin //30°
                       cos<=16'b0011_0111_0110_1100;//根号3/2
                       sin<=16'b0010_0000_0000_0000;//1/2                    
                   end    
             5'd6: begin //-30°
                       cos<=16'b1100_1000_1001_0100;
                       sin<=16'b0010_0000_0000_0000;                    
                   end 
             5'd7: begin //-45°
                       cos<=16'b1101_0010_1011_1111;
                       sin<=16'b0010_1101_0100_0001;                    
                   end 
             5'd8: begin //-60°
                       cos<=16'b1110_0000_0000_0000;
                       sin<=16'b0011_0111_0110_1100;                    
                   end                                                                           
             5'd9: begin //-60°
                       cos<=16'b1110_0000_0000_0000;
                       sin<=16'b0011_0111_0110_1100;                    
                   end                        
             5'd10: begin //-45°
                       cos<=16'b1101_0010_1011_1111;
                       sin<=16'b0010_1101_0100_0001;                    
                    end                  
             5'd11: begin //-30°
                       cos<=16'b1100_1000_1001_0100;
                       sin<=16'b0010_0000_0000_0000;                    
                    end           
             5'd12: begin //30°
                      cos<=16'b0011_0111_0110_1100;//根号3/2
                      sin<=16'b0010_0000_0000_0000;//1/2                    
                  end    
            5'd13: begin //45°
                      cos<=16'b0010_1101_0100_0001;
                      sin<=16'b0010_1101_0100_0001;                   
                  end  
            5'd14: begin //60°
                      cos<=16'b0010_0000_0000_0000;
                      sin<=16'b0011_0111_0110_1100;                  
                  end 
            5'd15: begin //60°
                      cos<=16'b0010_0000_0000_0000;
                      sin<=16'b0011_0111_0110_1100;                  
                  end  
            5'd16: begin //45°
                      cos<=16'b0010_1101_0100_0001;
                      sin<=16'b0010_1101_0100_0001;                   
                  end                                                         
            5'd17: begin //30°
                      cos<=16'b0011_0111_0110_1100;//根号3/2
                      sin<=16'b0010_0000_0000_0000;//1/2                    
                  end    
            5'd18: begin //-30°
                      cos<=16'b1100_1000_1001_0100;
                      sin<=16'b0010_0000_0000_0000;                    
                  end 
            5'd19: begin //-45°
                      cos<=16'b1101_0010_1011_1111;
                      sin<=16'b0010_1101_0100_0001;                    
                  end 
            5'd20: begin //-60°
                      cos<=16'b1110_0000_0000_0000;
                      sin<=16'b0011_0111_0110_1100;                    
                  end                                                                           
            5'd21: begin //-60°
                      cos<=16'b1110_0000_0000_0000;
                      sin<=16'b0011_0111_0110_1100;                    
                  end                        
            5'd22: begin //-45°
                      cos<=16'b1101_0010_1011_1111;
                      sin<=16'b0010_1101_0100_0001;                    
                   end                  
            5'd23: begin //-30°
                      cos<=16'b1100_1000_1001_0100;
                      sin<=16'b0010_0000_0000_0000;                    
                   end                                                                                                                                                 
             default: 
                   begin
                       cos<=16'b0000_0000_0000_0000;
                       sin<=16'b0000_0000_0000_0000;
                   end
     endcase
 end

 cordic_1 arctan (
      .aclk(clk),                                        // input wire aclk
      .s_axis_cartesian_tvalid(arctan_in_tvalid),  // input wire s_axis_cartesian_tvalid
      .s_axis_cartesian_tdata(arctan_in_tdata),    // input wire [31 : 0] s_axis_cartesian_tdata
      .m_axis_dout_tvalid(arctan_out_tvalid),            // output wire m_axis_dout_tvalid
      .m_axis_dout_tdata(arctan_dout)              // output wire [31 : 0] m_axis_dout_tdata
    );   
    
endmodule

问题分析

调试过程中有遇到一个问题。

  1. 当输入的正弦值为正，余弦值为正时，反正切得到的结果是对的，例如输入sin30°，cos30°时，输出结果为pi/6。
  2. 当输入的正弦值为负，余弦值为正时，反正切得到的结果也是对的，例如输入-sin30°，cos30°时，输出结果为-pi/6
  3. 但是当输入的正弦值为正，余弦值为负时，反正切得到的结果为正值，例如输入sin30°，-cos30°时，输出结果为6pi/5，不符号用户手册中对输出结果为-pi~pi的描述。

为了使第三种情况输出结果也为-pi/6，查了很多资料，最后在xilinx社区找到了解决办法。
在这里插入图片描述
因此，在代码中加了一段对cos值符号位的判断，由此来决定是否对输出结果进行“减去2pi”的操作。

 always@(posedge clk or negedge rst_n)
 begin
     if(rst_n==1'b0)
         arctan_out_tdata<=0;
     else if(next_state==S_STOP)
          begin
              if(cos[15]==1'b0)
                  arctan_out_tdata<=arctan_dout;
              else
                  arctan_out_tdata<=arctan_dout-16'b0110_0100_1000_0111;
          end
     else
         arctan_out_tdata<=arctan_out_tdata;
 end

最终输出正常，仿真结果如下。（图中为输入输出均为32bit的情况）
在这里插入图片描述

参考链接

xilinx社区的完整问题描述和解决方案： link.

进击的小欧
关注
  • 12
    点赞
  • 踩
  • 74
    收藏
    觉得还不错? 一键收藏
  • 5
    评论
专栏目录
Xilinx Cordic IP核的用法( 以arctan为例）
01-19
在ISE14.6开发环境中，调用Xilinx的Cordic IP核实现arctan算法。
用FPGA CORDIC IP核实现信号的相位检测，计算相位角
www_haha__的博客
03-04 2436
用FPGA CORDIC IP核实现信号的相位检测，计算相位角
FPGA Verilog Cordic算法实现三角函数计算，精度达到10e-5，附完整资料说明，包含串口收发功能
最新发布
mJXhltZhJl的博客
04-19 958
本文详细介绍了Cordic算法的原理和Verilog实现方法，并提供了完整的资料说明。同时，本文在altera板子上进行了完整的实验验证，证明了Cordic算法的正确性和性能。本文提供了完整的资料说明，包括FPGA Verilog Cordic算法的详细原理和实现方法。旋转因子的选择：Cordic算法的核心是选择合适的旋转因子，以确保计算的精度和效率。FPGA Verilog Cordic算法实现三角函数计算，可计算sin cos arctan，精度达到，10e-5，有完整资料说明。
FPGA数字信号处理（十四）Vivado Cordic IP核计算arctan
热门推荐
FPGADesigner的博客
06-23 2万+
在数字信号处理系统中经常需要计算arctan函数，比如在解调系统中由DDC得到复基带信号q和i支路计算arctan(q/i)即可得到基带信号的相位。在FPGA设计中可以使用CORDIC算法来实现arctan。 本文将介绍在Vivado开发环境下如何使用Xilinx提供的Cordic(6.0) IP核计算arctan。该IP核还可以实现其它CORDIC算法可实现的功能，将在后面的文章中介绍。 ...
Vivado Cordic IP核计算arctan
qq_33269442的博客
08-27 1462
Vivado Cordic IP核计算arctan 参考资料 IP核配置 Functional Selection：Arctan input wire [31 : 0] s_axis_cartesian_tdata→设置width为16 round mode：Nearest Even 参数配置(具体参照手册例子) 要设定好输入数据的整数和小数 CORDIC的输入数据范围应限制在-11（否则会出错），输出数据范围为（-pipi，选中Coarse Rotation）。且数据格式为带符号小数，输入数据整数位有2
vivado CORDIC ip核计算arctan记录
qq_41412445的博客
01-17 2758
本文主要记录自己使用vivado CORDIC ip 核计算arctan的过程，给出了具体过程进行参考。
时序模型：循环神经网络（RNN）
xunyishuai5020的博客
02-22 1万+
1. 循环神经网络（RNN）定义 循环神经网络（recurrent neural network, RNN）是一类专门设计处理不定长序列数据的神经网络。 与使用一种新计算1作为核心的卷积神经网络不同，循环神经网络仍使用特征的线性组合作为计算核心，并使用共享参数策略使模型能泛化不同长度的序列数据。 2. 循环神经网络的由来：从全连接神经网络（DNN）到循环神经网络（RNN） 2.1 全连接神经网络的不足 最初，研究者直接将时间序列看作特征向量，将序列每一时刻的观测作为一种特征输入全连接神经网络，来预测序列的标
循环神经网络RNN模型
个人学术整理
06-05 1732
rnn模型
RNN, LSTM, GRU模型的作用, 构建, 优劣势比较,attention机制
码python的Vinsmoke
10-29 2万+
1. RNN架构解析 1.1 认识RNN模型 学习目标 了解什么是RNN模型. 了解RNN模型的作用. 了解RNN模型的分类. 什么是RNN模型 RNN(Recurrent Neural Network), 中文称作循环神经网络, 它一般以序列数据为输入, 通过网络内部的结构设计有效捕捉序列之间的关系特征, 一般也是以序列形式进行输出. 一般单层神经网络结构: RNN单层网络结构: 以时间步对RNN进行展开后的单层网络结构: RNN的循环机制使模型隐层上一时间步产..
xilinx cordic IP核的用法- arctan的算法
桔色小西瓜
05-10 2万+
本文介绍如何使用xilinx的CORDIC核计算一个arctan的值。 此方法对于梯度运算有着重要的意义。 原理略 环境：xilinx ISE 14.7  cordic 4.0 ISim 首先是IP核的选项设置，如图： 标注1：选择函数的类型，我们选择计算arctan的值。 标注2：选择cordic的结构，是字串行还是并行，我们选择并行。 标注3：
CORDIC IP核 使用参考
12-21
CORDIC IP 核是 Field-Programmable Gate Array（FPGA）中的一种数字信号处理器件，可以实现CORDIC算法，用于计算三角函数值、指数函数值和对数函数值等。CORDIC IP 核使用参考主要介绍了在FPGA开发板上使用CORDIC ...
基于Cordic算法求反正切的Verilog及Matlab实现.zip
05-23
本文是根据Cordic算法原理，利用二分法思想，完成对任意输入值求反正切的Verilog编写及Matlab程序仿真，并且自创扩展到四象限，可以求（-180~180）全相位的反正切值，是求反正切非常实用的手段。
基于CORDIC算法的数字下变频技术设计与实现
01-31
为了解决这一问题，提出了一种基于CORDIC（coordinate rotation digital compute，坐标旋转数值计算）算法的数控振荡器的设计方法，而且该算法使数控本振和数字混频两个功能合在一起完成，省去了两个乘法器。...
cordic IP vhdl实现
05-28
提供的“cordic_IP_vhdl”文件应该包含了上述描述的Cordic算法的VHDL实现。这个IP核可以直接在FPGA或ASIC设计中复用，简化开发流程。使用者需要理解并适配该IP核的输入输出接口，确保与系统其他部分的正确连接。 ...
基于CORDIC的高速Sobel算法实现
10-15
为提高图像边缘检测的处理速度，提出一种基于CORDIC的高速Sobel算法实现。在FPGA平台上，在并行处理数据和流水线操作的基础上，使用扩展数据位和覆盖所有角度的流水线型CORDIC，提高Sobel的运算效率。实验结果表明，...
vivado中cordic IP 6.0 arctan 的使用
Ocean的机器学习之路
01-23 4888
有个项目需要用到arccos，然而，cordic中只有arctan，所以只能将就用一下了。哈哈 之前，有写过cordic中sin/cos 的使用。cordic IP 6.0 sin/cos 的使用  所以，也算对这个ip核比较熟悉了，这里只想记录一下  在使用arctan时，一些相关配置和输入输出的位数问题。如下图：   基本就需要注意一下这5个位置： 1，选择arctan模式。 2...
FPGA实现Cordic算法求解arctan和sqr(x*2 + y* 2)
weixin_44678052的博客
05-13 1777
由于在项目中需要使用的MPU6050，进行姿态解算，计算中设计到**arctan 和 sqr(x2 + y2),**这两部分的计算，在了解了一番之后，发现Cordic算法可以很方便的一次性求出这两个这两部分的计算。另外也可以一次性求出sin和cos的值。另外该算法还可以计算其他的一些公式(没做过多的了解)。
RNN神经网络
weixin_43391596的博客
11-26 1310
RNN神经网络---学习笔记
cordic ip核 sin/cos testbench测试
07-18
### 回答1： Cordic IP核是一种用于计算正弦和余弦函数的数学运算IP核。在进行Cordic IP核的sin/cos测试之前，我们需要编写一个测试台（testbench）来进行测试。 测试台的设计目标是在输入一组角度值后，计算出对应的正弦和余弦值，并与预期值进行对比，以确保Cordic IP核的功能正确性。 首先，我们需要确定测试的角度值范围。我们可以选择一个角度的正弦和余弦值是已知的，并以这个角度为中心，测试一定范围内的角度值。例如，我们可以选择0°到360°之间的某个角度值进行测试。 测试台的流程如下： 1. 定义输入信号和预期输出信号。 2. 将输入信号传递给Cordic IP核，并接收其计算得到的结果。 3. 将Cordic IP核计算的结果与预期输出进行对比，检查它们之间的误差是否在可接受范围内。 4. 如果误差在可接受范围内，则测试通过，输出测试结果。 5. 如果误差不在可接受范围内，则测试失败，输出测试结果。 在测试台中，我们可以通过生成随机角度值的方法，对Cordic IP核进行多组测试，以覆盖更广泛的角度范围。 测试台设计中需要考虑的几个重要因素包括输入输出数据的精度、测试角度的覆盖范围和测试结果的判定标准。通过合理设置这些参数，我们可以编写一个有效的Cordic IP核sin/cos测试台，以确保IP核的正确性。 ### 回答2： Cordic IP核是一种用于计算三角函数的硬件核心，通过使用迭代方法来实现高效的计算。为了测试Cordic IP核的sin/cos功能，我们可以编写一个Testbench。 Testbench的目的是验证Cordic IP核的输出结果是否正确。在Testbench中，我们首先需要生成一系列的输入值作为测试用例，这些输入值可以按照一定的步长递增或随机生成。然后，将这些输入值送入Cordic IP核进行计算，得到对应的sin和cos值。 接下来，我们需要知道这些测试用例的预期输出值。我们可以使用软件或数学库函数计算标准的sin和cos值作为参考。然后，将Cordic IP核的输出结果与预期输出值进行比较。如果两者非常接近，即可认为Cordic IP核的sin/cos功能正确。 在Testbench中，我们可以使用assert语句对比Cordic IP核的输出结果和预期输出值。如果两者不匹配，assert语句会给出错误提示。此外，我们还可以输出测试结果是否通过的信息，以便于判断Cordic IP核的性能和准确性。 最后，在Testbench中我们还可以记录测试的运行时间以及任何其他有用的信息，以便于分析和评估Cordic IP核的性能。 总结起来，通过编写一个完善的Testbench，我们可以对Cordic IP核的sin/cos功能进行全面的测试和评估，确保其在实际应用中的正确性和可靠性。 ### 回答3： Cordic是一种数字算法，通常用于计算三角函数（如sin和cos函数）。Cordic IP核是一种可以在FPGA中实现Cordic算法的IP核。 进行Cordic IP核的sin/cos测试需要编写一个测试平台，也就是所谓的testbench。测试平台主要用于验证Cordic IP核的正确性和性能。 测试平台的第一步是生成输入信号，即角度（或弧度）的值。可以使用随机数生成器来生成不同的角度输入。然后将这些角度值提供给Cordic IP核进行计算。 下一步是为Cordic IP核建立期望输出。可以使用系统级建模（System-Level Modeling）来计算期望的sin和cos值。使用已知的数学公式，将输入的角度值代入，计算出对应的sin和cos值。 然后，将输入信号提供给Cordic IP核，将计算结果与期望输出进行比较。可以使用assert语句来进行比较。如果计算结果与期望输出相差很小，可以认为测试通过。否则，可能需要检查IP核的实现是否存在问题。 在测试过程中，还可以记录每个输入角度和对应的计算结果，以便后续分析和调试。可以将这些记录在一个文件中，或者将其显示在仿真工具的波形窗口中。 测试完所有的输入角度后，可以进行一些统计分析，如计算平均误差、最大误差等。这些分析可以用于评估Cordic IP核的性能和准确性。 总之，Cordic IP核的sin/cos测试需要编写一个测试平台，生成输入信号，计算期望输出，并将计算结果与期望输出进行比较。通过分析比较结果，可以评估IP核的性能和准确性。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 5
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值