FIR Filter Verstärkung größer 1?

[TM (Gast)]

Hallo,

Bin grad dabei einen FIR Filter in Octave zu modellieren.

fir.m

function [fw,h]=fir(M,h,fs,step)
%Autor: Tobias Markus
%FIR Filter with M koeffizients
%fs is the sampling frequency and 
%h is the koefficent vector must have the size
%of M. And step defines the steps between 0Hz and fs
n=0;
m=0;
sum =0;
f=[0:step:fs/2];
hz=[];
% Calculating the transfer function
while m < (fs/(2*step)+1)
while n < M
  sum = sum +h(n+1)*exp(-j*n*(2*pi*f(m+1))/fs);
n=n+1;
end
hz(m+1)= sum;
sum=0;
m=m+1;
n=0;
end
h=hz;
fw=f;
endfunction

Hauptprogramm:

h=[];
h(1)=-0.045;
h(2)=0;
h(3)=0.075;
h(4)=0.159;
h(5)=0.225;
h(6)=0.25;
h(7)=0.225;
h(8)=0.159;
h(9)=0.075
h(10)=0
h(11)=-0.045
fs=16000;
M=11;
step=10;
f=[];
w=[];
[f,w]=fir(M,h,fs,step);
plot(f,abs(w));

Wenn ich nun den Filter plotten lasse habe ich Verstärkungen größer 
eins.
Wie kann das sein hab ich einen Fehler gemacht?

[N. M. (mani)]

Ist zwar schon eine Weile her, aber ist die Gleichspannungsverstärkung 
(also bei f(0) ) nicht die Summe der Filterkoeffizienten!?
Wenn ich mich gerade nicht verrechnet habe, dann müsste es A(0)=1,078 
sein.
Das könnte doch einigermaßen mit deinem Bild zusammenpassen.

Grüße

[TM (Gast)]

Ahh ok das kann sein.
Werd mich wohl nen bisschen genauer reinlesen müssen.

Gruß Tobias

[z^-1 (Gast)]

TM schrieb:
> Wie kann das sein hab ich einen Fehler gemacht?

Du hast deine Filterkoeffizienten nicht richtig normiert.

[TM (Gast)]

Um die gehts nicht, geht nur um den Algorithmus.