ИЭ (станции+реле лектор Горячевский) / Лабы / лаба 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