Matlab function Eulero

[jayeffe]

Salve a tutti. In un esercizio riguardo le discretizzazioni con metodo di eulero mi viene fornito questo codice .
Volevo sapere cosa fa a a cosa si riferisce.
Ho letto "metodo computazionale"..

Codice: Seleziona tutto

%Eulero
% global u1_e u2_e  e1_e e2_e num_e den_e;
u1_e=0;
u2_e=0;
e1_e=0;
e2_e=0;
[num_e,den_e] = tfdata(Cdeulero,'v');

%% matlab function 2 polo 2 zero eulero in avanti
function ud2 = Funzione_eulero( ee )

global u1_e u2_e e1_e e2_e num_e den_e;

app2 = -(den_e(2)/den_e(1))*u1_e -(den_e(3)/den_e(1))*u2_e +ee*(num_e(1)/den_e(1)) +e1_e*(num_e(2)/den_e(1)) +e2_e*(num_e(3)/den_e(1));

u2_e=u1_e;
u1_e=app2;
e2_e=e1_e;
e1_e=ee;

ud2=app2;


end

[SandroCalligaro]

Dando un'occhiata, si intuisce che il codice implementa (probabilmente) un filtro IIR del 2o ordine, la cui funzione di trasferimento, però, è definita altrove.

In particolare, ee sembra essere l'ingresso, app o ud_2 l'uscita, num_e e den_e i vettori che contengono i coefficienti di numeratore e denominatore, e1_e,e2_e dovrebbero essere gli ingressi agli istanti precedenti, mentre u2_e, u1_e dovrebbero essere le uscite ritardate.

[jayeffe]

Si , è scritto metodo computazionale, potrebbe essere un equazione alle differenze

[SandroCalligaro]

Scusami, ma non mi è chiaro quanto sai di sistemi a tempo discreto...
Dovresti almeno dire se conosci la trasformata zeta...

[jayeffe]

Si Conosco la trasformata zeta. Il problema che me la cavo poco con matlab.