indata_plot.m

来自「FPGA开发光盘各章节实例的设计工程与源码」· M 代码 · 共 33 行

%---------------将仿真输入数据从文件中读出------------------%
%           f1=150;                                                               %
%           f2=330;                                                               %
%           fs=1000;                                                             %
%           fir_coef=0.5;                                                       %
%           i=3.1415926;                                                      %    
%           dets=1/fs;                                                           %
%           fp=fir_coef*fs/2;                                                  %    
%---------------------------------------------------------------------%
fid=fopen('filter.vec','r');
[data_in]=fscanf(fid,'%d',[1,inf]);
fclose(fid);

fs=1000;
dets=1/fs;

%----------绘制输入数据的时域和频域图------------%
t=0:160;
t=t*dets;
subplot(2,2,1);
plot(t,data_in/1024,'-k');
axis([0,0.1,-2,2])
title  ('输入数据的时域图','fontsize',8);
xlabel ('时间t（秒）','fontsize',8);
ylabel ('幅度y','fontsize',8);

subplot(2,2,2);
stem(abs(fft(data_in/1024)),'*b');
axis([0,80,0,80])
title  ('输入数据的频谱图','fontsize',8);
xlabel ('频率（赫兹）','fontsize',8);
ylabel ('幅度Y','fontsize',8);