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

%Этот файл отличается тем, что здесь используется расчет в несколько шагов
%в цикле for
function []=lab_2_for()
format short
clc,clear;%очистить все предыдущие результаты
%исходные данные не меняются
S = [
     0.6-0.9i   ;
     0.8+0.3i   ;
     -1+0i      ;
     -0.1+0.7i  ;
     -0.9+1i    ;
     0+0i       ];

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';

%конец исходных данных
    for i = 1:100
        J = conj ( S ./ (deltaU+U0) /sqrt(3) ); 
        deltaU = sqrt(3) * inv (Y) * J ;
    end
Iv = A' * deltaU .* Yv /sqrt(3);
%вывод результатов
disp('Напряжения в узлах U, кВ')
abs(deltaU+U0)
disp('токи в ветвях I, кА')
abs(Iv)

end