Matlab产生波形数据文件

reference: https://blog.csdn.net/HOOKNET/article/details/79254351

前言

在FPGA开发过程中几乎都要用到仿真的功能，对于一些简单的外部激励（如时钟、复位、简单数据或者信号等）直接在testbench中编写产生就行了，但对于复杂的外部激励数据，很难在testbench中产生，这时就要通过读取外部文件里的数据来实现。通过和matlab的配合使用，基本上可以模拟各种外部激励。
举例来说：输入信号是三个不同频率的正弦波的相加，经过FIR低通滤波器滤除高频分量，输出频率最低的那个正弦信号。这种情况下测试用的输入信号不能通过testbench编写产生。
简单来说有以下两种方法可以模拟输入信号：

在FPGA内部通过DDS产生三个正弦波，然后将三个波形相加作为输入信号。
利用matlab产生输入信号，将数据导出为.txt文件，在仿真时读取文件内的数据作为外部激励。

显然第二种方法更加灵活和便捷。下面，具体介绍一下这种方法的使用。

平台：

Vivado 16.4
Matlab R2017b

Matlab程序编写：

代码如下：

%=============设置系统参数==============%

f1=1e6;        %设置波形频率

f2=500e3;

f3=800e3;

Fs=20e6;        %设置采样频率

L=1024;         %数据长度

N=14;           %数据位宽

%=============产生输入信号==============%

t=0:1/Fs:(1/Fs)*(L-1);

y1=sin(2*pi*f1*t);

y2=sin(2*pi*f2*t);

y3=sin(2*pi*f3*t);

y4=y1+y2+y3;

y_n=round(y4*(2^(N-3)-1));      %N比特量化;如果有n个信号相加，则设置（N-n）

%=================画图==================%

a=10;           %改变系数可以调整显示周期

stem(t,y_n);

axis([0 L/Fs/a -2^N 2^N]);      %显示

%=============写入外部文件==============%

fid=fopen('E:\Workspace\Vivado_16.4\TEST\TestBench\sin_data.txt','w');    %把数据写入sin_data.txt文件中，如果没有就创建该文件

for k=1:length(y_n)

    B_s=dec2bin(y_n(k)+((y_n(k))<0)*2^N,N);

    for j=1:N

        if B_s(j)=='1'

            tb=1;

        else

            tb=0;

        end

        fprintf(fid,'%d',tb);

    end

    fprintf(fid,'\r\n');

end

fprintf(fid,';');

fclose(fid);

此程序中设置了三个频率分别为1M、500k和800k的正弦波，然后将三个波形相加并且量化后作为输出。最后将路径设置为相应文件所在路径即可，需要注意的是如果对应路径下没有相应文件，则会自动新建文件并写入数据。

运行程序：
打开对应的文件目录

可以看到二进制数据文件已经生成。

接下来就可以进行在testbench中读取外部数据的操作了。

Testbench的编写：

代码如下：

`timescale 1ns/1ps

module TB_readfile();

reg     SCLK;

reg    [:]  data_out;

//--------------时钟部分----------------//

initial     SCLK = ;

always      #     SCLK = ~SCLK;

//-------------------------------------//

parameter data_num = 'd1024;

integer   i = ;

reg [:]  data_men[:data_num];

reg [:]  data_reg = ;

initial begin

    $readmemb("E:/Workspace/Vivado_16.4/2017_8_28_TEST/TestBench/sin_data.txt",data_men);   //注意斜杠的方向，不能反<<<<<<<

end

always @(posedge SCLK) begin

    data_out <= data_men[i];

    i <= i + 'd1;

end

//------------------------------------//

endmodule

里只需要读取外部数据，所以Vivado工程里只需要添加仿真文件就行了。

运行仿真：

可以看到，仿真的波形和matlab中显示的波形一致，说明结果正确。