本系列的10~13篇介绍了数字通信中ASK调制解调技术的实现。本文将介绍数字通信系统中FSK调制技术的FPGA实现，后面的文章将介绍FSK解调。

FSK调制

FSK是利用载波的频率表示基带信息，比如在2FSK系统中，采用单频信号f1表示信息0、单频信号f2表示信息1。根据码元转换时载波的相位是否连续，分为非连续相位FSK和连续相位FSK，如下图所示：

第三幅图为连续相位FSK调制，也称作CPFSK，可视作振荡频率随基带信号线性变化；第四幅图为非连续相位FSK，在码元转换时有一个相位的跳变，可视作两个ASK调制信号的叠加。

FSK的一个重要参数是频移指数h=(f2-f1)*T，T为码元周期。该参数影响着两个单频信号之间的相关系数ρ。ρ越小，越容易将两个信号区分开，即解调系统的解调性能越好。h的几个特殊取值如下：

当h=0.715时，ρ可以取到最小值-2π/3；
当h=kπ(k≥1)时，ρ为0，即两信号正交；
当h=N(整数)时，FSK总是会体现出连续相位特性。

调制信号的频谱与相位变化息息相关，如果存在相位突变，信号频谱宽度会扩展（产生旁瓣）。CPFSK是一种恒包络调制技术，在码元速率和h相同时，CPFSK比非连续相位的FSK主瓣宽度更窄，且旁瓣分量更小，通过限带滤波后信号包络起伏小，适合应用于非线性信道。

MATLAB设计

MATLAB中可以按照FSK原理完成FSK调制，但是只能产生非连续相位的FSK调制信号（除非h=1）。如果想要产生连续相位的FSK调制信号，必须使用MATLAB提供的fskmod函数，示例代码如下：

M = 4;        % Modulation order
freqsep = 8;  % Frequency separation (Hz)
nsamp = 8;    % Number of samples per symbol
Fs = 32;      % Sample rate (Hz)% Generate random M-ary symbols.
x = randi([0 M-1],1000,1);
% Apply FSK modulation.
y = fskmod(x,M,freqsep,nsamp,Fs);

上面将代表0和1两个频率之间的间隔直接设置为8，也可以根据频移指数h的计算公式来设定。fskmod最后还可以增加一个参数，默认值为’cont’表示连续相位；’discont’表示非连续相位。

需要注意的是，老版本的MATLAB可能会用到randint函数产生随机序列；而2017以后的版本只能使用randi函数。

FSK的FPGA设计

在Vivado开发环境下完成2FSK（非连续相位）调制技术的设计并进行仿真，系统时钟32Mhz，码元速率1Mhz，载波频率6Mhz，频移指数h=3.5。模块接口如下：

`timescale 1ns / 1ps
//-----------------------------------------------------
//   FSK调制，非连续相位，频率选择法，h=3.5
//-----------------------------------------------------
module FSK_liuqi
(input clk,     //32MHz采样时钟input rst,  input din,     //码元速率1MHzoutput [15:0] dout
);

配置两个DDS产生两个代表1和0的单频信号。f2-f1=h/T=3.5/1ns=3.5Mhz，则在6Mhz载波频率的基础上选择4.25Mhz和7.75MHz。DDS Compiler IP核使用方法可参考https://blog.csdn.net/fpgadesigner/article/details/80512067。本设计将DDS的无杂散动态范围设置为95dB，对应的位宽为16bit。实例化代码如下：

//-----------------------------------------------------
//   载波Fc=6MHz，h=3.5，则f1=4.25MHz、f2=7.75Mhz
//-----------------------------------------------------
wire [15:0] f1, f2;dds_f1 dds_f1 (            //4.25MHz.aclk(clk),                             .m_axis_data_tvalid(), .m_axis_data_tdata(f1)
);
dds_f2 dds_f2 (            //7.75MHz.aclk(clk),                             .m_axis_data_tvalid(), .m_axis_data_tdata(f2)
);

FSK调制信号的产生有两种方法。第一种方法将FSK视作两路ASK信号的合成（一路ASK有载波代表1，一路ASK无载波代表1），代码如下：

//-----------------------------------------------------
//   方法一：两路ASK合成法
//-----------------------------------------------------
wire din_inv = ~din;
reg [15:0] dout_reg;always @ (posedge clk or posedge rst)if (rst) dout_reg <= 0;else dout_reg <= din*f1 + din_inv*f2;assign dout =  dout_reg;

采用此方法的系统原理图如下：

可以看到一路ASK信号为din与f1的乘法，另一路信号为din取反后与f2的乘法，两路加起来便得到FSK调制信号。另一种方法直接使用多路选择器对输入的基带数据进行判决输出：

//-----------------------------------------------------
//   方法二：多路选择器（Mux）判决
//-----------------------------------------------------
reg [15:0] dout_reg;always @ (posedge clk or posedge rst)if (rst) dout_reg <= 0;else case (din)1'b1 : dout_reg <= f1;1'b0 : dout_reg <= f2;default : dout_reg <= 'd0;endcaseassign dout =  dout_reg;

采用此方法的系统原理图如下：

RTL_MUX即为多路选择器，明显第二种方法更加节省资源。

CPFSK的FPGA设计

连续相位的FSK调制信号产生可以借助DDS的特性，DDS的实现原理就使其输出信号具有连续相位特性，即使频率突然改变。配置一个DDS，Implementation->Phase Increment Programmability设置为Streaming，可以通过S_AXIS_PHASE接口动态配置DDS的输出频率。完整代码如下：

`timescale 1ns / 1ps
//-----------------------------------------------------
//   CPFSK调制，连续相位，h=3.5
//-----------------------------------------------------
module CPFSK_liuqi
(input clk,     //32MHz采样时钟input rst,  input din,     //码元速率1MHzoutput [15:0] dout
);//-----------------------------------------------------
//   载波Fc=6MHz，h=3.5，则f1=4.25MHz、f2=7.75Mhz
//-----------------------------------------------------
wire [14:0] phase;dds_compiler_0 dds (.aclk(clk),                                .s_axis_phase_tvalid(1'b1), .s_axis_phase_tdata({1'b0,phase}),    .m_axis_data_tvalid(),   .m_axis_data_tdata(dout)
);assign phase = din ? 'd4351 : 'd7935;endmodule

