Молодой учёный
Молодой учёный Молодой учёный
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);
“Young Scientist” . #3 (50) . March 2013 Computer Science (X(4)-X(2))*Ir(3)/(2*tz); (X(5)-X(3))*Ir(4)/(2*tz); (X(6)-X(4))*Ir(5)/(2*tz); (X(7)-X(5))*Ir(6)/(2*tz); (X(8)-X(6))*Ir(7)/(2*tz); (X(9)-X(7))*Ir(8)/(2*tz); (X(10)-X(8))*Ir(9)/(2*tz); (X(11)-X(9))*Ir(10)/(2*tz); (X(12)-X(10))*Ir(11)/(2*tz); (X(13)-X(11))*Ir(12)/(2*tz); (X(14)-X(12))*Ir(13)/(2*tz); (X(15)-X(13))*Ir(14)/(2*tz); (X(16)-X(14))*Ir(15)/(2*tz); (X(17)-X(15))*Ir(16)/(2*tz); (X(18)-X(16))*Ir(17)/(2*tz); (X(19)-X(17))*Ir(18)/(2*tz); (X(20)-X(18))*Ir(19)/(2*tz); (-X(19))*Ir(20)/(2*tz)]; % Скорость v0=v0+(sum(F)/m)*dt; end; % Построение графиков k=0:K; subplot(2,1,1); plot(k*dt,v); title('Скорость'); xlabel('t, с'); ylabel('v, м/с'); grid on subplot(2,1,2); plot(k*dt,f); title('Электромагнитное усилие'); xlabel('t, с'); ylabel('F, Н'); grid on end Результаты моделирования представлены в таблицах 1, 2 и 3, а также на рис. 3. Таблица 1. Значения элементов матрицы А при k = 1 a1,1 1,313E+05 a1,2 -6,514E+04 a1,3 -7,370E-04 a2,1 -6,514E+04 a2,2 7,265E+04 a2,3 -6,514E+03 a2,4 -7,370E-05 a3,1 7,370E-03 a3,2 -6,514E+03 a3,3 8,165E+03 a3,4 -6,514E+02 a3,5 -1,474E-05 a4,2 7,370E-04 a4,3 -6,514E+02 a4,4 1,782E+03 a4,5 -1,303E+02 a4,6 -1,474E-05 a5,3 7,370E-05 a5,4 -1,303E+02 a5,6 1,261E+03 a5,7 -1,303E+02 a5,8 -1,474E-05 a16,14 1,474E-05 a16,15 -1,303E+02 a16,16 1,261E+03 a16,17 -1,303E+02 a16,18 -7,370E-05 a17,15 1,474E-05 139
- Page 94 and 95: 88 Технические наук
- Page 96 and 97: 90 Технические наук
- Page 98 and 99: 92 Технические наук
- Page 100 and 101: 94 Технические наук
- Page 102 and 103: 96 Технические наук
- Page 104 and 105: 98 Технические наук
- Page 106 and 107: 100 Технические наук
- Page 108 and 109: 102 Технические наук
- Page 110 and 111: 104 Технические наук
- Page 112 and 113: 106 Технические наук
- Page 114 and 115: 108 Технические наук
- Page 116 and 117: 110 Технические наук
- Page 118 and 119: 112 Технические наук
- Page 120 and 121: 114 Технические наук
- Page 122 and 123: 116 Технические наук
- Page 124 and 125: 118 Технические наук
- Page 126 and 127: 120 Технические наук
- Page 128 and 129: 122 Технические наук
- Page 130 and 131: 124 Технические наук
- Page 132 and 133: 126 Технические наук
- Page 134 and 135: 128 Технические наук
- Page 136 and 137: 130 Информатика «Мол
- Page 138 and 139: 132 Информатика «Мол
- Page 140 and 141: 134 Информатика «Мол
- Page 142 and 143: 136 Информатика «Мол
- Page 146 and 147: 140 Информатика «Мол
- Page 148 and 149: 142 Информатика «Мол
- Page 150 and 151: 144 Информатика «Мол
- Page 152 and 153: 146 Информатика «Мол
- Page 154 and 155: 148 Информатика «Мол
- Page 156 and 157: 150 Химия «Молодой у
- Page 158 and 159: 152 Биология «Молодо
- Page 160 and 161: 154 Биология «Молодо
- Page 162 and 163: 156 Биология «Молодо
- Page 164 and 165: 158 Биология «Молодо
- Page 166 and 167: 160 Биология «Молодо
- Page 168 and 169: 162 Экология «Молодо
- Page 170 and 171: 164 Экология «Молодо
- Page 172 and 173: 166 Экология «Молодо
- Page 174 and 175: 168 Экология «Молодо
- Page 176 and 177: 170 Экология «Молодо
- Page 178 and 179: 172 Гeография «Молод
- Page 180 and 181: 174 Гeография «Молод
- Page 182 and 183: 176 Гeография «Молод
- Page 184 and 185: 178 Гeография «Молод
- Page 186 and 187: 180 Гeография «Молод
- Page 188 and 189: 182 Гeография «Молод
- Page 190 and 191: 184 Гeография «Молод
- Page 192: Молодой ученый Еже
“Young Scientist” . #3 (50) . March 2013 Computer Science<br />
(X(4)-X(2))*Ir(3)/(2*tz);<br />
(X(5)-X(3))*Ir(4)/(2*tz);<br />
(X(6)-X(4))*Ir(5)/(2*tz);<br />
(X(7)-X(5))*Ir(6)/(2*tz);<br />
(X(8)-X(6))*Ir(7)/(2*tz);<br />
(X(9)-X(7))*Ir(8)/(2*tz);<br />
(X(10)-X(8))*Ir(9)/(2*tz);<br />
(X(11)-X(9))*Ir(10)/(2*tz);<br />
(X(12)-X(10))*Ir(11)/(2*tz);<br />
(X(13)-X(11))*Ir(12)/(2*tz);<br />
(X(14)-X(12))*Ir(13)/(2*tz);<br />
(X(15)-X(13))*Ir(14)/(2*tz);<br />
(X(16)-X(14))*Ir(15)/(2*tz);<br />
(X(17)-X(15))*Ir(16)/(2*tz);<br />
(X(18)-X(16))*Ir(17)/(2*tz);<br />
(X(19)-X(17))*Ir(18)/(2*tz);<br />
(X(20)-X(18))*Ir(19)/(2*tz);<br />
(-X(19))*Ir(20)/(2*tz)];<br />
% Скорость<br />
v0=v0+(sum(F)/m)*dt;<br />
end;<br />
% Построение графиков<br />
k=0:K;<br />
subplot(2,1,1);<br />
plot(k*dt,v);<br />
title('Скорость');<br />
xlabel('t, с');<br />
ylabel('v, м/с');<br />
grid on<br />
subplot(2,1,2);<br />
plot(k*dt,f);<br />
title('Электромагнитное усилие');<br />
xlabel('t, с');<br />
ylabel('F, Н');<br />
grid on<br />
end<br />
Результаты моделирования представлены в таблицах 1, 2 и 3, а также на рис. 3.<br />
Таблица 1. Значения элементов матрицы А при k = 1<br />
a1,1 1,313E+05 a1,2 -6,514E+04 a1,3 -7,370E-04 a2,1 -6,514E+04<br />
a2,2 7,265E+04 a2,3 -6,514E+03 a2,4 -7,370E-05 a3,1 7,370E-03<br />
a3,2 -6,514E+03 a3,3 8,165E+03 a3,4 -6,514E+02 a3,5 -1,474E-05<br />
a4,2 7,370E-04 a4,3 -6,514E+02 a4,4 1,782E+03 a4,5 -1,303E+02<br />
a4,6 -1,474E-05 a5,3 7,370E-05 a5,4 -1,303E+02 a5,6 1,261E+03<br />
a5,7 -1,303E+02 a5,8 -1,474E-05 a16,14 1,474E-05 a16,15 -1,303E+02<br />
a16,16 1,261E+03 a16,17 -1,303E+02 a16,18 -7,370E-05 a17,15 1,474E-05<br />
139
Change language