找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 2769|回复: 0
打印 上一主题 下一主题
收起左侧

cordic算法详细Verilog的代码

[复制链接]
跳转到指定楼层
楼主
cordic算法的详细代码,包括tb代码


源程序如下:
  1. `timescale 1ns / 1ps
  2. //////////////////////////////////////////////////////////////////////////////////
  3. // Company:
  4. // Engineer:
  5. // Design Name:
  6. // Module Name:    Cordic
  7. // Project Name:
  8. // Target Devices:
  9. // Tool versions:
  10. // Description:
  11. //
  12. // Dependencies:
  13. //
  14. // Revision:
  15. // Revision 0.01 - File Created
  16. // Additional Comments:
  17. //
  18. //////////////////////////////////////////////////////////////////////////////////

  19. module Cordic(clk, phi, cos, sin);
  20.         parameter W = 10,W_Z = 11;
  21.         parameter PIPELINE=10;
  22.    input clk;
  23.    input [W_Z-1:0] phi;
  24.    output[W-1:0] cos,sin;
  25.        
  26.         reg [W-1:0] cos,sin;
  27.         reg [W-1:0] x[10:0],y[10:0];
  28.         reg [W_Z-1:0] z[9:0];
  29.         reg  sign[PIPELINE:0];
  30.         integer i;
  31.        
  32.         always @(posedge clk)
  33.         begin
  34.                 x[0] <= 10'h137; // 修正CORDIC算法的比例因子,An的倒数
  35.                 y[0] <= 10'h00;       
  36.                 z[0] <= phi;
  37.                 // 旋转45度
  38.                 if(z[0][W_Z-1]) begin
  39.                         x[1] <= x[0] + y[0];       
  40.                         y[1] <= y[0] - x[0];       
  41.                         z[1] <= z[0] + 11'h192;
  42.                 end
  43.       else begin
  44.                         x[1] <= x[0] - y[0];       
  45.                         y[1] <= y[0] + x[0];       
  46.                         z[1] <= z[0] - 11'h192;
  47.                 end
  48.                 // 旋转26.57度
  49.                 if(z[1][W_Z-1]) begin
  50.                         x[2] <= x[1] + {{1{y[1][W-1]}}, y[1][W-1:1]};
  51.                         y[2] <= y[1] - {{1{x[1][W-1]}}, x[1][W-1:1]};
  52.                         z[2] <= z[1] + 11'hF1;
  53.                 end
  54.                 else begin
  55.                         x[2] <= x[1] - {{1{y[1][W-1]}}, y[1][W-1:1]};
  56.                         y[2] <= y[1] + {{1{x[1][W-1]}}, x[1][W-1:1]};
  57.                         z[2] <= z[1] - 11'hF1;
  58.                 end
  59.                 // 旋转14.04度
  60.                 if(z[2][W_Z-1]) begin
  61.                         x[3] <= x[2] + {{2{y[2][W-1]}}, y[2][W-1:2]};
  62.                         y[3] <= y[2] - {{2{x[2][W-1]}}, x[2][W-1:2]};
  63.                         z[3] <= z[2] + 11'h80;
  64.                 end
  65.                 else begin
  66.                         x[3] <= x[2] - {{2{y[2][W-1]}}, y[2][W-1:2]};
  67.                         y[3] <= y[2] + {{2{x[2][W-1]}}, x[2][W-1:2]};
  68.                         z[3] <= z[2] - 11'h80;
  69.                 end
  70.                 // 旋转7.13度
  71.                 if(z[3][W_Z-1]) begin
  72.                         x[4] <= x[3] + {{3{y[3][W-1]}}, y[3][W-1:3]};
  73.                         y[4] <= y[3] - {{3{x[3][W-1]}}, x[3][W-1:3]};
  74.                         z[4] <= z[3] + 11'h40;
  75.                 end
  76.                 else begin
  77.                         x[4] <= x[3] - {{3{y[3][W-1]}}, y[3][W-1:3]};
  78.                         y[4] <= y[3] + {{3{x[3][W-1]}}, x[3][W-1:3]};
  79.                         z[4] <= z[3] - 11'h40;
  80.                 end
  81.                 // 旋转3.58度
  82.                 if(z[4][W_Z-1]) begin
  83.                         x[5] <= x[4] + {{4{y[4][W-1]}}, y[4][W-1:4]};
  84.                         y[5] <= y[4] - {{4{x[4][W-1]}}, x[4][W-1:4]};
  85.                         z[5] <= z[4] + 11'h20;
  86.                 end
  87.                 else begin
  88.                         x[5] <= x[4] - {{4{y[4][W-1]}}, y[4][W-1:4]};
  89.                         y[5] <= y[4] + {{4{x[4][W-1]}}, x[4][W-1:4]};
  90.                         z[5] <= z[4] - 11'h20;
  91.                 end
  92.                 // 旋转1.79度
  93.                 if(z[5][W_Z-1]) begin
  94.                         x[6] <= x[5] + {{5{y[5][W-1]}}, y[5][W-1:5]};
  95.                         y[6] <= y[5] - {{5{x[5][W-1]}}, x[5][W-1:5]};
  96.                         z[6] <= z[5] + 11'h10;
  97.                 end
  98.                 else begin
  99.                         x[6] <= x[5] - {{5{y[5][W-1]}}, y[5][W-1:5]};
  100.                         y[6] <= y[5] + {{5{x[5][W-1]}}, x[5][W-1:5]};
  101.                         z[6] <= z[5] - 11'h10;
  102.                 end
  103.                 // 旋转0.90度
  104.                 if(z[6][W_Z-1]) begin
  105.                         x[7] <= x[6] + {{6{y[6][W-1]}}, y[6][W-1:6]};
  106.                         y[7] <= y[6] - {{6{x[6][W-1]}}, x[6][W-1:6]};
  107.                         z[7] <= z[6] + 11'h8;
  108.                 end
  109.                 else begin
  110.                         x[7] <= x[6] - {{6{y[6][W-1]}}, y[6][W-1:6]};
  111.                         y[7] <= y[6] + {{6{x[6][W-1]}}, x[6][W-1:6]};
  112.                         z[7] <= z[6] - 11'h8;
  113.                 end
  114.                 // 旋转0.45度
  115.                 if(z[7][W_Z-1]) begin
  116.                         x[8] <= x[7] + {{7{y[7][W-1]}}, y[7][W-1:7]};
  117.                         y[8] <= y[7] - {{7{x[7][W-1]}}, x[7][W-1:7]};
  118.                         z[8] <= z[7] + 11'h4;
  119.                 end
  120.                 else begin
  121.                         x[8] <= x[7] - {{7{y[7][W-1]}}, y[7][W-1:7]};
  122.                         y[8] <= y[7] + {{7{x[7][W-1]}}, x[7][W-1:7]};
  123.                         z[8] <= z[7] - 11'h4;
  124.                 end
  125.                 // 旋转0.22度
  126.                 if(z[8][W_Z-1]) begin
  127.                         x[9] <= x[8] + {{8{y[8][W-1]}}, y[8][W-1:8]};
  128.                         y[9] <= y[8] - {{8{x[8][W-1]}}, x[8][W-1:8]};
  129.                         z[9] <= z[8] + 11'h2;
  130.                 end
  131.                 else begin
  132.                         x[9] <= x[8] - {{8{y[8][W-1]}}, y[8][W-1:8]};
  133.                         y[9] <= y[8] + {{8{x[8][W-1]}}, x[8][W-1:8]};
  134.                         z[9] <= z[8] - 11'h2;
  135.                 end
  136.                 // 旋转0.11度
  137.            if(z[9][W_Z-1])
  138.                 begin
  139.                         x[10] <= x[9] + {{9{y[9][W-1]}}, y[9][W-1:9]};
  140.                         y[10] <= y[9] - {{9{x[9][W-1]}}, x[9][W-1:9]};
  141.                 //        z[10] <= z[9] + 11'h1;
  142.                 end
  143.                 else begin
  144.                         x[10] <= x[9] - {{9{y[9][W-1]}}, y[9][W-1:9]};
  145.                         y[10] <= y[9] + {{9{x[9][W-1]}}, x[9][W-1:9]};
  146.                 //        z[10] <= z[9] - 11'h1;
  147.                 end

  148.        
  149.                        
  150.                 cos <= x[10][9]?~x[10]:x[10];
  151.                 sin <= sign[PIPELINE]? ~(y[10][9]?~y[10]:y[10]) + 1 : y[10][9]?~y[10]:y[10];
  152.         end
  153.        
  154.         always @(posedge clk)
  155.                 begin
  156.                   for(i=0;i<PIPELINE;i=i+1)
  157.          sign[i+1] <= sign[i];
  158.                 sign[0] <= phi[W_Z-1];
  159.                 end
  160.        
  161. endmodule
复制代码

所有资料51hei提供下载:
7660161.zip (1.36 MB, 下载次数: 17)

评分

参与人数 1黑币 +50 收起 理由
admin + 50 共享资料的黑币奖励!

查看全部评分

分享到:  QQ好友和群QQ好友和群 QQ空间QQ空间 腾讯微博腾讯微博 腾讯朋友腾讯朋友
收藏收藏 分享淘帖 顶 踩
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

手机版|小黑屋|51黑电子论坛 |51黑电子论坛6群 QQ 管理员QQ:125739409;技术交流QQ群281945664

Powered by 单片机教程网

快速回复 返回顶部 返回列表