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

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

magz.elibraries.eu
from magz.elibraries.eu More from this publisher
20.07.2013 Views

136 Информатика «Молодой учёный» . № 3 (50) . Март, 2013 г. A(12,10)=R(11)*Lr*v0/(2*tz); A(12,11)=-(rr+Lr/dt)*R(12)-(Lr*(R(12)+R(11))+1)*v0/(2*tz); A(12,12)=(rr+Lr/dt)*(R(13)+R(12))+(R(12)-R(13))*Lr*v0/(2*tz)+1/dt; A(12,13)=-(rr+Lr/dt)*R(13)+(Lr*(R(14)+R(13))+1)*v0/(2*tz); A(12,14)=-R(14)*Lr*v0/(2*tz); % n=13 A(13,11)=R(12)*Lr*v0/(2*tz); A(13,12)=-(rr+Lr/dt)*R(13)-(Lr*(R(13)+R(12))+1)*v0/(2*tz); A(13,13)=(rr+Lr/dt)*(R(14)+R(13))+(R(13)-R(14))*Lr*v0/(2*tz)+1/dt; A(13,14)=-(rr+Lr/dt)*R(14)+(Lr*(R(15)+R(14))+1)*v0/(2*tz); A(13,15)=-R(15)*Lr*v0/(2*tz); % n=14 A(14,12)=R(13)*Lr*v0/(2*tz); A(14,13)=-(rr+Lr/dt)*R(14)-(Lr*(R(14)+R(13))+1)*v0/(2*tz); A(14,14)=(rr+Lr/dt)*(R(15)+R(14))+(R(14)-R(15))*Lr*v0/(2*tz)+1/dt; A(14,15)=-(rr+Lr/dt)*R(15)+(Lr*(R(16)+R(15))+1)*v0/(2*tz); A(14,16)=-R(16)*Lr*v0/(2*tz); % n=15 A(15,13)=R(14)*Lr*v0/(2*tz); A(15,14)=-(rr+Lr/dt)*R(15)-(Lr*(R(15)+R(14))+1)*v0/(2*tz); A(15,15)=(rr+Lr/dt)*(R(16)+R(15))+(R(15)-R(16))*Lr*v0/(2*tz)+1/dt; A(15,16)=-(rr+Lr/dt)*R(16)+(Lr*(R(17)+R(16))+1)*v0/(2*tz); A(15,17)=-R(17)*Lr*v0/(2*tz); % n=16 A(16,14)=R(15)*Lr*v0/(2*tz); A(16,15)=-(rr+Lr/dt)*R(16)-(Lr*(R(16)+R(15))+1)*v0/(2*tz); A(16,16)=(rr+Lr/dt)*(R(17)+R(16))+(R(16)-R(17))*Lr*v0/(2*tz)+1/dt; A(16,17)=-(rr+Lr/dt)*R(17)+(Lr*(R(18)+R(17))+1)*v0/(2*tz); A(16,18)=-R(18)*Lr*v0/(2*tz); % n=17 A(17,15)=R(16)*Lr*v0/(2*tz); A(17,16)=-(rr+Lr/dt)*R(17)-(Lr*(R(17)+R(16))+1)*v0/(2*tz); A(17,17)=(rr+Lr/dt)*(R(18)+R(17))+(R(17)-R(18))*Lr*v0/(2*tz)+1/dt; A(17,18)=-(rr+Lr/dt)*R(18)+(Lr*(R(19)+R(18))+1)*v0/(2*tz); A(17,19)=-R(19)*Lr*v0/(2*tz); % n=18 A(18,16)=R(17)*Lr*v0/(2*tz); A(18,17)=-(rr+Lr/dt)*R(18)-(Lr*(R(18)+R(17))+1)*v0/(2*tz); A(18,18)=(rr+Lr/dt)*(R(19)+R(18))+(R(18)-R(19))*Lr*v0/(2*tz)+1/dt; A(18,19)=-(rr+Lr/dt)*R(19)+(Lr*(R(20)+R(19))+1)*v0/(2*tz); A(18,20)=-R(20)*Lr*v0/(2*tz); % n=19 A(19,17)=R(18)*Lr*v0/(2*tz); A(19,18)=-(rr+Lr/dt)*R(19)-(Lr*(R(19)+R(18))+1)*v0/(2*tz); A(19,19)=(rr+Lr/dt)*(R(20)+R(19))+(R(19)-R(20))*Lr*v0/(2*tz)+1/dt; A(19,20)=-(rr+Lr/dt)*R(20)+(Lr*(R(21)+R(20))+1)*v0/(2*tz); % n=20 A(20,18)=R(19)*Lr*v0/(2*tz); A(20,19)=-(rr+Lr/dt)*R(20)-(Lr*(R(20)+R(19))+1)*v0/(2*tz); A(20,20)=(rr+Lr/dt)*(R(21)+R(20))+(R(20)-R(21))*Lr*v0/(2*tz)+1/dt; % Напряжения A(21,5)=U; A(21,6)=U; A(21,9)=-U; A(21,10)=-U; A(21,11)=-U;

