Simulazione matlab motore asincrono

[oressac]

Salve
Ho provato a realizzare un modello matlab (non tramite simulink) del motore asincrono, ma non riesco farlo girare correttamente:
Ecco il codice:

Codice: Seleziona tutto

f=50;
ampiezza=230*sqrt(2);
w=2*pi*f;
durata_secondi=0.4;


%caratteristiche motore

Rs=1.95;
Ls=0.255;
Rr=2.65;
Lr=0.254;
Lm=0.214;
p=2;

Jm=0.0054;
Fm=0.003;

%caratteristiche carico

Cm=0;
Jl=0.01;
Fl=0.001;


J=Jl+Jm;
F=Fl+Fm;

    function di = equaz(t,i)
               
                v1=ampiezza*sin(w*t);
                v2=ampiezza*sin(w*t-2/3*pi);
                v3=ampiezza*sin(w*t+2/3*pi);
               
                v=park([v1,v2,v3]);
                vas=v(1);
                vbs=v(2);
             
               di = zeros(5,1);   
               
               di(1) = (-Rs*i(3)+vas-Ls*di(3))/Lm;  %diar
               di(2) = (-Rs*i(4)+vbs-Ls*di(4))/Lm;  %dibr
               di(3) = -(Rr*i(1)+Lr*i(5)*p*i(2)+Lm*i(5)*p*i(4)+Lr*di(1))/Lm; %dias   
               di(4) = (Lr*i(5)*p*i(1)+Lm*i(5)*p*i(3)-Rr*i(2)-Lr*di(2))/Lm;  %dibs   
               
               Ce    = 1.5*p*Lm*(i(3)*i(2)-i(4)*i(1));  %Ce
               di(5) = (Ce-F*i(5)-Cm)/J;   %dwm
               
    end

            [T1,I] = ode45(@equaz,[0 , durata_secondi],[0 0 0 0 0]);

Le correnti mi divergono e non si capisce più nulla. E strano perche faccio risolvere tutto al matlab, io gli passo solo il sistema di eq. differenziali che caratterizzano il motore, le ho controllate e ricontrollate, sembreno scritte bene. Qualcuno puo darmi qualche dritta?

[rini]

Ho visto che hai fatto la trasformata di Park.

Codice: Seleziona tutto

v=park([v1,v2,v3]);

Controlla come è definita questa trasformata oppure postala qui che proviamo a controllarla noi.

[oressac]

L'ho definita cosi:

To=sqrt(2/3)*[1 -1/2 -1/2 ; 0 sqrt(3)/2 -sqrt(3)/2; 1/sqrt(2) 1/sqrt(2) 1/sqrt(2)];
v=[v1,v2,v3]*To';
vas=v(1);
vbs=v(2);

Non credo ci siano problemi,proprio non riesco a capire perche diverge!

[oressac]

Non so proprio dove sbatterci!!!!!