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

基于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的设置就可以达到改变地址的变化频率,其实就是输出的频率。不知道这样说看者能否理解。

这里相位累加器位数为 N 位(N 的取值范围实际应用中一般为 24~32),相当于把正弦信号在相位上的精度定义为 N 位,所以其分辨率为1 /2𝑁 。 若 DDS 的时钟频率为𝐹𝑐𝑙𝑘,频率控制字 fword 为 1,则输出频率为𝐹𝑜𝑢𝑡 = 𝐹𝑐𝑙𝑘/2𝑁 ,这个频率相当于“基频”。若 fword 为 B,则输出频率为𝐹𝑜𝑢𝑡 = 𝐵 × 𝐹𝑐𝑙𝑘/2𝑁 。因此理论上由以上三个参数就可以得出任意的𝑓𝑜输出频率。且可得出频率分辨率由时钟频率和累加器的位数决定。当参考时钟频率越高,累加器位数越高,输出频率分辨率就越高。
代码:
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值太慢了。mif_malker2010.exe下载地址:http://pan.baidu.com/s/1bCqAp4

评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值