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

134 Информатика
«Молодой учёный» . № 3 (50) . Март, 2013 г.
В матричной форме это уравнение можно получить следующим образом. Так как в матрице-столбце токи ia, ib и ic занимают
положения 21, 22, и 23 соответственно, то в матрице А в 23-ей строке необходимо приравнять элементы a23,1, a23,2, …, a23,20 к нулю, а элементы a23,20, a23,21, a23,23 к единице. То есть 23,21 23,22 23,23 1
a a a = = = и s 23 = 0 .
Произведем построение математической модели расчета линейного асинхронного двигателя метода Гаусса-Жордана
с использованием языка программирования MatLab. Ниже приведен пример кода:
function Linear_AD1
% Исходные данные линейного асинхронного двигателя
Rb=0.1003*10^7;
rs=19;
Ls=0.074;
rr=9.269*10^-5;
Lr=0.0372*10^-5;
dt=0.001;
tz=9.769*10^-3;
m=1.9;
v0=0;
wn=200;
f=50;
w=2*pi*f;
U=wn/dt;
Um=310;
X=zeros(23,1);
F=0;
K=input('длительность цикла k=');
% Создание матрицы сопротивлений
R=zeros(21,1);
R(1,1)=500*Rb;
R(2,1)=500*Rb;
R(3,1)=50*Rb;
R(4,1)=5*Rb;
for i=5:17
R(i,1)=Rb;
end;
R(18,1)=5*Rb;
R(19,1)=50*Rb;
R(20,1)=500*Rb;
R(21,1)=500*Rb;
for k=1:(K+1)
v(1,k)=v0; %создание вектор-строки для графика скорости
f(1,k)=sum(F); %создание вектор-строки для графика усилия
Uab=Um*cos(w*(k-1)*dt+2*pi/3);
Ubc=Um*cos(w*(k-1)*dt);
% Матрица токов статора
IS=[0,0,0,0,X(21)X(21),-X(22),-X(22),X(23),X(23),-X(21),-X(21),X(22),
X(22),X(23),-X(23),0,0,0,0];
% Создание матрицы А
A=zeros(23);
% n=1
A(1,1)=(rr+Lr/dt)*(R(2)+R(1))+(R(1)-R(2))*Lr*v0/(2*tz)+1/dt;
A(1,2)=-(rr+Lr/dt)*R(2)+(Lr*(R(3)+R(2))+1)*v0/(2*tz);
A(1,3)=-R(3)*Lr*v0/(2*tz);
% n=2
A(2,1)=-(rr+Lr/dt)*R(2)-(Lr*(R(2)+R(1))+1)*v0/(2*tz);
A(2,2)=(rr+Lr/dt)*(R(3)+R(2))+(R(2)-R(3))*Lr*v0/(2*tz)+1/dt;
A(2,3)=-(rr+Lr/dt)*R(3)+(Lr*(R(4)+R(3))+1)*v0/(2*tz);