ИЭ (станции+реле лектор Горячевский) / Лабы / лаба 3 / lab_3

function []=lab_3()
format short
clc;
S = [
     0.6-0.9i   ;
     0.8+0.3i   ;
     -1+0i      ;
     -0.1+0.7i  ;
     -0.9+1i    ;
     0          ];
U0 = 10;
deltaU  = [0;0;0;0;0;0];
Zv = [
    0.8+1.6i    ;
    1+2i        ;
    0.2+0.4i    ;
    0.6+1.2i    ;
    0.8+1.6i    ;
    0.8+1.6i    ;
    0.4+0.8i   ];
Yv = 1./Zv;
A = [
    0   0    0    1   0   -1  0;
    1   0    -1   0   0   0   0;
    0   1    0    -1  0   0   0;
    0   0    0    0   0   1   -1;
    0   -1   1    0   -1  0   0;
    0   0    0    0   1   0   0];
Y = A * diag (Yv) * A';
J = conj (S ./ (deltaU+U0) /sqrt(3));
Y = [Y J];
for i = 1 : size(Y,1)
    Y(i,:) = Y(i,:) / Y(i,i);
    for j = i+1 : size(Y,1)
        if (Y(j,i)==0)
            continue;
        end
        Y(j,:) = Y(j,:) / Y(j,i) - Y(i,:);
    end
end
for i = size(Y,1) : -1 : 1
    Y(i,:) = Y(i,:) / Y(i,i);
    for j = i-1 : -1 : 1
        if (Y(j,i)==0)
            continue;
        end
        Y(j,:) = Y(j,:) / Y(j,i) - Y(i,:);
    end
end
deltaU = Y(:,size(Y,2)) * sqrt(3);
U=deltaU+U0;
disp('Напряжения в узлах, кВ');
disp(abs(U));
disp('Токи в ветвях, кА');
Iv = A' * deltaU .* Yv /sqrt(3);
disp(abs(Iv));
end