基于LUT的DDS的设计
DDS(Direct Digital Synthesis)直接数字频率合成技术由USA J.Tierncy首先提出。是以数字信号处理理论为基础，从相位概念出发直接合成所需波形的新的全数字技术的频率合成方法。DDS主要出现在数字混频系统中。在数字混频中，通过DDS产生正交的本地振荡信号即正、余弦信号与输入信号相乘实现频谱搬移，如通信系统的调制、解调。目前FPGA实现DDS有三种途径：基于IIR滤波器的实现方法、基于查找表LUT的实现方法以及基于CORDIC算法的实现方法。其中采用LUT的方法较为通用且比较容易实现。
1 基于LUT的DDS理论
一个典型的基于LUT的DDS系统由相位累加器和波形存储器两部分构成，如图1所示。图中相位累加器的位宽为nbit,步进值为μ,LUT的深度N为2n,宽度为Lbit。LUT中依相位顺序存储一个周期的波形数据。
图1 基于LUT的DDS架构
2 matlab 仿真
图2 相位累加器位宽为4bit时LUT地址与存储数据的对应关系
如图2所示，我们要制作地址位宽为10bit数据位宽为32bit的LUT查找表。Matlab源码如下：
clear all
close all
clc
FS = 44100;%采样率
fc = 1000;
N = 1024;%1024点一个正玄周期
t =0:2*pi/(N-1):2*pi;
sin_dds = sin(t);
figure,
hold on
plot(t,sin_dds,'*');
legend('sin');
grid;
hold off
图3 matlab 1024点1个周期的正玄波
3 FPGA仿真
整个仿真结构如图1所示，由相位累加控制器和sin波形存储器组成。仿真生成采样率为44.1KHZ @1KHZ正玄波和余弦波(相位相差90度)。
tb_dds源码：
1. `timescale 1ns / 1ps
2. `define NULL 0
3. Module tb_dds();
4.
5. paRAMeter KHZ1 =44;
6. paRAMeter KHZ2 =22;
7. parameter KHZ4 =11;
8. reg clk;
9. initial begin
10. clk = 0;
11. #313333 clk = 0;
12. forever #11.072 clk = ~clk;
13. end//45.1584MHZ
14.
15. reg reset_n;
16. initial begin
17. reset_n = 1;
18. #10000 reset_n = 0;
19. #100000 reset_n = 1;
20. end
21.
22. reg [9:0] addr;
23. reg [9:0] cnt;
24. reg [31:0] k;
25. integer fb;
26.
27. wire signed [31:0] sin;
28. wire signed [31:0] cos;
29.
30. //Phase accumulation controller
31. always @(posedge clk or negedge reset_n) begin
32. if(!reset_n) begin
33. addr <= 0;
34. cnt <= 0;
35. end
36. else if(cnt == KHZ1-1) begin
37. cnt <= 0;
38. addr <= addr +10'd1;
39. end
40. else begin
41. cnt <= cnt + 1;
42. addr <= addr;
43. end
44. end
45.
46. //reg [8:0] i;//88.2khz
47. reg [9:0] i;//44.1khz Sample rate
48. reg signed [31:0] sin_slow;
49. reg signed [31:0] cos_slow;
50.
51. always @(posedge clk or negedge reset_n) begin
52. if(!reset_n) begin
53. i <= 0;
54. sin_slow <= 0;
55. cos_slow <= 0;
56. k <=0;
57. fb = $fopen("sin.dat","w");
58. if(fb == `NULL) begin
59. $display("Can not open sin.dat");
60. $finish;
61. end
62. end
63. else begin
64. i <= i+1;
65. if(i == 0) begin
66. sin_slow <= sin;
67. cos_slow <= cos;
68. k <= k +1;
69. $display("time=[%d],%d",$realtime,sin);
70. $fwrite(fb,"%d ",sin);
71. if(k == 4096) begin
72. $fclose(fb);
73. $stop;
74. end
75. end
76. end
77. end
78.
79.
80.
82. .clk(clk),
83. .addr(addr),//0-1023 1T
84. .sin(sin),
85. .cos(cos)
86. );
87. endModule
仿真结果：
生成Sample rate 44.1KHZ@1KHZ的正玄和余弦波相位相差90°。
图4 modelsim 仿真结果
Matlab FFT分析:
图5 FFT分析结果
由图5可知产生的正玄波频率为1KHZ。基于LUT的DDS设计完成，此节将是后期基于FPGA数字信号处理的基础和关键。
编辑：黄飞
