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.当然也可以试试其他角度。