%ogniwo

close all,clear all

Y=zeros(4);r1=10e3;r2=10e3;r3=479.7;r4=1.426e3;r5=10e3;

c1=7.019e-6;c2=1.125e-3;

w=linspace(0.01,10,300);

Ip=zeros(4,1);Ip(1,1)=1;

for i=1:300

s=j*w(i);

Y(1,1)=1/r4+s*c2;Y(1,3)=-1/r4;Y(1,4)=-s*c2;

Y(2,2)=-s*c1;Y(2,3)=1/r1+s*c1;Y(2,4)=-1/r1;

Y(3,1)=-1/r4;Y(3,3)=1/r4+1/r5;

Y(4,1)=-s*c2;Y(4,2)=-1/r2;Y(4,3)=-1/r1;Y(4,4)=1/r1+1/r2+1/r3+s*c2;

V=Y\Ip;

uwy(i)=20*log10(abs(V(2,1)/V(1,1)));

end

plot(w,uwy);grid