Digitaler Tiefpass VHDL

[Tobias (Gast)]

Für ein Projekt benötige ich einen einfachen digitalen Tiefpassfilter 
z.B. 1. Ordnung. Hierfür bin ich im Netz auf folgenden Code gestoßen den 
ich nicht ganz verstehe! Ich weis das es sich um ein TP 1. Ordnung 
handelt, sehe allerdings nicht wo ich die Werte für die Grenzfrequenz 
einstellen kann/ bzw. wie der Zusammenhang mit der Übertragungsfunktion 
des äuivalenten analogen Filters aussieht!?
Wer kann mir weiterhelfen?

  LPF: process (Clk, Rst)
    variable sumI, sumQ: signed(baseband_bits-1 downto 0);
    
    function "sra"(x: signed; n: natural) return signed is
      variable xx: signed(x'length-1 downto 0) := x;
    begin
      for i in 1 to n loop
        xx := xx(xx'left) & xx(xx'left downto 1);
      end loop;
      return xx;
    end;

  begin

    if Rst = '1' then
      filterI <= (others => '0');
      filterq <= (others => '0');
      sumI    := (others => '0');
      sumQ    := (others => '0');
    idumpI  <= (others => '0');
    idumpQ  <= (others => '0');
    elsif rising_edge(Clk) then
      if dump = '1' then -- Dump integrator, generate output
        -- Output lowpass
        filterI <= resize(
                     filterI  + ((signed'((sumI - filterI) * filter_rate)) sra filter_shift),
                     filterI'length);
        filterQ <= resize(
                     filterQ  + ((signed'((sumQ - filterQ) * filter_rate)) sra filter_shift),
                     filterQ'length);
        -- dump the integrator
        sumI := (others => '0');
        sumQ := (others => '0');
      end if;
      -- Integrate (or initialise, if it's the dump cycle)
      sumI := sumI + prodI;
      sumQ := sumQ + prodQ;
      idumpI <= sumI;
      idumpQ <= sumQ;
    end if;

  end process;


  -- Produce the output
  --
  I <= filterI;
  Q <= filterQ;

end;

[Armin (Gast)]

Du suchst nach der Tustin-Transformation
http://de.wikipedia.org/wiki/Bilineare_Transformation_(Signalverarbeitung)


damit kannst du eine bekannte Übertragungsfunktion G(s), z.B. für 
besagten Tiefpass, in eine zeitdiskrete z-Übertragungsfunktion G(z) 
transformieren (Näherung).
Das Ergebnis lässt sich dann direkt mit "Zeitverzögerungen"=Register in 
PLD implementieren.

Was ich aber so bisher gelesen habe, macht man das heutzutage nicht mehr 
per Hand. Dafür gibt's in vielen Toolchains spezielle Filterrechner.

[Tobias (Gast)]

Ahh ok..  danke für die schnelle Antwort
Kannst du mir da einen speziellen Filterrechner empfehlen?
(Im Falle das du schon mal mit einem gearbeitet hast)

[Christoph db1uq K. (christoph_kessler)]

http://www.mikrocontroller.net/articles/Digitalfilter_mit_ATmega
Filterberechnung mit Scilab und Octave