在FPGA中实现DDS方案详解(频率,幅度,波形可调)

在FPGA中实现DDS方案详解(频率,幅度,波形可调)

目录

1. DDS原理简介

2. 代码实现

2.1 控制模块

2.2 DDS实现模块

2.3 ip核ROM的实现

3. 实验结果

3.1 仿真结果

3.2 实际输出结果

3.2.1波形选择

3.2.2 频率控制:

3.2.3 幅度控制


1. DDS原理简介

        DDS技术是从波形相位概念出发,直接对相应的波形数据进行抽样,得到不同的相位,通过DAC转换成模拟波形,最后通过低通滤波器平滑输出所需频率的波形。下图为DDS的示意流程图。


两个重要参数解释:      

         频率控制字FTW:一般为整数,数值大小控制输出信号的频率大小,数值越大输出信号频率越高,实际输出的信号的频率与频率控制字关系为:

频率控制字 = 脉冲频率 * 2^N / 时钟频率

        相位控制字P_WORD:一般为整数,数值大小控制输出信号的相位偏移,主要用于相位的信号调制


一个普通的DDS设计应该由:相位累加器,相位调制器,波形数据表,DAC转换组成。

        相位累加器:是整个DDS的核心,用于控制波形的相位累加,组成完整的波形显示。

相位累加器的输入为频率控制字,表示相位增量,其位宽为N。

        相位调制器:接收相位累加器输出的相位码,在这里加上一个相位偏移值P,主要用于信号的相位调制。

        波形数据表:波形数据表ROM中存有一个完整周期的正弦波信号。举个例子:假设波形数据ROM的地址位宽为12位,储存数据位宽为8位,即ROM有

2^{12}=4096

个存储空间,每个存储空间可以存储1字节数据。将一个周期的正弦波信号,沿横轴等间隔采样4096次,每次采集的信号幅度用1字节数据表示,最大值为255,最小值为0。将4096次采样结果按顺序写入ROM的4096个存储单元,一个完整周期正弦波的数字幅度信号就写入了波形数据表ROM中。波形数据表ROM以相位调制器传入的相位码为ROM读地址,将地址对应存储单元的电压幅值数字量进行输出。

        DAC转换:将输入的电压幅值数字量转换为模拟量输出。由上个公式可得,当频率控制字FTW为1时,此时输出信号频率最低。根据奈奎斯特采样定理,频率控制字的最大值应小于 2^N / 2 

运行流程为:


相位累加器得到的相位码是如何实现ROM寻址的呢?

        对于N位的相位累加器,其对应的相位累加值为2^N,如果正弦ROM中存储单元的个数也是2^N的话,那这个问题就很容易解决,但是这对ROM的存储容量要求较高。

        那么在实际中,我们使用相位累加器的高几位来对ROM进行寻址,也就是说并不是每个系统时钟都对ROM进行数据读取,而是多个时钟读取一次,因为这样能保证相位累加器溢出时,从正弦ROM表中取出正好一个正弦周期的样点。

        因此,相位累加器每计数2^N次,对应一个正弦周期。而相位累加器1秒钟计数fclk次(工作时钟频率),当频率控制字FTW为1时,DDS输出的时钟频率就是频率分辨率。当FTW增加时,相位累加器溢出的频率增加,对应DDS输出信号频率就变成了FTW倍的DDS频率分辨率。

        举个例子:

        假设ROM的存储单元个数为4096,每个存储数据用8位二进制表示。即ROM地址线宽为12,数据线宽为8;相位累加器位宽设置为N=32。

        由于ROM地址线位宽为12,所以相位调制器的输出位宽也应该为12。(因为DDS原理就是通过相位调制器的输出去ROM寻址)在相位调制器中与相位控制字进行累加时,用相位累加器的高12位累加,而低位只与频率控制字累加

        以频率控制字FTW=1为例,相位累加器的低20位一直加1,直到低20位溢出而向高12位进位,此时ROM为0,也就是说,ROM的0地址中的数据被读了2^20次,继续往下,ROM的4096个点全部都会读2^20次,最终输出的波形频率应该是参考时钟频率的1/2^20,周期被扩大了2^20倍。同样当频率控制字为100时,所以最终输出频率是上述的100倍。


