summaryrefslogtreecommitdiff
path: root/macros/latex/contrib/withesis/derivs.m
blob: c723cdbbe4fda29ba9e0958d84558b3f13377b93 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
function dPsi = derivs(t,Psi); 
% derivs computes the state derivatives for an induction machine 
% The machine is in the synchronous reference frame. 
% Eric Benedict, Spring 1996 
 
% Machine parameters... 
global Rs Rqr Rdr Xls Xlr Xm Xms J P We Wb Vs 
 
% rename state vector... 
Phqs=Psi(1); 
Phds=Psi(2); 
Phqr=Psi(3); 
Phdr=Psi(4); 
Wr=Psi(5); 
Theta=Psi(6); 
 
 
% Select Reference Frame... 
% Assume no angle differnce between a-axis and q-axis at t=0 

% Synchronous
W=We; 
ThetaR=W*t; 
  
% Compute Machine Currents... 
Phmq=Xms*Phqs/Xls + Xms*Phqr/Xlr; 
Phmd=Xms*Phds/Xls + Xms*Phdr/Xlr; 
 
Iqs=(Phqs - Phmq)/Xls; 
Ids=(Phds - Phmd)/Xls; 
Iqr=(Phqr - Phmq)/Xlr; 
                                    Idr=(Phdr - Phmd)/Xlr; 
 
 
% Compute Torques... 
Te=3*P*(Phds*Iqs-Phqs*Ids)/(4*Wb);    %  Electrical Torque 
Tl=0;                                 %  Load Torque 
  
 
% Compute Voltages... 
Vqs=Vs*cos(We*t - ThetaR); 
Vds=-1*Vs*sin(We*t - ThetaR); 
Vqr=0; 
Vdr=0; 
  
% Compute new derivative of state vector... 
dPsi(1) = Vqs - Rs*(Phqs - Phmq)/Xls - W*Phds/Wb;          % PHqs/Wb 
dPsi(2) = Vds - Rs*(Phds - Phmd)/Xls + W*Phqs/Wb;          % PHds/Wb 
dPsi(3) = Vqr - Rqr*(Phqr - Phmq)/Xlr - (W-Wr)*Phdr/Wb;    % PHqr/Wb 
dPsi(4) = Vdr - Rdr*(Phdr - Phmd)/Xlr + (W-Wr)*Phqr/Wb;    % PHdr/Wb 
 
dPsi(5) = (Te - Tl)*P/(2*J*Wb);           % Wr/Wb 
 
dPsi(6) = Wr/Wb;                          % Theta/Wb 

dPsi=Wb*dPsi;