reg [7:0]data; initial begin #10 clk =0; forever #4 clk = ~clk; end initial begin #20 rst=1; #20 rst=0; end reg [7:0]data_sin[160000:0]; ////改動點數據矩陣長度設置 integer i; initial begin i=0; begin $readmemb("D:/Chapter_8/E8_1_QAMModem/QAM.txt",data_sin,0,160000); ///改動點數據矩陣長度 注意是“/” 而不是“\” end forever begin @(posedge clk) begin i <= i+1; din <= data_sin[i]; end end end endmodule
matlab 寫txt文本的代碼
fid = fopen('data.txt','w');
for oo=1:1:i
if mod(oo,10) == 0
fprintf(fid,'%f,%f,\n',sI1(oo),sQ1(oo));
else
fprintf(fid,'%f,%f,',sI1(oo),sQ1(oo));
end
end
fclose(fid);
verilog 對應的寫文件,寫入IQ數據
integer file_out; initial begin file_out = $fopen("mI.txt"); if (!file_out) begin $finish; end end wire signed [23:0] dout_s = fifo_fft_data[47:24]; always @ (posedge clk) begin if(fifo_fft_valid) $fdisplay(file_out, "%d", dout_s); end integer file_out_Q; initial begin file_out_Q = $fopen("mQ.txt"); if (!file_out_Q) begin $finish; end end wire signed [23:0] dout_s_Q = fifo_fft_data[23:0]; always @ (posedge clk) begin if(fifo_fft_valid) $fdisplay(file_out_Q, "%d", dout_s_Q); end
對應以上文件的matlab 讀取數據:
%讀取FPGA仿真出的數據
clc;
clear;close all;
fid=fopen('mI.txt','r');
[di,N]=fscanf(fid,'%lg',inf);
fclose(fid);
fid=fopen('mQ.txt','r');
[dq,N]=fscanf(fid,'%lg',inf);
fclose(fid);
exp1=di+dq*1i;
% exp2 = exp1(25000:30000);
% exp2 = exp1(1024:8192);
exp2 = exp1;
plot(di);
figure;
plot(dq);
figure;
plot(20*log10(abs(fft((exp2).* window(@gausswin,length(exp2),4)))));
