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

function []=lab_6()

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

%конец исходных данных


%параметры расчета
%здесь попробовать разный шаг от 0.0001 до 1 и разное время расчета
h = 0.001; %шаг расчета
t = 2; %время расчета
n_max = ceil(t/h); %кол-во шагов расчета

%далее делаем n_max шагов методом обратной матрицы
for n = 0:n_max

%начальные условия
%это будет вектор токов Iv
%начальным условием будет режим, который считали в предыдущих лабах
J = conj ( S ./ (deltaU+U0) /sqrt(3) ); %правая часть МУН
deltaU = sqrt(3) * inv(Y) * J;
Iv = A' * deltaU .* Yv /sqrt(3);

%аварийное возмущение в момент времени 1 сек
if (h*n>1)
Yv(1)=0;
Zv=1./Yv;
Y = A * diag (Yv) * A';
deltaU = sqrt(3) * inv(Y) * J;
Iv = A' * deltaU .* Yv /sqrt(3);
end
%промежуточное сохранение для вывода на график
Jv(:,n+1) = abs(Iv);
end
disp('Напряжения в конце расчёта:');
disp(abs(deltaU+U0));
disp('Токи в ветвях:');
disp(abs(Iv*1000));
%вывод на график токов в ветвях
plot( ...
[0:h:t],Jv(1,:)*1000,'r',...
[0:h:t],Jv(2,:)*1000,'g',...
[0:h:t],Jv(3,:)*1000,'b',...
[0:h:t],Jv(4,:)*1000,'k',...
[0:h:t],Jv(5,:)*1000,'r--',...
[0:h:t],Jv(6,:)*1000,'g--',...
[0:h:t],Jv(7,:)*1000,'b--');
xlabel('t, c');
ylabel('I_{v}, А');
legend('I_{v1}', 'I_{v2}', 'I_{v3}', 'I_{v4}', 'I_{v5}', 'I_{v6}', 'I_{v7}');
legend('FontSize', 12, 'Location','nw');

end

%В отчёте показать графики для шага расчёта h=1; 0.1; 0.01; 0.001
%Для каждого шага h подобрать оптимальное время расчёта t для лучшей наглядности
%В выводе написать на что влияет шаг расчёта h и как от него зависит точность расчёта? Как отключение линии влияет на напряжения в узлах?
%Отчёт по этой лабораторной выложить до 17:00 12.11.2021 в группу в Тимс Путь:Файлы/Гр.90201/Л.Р. №6
%Отчёты по остальным лабораторным выкладываем в соответствующую папку в Тимс Путь:Файлы/Гр.90201/