%48階FIR濾波器 9BIT量化
b=[-6 -9 -9 -4 6 17 23 21 9 -6 -16 -13 2 22 32 24 -3 -37 -54 -36 24 113 201 255 ...
    255 201 113 24 -36 -54 -37 -3 24 32 22 2 -13 -16 -6 9 21 23 17 6 -4 -9 -9 -6];
[h,w]=freqz(b,47);
plot(w,abs(h));

a=[1 0 0 1 0 0 1 0 0 1 ...
   0 1 0 1 0 0 1 0 1 0 ...
   1 0 1 0 1 0 1 0 1 0 ...
   0 1 1 1 1 0 1 0 0 1 ...
   1 0 0 1 1 0 1 0 0 1 ...
   0 1 0 1 0 1 0 1 0 0 ...
   0 1 1 1 1 0 0 1 1 0 ...
   0 1 0 1 0 1 1 0 1 0 ...
   0 1 1 1 0 0 0 1 0 1 ...
   0 0 0 1 1 1 1 0 0 1 ...
   1 0 0 1 ...
   1 0 0 1 0 0 1 0 0 1 ...
   0 1 0 1 0 0 1 0 1 0 ...
   1 0 1 0 1 0 1 0 1 0 ...
   0 1 1 1 1 0 1 0 0 1 ...
   1 0 0 1 1 0 1 0 0 1 ...
   0 1 0 1 0 1 0 1 0 0 ...
   0 1 1 1 1 0 0 1 1 0 ...
   0 1 0 1 0 1 1 0 1 0 ...
   0 1 1 1 0 0 0 1 0 1 ...
   0 0 0 1 1 1 1 0 0 1 ...
   1 0 0 1 0 0 1 0 0 1 ...
   0 1 0 1 0 0 1 0 1 0 ...
   1 0 1 0 1 0 1 0 1 0 ...
   0 1 1 1 1 0 1 0 0 1 ...
   1 0 0 1 1 0 1 0 0 1 ...
   0 1 0 1 0 1 0 1 0 0 ...
   0 1 1 1 1 0 0 1 1 0 ...
   0 1 0 1 0 1 1 0 1 0 ...
   0 1 1 1 0 0 0 1 0 1 ...
   0 0 0 1 1 1 1 0 0 1 ...
   ];
y=conv(a,b);

figure
plot(y);