Discrete Cosine Transform Filter in Simulink

DSP Nubi Gast 07.11.2011 22:36 #2411451

Lesenswert?

•

Hi :o)

Ich möchte einen Filter in Simulink implementieren der folgendermaßen 
aussieht:

2/N * Summe von i=0 bis N-1 von {(cos [2*pi/N*i]) * z^-i}

N ist die Anzahl meiner Samples
z^-i ist mein Delay

Es ist ein FIR Passbandfilter mit N Taps, eine Form des DCT filters.

Später möchte ich das Ding folgendermaßen erweitern:

2/N *
Summe von i=0 bis N-1 von {
Summe aller h von {
(cos [2*pi/N*h*i]) * z^-i
}
}

h sind in dem Fall die Vielfachen meiner Grundfrequenz.

Ich bin eher der Regelungsfreak und weniger der DSP-Profi. Ich hoffe ihr 
könnt mir bei der Aufgabe auf die Sprünge helfen :)

Antwort Zitieren Bearbeiten Thread verschieben Thread sperren Anmeldepflicht aktivieren Anpinnen Thread löschen Thread mit anderem zusammenführen Beitrag melden

DSP Nubi Gast 07.11.2011 22:41 #2411460

Lesenswert?

•

▲
▼

Als Nachtrag:

Ich habe mal einen simplen MAF mit folgender Syntax aufgebaut:

[1/N .*ones(1,N)]

Also 1/N und N Koeffizienten von z mit 1 belegt.

Vielleicht kann man ähnlich den obigen Filter aufbauen? Was mir nur 
etwas Kopfschmerzen bereitet:

Ich fülle die Koeffizienten nicht nur mit 1 und wie bringe ich das 
mitlaufende z^-1 unter?

Antwort Zitieren Bearbeiten Löschen Beitrag melden

DSP Nubi Gast 08.11.2011 20:32 #2412715

Lesenswert?

•

▲
▼

Ich bins wieder. Ich will das Problem mal herunterbrechen und garnicht 
explizit auf den DCT Filter beziehen.

Letztendlich brauche ich nur ein Array mit folgender Reihe:

[cos(2*pi/N*1 , cos(2*pi/N*2 , cos(2*pi/N*3 , cos(2*pi/N*4 ... 
cos(2*pi/N*i]

Da i bei mir 200 ist brauche ich eine Möglichkeit ein solches Feld zu 
erstellen ohne mich Wund zu tippen.

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Schnuffikus L. (squizzel) 09.11.2011 22:30 #2414344

Lesenswert?

•

▲
▼

N = abtastrate / grundfrequenz

for i = 0:(N-1)
 num = cos(2*pi/N*i)
end

In den Numerator eines FIR trägt man dann 2/N*num ein.

Erweitert auf mehrere Oberschwingungen:

for i = 0:(N-1)
 num = cos(2*pi/N*1*i) + cos(2*pi/N*2*i) + ... + cos(2*pi/N*k*i)
end

Das sollte klappen

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Max G. (l0wside) 10.11.2011 17:01 #2415352

Lesenswert?

•

▲
▼