DDS的频率分辨率设置为976.5625Hz，正好对应相位位宽为15Bits（Fs/2^B=32MHz/2^15=976.5625Hz）。因此f1=4.25MHz对应的频率控制字为4351（f1/976.5625-1=4351）；f2=7.75MHz对应的频率控制字为7935（f2/976.5625-1=7935）。

代码中根据输入基带的情况设定对应的频率控制字。可见CPFSK只需要使用一个DDS Compiler IP核即可，比FSK要节省不少资源。

仿真测试

编写testbench，产生一个周期信号模拟二进制基带数据。首先仿真非连续相位的FSK：

可以明显看到在码元的跳变时刻FSK信号的相位有明显的突变。再仿真CPFSK：

看到码元跳变时刻FSK信号的相位仍保持连续变化。注意由于DDS有一定的Latency，导致上图FSK信号相比基带数据有一定的延迟。

Quartus环境下实现

Quartus中对应的IP核为NCO，在产生FSK信号时与DDS用法基本相同；但在产生CPFSK信号时，NCO IP核有独特的配置方式。

NCO的输出频率由phi_inc_i接口控制。需要使用FSK调制功能时，在Implementation标签下选中“Frequency Modulation Input”：

IP核会多出一个freq_mod_i接口，该接口的作用时在phi_inc_i输出频率的基础上进行调整。比如：频率控制字位宽为25Bits，phi_inc_i输出频率6Mhz，系统时钟32MHz，需要输出f1=4.25Mhz，则freq_mod_i对应的控制字为(4.25Mhz-6MHz)/(32Mhz/2^25)=-1835008。

根据基带数据改变freq_mod_i接口的控制字即可实现CPFSK信号的产生。

FPGA数字信号处理（22）FSK调制技术相关推荐

FPGA 数字信号处理之 FSK 调制、解调的实现与仿真基于 verilog + ise + modelsim + matlab （保姆级）
试验说明: 笔者复现FSK调制.解调的过程中,在网上查了很多很多相关文章,至少在本文之前,没发现能在参考一篇博文的情况下能完整复现实验结果的.总之,全乎的基本没有. 本笔记(文)的目标是让入门者或者说 ...
FPGA数字信号处理（23）FSK解调技术（包络检波法）
上一篇介绍了数字通信系统中FSK调制技术的FPGA实现.调制信号经过DAC.可选的带通滤波器.功率放大器.天线发送出去后,在接收端收到FSK信号后需要对其解调,提取出包含的信息(基带信号). FSK解 ...
FPGA数字信号处理（十一）ASK解调技术
上一篇介绍了数字通信系统中ASK调制技术的FPGA实现.调制信号经过DAC.可选的带通滤波器.功率放大器.天线发送出去后,在接收端收到ASK信号后需要对其解调,提取出包含的信息(基带信号).无论在哪种 ...
FPGA数字信号处理（1)- AM调制的FPGA实现
FPGA数字信号处理(1)-AM调制的实现一:前言本内容分享为本人自学经历.受限于作者水平可能有不准确的地方.欢迎诸位批评指正. 分享的文章需要一些基本的FPGA开发基础二:概述这部分简单,但 ...
FPGADesigner《FPGA数字信号处理系列》目录与传送门
FPGA数字信号处理(1)数字混频(NCO与DDS的使用): https://blog.csdn.net/fpgadesigner/article/details/80512067 FPGA数字信号处 ...
新书预告：Xilix FPGA数字信号处理设计——基础版
掌握FPGA数字信号处理设计需满足三个条件:熟悉FPGA设计方法.理解数字信号处理理论.掌握理论的工程实现方法.对初学者来讲,每个条件看似都难以逾越.杜勇老师完美融合课程教学与工程设计的需求特点,以独 ...
FPGA数字信号处理（一）数字混频（NCO与DDS的使用）
这是数字信号处理系列的第一篇,以简单的数字混频为例,介绍在FPGA程序设计中很重要的二进制原码.补码:有符号数.无符号数的问题.本文不是像课本那样介绍这些基础概念,而是介绍很实际的设计方法. 借助于数 ...
FPGA数字信号处理（二）并行FIR滤波器Verilog设计
该篇是FPGA数字信号处理的第二篇,选题为DSP系统中极其常用的FIR滤波器.本文将简单介绍FIR滤波器的原理,详细介绍使用Verilog HDL设计并行FIR滤波器的流程和方法.接下来几篇会介绍串行 ...
FPGA数字信号处理（六）直接型IIR滤波器Verilog设计
该篇是FPGA数字信号处理的第六篇,2-5篇介绍了DSP系统中极其常用的FIR滤波器.本文将简单介绍另一种数字滤波器--IIR滤波器的原理,详细介绍使用Verilog HDL设计直接型IIR滤波器的方 ...

FPGA数字信号处理（22）FSK调制技术