2. 代码实现

        本次实验的硬件平台为小梅哥AC101-EDA,开发软件为Quartus13.1。

        由上节对DDS原理的介绍,可以分为以下几个模块。        


2.1 控制模块

频率控制模块、相位控制模块、幅度控制模块这三个模块的逻辑一模一样,即:指定一个按键,通过按键的值将对应的参数进行改变。下面以幅度控制模块进行举例说明。

当系统检测到按键被按下后,将变量cnt进行计数:

即用cnt来表示按键的按下的次数,这里是当cnt按下5次后,被清零。其中key_flag表示按键按下,由按键消抖模块输出。

        amplitude的值由cnt的值决定,用一个case语句即可实现。

 其他两个模块几乎一模一样,这里不再赘述。


2.2 DDS实现模块

还需要设计一个波形选择的功能,可按如下实现。

其中dac_data0、dac_data1、dac_data2、dac_data3分别对应四个rom的输出。


2.3 ip核ROM的实现

在Quatrtus13.1中,rom使用ip核进行实现:由于实验所用的开发板的DAC为8位,所以ROM的数据位宽这里也设置为8。

总共使用4个ROM,分别对应正弦波、三角波、方波、锯齿波,ROM所需的.mif文件由matlab实现。如正弦波文件的matlab代码如下:

clc; %清除命令行命令 clear all; %清除工作区变量,释放内存空间 F1=1; %信号频率 Fs=2^12; %采样频率 P1=0; %信号初始相位 N=2^12; %采样点数 t=[0:1/Fs:(N-1)/Fs]; %采样时刻 ADC=2^7 - 1; %直流分量 A=2^7; %信号幅度 %生成正弦信号 s=A*sin(2*pi*F1*t + pi*P1/180) + ADC; plot(s); %绘制图形 %创建mif文件 fild = fopen('sin_wave_4096x8.mif','wt'); %写入mif文件头 fprintf(fild, '%s\n','WIDTH=8;'); %位宽 fprintf(fild, '%s\n\n','DEPTH=4096;'); %深度 fprintf(fild, '%s\n','ADDRESS_RADIX=UNS;'); %地址格式 fprintf(fild, '%s\n\n','DATA_RADIX=UNS;'); %数据格式 fprintf(fild, '%s\t','CONTENT'); %地址 fprintf(fild, '%s\n','BEGIN'); %开始 for i = 1:N s0(i) = round(s(i)); %对小数四舍五入以取整 if s0(i) <0 %负1强制置零 s0(i) = 0 end fprintf(fild, '\t%g\t',i-1); %地址编码 fprintf(fild, '%s\t',':'); %冒号 fprintf(fild, '%d',s0(i)); %数据写入 fprintf(fild, '%s\n',';'); %分号,换行 end fprintf(fild, '%s\n','END;'); %结束 fclose(fild);

这里需要注意一点的是,如果在modelsim进行仿真时,rom输出一直为0,是因为没有把生成的.mif文件放在quartus的工程目录下。具体可以看一下我的这篇文章modelsim无法读取.mif文件

但是是否对实际输出波形有影响不知道,我没有实践过。

DAC驱动模块:

        开发板上的DAC是MCP4802,采用的SPI通信。代码部分参考这篇文章基于FPGA的目标反射回波检测算法及其实现(准备篇)

在Quartus中查看RTL视图:


3. 实验结果

3.1 仿真结果

3.2 实际输出结果

3.2.1波形选择

3.2.2 频率控制

3.2.3 幅度控制


若有不懂,不清楚的地方或有相关建议,欢迎留言!


Read more

从零开始打造高性能数据结构——手把手教你实现环形缓冲

从零开始打造高性能数据结构——手把手教你实现环形缓冲

