基于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<<
