A7TZrRHEJI

r2a)+i20*...

(r1b+2*r2a)-i30*(2*r1b+r2a)+i40*(2*r1b+r2a)+i50*(r2a-r1b)- i60*(r1b+2*r2a))+...

Ld*(ea-ec+r2a*(i10+i20+i30+i40+i50+i60))); ni20=R*(ea-ec)+ea*r1b+eb*r2a-ec*(r1b+r2a); %"deni2" di22=(2*L*(R+r1b+r2a+r2c)+Ld*(r2a+r2c)); di21=(R*(r2a+r2c)+r1b*(r2a+r2c)+r2a*r2c); end

%i3-i5

if s<=ss6; %"numi3"

ni32=L*(Ld*(2*i10-i20+2*i30+3*i60)-3*L*(i20-i30+i40-i60)); ni31=(L*(R*(i10-2*i20+i30-i40-i50+2*i60)-ea-eb+2*ec+i10*(r2a-r1b)- i20*(r1b+...

2*r2a)+i30*(2*r1b+r2a)-i40*(2*r1b+r2a)+i50*(r1b-r2a)+i60*(r1b+2*r2a))- Ld*(eb-ec-r1b*(i10+i20+i30+i40+i50+i60))); ni30=R*(ec-eb)-ea*r1b-eb*r2a+ec*(r1b+r2a);

%"deni3"

di32=(2*L*(R+r1b+r1c+r2a)+Ld*(r1b+r1c));

di31=(R*(r1b+r1c)+r1b*(r1c+r2a)+r1c*r2a); end

if s>ss5; %"numi5"

ni52=L*(Ld*(3*i20+2*i40+2*i50-i60)-3*L*(i10-i20-i50+i60)); ni51=(Ld*(eb-ec+r1c*(i10+i20+i30+i40+i50+i60))-L*(R*(i10-2*i20+i30- i40-...

i50+2*i60)+ea-2*eb+ec+i10*(2*r1c+r2a)-i20*(r1c+2*r2a)+i30*(r2a-r1c)+...

i40*(r1c-r2a)-i50*(2*r1c+r2a)+i60*(r1c+2*r2a))); ni50=R*(eb-ec)-ea*r1c+eb*(r1c+r2a)-ec*r2a; %"deni5" di52=(2*L*(R+r1b+r1c+r2a)+Ld*(r1b+r1c)); di51=(R*(r1b+r1c)+r1b*(r1c+r2a)+r1c*r2a);

end

if s>ss6; ni32=L*(Ld*(i10+2*i20+2*i30+3*i50)-L*(i10-3*(i30-i40+i50)));

21

ni31=(L*(R*(i20+i30-i40+2*i50-i60)-2*ea+eb+ec-i10*r2b+i20*(r1c- r2b)+...

i30*(r1c+2*r2b)-i40*(r1c+2*r2b)+i50*(2*r1c+r2b)+i60*(r2b-r1c))-Ld*(ea-...

eb-r2b*(i10+i20+i30+i40+i50+i60))); ni30=R*(eb-ea)-ea*(r1c+r2b)+eb*r1c+ec*r2b; %"deni3" di32=(2*L*(R+r1c+r2a+r2b)+Ld*(r2a+r2b)); di31=(R*(r2a+r2b)+r1c*(r2a+r2b)+r2a*r2b); end

%"numi6"

if (s<=ss7)&(s<per);NN=1; ni62=L*(Ld*(i10+2*i40-i50+2*i60)-L*(i10+3*(i20+i50-i60))); ni61=(Ld*(ea-eb+r2a*(i10+i20+i30+i40+i50+i60))-L*(R*(i20+i30- i40+2*i50-...

i60)-ea+2*eb-ec+i10*r2a+i20*(r1c+2*r2a)+i30*(r1c-r2a)+i40*(r2a-r1c)+...

i50*(2*r1c+r2a)-i60*(r1c+2*r2a))); ni60=R*(ea-eb)+ea*r1c-eb*(r1c+r2a)+ec*r2a; %"deni6" di62=(2*L*(R+r1c+r2a+r2b)+Ld*(r2a+r2b)); di61=(R*(r2a+r2b)+r1c*(r2a+r2b)+r2a*r2b); end

