HIlfe bei Übertragungsfunktion von CIC-FIlter

[gerbej l. (lab)]

Hallo zusammen,

ich versuche seit gestern das Verhalten eines CIC-Filters im Zeit und 
Frequenzbreich mit Matlabzu untersuchen.daher wollte ich die 
Impulsantwort und seine Reaktion auf eine beliebige Signalfolge im 
Zeitbereich simulieren. hier ist seine Übertragungsfunktion:

h(z)=[(1-z^(-R*M))/(1-z^-1)]^N

R=10,M=1, N=8 als Beispiel
mein problem ist die Argumente für die funktion tf() zu bestimmen, weis 
Jemand, wie man die Koeffizienten von der Zähler und Nennerploynom 
bestimmen kann.

ich bin für jede Hilfe sehr dankbar

Gruß
Lab

[NopNop (Gast)]

Kenne mich in Matlab nicht sonderlich aus.
Könnte dass so passen?

%CIC-Filter

clear all;
close all;

R=10;
M=1;
N=8;

%sys = tf(num,den,Ts)
z=tf('z'); %Sampletime noch angeben?

h= (1-z^(-R*M))/(1-z^(-1))^N;

num= h.num{1};
denum= h.den{1};

fvtool(num,denum);

Wenn das so passt könnte man noch kurz ne Funktion oder so stricken.
Grüße

[NopNop (Gast)]

Gerade ausprobiert mit nem kleinen FIR-Filter.
Händisch berechnete & gezeichnete Impulsantwort stimmt mit der von 
MATLAB überein.

G(z) = z * z/(z-1)^2 - 2z/(z-1)^2 * z^-4 + z/(z-1)^2 * z^-9
        Gerade m.        Gerade m.              Gerade m.
        Steigung 1       Steigung -2            Steigung 1

Jeweils bei z              z^-4                   z^-9

Also ein Dreieck:
                 +
                 |
                 |
                 |
                 |
                5|      /´ \
                 |    ,´    `\
                 |  ,´        `.
                1| /            `.
                 |
                 +------------------+
                 0        4      8

Grüße

[gerbej l. (lab)]

Hallo NopNop,

danke sehr Ihre Hilfe; ich habe es probiert und es klappt bei  mir 
immer noch nicht ;ich habe so Ihr code probiert

R=10;
M=1;
N=8;
Ts=0:001:1;
sys = tf(num,den,Ts);

h= (1-z^(-R*M))/(1-z^(-1))^N;
num= h.num{1}
denum= h.den{1};
fvtool(num,denum);

bokomme folgende fehler

error: `num' undefined near line 5 column 10
error: evaluating argument list element number 1
>>>error: evaluating argument list element number 1
error: called from:

konnten Sie mir bitte Ihr code bisschen erklären
danke im voraus


PS
sorry ich arbeite nicht mit matlab sondern Octave
aber bis zu fvtool sollte es eigentlich gehen

[NopNop (Gast)]

gerbej lab schrieb:
> sorry ich arbeite nicht mit matlab sondern Octave

Dann sollte man in seiner Frage auch nicht schreiben dass man Matlab 
verwendet ;-)
Bei der Umsetzung in Octave kann ich Dir nicht helfen.
Sollte aber gleich bzw. ähnlich gehen.

gerbej lab schrieb:
> ich habe es probiert und es klappt bei  mir
> immer noch nicht ;ich habe so Ihr code probiert

Das wundert mich nicht. So wie Du das geschrieben hast KANN es garnicht 
funktionieren. Wo definierst Du denn "num" und "den" ?

R=10;
M=1;
N=8;
Ts=0:001:1;
sys = tf(num,den,Ts); % <== Fehler, da num&den nicht definiert

h= (1-z^(-R*M))/(1-z^(-1))^N;
num= h.num{1}         % <== Erst hier wird num definiert
denum= h.den{1};      % <== den gibt es garnicht, sondern nur "denum"
fvtool(num,denum);

gerbej lab schrieb:
> konnten Sie mir bitte Ihr code bisschen erklären

%CIC-Filter
clear all;  %Alle Variablen löschen
close all;  %Alle Fenster schließen

%===================================
% Parameter
%===================================
R=10;
M=1;
N=8;

%z = tf('z',Ts) to specify a TF model with sample time Ts using a rational
%function in the discrete-time variable, z. Set Ts = -1 or Ts = [] to
%leave the sample time unspecified.

z=tf('z', -1);

%Erzeugen der Impulsantwort
h= (1-z^(-R*M))/(1-z^(-1))^N;

%Extrahieren der Koeffizienten aus der Transferfunktion h
nume= h.num{1};
denum= h.den{1};

%Aufruf des "Filter Visualization Tool" um Amplitudengang, Phase,
%Pol-Nullstellendiagramm anzuzeigen
fvtool(nume,denum);

Grüße