基于FPGA的DDS任意波形发生器设计

一、简介

DDS技术最初是作为频率合成技术提出的，由于其易于控制，相位连续，输出频率稳定度高，分辨率高, 频率转换速度快等优点，现在被广泛应用于任意波形发生器(AWG)。基于DDS技术的任意波形发生器用高速存储器作为查找表，通过高速D/A转换器来合成出存储在存储器内的波形。所以它不仅能产生正弦、余弦、方波、三角波和锯齿波等常见波形，而且还可以利用各种编辑手段，产生传统函数发生器所不能产生的真正意义上的任意波形。

二、原理

根据傅立叶变换定理可知，任何周期信号都可以分解为一系列正弦或余弦信号之和，不失一般性，以正弦信号的产生为例详细说明直接数字频率合成技术的原理。比如一个频率为f_c的正弦信号，其时域表达式为：

其相位表达式为：

从两式可以看出，正弦信号是关于相位的一个周期函数，下图更加直观的描述相位与幅度的关系，16个相位与16个幅度值相对应，即每一个相位值对应一个幅度值，比如1100对应的相位为3π/2,对应的幅度值为-1.

相位和幅值的一一对应关系就好比存储器中地址和存储内容的关系，如果把一个周期内每个相位对应的幅度值存入存储器当中，那么对于任意频率的正弦信号，在任意时刻，只要已知相位Φ(t)，也就知道地址，就可通过查表得到s(t)。下图是DDS的基本结构框图：

相位累加器在每个时钟脉冲输入时，把频率控制字累加一次，相位累加器的输出数据就是信号的相位，用输出的数据作为波形存储器(ROM)的相位取样地址，这样就可以把存取在波形存储器内的波形抽样值经查找表查处，完成相位到幅值的转换。频率控制字相当于 Φ(t)中的2πf_c,相位控制字相当于Φ(t)中的θ_0。

由于相位累加器字长的限制，相位累加器累加到一定值后，其输出将会溢出，这样波形存储器的地址就会循环一次，即意味着输出波形循环一周。故改变频率控制字即相位增量，就可以改变相位累加器的溢出时间，在时钟频率不变的条件下就可以改变输出频率。改变查表寻址的时钟频率，同样也可以改变输出波形的频率。

为了获得较高的频率分辨率，则只有增加相位累加器的字长N，故一般N都取值较大。但是受存储器容量的限制，存储器地址线的位数w不可能很大，一般都要小于N。这样存储器的地址线一般都只能接在相位累加器输出的高w位，而相位累加器输出余下的(N-W)个低位则只能被舍弃，这就是相位截断误差的来源。

DDS模块的输出频率f _out是系统工作频率f _c、相位累加器位数N及频率控制字K满足如下关系

频率分辨率，即频率的变化间隔

三、实现代码