if (s<ds)&(s>ss77+1); ni12=L*(3*L*(i10-i30+i40-i50)+Ld*(2*i10+2*i20-i30+3*i40)); ni11=(L*(R*(i10+i20-2*i30+2*i40-i50-i60)+2*ea-eb- ec+i10*(2*r1c+r2b)+i20*(r2b-r1c)-i30*(r1c+2*r2b)+...

i40*(r1c+2*r2b)-i50*(2*r1c+r2b)+i60*(r1c-r2b))+Ld*(ea- ec+r1c*(i10+i20+i30+i40+i50+i60)));

ni10=-ec*r2b; di13=L*(3*L+2*Ld);

di12=(2*L*(R+r1a+r1c+r2b)+Ld*(r1a+r1c));

di11=(R*(r1a+r1c)+r1a*(r1c+r2b)+r1c*r2b);

di10=0;

…

ni62=-(L*(3*L*(i20-i30+i40-i60)-Ld*(3*i30-i40+2*(i50+i60)))); ni61=-((L*(R*(i10+i20-2*i30+2*i40-i50-i60)+ea+eb-2*ec+i10*(r2b- r1a)+i20*(2*r1a+r2b)-i30*(r1a+2*r2b)+...

22

i40*(r1a+2*r2b)+i50*(r1a-r2b)-i60*(2*r1a+r2b))+Ld*(ea-ec-i10*r1a- i20*r1a-i30*r1a-i40*r1a-i50*r1a-i60*r1a)));

ni60=ec*(r1a+r2b);

di63=L*(3*L+2*Ld);

di62=(2*L*(R+r1a+r1c+r2b)+Ld*(r1a+r1c));

di61=(R*(r1a+r1c)+r1a*(r1c+r2b)+r1c*r2b);

di60=0;

end NN=0;

%-------------------------------------------------------

ni1=[ni12 ni11 ni10]; di1=[di13 di12 di11 0];

ni4=[ni42 ni41 ni40]; di4=[di13 di42 di41 0];

ni2=[ni22 ni21 ni20]; di2=[di13 di22 di21 0];

ni5=[ni52 ni51 ni50]; di5=[di13 di52 di51 0];

ni3=[ni32 ni31 ni30]; di3=[di13 di32 di31 0];

ni6=[ni62 ni61 ni60]; di6=[di13 di62 di61 0];

y1=impulse(ni1,di1,t);y2=impulse(ni2,di2,t);y3=impulse(ni3,di3,t);

y4=impulse(ni4,di4,t);y5=impulse(ni5,di5,t);y6=impulse(ni6,di6,t); g1=y1(2,1); g2=y2(2,1); g3=y3(2,1); g4=y4(2,1); g5=y5(2,1);g6=y6(2,1); yy6(s,1)=g6;yy1(s,1)=g1;

yy2(s,1)=g2;yy3(s,1)=g3;

yy4(s,1)=g4;yy5(s,1)=g5;

z1=ss1+z11;z2=ss2+z11;z3=ss3+z11;z4=ss4+z11;z5=ss5+z11;z6=ss6+z11;z

23

7=ss7+z11;

if((s>=ss1)&(s<ss2))&(s<per);

g66=g6;g11=g1;g22=g2;g33=g3;g44=g4;g55=g5; end

if (s>=ss2 & s<ss3);

if g1<a; g11=0;g44=g4+g1; else g11=g1;g22=g2;g33=g3;g44=g4;g55=g5;g66=g6; end end

if (s>=ss3 & s< ss4);

if g4<a; g44=0;g22=g2+g4; else g11=g1;g33=g3;g55=g5;g66=g6;g44=g4;g22=g2; end end

if (s>=ss4 & s < ss5);

if g2<a; g22=0;g55=g2+g5; else g55=g5;g22=g2;g11=g1;g33=g3;g44=g4;g66=g6; end end

if (s >=ss5 & s<ss6);

