RLC-calculations.rar_RLC circuit_RLC circuit matlab

%lnjV.m - liniju modulis if start==1 %!!!!!!!!!!!!!!!!!!!! Startup !!!!!!!!!!!!!!!!!!!! %parengia liniju failus, kai start=1 LN=round(lL/(0.3*h));% velinimosi funkciju vektoriu gabaritai % gesimo koeficientai for j=1:nL kg(j)=exp(-RLo(j)*lL(j)*0.5/WL(j)); end % didziausias velinimosi funkciju gabaritas Lnm=0;for j=1:nL Lnm=Lnm+LN(j);end % pradiniai SFF SBB indeksai Si(1)=1; for j=2:nL Si(j)=Si(j-1)+LN(j-1); end %linijos uzkraunamos pradine Uo itampa a1=1; KK=0; for j=1:nL KK=KK+LN(j); SFO=Uo/WL(j); SBO=Uo/WL(j); for k=a1:KK SFF(k)=SFO; SBB(k)=SBO; end a1=1+LN(j); end U=Uo*ones(1,N); %skaiciuojamas mazgo ekv. banginis laidumas Y=zeros(1,N); for j=1:nL a1=L(j,1);Y(a1)=Y(a1)+1/WL(j); a1=L(j,2);Y(a1)=Y(a1)+1/WL(j); end %!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! else % 1. Atsispindejusios bangos for j=1:nL a1=Si(j); SBO=(2*U(L(j,2))/WL(j))-SFF(a1); SFO=(2*U(L(j,1))/WL(j))-SBB(a1); SBB(a1)=SBO*kg(j);%SBO - atsispindejusiu nuo linijos pradzios bangos reiksmiu m-ca SFF(a1)=SFO*kg(j);%SFO - atsispindejusiu nuo linijos galo bangos reiksmiu m-ca end % 2.Modeliuojamas elektromagnetiniu bangu plitimas %surandamas SFF indeksas a1=0; for j=1:nL % nL - liniju skaicius a1=a1+LN(j); if Si(j)<a1 Si(j)=Si(j)+1;else Si(j)=a1-LN(j)+1;end end % 3.Surandamos mazgu ekvivalentiniu sroves versmiu reiksmes % sumines SFF SBB reiksmes I=zeros(1,N); % I - mazgu ekvivalentiniu sroves versmiu m-ca for j=1:nL a1=Si(j); I(L(j,1))=I(L(j,1))+SBB(a1); I(L(j,2))=I(L(j,2))+SFF(a1); end end