利用matlab或者Guagle_wave工具生成波形文件，存入ROM。基于Quartus II平台，并且调用了ROM核。

   
    
   
   module DDS ( 
    
     
    
                  sys_clk,
    
     
    
                  sys_rst_n,
    
     
    
                  freq_word,
    
     
    
                  phase_word,
    
     
    
                  wave_out
    
     
    
                  );
    
     
    
    
    
     
    
    		input        sys_clk        ;   //系统工作时钟
    
     
    
    		input        sys_rst_n      ;   //复位，低有效
    
     
    
    		input [31:0] freq_word      ;   //频率控制字，控制输出波形频率
    
     
    
    		input [11:0] phase_word     ;   //相位控制字，控制初始相位
    
     
    
    		                         
    
     
    
    		output[9:0]  wave_out       ;   //输出波形,位宽10bit
    
     
    
    		                     
    
     
    
    		reg   [31:0] freq_word_reg  ;
    
     
    
    		reg   [11:0] phase_word_reg ;
    
     
    
    
    
     
    
    		reg   [31:0] phase_adder    ;   //相位累加器
    
     
    
    		reg   [9:0]  rom_address    ;   //存储深度2^10
    
     
    
    
    
     
    
    		always @(posedge sys_clk or negedge sys_rst_n) begin 
    
     
    
    				  if (sys_rst_n ==1'b0) begin 
    
     
    
    						freq_word_reg <= 32'h0000;
    
     
    
    				  end
    
     
    
    				  else  
    
     
    
    						freq_word_reg <= freq_word;
    
     
    
    		end
    
     
    
    
    
     
    
    		always @(posedge sys_clk or negedge sys_rst_n) begin 
    
     
    
    				  if (sys_rst_n ==1'b0) begin 
    
     
    
    						phase_word_reg <= 12'h0000;
    
     
    
    				  end
    
     
    
    				  else  
    
     
    
    						phase_word_reg <= phase_word;
    
     
    
    		end
    
     
    
    
    
     
    
    		always @(posedge sys_clk or negedge sys_rst_n) begin 
    
     
    
    				  if (sys_rst_n ==1'b0) begin 
    
     
    
    						phase_adder <= 32'h0000;
    
     
    
    				  end
    
     
    
    				  else  
    
     
    
    						phase_adder <= phase_adder + freq_word_reg;        //对频率控制字进行累加
    
     
    
    		end
    
     
    
    
    
     
    
    		always @(posedge sys_clk or negedge sys_rst_n) begin 
    
     
    
    				  if (sys_rst_n ==1'b0) begin 
    
     
    
    						rom_address <= 32'h0000;
    
     
    
    				  end
    
     
    
    				  else  
    
     
    
    						rom_address <= phase_adder[31:20] + phase_word_reg;       
    
     
    
    		end
    
     
    
    
    
     
    
    		sin_rom DDS_ROM_U0 (
    
     
    
    							.address (rom_address) ,
    
     
    
    							.clock   (sys_clk)  ,  
    
     
    
    							.q       (wave_out)   
    
     
    
    							);
    
     
    
    
    
     
    
    endmodule

测试文件代码

   
    
   
   `timescale 1 ns/ 1 ns
    
     
    
    module DDS_tb();
    
     
    
    	reg [31:0] freq_word  ;
    
     
    
    	reg [11:0] phase_word ;
    
     
    
    	reg        sys_clk    ;
    
     
    
    	reg        sys_rst_n  ;
    
     
    
    																
    
     
    
    	wire [9:0]  wave_out  ;
    
     
    
    									
    
     
    
    	DDS i1 (
    
     
    
    		.freq_word(freq_word),
    
     
    
    		.phase_word(phase_word),
    
     
    
    		.sys_clk(sys_clk),
    
     
    
    		.sys_rst_n(sys_rst_n),
    
     
    
    		.wave_out(wave_out)
    
     
    
    	       );
    
     
    
    	
    
     
    
    	initial                                                
    
     
    
    		begin                                                  
    
     
    
    			sys_clk    = 0;
    
     
    
    			sys_rst_n  = 0;
    
     
    
    			freq_word  = 0;
    
     
    
    			phase_word = 0;
    
     
    
    			#2000;
    
     
    
    			sys_rst_n  = 1; 
    
     
    
    			freq_word  = 32'd1024;   //25M
    
     
    
    			phase_word = 32'd0;      //相位0度                                  
    
     
    
    			#200000000;
    
     
    
    			freq_word  = 32'd2048;   //50M
    
     
    
    			phase_word = 32'd512;    //相位90度
    
     
    
    			#200000000;
    
     
    
    			$stop;                  
    
     
    
    		end  
    
     
    
    	always sys_clk = #20 ~sys_clk;   //系统时钟25M                                        
    
     
    
    endmodule

仿真波形如下：

改变ROM中的波形数据就可以输出相对应的波形。

参考资料： A Technical Tutorialon Digital Signal Synthesis

无线通信FPGA设计

来自为知笔记(Wiz)

来源：oschina

链接：https://my.oschina.net/u/4360442/blog/4313902

标签

matlab

Quartus II