if g5<a; g55=0;g33=g3+g5; else g33=g3;g55=g5;g11=g1;g22=g2;g44=g4;g66=g6; end end

if (s >=ss6 & s <ss7);

if g3<a; g33=0;g66=g6+g3; else g33=g3;g66=g6;g11=g1;g22=g2;g44=g4;g55=g5; end end

if((s>=ss1)&(s<ss2))&(s>=per);

if g6<a; g66=0; g11=g1+g6; g22=g2; g33=g3; g44=g4; g55=g5; else g66=g6;g11=g1;g22=g2;g33=g3;g44=g4;g55=g5; end

end

%Выбор последнего значения тока и присоединение к вектору вычис-

%ленных

i1=[i1;g11]; i2=[i2;g22]; i3=[i3;g33]; i4=[i4;g44]; i5=[i5;g55]; i6=[i6;g66]; % Моделирование алгоритма переключений

if s==ss1; end

if (s >= ss1)&(s<ss1+3); r1a=rl;r2b=rl;

24

end % 1-й интервал

if (s>z1 & g1+g6 <0.1);r2b=rm; end

if s == ss2; r2c=rl;r1a=rl; end %2-й интервал

if (g4+g1<0.1 & s>z2); r1a=rm; end

if s == ss3; r1b=rl;r2c=rl; end % 3-й интервал

if (s>z3 & g2+g4<0.1);r2c=rm; end

if s == ss4; r2a=rl; r1b=rl; end % 4-й интервал

if (s>z4 & g5+g2<0.1); r1b=rm; end

if s == ss5; r1c=rl; r2a=rl; end % 5-й интервал

if (s>z5 & g3+g5<0.1); r2a=rm; end

if s == ss6; r2b=rl; r1c=rl; end% 6-й интервал

if (s>ss66 & g6+g3<0.1); r1c=rm; ds=s; end

i10=g1;i20=g2;i30=g3;i40=g4;i50=g5;i60=g6;

%Накопление значений сопротивлений вентилей

R1a(1,s)=r1a;R2a(1,s)=r2a;R1b(1,s)=r1b;R2b(1,s)=r2b;R1c(1,s)=r1c;R2c(1,

s)=r2c;

%Построение графиков токов

axis([0,S,-20,200]); grid on

plot(s,g11,'or','LineWidth',2),title('Switch currents'),xlabel('Number of step'),ylabel('IVS, A')

hold on plot(s,g22,'og','LineWidth',2),plot(s,g33,'ob','LineWidth',2), plot(s,g44,'xr','LineWidth',2) plot(s,g55,'xg','LineWidth',2),plot(s,g66,'xb','LineWidth',2) pause(0.1)

25

s=s+1;

end % Конец цикла hold off

shg pause

I=(i1+i2+i3+i4+i5+i6);

plot(I,'k','LineWidth',2),grid,title('Id') pause

%Напряжение на RL нагрузке и на активном сопротивлении

ULd=Ld*gradient(I',dt); UR=I'.*RR; Ud=ULd+UR;

plot(Ud,'m','LineWidth',3),grid,title('U load'),pause Sco1=ceil(S/6/N);

Ud1=Ud(Sco1:S);

%Результат: матрица основных сигналов

sv=[i1 i2 i3 i4 i5 i6 I Ud' Ea Eb Ec]; Ud2=sqrt(sum(Ud1.^2)/length(Ud1)); % Действующее значение Ud % Конец программы

26

27

Кошелев Петр Александрович Парамонов Сергей Владимирович

Компьютерно-информационные технологии Анализа и синтеза электротехнических комплексов (ЭТК)

Часть 2

Нелинейные и нестационарные динамические системы Учебное пособие

Редактор Э. К. Долгатов

————————————————————————————————

Подписано в печать 23.08.16. Формат 60×84 1/16. Бумага офсетная. Печать цифровая. Печ. л. 1,75.

Гарнитура «Times New Roman». Тираж 45 экз. Заказ 218.

————————————————————————————————

Издательство СПбГЭТУ «ЛЭТИ» 197376, С.-Петербург, ул. Проф. Попова, 5

28