лабораторная работа / мэмс / cornus_text_4lw

R1=20;

R2=10;

R3=30;

L3=0.002;

C2=3e-04;

f=50;

ph1=60;

omeg=2*pi*f;

E1=220*exp(ph1*(pi/180)*i);

h=0.04/200;

%impedance of before commutation contour

XLb=omeg*L3;

Zb=R1+R3+i*XLb;

%N=input('number of iterations N=')

%initial conditions

I3L=E1/(Zb);

Uc(1)=0;

t(1)=0;

I3LM=abs(I3L);

phI3L=(angle(I3L)/pi*180);

I3(1)=current(I3LM,f,phI3L,t(1));

%definition of conductances

G11=1/R1;G33=(1/R3)+(h/L3);G22=(1/R2)+(C2/h);

G13=1/R3;G31=G13;

G12=1/R2;G21=G12;G32=0;G23=G32;

G=[G11 -G12 -G13;

-G21 G22 -G21;

-G31 -G32 G33];

for k=2:201

t(k)=t(k-1)+h;

Emf(k)=voltas(220,f,ph1,t(k));

I11=Emf(k)/R1;

I22=Uc(k-1)*(C2/h);

I33=I3(k-1);

%I3(k)=current(I3LM,f,phI3L,t(k));

I=[I11;I22;I33];

Pot=G\I;

I1(k)=(Emf(k)-Pot(1))/R1;

I2(k)=(Pot(1)-Pot(2))/R2;

I3(k)=(Pot(1)-Pot(3))/R3;

Uc(k)=Pot(2);

end

plot(t,I1,'.-',t,I2,'--',t,I3,t,Uc,'y');

grid on;

+++++++++++++++++

function em=voltas(EM,f,ph1,t)

em=sqrt(2)*EM*sin(2*pi*f*t+ph1);

++++++++++++++

function im=current(IM,f,phI,t)

im=sqrt(2)*IM*sin(2*pi*f*t+phI);