基于Verilog的DDS波形发生器的分析与实现（三角波、正弦波）

　　最近学习了一下关于DDS的相关知识，本篇概要记录一下自己的理解与实现。

　　DDS信号发生器采用直接数字频率合成(Direct Digital Synthesis，简称DDS)技术，把信号发生器的频率稳定度、准确度提高到与基准频率相同的水平，并且可以在很宽的频率范围内进行精细的频率调节。采用这种方法设计的信号源可工作于调制状态，可对输出电平进行调节，也可输出各种波形。

　　下图为DDS 的基本结构图

　　由上图 可以看出，DDS 主要由相位累加器、相位调制器、波形数据表以及 D/A 转换器构成。

　　相位累加部分控制输出波形频率，相位字输入部分来改变相位，ROM表中存储一个周期波形的幅度值。

　　其中相位累加器由 N 位加法器与 N 位寄存器构成。每来一个时钟，加法器就将频率控制字与累加寄存器输出的相位数据相加，相加的结果又反馈至累加寄存器的数据输入端，以使加法器在下一个时钟脉冲的作用下继续与频率控制字相加。这样，相位累加器在时钟作用下，不断对频率控制字进行线性相位累加。即在每一个时钟脉冲输入时，相位累加器便把频率控制字累加一次。相位累加器输出的数据就是合成信号的相位。相位累加器的溢出频率，就是 DDS 输出的信号频率。（解释：定义一个N位寄存器，一般为32位，如果来一个时钟计一次，那就要计2^32次才满，这样太慢，因此引入频率控制字设为A，以前以1为单位，现在以A为单位计数累加，可以控制计数更新的速度）

　　用相位累加器输出的数据，作为波形存储器的相位采样地址，这样就可以把存储在波形存储器里的波形采样值经查表找出，完成相位到幅度的转换。波形存储器的输出送到 D/A 转换器，由 D/A 转换器将数字信号转换成模拟信号输出。

　　一般32位累加器不会全用来输出作为ROM地址，会根据ROM深度来适当截取高位作为地址，其余位可以作为控制频率。例如现在ROM中存储波形一个周期数据每个数据位宽8位，则数据范围为0-2^8（256），但是要产生一个周期波形需要512个，因为0-256一般是上升期，256-0处于下降期，那么现在ROM深度为512，则地址位宽应为9位，2的9次方=512，则32位累加器只需高9位即可[31：23]寻址，其余位用来控制地址改变的频率。假如现在每来一个时钟地址变一次，那么其余位(即频率控制字A)应设置为32‘h800000即32’b0000_0000_10000000_0000_0000_0000_0000，最高位1其实就是地址的最低位，累加器初始为0，来个时钟沿加一次A，高9位地址变化一次。要想两个时钟变化一次，那么A就是32‘b0000_00000100_0000_0000_0000_0000_0000，两个时钟之后地址的最低位才会变化。这样通过对频率控制字A的设置就可以达到改变地址的变化频率，其实就是输出的频率。不知道这样说看者能否理解。

　　代码：

　　module dds(

　　clk,

　　rst,

　　Fword, //频率控制字A

　　Pword, //相位控制字

　　data

　　);

　　input clk,rst;

　　input [31:0] Fword;

　　input [8:0] Pword;

　　output [7:0] data; //8位数据

　　reg [31:0] r_Fword;

　　reg [8:0] r_Pword;

　　reg [31:0] cnt;

　　wire [8:0] rom_adder;

　　always@(posedge clk)

　　begin

　　r_Fword<=Fword;

　　r_Pword<=Pword;

　　end

　　always@(posedge clk or negedge rst)    //累加器部分

　　begin

　　if(!rst)

　　cnt<=32'd0;

　　else

　　cnt<=cnt+r_Fword;

　　end

　　assign rom_adder=cnt[31:23]+r_Pword;     //ROM地址

　　rom  rom(                              //例化一个ROM ip核

　　.address(rom_adder),

　　.clock(clk),

　　.q(data)

　　);

　　endmodule

　　tb:

　　`timescale 1ns/1ns

　　module dds_tb;

　　reg clk,rst;

　　reg [31:0] Fword;

　　reg [8:0] Pword;

　　wire [7:0] data;

　　dds u0(

　　.clk(clk),

　　.rst(rst),

　　.Fword(Fword),

　　.Pword(Pword),

　　.data(data)

　　);

　　initial clk=1;

　　always #10 clk=~clk;

　　initial begin

　　rst=0;

　　Fword=32'h800000;

　　Pword=9'd0;

　　#101;

　　rst=1;

　　#500;

　　$stop;

　　end

　　endmodule

　　ROM中值，用来比对，在例化ip是输出加了一个寄存器，因此会延时一拍输出。主时钟50M,下图把data转换为模拟值，两个黄线之间为频率97.65khz，A为32‘h800000=32’d8388608，50_000_000乘以A等于419430400000000，除以2^32等于97656.25khz验证正确。

　　更改ROM中mif文件，添加三角波，ROM中三角波这里一共有256个数值，0-127,127-0，地址位宽为8位，那么[31:24]作为地址，其余位作为频率控制字。只需改动程序中地址位宽就可以了，A暂时设为32‘h01000000，下图为相位控制字为8‘h0,8’h5；可以看出相位控制字不会改变输出频率.

　　附加：如果此时知道主时钟50M，想生成一个25M的波形，那么首先根据公式计算出A=2147483648即二进制1000_0000_0000_0000_0000_0000_0000_0000.仿真如下图为方波，因为根据地址变化只取0、127.

　　总结：本文内容也是在看相关视频后自己更改ROM深度，更换波形仿真后得出，由于叙述水平有限，其中原理可以自行查看其他文章内容了解，关于本文如有不懂之处可以联系我，共同再探讨。同时其中主要思想可以用来作为任意分频设计，后续再做。

　　文末附上一个mif文件生成器，不然手动输入ROM值太慢了

　更多信息可以来这里获取==>>电子技术应用-AET<<