A Quadrature Signals Tutorial: Complex, But Not Complicated

Understanding the 'Phasing Method' of Single Sideband Demodulation

Complex Digital Signal Processing in Telecommunications

Introduction to Sound Processing

Introduction of C Programming for DSP Applications

There are **9** messages in this thread.

You are currently looking at messages 1 to .

**Is this discussion worth a thumbs up?**

I need to implement a decimation filter( using CIC structure )in C. I have created integrator block,another block for decimation and finally a third block for comb structure and cascaded them in that order.But the final frequency response doesnot match the theoretical graph.I have used 2's compliment notation,Q15 format. It shall be of great help if anyone has sample codes for this filter. Regards, Sowmini

"sowmini" <s...@gmail.com> wrote in message news:Q...@giganews.com... > I need to implement a decimation filter( using CIC structure )in C. > I have created integrator block,another block for decimation and finally a > third block for comb structure and cascaded them in that order.But the > final frequency response doesnot match the theoretical graph.I have used > 2's compliment notation,Q15 format. > > It shall be of great help if anyone has sample codes for this filter. > > Regards, > Sowmini > > How wide are your integrator and comb registers? They have to be very large if you are doing a significant amount decimation, something like input_width + N*log2(R) where R is decimation rate and N is the number of stages. In my vhdl CIC they are about 80 bits wide. So unless you are using floating point, I doubt it will work in straight C code. -Clark

```
... and you can't use floating point for CIC filters.
Dirk
On Oct 5, 8:36 am, "cpope" <cep...@nc.rr.com> wrote:
So unless you are using floating
> point, I doubt it will work in straight C code.
>
> -Clark
```

```
On Oct 5, 8:56 am, dbell <bellda2...@cox.net> wrote:
> ... and you can't use floating point for CIC filters.
>
> Dirk
>
> On Oct 5, 8:36 am, "cpope" <cep...@nc.rr.com> wrote:
> So unless you are using floating
>
>
>
> > point, I doubt it will work in straight C code.
>
> > -Clark- Hide quoted text -
>
> - Show quoted text -
Forgot to say...
I have done them successfully in straight C code with (long) integer
calulations. I don't remember how many stages though.
Dirk
Clark,
How is the world after W-J?
Dirk
```

As per the relation,RegisterWidth = k Ã log2(Mc)+ Bin for k=4,Mc=4 and Bin , the registerwidth$ and i have used 32 bit wide registers instead.

```
sowmini wrote:
> I need to implement a decimation filter( using CIC structure )in C.
> I have created integrator block,another block for decimation and finally a
> third block for comb structure and cascaded them in that order.But the
> final frequency response doesnot match the theoretical graph.I have used
> 2's compliment notation,Q15 format.
>
> It shall be of great help if anyone has sample codes for this filter.
>
> Regards,
> Sowmini
>
The CIC only works with fixed point arithmetic, and you cannot
round/truncate bits in the integrator stages. The integrators must have
enough bits to represent the full scale input multiplied by the CIC gain
without overflow. The CIC gain is (N*R)^M where N is the number of
samples delay in the comb section (typically 1), R is the decimation
ratio, and M is the order of the filter (i.e. the number of cascaded
integrator and comb stages). Typical CIC filters have order (M) of
around 4 or 5, so the gain is related to R^5. The number of bits
required in the accumulator is then I+M*log2(N*R), which can easily
exceed 50 bits. There are not many CIC filters with only 16 bits.
You can drop LSBs in the decimation stage before the comb sections, so
the comb sections don't usually need nearly as many bits as the
integrators in a decimating CIC.
```