“Young Scientist” . #3 (50) . March 2013 Computer Science A(21,12)=-U; A(21,15)=U; A(21,16)=U; A(22,7)=U; A(22,8)=U; A(22,8)=U; A(22,9)=U; A(22,10)=U; A(22,13)=-U; A(22,14)=-U; A(22,15)=-U; A(22,16)=-U; % Токи A(4,21)=-wn*Lr*v0/(2*tz); A(5,21)=-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz); A(6,21)=wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn; A(6,22)=(-wn*Lr*v0/(2*tz))*(-1); A(7,21)=wn*Lr*v0/(2*tz); A(7,22)=(-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz))*(-1); A(8,22)=(wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn)*(-1); A(8,23)=-wn*Lr*v0/(2*tz); A(9,22)=(wn*Lr*v0/(2*tz))*(-1); A(9,23)=-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz); A(10,21)=(-wn*Lr*v0/(2*tz))*(-1); A(10,23)=wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn; A(11,21)=(-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz))*(-1); A(11,23)=wn*Lr*v0/(2*tz); A(12,21)=(wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn)*(-1); A(12,22)=-wn*Lr*v0/(2*tz); A(13,21)=(wn*Lr*v0/(2*tz))*(-1); A(13,22)=-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz); A(14,22)=wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn; A(14,23)=(-wn*Lr*v0/(2*tz))*(-1); A(15,22)=wn*Lr*v0/(2*tz); A(15,23)=(-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz))*(-1); A(16,23)=(wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn)*(-1); A(17,23)=wn*Lr*v0/(2*tz); A(21,21)=rs+Ls/dt; A(21,23)=-(rs+Ls/dt); A(22,22)=-(rs+Ls/dt); A(22,23)=rs+Ls/dt; A(23,21)=1; A(23,22)=1; A(23,23)=1; % Матрица свободных членов S=[((Lr*(R(2)+R(1))+1)/dt)*X(1)+(-Lr*R(2)/dt)*X(2); (-Lr*R(2)/dt)*X(1)+((Lr*(R(3)+R(2))+1)/dt)*X(2)+(-Lr*R(3)/dt)*X(3); (-Lr*R(3)/dt)*X(2)+((Lr*(R(4)+R(3))+1)/dt)*X(3)+(-Lr*R(4)/dt)*X(4); (-Lr*R(4)/dt)*X(3)+((Lr*(R(5)+R(4))+1)/dt)*X(4)+(-Lr*R(5)/dt)*X(5); -Lr*wn*IS(5)/dt+(-Lr*R(5)/dt)*X(4)+((Lr*(R(6)+R(5))+1)/dt)*X(5)+ (-Lr*R(6)/dt)*X(6); -Lr*wn*IS(6)/dt+(-Lr*R(6)/dt)*X(5)+((Lr*(R(7)+R(6))+1)/dt)*X(6)+ (-Lr*R(7)/dt)*X(7); -Lr*wn*IS(7)/dt+(-Lr*R(7)/dt)*X(6)+((Lr*(R(8)+R(7))+1)/dt)*X(7)+ (-Lr*R(8)/dt)*X(8); 137