◆ 博主名称: 小此方-ZEEKLOG博客 大家好,欢迎来到小此方的博客。 ⭐️个人专栏:《C语言》_小此方的博客-ZEEKLOG博客 算法_小此方的博客-ZEEKLOG博客  ⭐️踏破千山志未空,拨开云雾见晴虹。 人生何必叹萧瑟,心在凌霄第一峰。 目录 一,普通队列的劣势 1. 空间浪费严重(“假溢出”问题) 2. 需要频繁移动元素(若避免浪费) 3. 扩容成本高 4. 无法解决“假溢出”导致的提前扩容 二,环形缓冲结构分析  1. “循环”取模实现指针回绕  2.“循环”,轮流入座而不是排长队 三,实现环形缓冲 1,MyCircularQueue(k): 构造器   1,结构体搭建   2,初始化 3,为什么选择k+1块空间而不是k块空间?

By Ne0inhk
【动态规划篇】专题(六):子序列问题——不连续的艺术

【动态规划篇】专题(六):子序列问题——不连续的艺术

文章目录 * LIS 模型及其衍生:回头看,全是风景 * 一、 前言:从 O(N) 到 O(N²) * 二、 最长递增子序列 (Medium) * 2.1 题目描述 * 2.2 核心思路:LIS 模型 * 2.3 代码实现 * 三、 摆动序列 (Medium) * 3.1 题目描述 * 3.2 状态定义:波峰与波谷 * 3.3 代码实现 * 四、 最长递增子序列的个数 (Medium) * 4.1 题目描述 * 4.2 双重状态 * 4.

By Ne0inhk
Python | 栅格数据Sen+MK长时间序列趋势分析+显著性检验保姆级教程

Python | 栅格数据Sen+MK长时间序列趋势分析+显著性检验保姆级教程

遥感代码星球的第002篇代码分享  今天给大家分享的是: Theil-Sen Median趋势分析 + Mann-Kendall显著性检验   基于遥感数据的植被碳水利用效率时空变化和归因分析[J].生态学报,2024,44(01):377-391.DOI:10.20103/j.stxb.202207232112. 在研究气候变化、环境变化、生态变化等领域时,如何准确分析时间序列数据的趋势是一个重要的课题。Theil-Sen Median趋势分析(简称Sen分析)和Mann-Kendall显著性检验(简称MK检验)是两种常用的非参数方法,它们在处理气候、环境、生态等领域的趋势分析中有着广泛的应用。本文将详细介绍这两种方法的基本原理,并探讨它们如何结合使用以提高趋势分析的准确性。同时结尾附完整代码及案例数据供大家学习使用。 01 原理及公式 1. Theil-Sen Median 趋势分析(Sen分析) Theil-Sen Median 趋势分析是一种稳健的非参数统计方法,适用于具有异常值或缺失值的时间序列数据。该方法通过计算数据的中位数斜率(即Sen斜率

By Ne0inhk
【数据结构和算法】面试必刷之随机链表复制:这三步让你彻底吃透 random 指针

【数据结构和算法】面试必刷之随机链表复制:这三步让你彻底吃透 random 指针

🔥小龙报:个人主页 🎬作者简介:C++研发,嵌入式,机器人等方向学习者 ❄️个人专栏:《C语言》《【初阶】数据结构与算法》 ✨ 永远相信美好的事情即将发生 文章目录 * 前言 * 一、随即链表的复制 * 1.1 题目 * 1.2 算法原理 * 1.3 代码 * 总结与每日励志 前言 随机链表的复制是数据结构中的经典难题,核心难点在于复制节点的random指针——其指向的节点可能尚未创建,也可能指向链表中的任意节点。本文采用“原地拷贝+拆分”的最优思路,分三步拆解解题逻辑,结合代码实现与原理分析,清晰讲解如何高效解决该问题,帮助读者吃透random指针的处理技巧,掌握链表操作的核心思维。 一、随即链表的复制 1.1 题目 链接:随机链表的复制 1.2 算法原理

By Ne0inhk