Notes

function La5Lewandowski_iRL(R,X,U,psi)
% funkcja iRL(R,X,U,psi)
% wyprowadza przebiegi czasowe pradu nieustalonego w obwodzie RL
% po zamknieciu wylacznika
imie = 'Przemyslaw';
nazwisko = 'Lewandowski';
album = '249857';
albumNum = str2num(album);
sumaAlbum = 0;
while(albumNum > 0)
sumaAlbum = sumaAlbum + mod(albumNum,10);
albumNum=(albumNum-mod(albumNum,10))/10;
end
a = length(imie) + length(nazwisko) + sumaAlbum;
disp("Nazwa: "+imie+" "+nazwisko);
disp("a = "+a);
if nargin < 4
psi = 0; % kat nap. psi w stopniach
end
if nargin < 3
U = 200+a; % wartosc skuteczna napiecia zasilajca obwod, V
end
if nargin < 2
X = a; % reaktancja obwodu, ohm
end
if nargin < 1
R = a; % rezystancja obwodu, ohm
end
Um = sqrt(2)*U; % amplituda nap. zasilajacego
Z = sqrt(R^2+X^2); % impedancja pozorna obwodu
fi = atan(X/R); % kat impedancji
psi=psi/180*psi; % przeliczenie kata poczatkowego napiecia na radiany
f = 50; % czestotliwosc napiecia , Hz
w = 2*pi*f; % czestosc napiecia
L = X/w; % indukcyjnosc obwodu
tau = L/R; %stala zanikania
disp('Parametry obwodu:');
disp(' R= '),disp(R);
disp(' X= '),disp(X);
disp(' Z= '),disp(Z);
disp(' tau='),disp(tau);
disp(' U= '),disp(U);
disp(' psi= '),disp(psi);
disp(' fi = '),disp(fi);
% WYZNACZANIE PRZEBIEGU NIEUSTALONEGO PRADU
Tmax = 0.2 ; % max czas analizy
T = 0:0.002:Tmax; % przedzial czasu analizy od 0 do 1 sekundy
nT=length(T); % liczba punktow czasowych
p0 = zeros(1,nT); % os zerowa
Im = Um/Z; % amplituda pradu
iAC = Um/Z * sin(w.*T + psi - fi); % skladowa okresowa
iDC = - Um/Z*exp(-T./tau)*sin(psi-fi); %skladowa nieokresowa
iRL = iAC + iDC;
plot(T,iAC,'-b', T,p0,'-r');
title("psi = "+psi+", U = "+U+", R = "+R+", X = "+X), title(imie+" "+nazwisko+", "+album+", "+a+", grupa A, godz. 14.30-16.00 wykres: iAC");
xlabel('t,s'); ylabel('A');
axis normal; grid on;
pause(2);close; % pauza/zamkniecie
plot(T,iDC,'--g', T, p0,'-r');
title("psi = "+psi+", U = "+U+", R = "+R+", X = "+X), title(imie+" "+nazwisko+", "+album+", "+a+", grupa A, godz. 14.30-16.00 wykres: iDC");
xlabel('t,s'); ylabel('A');
axis normal; grid on;
pause(2);close; % pauza/zamkniecie
plot(T,iRL,'-k', T, p0,'-r');
title("psi = "+psi+", U = "+U+", R = "+R+", X = "+X), title(imie+" "+nazwisko+", "+album+", "+a+", grupa A, godz. 14.30-16.00 wykres: iRL");
xlabel('t,s'); ylabel('A');
axis normal; grid on;
pause(2);close; % pauza/zamkniecie
plot(T,iRL,'-k',T,iAC,'-b',T,iDC,'--g', T,p0,'-r');
xlabel('t,s'); ylabel('A');
title(imie+" "+nazwisko+", "+album+", grupa A, godz. 14.30-16.00");
title("a = "+a);
title("psi = "+psi+", U = "+U+", R = "+R+", X = "+X);
title(imie+" "+nazwisko+", "+album+", "+a+", grupa A, godz. 14.30-16.00 wykres: iRL(t), iAC(t), iDC(t)");
axis normal; grid on;
pause(2);close; % pauza/zamkniecie
disp('KONIEC wyprowadzania przebiegow pradu w obwodzie RL');
end