cordic算法的fpga实现
cordic算法参考:http://wenku.baidu.com/view/6c623aa8910ef12d2bf9e732.html
这是百度文库的一个文档,详细介绍了cordic算法的基本内容。我就是读了这个文档后做的cordic算法实现sin cos的运算的。
先谈感悟:做这种迭代型的算法,我觉得用for语句实现是最容易理解的,所以我的设计就是基于for结构实现的,关于算法中重要的位数确定,上边的文档里也有说明,我的选择是z选17位,1位符号位,8位整数,8位小数,x,y选21位,20位小数,一位符号位。所有数都是照这个标准的。
贴代码(贴后格式改变了,没缩进了,凑合看)
module test_cordic(clk,rst_n,z0,en,sin_out,cos_out);//z0是输入角度值,en是运算使能,sin_out,cos_out是运算输出。
input clk,rst_n;
input [16:0] z0;
input en;
output [10:0] sin_out,cos_out;
reg en_buff;
always@(posedge clk or negedge rst_n)
begin
if(!rst_n)
en_buff<=0;
else
en_buff<=en;
end
reg [16:0] angel[0:9];
always@(posedge clk or negedge rst_n)
begin
if(!rst_n)
begin
angel[0]<=17'b00010110100000000;//45
angel[1]<=17'b00001101010011001;//26.6
angel[2]<=17'b00000111000000000;//14
angel[3]<=17'b00000011100011001;//7.1
angel[4]<=17'b00000001110011001;//3.6
angel[5]<=17'b00000000111001100;//1.8
angel[6]<=17'b00000000011100110;//0.9
angel[7]<=17'b00000000001100110;//0.4
angel[8]<=17'b00000000000110011;//0.2
angel[9]<=17'b00000000000011001;//0.1
end
end
reg [16:0]reg_z[0:9];//1 fuhao wei 8 zhengshu 8xiaoshu
reg [10:0] reg_x[0:9];
reg [10:0] reg_y[0:9];
reg [4:0] i;
always@(posedge clk or negedge rst_n)
begin
if(!rst_n)
begin
for(i=0;i<=9;i=i+1'b1)
begin
reg_x[i]<=0;
reg_y[i]<=0;
reg_z[i]<=0;
end
end
else if(en&&!en_buff)
begin
reg_x[0]<=11'b01001101101;//0.6073<<10 + 1wei fuhao wei
reg_y[0]<=0;
reg_z[0]<=z0;
end
else
begin
for(i=1;i<=9;i=i+1'b1)
begin
if(reg_z[i-1][16])
begin
reg_x[i]<=reg_x[i-1]+(reg_y[i-1]>>i-1);
reg_y[i]<=reg_y[i-1]-(reg_x[i-1]>>i-1);
reg_z[i]<=reg_z[i-1]+angel[i-1];
end
else
begin
reg_x[i]<=reg_x[i-1]-(reg_y[i-1]>>i-1);
reg_y[i]<=reg_y[i-1]+(reg_x[i-1]>>i-1);
reg_z[i]<=reg_z[i-1]-angel[i-1];
end
end
end
end
assign sin_out=reg_y[9];
assign cos_out=reg_x[9];
endmodule
测试代码较简单,不贴了,贴一个30度求值结果
结果是sin30=00111111111;什么意思?最高位符号位是0,代表正数,后10位是小数部分即0.0111111111;转为10进制就是0.499;同理cos30=01101110110转为10进制就是0.8652.当然也可以试试其他角度。