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

“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