Молодой учёный

138 Информатика
«Молодой учёный» . № 3 (50) . Март, 2013 г.
-Lr*wn*IS(8)/dt+(-Lr*R(8)/dt)*X(7)+((Lr*(R(9)+R(8))+1)/dt)*X(8)+
(-Lr*R(9)/dt)*X(9);
-Lr*wn*IS(9)/dt+(-Lr*R(9)/dt)*X(8)+((Lr*(R(10)+R(9))+1)/dt)*X(9)+
(-Lr*R(10)/dt)*X(10);
-Lr*wn*IS(10)/dt+(-Lr*R(10)/dt)*X(9)+((Lr*(R(11)+R(10))+1)/dt)*X(10)+
(-Lr*R(11)/dt)*X(11);
-Lr*wn*IS(11)/dt+(-Lr*R(11)/dt)*X(10)+((Lr*(R(12)+R(11))+1)/dt)*X(11)+
(-Lr*R(12)/dt)*X(12);
-Lr*wn*IS(12)/dt+(-Lr*R(12)/dt)*X(11)+((Lr*(R(13)+R(12))+1)/dt)*X(12)+
(-Lr*R(13)/dt)*X(13);
-Lr*wn*IS(13)/dt+(-Lr*R(13)/dt)*X(12)+((Lr*(R(14)+R(13))+1)/dt)*X(13)+
(-Lr*R(14)/dt)*X(14);
-Lr*wn*IS(14)/dt+(-Lr*R(14)/dt)*X(13)+((Lr*(R(15)+R(14))+1)/dt)*X(14)+
(-Lr*R(15)/dt)*X(15);
-Lr*wn*IS(15)/dt+(-Lr*R(15)/dt)*X(14)+((Lr*(R(16)+R(15))+1)/dt)*X(15)+
(-Lr*R(16)/dt)*X(16);
-Lr*wn*IS(16)/dt+(-Lr*R(16)/dt)*X(15)+((Lr*(R(17)+R(16))+1)/dt)*X(16)+
(-Lr*R(17)/dt)*X(17);
(-Lr*R(17)/dt)*X(16)+((Lr*(R(18)+R(17))+1)/dt)*X(17)+
(-Lr*R(18)/dt)*X(18);
(-Lr*R(18)/dt)*X(17)+((Lr*(R(19)+R(18))+1)/dt)*X(18)+(-Lr*R(19)/dt)*X(19)
(-Lr*R(19)/dt)*X(18)+((Lr*(R(20)+R(19))+1)/dt)*X(19)+(-Lr*R(20)/dt)*X(20);
(-Lr*R(20)/dt)*X(19)+((Lr*(R(21)+R(20))+1)/dt)*X(20);
(X(5)+X(6)-X(9)-X(10)-X(11)-X(12)+X(15)+X(16))*wn/dt+(Ls/dt)*(X(21)-X(23))+Uab;
(X(9)+X(10)-X(15)-X(16)-X(13)-X(14)+X(7)+X(8))*wn/dt+(Ls/dt)*(X(23)-X(22))+Ubc;
0];
% Решение методом Гаусса-Жордана
Z=rref([A S]);
X=Z(1:23,24:24); %Выделение последнего столбца из матрицы
% Матрица токов ротора
Ir=[X(1)*(R(1)+R(2))-X(2)*R(2);
-X(1)*R(2)+X(2)*(R(2)+R(3))-X(3)*R(3);
-X(2)*R(3)+X(3)*(R(3)+R(4))-X(4)*R(4);
-X(3)*R(4)+X(4)*(R(4)+R(5))-X(5)*R(5);
-wn*IS(5)-X(4)*R(5)+X(5)*(R(5)+R(6))-X(6)*R(6);
-wn*IS(6)-X(5)*R(6)+X(6)*(R(6)+R(7))-X(7)*R(7);
-wn*IS(7)-X(6)*R(7)+X(7)*(R(7)+R(8))-X(8)*R(8);
-wn*IS(8)-X(7)*R(8)+X(8)*(R(8)+R(9))-X(9)*R(9);
-wn*IS(9)-X(8)*R(9)+X(9)*(R(9)+R(10))-X(10)*R(10);
-wn*IS(10)-X(9)*R(10)+X(10)*(R(10)+R(11))-X(11)*R(11);
-wn*IS(11)-X(10)*R(11)+X(11)*(R(11)+R(12))-X(12)*R(12);
-wn*IS(12)-X(11)*R(12)+X(12)*(R(12)+R(13))-X(13)*R(13);
-wn*IS(13)-X(12)*R(13)+X(13)*(R(13)+R(14))-X(14)*R(14);
-wn*IS(14)-X(13)*R(14)+X(14)*(R(14)+R(15))-X(15)*R(15);
-wn*IS(15)-X(14)*R(15)+X(15)*(R(15)+R(16))-X(16)*R(16);
-wn*IS(16)-X(15)*R(16)+X(16)*(R(16)+R(17))-X(17)*R(17);
-X(16)*R(17)+X(17)*(R(17)+R(18))-X(18)*R(18);
-X(17)*R(18)+X(18)*(R(18)+R(19))-X(19)*R(19);
-X(18)*R(19)+X(19)*(R(19)+R(20))-X(20)*R(20);
-X(19)*R(20)+X(20)*(R(20)+R(21))];
% Электромагнитное усилие
F=[ (X(2))*Ir(1)/(2*tz);
(X(3)-X(1))*Ir(2)/(2*tz);