?????????????? ?????? ? ????????? ?????????? ??????? ?????

Все этапы можно свести в одну формулу для расчета по частям контурных
схем:
ik= (1 – ykk · A · Z · At)·ykk · ek
где Z = (At · ykk · A) -1 , матрица сопротивлений цепи пересечений, ykk – матрица
проводимостей контуров изолированных подсхем.
Формирование топологических матриц, расчет подсхем и расчет соединенной
схемы по приведенным выше формулам:
clear;
R1=1;R2=2;R3=3;R4=4;R5=5;R6=6;R7=7;R8=8;R9=9;R10=10;R11=11;
// ввод данных для 1-ой подсхемы
e1=1;
ev1=[e1,0,0,0,0,0];
ev2=[e1,0,0,0,0,0,0,0];
ev3=[e1,0,0,0,0,0,0,0,0,0,0];
z1=[R1,R2,R3,R4,R5,R6];
z1=diag(z1);
// контурные матрицы 1-й подсхемы
C1=[1, 0, 0, 0; 0, 0, 1, 0; 1, 0,-1, 0; -1, 1, 0, 0; 0, 0,-1, 1; 0,-1, 0, 1];
// расчет контурных сопротивлений и контурных э.д.с 1-й подсхемы
zk1=C1'*z1*C1; ek1=C1'*ev1';
// ввод данных для 2-й подсхемы
z2=[R1,R2,R3,R4,R5,R6,R7,R8];
z2=diag(z2);
// контурные матрицы 2-й подсхемы
C2=[0, 0, 1, 0; 1, 0, 0, 0; -1, 0, 1, 0; 1, 0, 0, 0; 0, 0,-1, 1; -1, 1, 0, 0; 0, -1, 0, 1; 0, 1, 0, 0];
// расчет контурных сопротивлений и контурных э.д.с 2-й подсхемы
zk2=C2'*z2*C2; ek2=C2'*ev2';
// ввод данных для 3-й подсхемы
z3=[R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11];
z3=diag(z3);
// контурные матрицы 3-й подсхемы
C3=[0, 0, 1, 0, 0; 0, 0,-1,-1, 0; 0, 0, 0, 1,-1; 0, 0, 0, 0, 1; 1, 0, 1, 0, 0; 0, 1, 0, 1, 0;
1, 0, 0, 0, 0; 1, 1, 0, 0, 0; 0, 1, 0, 0, 0; 1, 0, 0, 0, 0; 0, 1, 0, 0, 0];
// расчет контурных сопротивлений и контурных э.д.с 3-й подсхемы
zk3=C3'*z3*C3; ek3=C3'*ev3';
// ввод матриц в соотв с п2-п4
A11=[0,0; 0,0; 1,0; 0,1];
A12=[0,0,0; 1,0,0; 0,1,0; 0,0,1];
A13=[0,0,0; 0,0,0; 1,0,0; 0,1,0; 0,0,1];
A21=[1,0,0; 0,1,1];
A22=[0,0,0,1; 1,0,0,0; 0,1,1,0];
A23=[1,0,0; 0,1,0; 0,0,1];
A31=[1,0,0,0,0; 0,1,0,0,0; 0,0,1,0,0];
A32=[1,0,0,0,0; 0,1,0,0,0; 0,0,0,1,0; 0,0,0,0,1];
A33=[0,0,1,0,0; 0,0,0,1,0; 0,0,0,0,1];
A1=A11*A21*A31; A2=A12*A22*A32; A3=A13*A23*A33;
A=[A1;A2;A3]
// расчет граничных контурных проводимостей подсхем
yk1=zk1^(-1); yk2=zk2^(-1); yk3=zk3^(-1);
ik1=yk1*ek1; ik2=yk2*ek2; ik3=yk3*ek3;
114