Octave CIC Filter

[tm (Gast)]

Hallo,

Ich möchte mich etwas näher in DSPs einarbeiten und bin gerade dabei 
einen CIC Filter in Octave zu implementieren.
Ich arbeite mit der Übertragungsfunktion um daraus dann den Frequenzgang 
zu erschließen. Allerdings glaube das ich genau da ein Denkfehler 
gemacht habe.

Hier ist mein Octave Code:

R=32; % Dezimierungsfaktor
N=3; % Anzahl Stufen
M=1;
f=[0:10:2000000]; % Frequenz bis 2 MHz da Systemtakt 4MHz
w=f/4000000; % Frequenz bezogen auf den Systemtakt
z=exp(j*w); % Ist hier vllt mein Fehler 
hz=((1-z.^(-R*M))./(1-z.^(-1))).^N;
plot(w,abs(hz));

[tm (Gast)]

tm schrieb:
> Hallo,
>
> Ich möchte mich etwas näher in DSPs einarbeiten und bin gerade dabei
> einen CIC Filter in Octave zu implementieren.
> Ich arbeite mit der Übertragungsfunktion um daraus dann den Frequenzgang
> zu erschließen. Allerdings glaube das ich genau da ein Denkfehler
> gemacht habe.
>
> Hier ist mein Octave Code:
> R=32; % Dezimierungsfaktor
>
> N=3; % Anzahl Stufen
>
> M=1;
>
> f=[0:10:2000000]; % Frequenz bis 2 MHz da Systemtakt 4MHz
>
> w=f/4000000; % Frequenz bezogen auf den Systemtakt
>
> z=exp(j*w); % Ist hier vllt mein Fehler
>
> hz=((1-z.^(-R*M))./(1-z.^(-1))).^N;
>
> plot(w,abs(hz));

Hey,

Ich glaub fast ich hab alles richtig gemacht musste nur noch hz in dB 
umwandeln und jetzt siehts eigentlich recht ordentlich aus.

Falls doch noch jemand einen Fehler entdeckt bitte einfach raus damit ;)

Gruß Tobias

[tm (Gast)]

tm schrieb:
>> w=f/4000000; % Frequenz bezogen auf den Systemtakt


Ok hier ist mein ich ein Fehler gewesen die normierte frequenz ist 
soweit ich nun herrausgefunden habe mit

$$w=2*pi*f/4000000$$

 definiert wobei die 4MHz hier nun meine Samplefrequenz ist.

[Stefan (Gast)]

tm schrieb:
> 2 MHz da Systemtakt 4MHz

Ist das wirklich der Systemtakt oder das  was für die Berechung des CIC 
zur Verfügung steht?

[tm (Gast)]

Da hab ich mich wohl im Kommentar falsch ausgedrückt es ist natürlich 
die Samplefrequenz mit welcher Daten in den Filter kommen.

Gruß Tobias