>sowmini wrote: > >> I need to implement a decimation filter( using CIC structure )in C. >> I have created integrator block,another block for decimation and finally a >> third block for comb structure and cascaded them in that order.But the >> final frequency response doesnot match the theoretical graph.I have used >> 2's compliment notation,Q15 format. >> >> It shall be of great help if anyone has sample codes for this filter. >> >> Regards, >> Sowmini >> > >The CIC only works with fixed point arithmetic, and you cannot >round/truncate bits in the integrator stages. The integrators must have >enough bits to represent the full scale input multiplied by the CIC gain >without overflow. The CIC gain is (N*R)^M where N is the number of >samples delay in the comb section (typically 1), R is the decimation >ratio, and M is the order of the filter (i.e. the number of cascaded >integrator and comb stages). Typical CIC filters have order (M) of >around 4 or 5, so the gain is related to R^5. The number of bits >required in the accumulator is then I+M*log2(N*R), which can easily >exceed 50 bits. There are not many CIC filters with only 16 bits. >You can drop LSBs in the decimation stage before the comb sections, so >the comb sections don't usually need nearly as many bits as the >integrators in a decimating CIC. > Hi pals, Im now doing a CIC decimation filter design for decimation rate2, delay 2, and 6stages. For the bits overflow, I have try to use the Hogenauer Pruning method from the text book UMeyerBaese, by using the cic.exe i generate the code below. It works fine in VHDL. But for the matlab, when in short bit width it will have the same correct output with the VHDL, somehow it have "noise" when the width is large where it is more than 1000 and above. If can Ill upload a screenshot so that u have a better picture of my problem. Is there any other method to do the hogenauer pruning in Matlab. %Hogenauer pruning cic.exe -- -------------------------------------------------------- -- Program for the design of a CIC decimator. -- -------------------------------------------------------- -- Input bit width Bin = 32 -- Output bit width Bout = 32 -- Number of stages S = 6 -- Decimation factor R = 2 -- COMB delay D = 2 -- Frequency resolution DR = 4 -- Passband freq. ratio P = 8 -- -------------------------------------------------------- -- ----------------- Results of the Design ---------------- -- -------------------------------------------------------- -- -------- Computed bit width: -- -------- Maximum bit growth over all stages = 12 -- -------- Maximum bit width including sign Bmax+1 = 44 -- Stage 1 INTEGRATOR. Bit width : 44 -- Stage 2 INTEGRATOR. Bit width : 42 -- Stage 3 INTEGRATOR. Bit width : 41 -- Stage 4 INTEGRATOR. Bit width : 40 -- Stage 5 INTEGRATOR. Bit width : 39 -- Stage 6 INTEGRATOR. Bit width : 38 -- Stage 1 COMB. Bit width : 38 -- Stage 2 COMB. Bit width : 37 -- Stage 3 COMB. Bit width : 37 -- Stage 4 COMB. Bit width : 36 -- Stage 5 COMB. Bit width : 35 -- Stage 6 COMB. Bit width : 34 -- ------- Maximum aliasing component : 0.000034 = 89.30 dB -- ------- Amplitude distortion : 0.532563 = 5.47 dB --MATLAB-- D = 2; % Differential delay R = 2; % Downsampler delayBuffer1 = zeros(1,D); delayBuffer2 = zeros(1,D); delayBuffer3 = zeros(1,D); delayBuffer4 = zeros(1,D); delayBuffer5 = zeros(1,D); delayBuffer6 = zeros(1,D); intOut1 = 0; intOut2 = 0; intOut3 = 0; intOut4 = 0; intOut5 = 0; intOut6 = 0; input = [Pulse100nsSigned' zeros(1,20)]; % impulse input output = []; % output for i = 1:length(input) % integrators intOut1 = intOut1 + input(i)/4; intOut2 = intOut2 + intOut1/2; intOut3 = intOut3 + intOut2/2; intOut4 = intOut4 + intOut3/2; intOut5 = intOut5 + intOut4/2; intOut6 = intOut6 + intOut5; % downsample if mod(i,R)==1 % combs combOut1 = intOut6/2 - delayBuffer1(2); delayBuffer1(2) = delayBuffer1(1); delayBuffer1(1) = intOut6/2; combOut2 = combOut1 - delayBuffer2(2); delayBuffer2(2) = delayBuffer2(1); delayBuffer2(1) = combOut1; combOut3 = combOut2/2 - delayBuffer3(2); delayBuffer3(2) = delayBuffer3(1); delayBuffer3(1) = combOut2/2; combOut4 = combOut3/2 - delayBuffer4(2); delayBuffer4(2) = delayBuffer4(1); delayBuffer4(1) = combOut3/2; combOut5 = combOut4/2 - delayBuffer5(2); delayBuffer5(2) = delayBuffer5(1); delayBuffer5(1) = combOut4/2; combOut6 = combOut5/4 - delayBuffer6(2); delayBuffer6(2) = delayBuffer6(1); delayBuffer6(1) = combOut5/4; output = [output combOut6]; end end output

clairechan <c...@n_o_s_p_a_m.gmail.com> wrote: >Hi pals, Im now doing a CIC decimation filter design for decimation rate2, >delay 2, and 6stages. For the bits overflow, I have try to use the >Hogenauer Pruning method from the text book UMeyerBaese, by using the >cic.exe i generate the code below. It works fine in VHDL. But for the >matlab, when in short bit width it will have the same correct output with >the VHDL, somehow it have "noise" when the width is large where it is more >than 1000 and above. If can Ill upload a screenshot so that u have a better >picture of my problem. Is there any other method to do the hogenauer >pruning in Matlab. Since it works in VHDL, can you patch the VHDL code into your matlab simulation and run it, so that you can avoid evne modeling it in matlab. Steve (always looking for a quick fix)

>clairechan <c...@n_o_s_p_a_m.gmail.com> wrote: > >>Hi pals, Im now doing a CIC decimation filter design for decimation rate2, >>delay 2, and 6stages. For the bits overflow, I have try to use the >>Hogenauer Pruning method from the text book UMeyerBaese, by using the >>cic.exe i generate the code below. It works fine in VHDL. But for the >>matlab, when in short bit width it will have the same correct output with >>the VHDL, somehow it have "noise" when the width is large where it is more >>than 1000 and above. If can Ill upload a screenshot so that u have a better >>picture of my problem. Is there any other method to do the hogenauer >>pruning in Matlab. > >Since it works in VHDL, can you patch the VHDL code into your matlab >simulation and run it, so that you can avoid evne modeling it in matlab. > > >Steve (always looking for a quick fix) > > > Thanks the tips is very helpful, and the problem is solved. God bless.