CZW.

clear all; %do czyszczenia workspace
format short e %format wyswietlania liczb
m0=[0.270, 0.216;0.216, 0.481] * 1e-3; %wartosci podane w mS - sprowadzamy do j.podstawowej - S
fa=[39.2, -158.1; -158.1, 31.5]; %faza (w stopniach)
y=m0.*exp(j*fa*pi/180); %postać zespolona, *pi/180 zamienia ze stopni na radiany
disp ('|Y|') %wyświetla w konsoli to co w cudzysłowie
abs(y) %wartość bezwględna/modul l.zespolonej
disp('arg{Y}[deg]')
angle(y)*180/pi %radiany -> stopnie
r = 3.3e3; %podane w poleceniu do zadania (UWAGA NA JEDNOSTKI!)
c=15e-9; %j.w.
w=2*pi*1.8e3; %pulsacja (wzor podany w poleceniu w=2*pi*f)
z=1/(1/r+1i*w*c) %wzór na impedancje
a=-1/y(2,1)*[y(2,2),1;det(y),y(1,1)]; %macierz a
disp('|A|')
abs(a)
disp('arg{A} [deg]')
angle(a)*180/pi %radiany -> stopnie

%Potem wyliczamy z wzorow KI i KU z macierzy a czy y. Potem chyba jest
%takie porownanie macierzy KI i KU, aby byly sobie rowne. Bez wzgledu na a
%i y.

disp('na podstawie macierzy Y')
kuy=(-y(2,1)*z)/(1+y(2,2)*z);
disp('|Ku| [V/V]')
abs(kuy)
disp('arg{Ku}[deg]')
angle(kuy)*180/pi

kiy=y(2,1)/(y(1,1)+det(y)*z);
disp('|Ki| [A/A]')
abs(kiy)
disp('arg{Ki}[deg]')
angle(kiy)*180/pi

disp ('na podstawie macierzy A')
kua=z/(a(1,1)*z+a(1,2));
disp('|Ku| [V/V]')
abs(kua)
disp('arg{Ku}[deg]')
angle(kua)*180/pi

kia=-1/(a(2,1)*z+a(2,2));
disp('|Ki| [A/A]')
abs(kia)
disp('arg{Ki}[deg]')
angle(kia)*180/pi