我是靠谱客的博主 安详口红,最近开发中收集的这篇文章主要介绍FPGA学习积累之AM调制解调(解调部分没搞太明白)随便说说一、平台二、DDS IP的使用三、FIR滤波(用于后面AM解调)四、AM调制五、整体代码六、仿真波形总结,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

基于FPGA的AM调制与解调以及从中学到的东西(注:解调没搞太清楚,没写)

  • 随便说说
  • 一、平台
  • 二、DDS IP的使用
  • 三、FIR滤波(用于后面AM解调)
  • 四、AM调制
  • 五、整体代码
  • 六、仿真波形
  • 总结


随便说说

(这次小杨第一次写博客,有问题,见谅!)
前几天没事做就想学习学习信号类的知识,然后就想到了去年电赛折磨我的AM与FM的调制与解调。AM的调制与解调用FPGA来实现还是不是很难,但是我还是看了两位博主的博客(链接在这一段结尾)。本次AM调制与解调的信号都是有DDS输出的,载波为1MHz,调制波为10KHz。
AM调制与解调:(博主:子木呀)
AM调制与解调:(博主:黄子炫)


一、平台

1.软件:Vivado 2018.3(就只会这个。。。。)
2.硬件:EGO1


二、DDS IP的使用

(这里加个链接(DDS的配置))

注意左边的RTL

个人能力有限,就只能看懂这四个地方(/手动糊脸/)。

1.Configuration Options(配置选项):配置DDS的类型,下来有三个选项,这里选着“Phase Generator and SIN Cos LUT”。

2.System Clock(配置DDS的输入系统时钟):对应的引脚(看左边的那个RTL,由aclk输入),这里配置10.24MHz(这个值后面会解释)。

3.Number of Channels(通道数):因为只用到一个输出,所以通道数为1。(其实是自己不会多通道的使用。。。。)

4.Parameter Selection(参数选择):这里我把选项下拉出来了。为了后面使用的方便,这里我选择“Hardware Parameters”

在这里插入图片描述
在这里插入图片描述
这里我把Hardware Parameters和System Parameters的两个界面给截图出来对比了一下。System Parameters需要去计算频率分辨率,这个计算起来就。。。。(我算不明白)。Hardware Parameters就很方便,它的Phase Width默认值为16,Output Width默认值为8。Phase Width的值会影响到我们的输入频率,我使用的时候配的是10位宽。下面展示Phase Width与输入频率的关系。(为了省事,我把一位博主写的部分搬过来,感谢子木呀!)

DDS核频率分辨率计算公式如下:

其中Δθ就是我们的频率控制字取值后文用fre_word表示,Bθ(n)是相位累加器位宽后文用B表示,fclk为DDS工作时钟。频率分辨率就是当Δθ=1时的fout。先假如我们需要的信号频率范围:1M-10MHz,分辨率0.01MHz。
取Δθ=1,fout=0.01MHz,B=10bit,我们可得DDS工作时fclk=10.24MHz。

这里展示一下Configuration配置完的界面(注意左边的RTL):

在这里插入图片描述
在这里插入图片描述
Implementation界面需要配置这四点。

1.选择可编程阶段增量(我觉得就是频率控制字的可编程。。)

2.相位偏移就不需要

3.输出波形为余弦波 [ AM调制公式:P(t)=(A+M_a cos⁡(t_1 ) )×cos⁡(t_2 ) ]

4.不需要相位输出

在这里插入图片描述
这个地方一开始是报错的,需要我们手动修改。前面我们的Phase Width为10,所以这个地方我们要配置为十位的二进制,这里为0000_0001。当这个界面没有报错的时候,左边的RTL的S_AXIS_CONFIG才会高亮。
在这里插入图片描述
在这里插入图片描述这是两个Summary的结果。

到此,DDS IP的用户界面就配置完了。


三、FIR滤波(用于后面AM解调)

FIR滤波器参数的生成

FIR滤波器的参数我用的MATLAB生成的,这里也顺带学习了一下MATLAB的Filter Designer工具。

在这里插入图片描述
1.采样频率我用的100MHz(后面发现100MHz程序运行太慢了。。。。),采样定律规定的是大于采样信号的2倍

2.这里截止频率我用的20KHz。调制信号为10KHz,所以大于10KHz就行,
但是如果有杂波,杂波频率在10KHz~1MHz之间的话,截止频率的要求就有点多了

3.滤波阶数为127

最后记得点击Design Filter,要不然那个波形不是这个滤波器的波形。
接下来生成.COE文件

在这里插入图片描述
FIR IP
在这里插入图片描述1.选择COE File

2.加载刚刚生成的.COE文件

在这里插入图片描述
1.输入的采样频率为10MHz(上面的100MHz改了,嘿嘿,10MHz我都嫌慢)

