xk_cmpr.m

fft c源程序

%function fft_cmpr;
ax_font_size = 22;
ax_lab_size = 24;
fmt = '-tiff -deps';
print_flag=0;

n = 16;

fname_c_model = input('C model input data: ','s');
fname_vhdl = input('testbench output data: ','s');

fp=fopen(fname_c_model,'r');
x=fscanf(fp,'%f');
fclose(fp);

fp=fopen(fname_vhdl,'r');
sim_data=fscanf(fp,'%f');
fclose(fp);

% asemble complex transform vector - C model

u = (x(1:2:2*n) + j*x(2:2:2*n))/32768.0;

% form complex series from VHDL simulation results

xsim_data = sim_data(1:2:length(sim_data)) + j*sim_data(2:2:length(sim_data));
xsim_data = xsim_data/32768.0;

figure(1)

b = 0:length(xsim_data)-1;

plot(b,real(xsim_data),b,real(u),'m')
grid
ax=axis; ax(1)=0; ax(2)=n-1; axis(ax);
title('TRANSFORM OVERLAY - REAL');

figure(2)

plot(b,imag(xsim_data),b,imag(u),'g')
grid
ax=axis; ax(1)=0; ax(2)=n-1; axis(ax);
title('TRANSFORM OVERLAY - IMAG.');

% error calculations for the real and imaginary components

er = real(u)-real(xsim_data);
ei = imag(u)-imag(xsim_data);

n = length(er);
q = 0:length(er)-1;
figure(3)
plot(q,er,'g')
grid
ax=axis; ax(1)=0;; ax(2)=n-1; axis(ax);
title('ERROR SIGNAL - REAL')

figure(4)
plot(q,ei,'b')
grid
ax=axis; ax(1)=0;; ax(2)=n-1;; axis(ax);
title('ERROR SIGNAL - IMAG');