“Young Scientist” . #3 (50) . March 2013 Computer Science<br />

A(21,12)=-U;<br />

A(21,15)=U;<br />

A(21,16)=U;<br />

A(22,7)=U;<br />

A(22,8)=U;<br />

A(22,8)=U;<br />

A(22,9)=U;<br />

A(22,10)=U;<br />

A(22,13)=-U;<br />

A(22,14)=-U;<br />

A(22,15)=-U;<br />

A(22,16)=-U;<br />

% Токи<br />

A(4,21)=-wn*Lr*v0/(2*tz);<br />

A(5,21)=-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz);<br />

A(6,21)=wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn;<br />

A(6,22)=(-wn*Lr*v0/(2*tz))*(-1);<br />

A(7,21)=wn*Lr*v0/(2*tz);<br />

A(7,22)=(-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz))*(-1);<br />

A(8,22)=(wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn)*(-1);<br />

A(8,23)=-wn*Lr*v0/(2*tz);<br />

A(9,22)=(wn*Lr*v0/(2*tz))*(-1);<br />

A(9,23)=-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz);<br />

A(10,21)=(-wn*Lr*v0/(2*tz))*(-1);<br />

A(10,23)=wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn;<br />

A(11,21)=(-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz))*(-1);<br />

A(11,23)=wn*Lr*v0/(2*tz);<br />

A(12,21)=(wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn)*(-1);<br />

A(12,22)=-wn*Lr*v0/(2*tz);<br />

A(13,21)=(wn*Lr*v0/(2*tz))*(-1);<br />

A(13,22)=-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz);<br />

A(14,22)=wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn;<br />

A(14,23)=(-wn*Lr*v0/(2*tz))*(-1);<br />

A(15,22)=wn*Lr*v0/(2*tz);<br />

A(15,23)=(-(rr+Lr/dt)*wn-wn*Lr*v0/(2*tz))*(-1);<br />

A(16,23)=(wn*Lr*v0/(2*tz)-(rr+Lr/dt)*wn)*(-1);<br />

A(17,23)=wn*Lr*v0/(2*tz);<br />

A(21,21)=rs+Ls/dt;<br />

A(21,23)=-(rs+Ls/dt);<br />

A(22,22)=-(rs+Ls/dt);<br />

A(22,23)=rs+Ls/dt;<br />

A(23,21)=1;<br />

A(23,22)=1;<br />

A(23,23)=1;<br />

% Матрица свободных членов<br />

S=[((Lr*(R(2)+R(1))+1)/dt)*X(1)+(-Lr*R(2)/dt)*X(2);<br />

(-Lr*R(2)/dt)*X(1)+((Lr*(R(3)+R(2))+1)/dt)*X(2)+(-Lr*R(3)/dt)*X(3);<br />

(-Lr*R(3)/dt)*X(2)+((Lr*(R(4)+R(3))+1)/dt)*X(3)+(-Lr*R(4)/dt)*X(4);<br />

(-Lr*R(4)/dt)*X(3)+((Lr*(R(5)+R(4))+1)/dt)*X(4)+(-Lr*R(5)/dt)*X(5);<br />

-Lr*wn*IS(5)/dt+(-Lr*R(5)/dt)*X(4)+((Lr*(R(6)+R(5))+1)/dt)*X(5)+<br />

(-Lr*R(6)/dt)*X(6);<br />

-Lr*wn*IS(6)/dt+(-Lr*R(6)/dt)*X(5)+((Lr*(R(7)+R(6))+1)/dt)*X(6)+<br />

(-Lr*R(7)/dt)*X(7);<br />

-Lr*wn*IS(7)/dt+(-Lr*R(7)/dt)*X(6)+((Lr*(R(8)+R(7))+1)/dt)*X(7)+<br />

(-Lr*R(8)/dt)*X(8);<br />

137

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!