2.系统输入时钟

到此FIR的配置就完了。


四、AM调制

AM调制公式:P(t)=(A + M_a × cos⁡(t_1) ) × cos⁡(t_2)
其中当A=1时,调制深度就是M_a;cos⁡(t_1 )为调制信号,cos⁡(t_2)为载波。
为了好解释代码,上式变为(D + depth_con × cos⁡(t_1) ) × cos⁡(t_2)

代码如下(还是嫖别人代码香,嘿嘿,虽然不道德,但是掌握了知识就是自己的了,再次感谢开篇链接的两位博主):

reg  [7:0] D = 127;          //AM波直流分量 
    reg signed[7:0] depth_con;  
 
    always@(posedge clk_out)         //设置调制深度
    begin
       case (depth)
            0:  depth_con <= 0 ;       //调制深度为0,直流分量对应值
            1:  depth_con <= 13;       //调制深度为0.1,直流分量对应值
            2:  depth_con <= 28;       //调制深度为0.2,直流分量对应值
            3:  depth_con <= 45;       //调制深度为0.3,直流分量对应值
            4:  depth_con <= 64;       //调制深度为0.4,直流分量对应值
            5:  depth_con <= 85;       //调制深度为0.5,直流分量对应值
            6:  depth_con <= 110;       //调制深度为0.6,直流分量对应值
            7:  depth_con <= 138;       //调制深度为0.7,直流分量对应值  
            8:  depth_con <= 171;       //调制深度为0.8,直流分量对应值               
            9:  depth_con <= 209;       //调制深度为0.9.,直流分量对应值
            10:  depth_con <= 255;       //调制深度为1,直流分量对应值       
       endcase
    end
    
    wire signed[15:0] modulate_mul8ma;
    wire signed[7:0] modulate_mulma;
    mult_genx mult_floating_point (        
      .CLK(clk_out),  // input wire CLK
      .A(modulate_out),      // input wire [7 : 0] A
      .B(depth_con),      // input wire [7 : 0] B
      .P(modulate_mul8ma)      // output wire [15 : 0] P
    );
    assign modulate_mulma = modulate_mul8ma>>8;
    reg [7:0] modulate_withdc;
    always@(posedge clk_out)
    begin
        modulate_withdc <=  modulate_mulma + D;
    end
    

    mult_gen_0 mult_modulation (        
      .CLK(clk_out),  // input wire CLK
      .A(modulate_withdc),      // input wire [7 : 0] A
      .B(carrier_out),      // input wire [7 : 0] B
      .P(AM_out)      // output wire [15 : 0] P
    );

因为Verilog语言不能像C语言、python等等其他编程语言一样可以直接用浮点数,我们要做一些必要的处理完成浮点数计算的问题,所以depth = depth_con>>8(depth_con / 256)。但是代码并没有直接做这一步就计算,而是先将载波与depth_con相乘之后在>>8送给modulate_mulma。为了保证输出的数据都是16位宽的,这里modulate_mulma还加了一个127。这里的两个乘法器要注意一下他们输入参数是否有无符号位,否则计算出的已调制波有问题。

到此,调制就完了。

五、整体代码

代码如下(还是嫖别人代码香,嘿嘿,虽然不道德,但是掌握了知识就是自己的了,再次感谢开篇链接的两位博主):

