Subversion Repositories wdsp

%Header file generation for Wishbone compliant FIR Core test

%Cleans work environment

clc;

clear all;

close all;

%Low-pass filter design

%Cut-off frequencies

fs=8000;

fc1=1500;

fc2=2000;

%Length of impulse response

Nh=50;

order=Nh-1;

%Designs filter

hn = firpm(order,[0 fc1 fc2 fs/2]/(fs/2),[1 1 0 0]);

%Bits to use for each coefficient

Nbits=16;

%Fractional part

Q=floor(Nbits-1-log2(max(abs(hn))))-1;

%Quantized impulse response

hn_Q=(round(hn*2^Q));

%Header file generation

handle = fopen('../../FIR50/coefs.h', 'wt');

fprintf(handle, '/*Impulse response for FIR filter automatically generated header file*/\n');

fprintf(handle, '/*M.Eng. Alexander López Parrado*/\n\n');

fprintf(handle, '#include <stdint.h>\n');

fprintf(handle, '/*Number of coeffcients*/\n');

fprintf(handle, '#define Nh %d\n\n',Nh);

fprintf(handle, '/*Number of bits in fractional part of coeffcients*/\n');

fprintf(handle, '/*Fixed point format with %d bits ([%d].[%d])*/\n',Nbits,Nbits-Q,Q);

fprintf(handle, '#define Q %d\n\n',Q);

fprintf(handle, '/*Filter Coefficients */');

fprintf(handle, '\nconst int16_t hn[Nh] = {\n');

hn_printf= [hn_Q; hn];

fprintf(handle, '%d, //%f\n',hn_printf);

fprintf(handle, '};\n');

fclose(handle);