`timescale 1ns / 1ps

module AM_generate_top(
    input sysclk
    );
    wire [15:0] AM_out;
    wire [3:0] depth= 9;
    wire [15:0] carrier_con = 10;
    wire [15:0] modulate_con = 1;
    wire [7:0] carrier_out;
    wire [7:0] modulate_out;
    wire clk_out;
    
    clk_div clk_1024k
       (
        // Clock out ports
        .clk_out1(clk_out),     // output clk_out1
       // Clock in ports
        .clk_in1(sysclk));      // input clk_in1
    
    carrier_dds DDS_carrier(
      .aclk(clk_out),                                  // input wire aclk
      .s_axis_config_tvalid(1),  // input wire s_axis_config_tvalid
      .s_axis_config_tdata(carrier_con),    // input wire [15 : 0] s_axis_config_tdata
      .m_axis_data_tvalid(),      // output wire m_axis_data_tvalid
      .m_axis_data_tdata(carrier_out)        // output wire [7 : 0] m_axis_data_tdata
    );
    
    modulation_dds DDS_modulation(
      .aclk(clk_out),                                  // input wire aclk
      .s_axis_config_tvalid(1),  // input wire s_axis_config_tvalid
      .s_axis_config_tdata(modulate_con),    // input wire [15 : 0] s_axis_config_tdata
      .m_axis_data_tvalid(),      // output wire m_axis_data_tvalid
      .m_axis_data_tdata(modulate_out)        // output wire [7 : 0] m_axis_data_tdata
    );
    
    
    reg  [7:0] A= 127;          //AM波直流分量 
    reg signed[7:0] depth_con;  
 
    always@(posedge clk_out)         //设置调制深度
    begin
       case (depth)
            0:  depth_con <= 0 ;       //调制深度为0,直流分量对应值
            1:  depth_con <= 13;       //调制深度为0.1,直流分量对应值
            2:  depth_con <= 28;       //调制深度为0.2,直流分量对应值
            3:  depth_con <= 45;       //调制深度为0.3,直流分量对应值
            4:  depth_con <= 64;       //调制深度为0.4,直流分量对应值
            5:  depth_con <= 85;       //调制深度为0.5,直流分量对应值
            6:  depth_con <= 110;       //调制深度为0.6,直流分量对应值
            7:  depth_con <= 138;       //调制深度为0.7,直流分量对应值  
            8:  depth_con <= 171;       //调制深度为0.8,直流分量对应值               
            9:  depth_con <= 209;       //调制深度为0.9.,直流分量对应值
            10:  depth_con <= 255;       //调制深度为1,直流分量对应值       
       endcase
    end
    (* use_dsp48 = "yes" *)    //这个没太大关系
    
    wire signed[15:0] modulate_mul8ma;
    wire signed[7:0] modulate_mulma;
    mult_genx mult_floating_point (        //调用乘法IP核,就是进行前面说的得到小数的那一部分操作
      .CLK(clk_out),  // input wire CLK
      .A(modulate_out),      // input wire [7 : 0] A
      .B(depth_con),      // input wire [8 : 0] B
      .P(modulate_mul8ma)      // output wire [15 : 0] P
    );
    assign modulate_mulma = modulate_mul8ma>>8;
    reg [7:0] modulate_withdc;
    always@(posedge clk_out)
    begin
        modulate_withdc <=  modulate_mulma + A;
    end
    

    mult_gen_0 mult_modulation (         //与载波相乘
      .CLK(clk_out),  // input wire CLK
      .A(modulate_withdc),      // input wire [7 : 0] A
      .B(carrier_out),      // input wire [7 : 0] B
      .P(AM_out)      // output wire [15 : 0] P
    );

reg [15:0] AM_abs;
always @(posedge clk_out) begin

    if(AM_out[15] == 1)    begin        //全波整流
        AM_abs <= -{AM_out};        //如果符号位是1,对数据取反
    end
    else if(AM_out[15] == 0)    begin
        AM_abs <= AM_out;           //如果符号位是0,数据不变
    end
end

    wire [39:0] demolate_signal;
    fir_my fir (             //低通滤波
      .aclk(sysclk),                              // input wire aclk
      .s_axis_data_tvalid(1),  // input wire s_axis_data_tvalid
      .s_axis_data_tready(),  // output wire s_axis_data_tready
      .s_axis_data_tdata(AM_abs),    // input wire [15 : 0] s_axis_data_tdata
      .m_axis_data_tvalid(),  // output wire m_axis_data_tvalid
      .m_axis_data_tdata(demolate_signal)    // output wire [39 : 0] m_axis_data_tdata
    );

    wire [7:0] demolate_final; 
    assign demolate_final[7:0] = demolate_signal[35:28];  //截位
    
endmodule

六、仿真波形

在这里插入图片描述
这里有两张有意思的:

在这里插入图片描述
在这里插入图片描述其实都是没问题的,就是有个小地方没搞对。

总结

这次AM调制与解调的学习是建立在我没有相关专业基础知识上,所以我最多只能理解到AM调制和解调的算法,但是并没有深入理解到它的内涵。不过在这个过程中学到了如何配置DDS和FIR的IP我觉得也挺值得的,万一明年电赛用到了呢?对吧。看着仿真的波形就得很美妙(特别是那个AM_out),我体会到了FPGA功能的另一种美以及信号处理的乐趣。(个人能力有限,有问题的话请大佬指明出来,我去学习)
(我还是想说一下那个解调:我的理解就是全波整流之后,做低通把那个10KHz的信号给留下来了,不知道是不是这个意思。)

最后

以上就是安详口红为你收集整理的FPGA学习积累之AM调制解调(解调部分没搞太明白)随便说说一、平台二、DDS IP的使用三、FIR滤波(用于后面AM解调)四、AM调制五、整体代码六、仿真波形总结的全部内容,希望文章能够帮你解决FPGA学习积累之AM调制解调(解调部分没搞太明白)随便说说一、平台二、DDS IP的使用三、FIR滤波(用于后面AM解调)四、AM调制五、整体代码六、仿真波形总结所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(37)

评论列表共有 0 条评论

立即
投稿
返回
顶部