Все этапы можно свести в одну формулу для расчета по частям контурных схем: 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
y1=A1'*yk1*A1; y2=A2'*yk2*A2; y3=A3'*yk3*A3; y=y1+y2+y3; for i=1:5,y(i,i)=y(i,i)+1;end; I1=-A1'*ik1; I2=-A2'*ik2; I3=-A3'*ik3; I=I1+I2+I3; E=y^(-1)*I; e=A*E; dik1=yk1*e(1:4); dik2=yk2*e(5:8); dik3=yk3*e(9:13); i1=ik1+dik1; i2=ik2+dik2; i3=ik3+dik3; 115
- Page 1:
Федеральное агентс
- Page 4 and 5:
УДК 62-83+ 621.3.01 ББК 31.29
- Page 6 and 7:
3.5. Алгоритм расчет
- Page 8 and 9:
Введение. При тензо
- Page 10 and 11:
1. Базовые понятия и
- Page 12 and 13:
1.2. Тензорные объек
- Page 14 and 15:
пендикулярно некот
- Page 16 and 17:
странства, разъеди
- Page 18 and 19:
Мощность электриче
- Page 20 and 21:
Кроме разностных а
- Page 22 and 23:
2. Моделирование це
- Page 24 and 25:
2.1.1. Разностные схе
- Page 26 and 27:
видно, для первого
- Page 28 and 29:
Можно заметить, что
- Page 30 and 31:
13. Записать выражен
- Page 32 and 33:
независимые источн
- Page 34 and 35:
Ucg(1)=Uc0; iLg(1)=iL0; t(1)=0. do
- Page 36 and 37:
Формуле (2.1.16) соотв
- Page 38 and 39:
uL=(iL+iL0) *L / d t; else // ин
- Page 40 and 41:
1k e1-e2 ek= 2k e2-e3-e4-e5 3k e5-e
- Page 42 and 43:
2.1.1.3. Формула диффе
- Page 44 and 45:
2.1.2. Разностно-итер
- Page 46 and 47:
i= E−u o Rd⋅i o RRd Динам
- Page 48 and 49:
i n+ 1 Zc(U n+1 ) U n+1 На при
- Page 50 and 51:
Промежуточные данн
- Page 52 and 53:
2.1.2.3. Разностно-ите
- Page 54 and 55:
Схема включается н
- Page 56 and 57:
2.1.3. Управляемые ис
- Page 58 and 59:
2.2. Топологические
- Page 60 and 61:
В этом выражении ма
- Page 62 and 63:
пару, то в контурны
- Page 64 and 65:
o io o 0 ir= k ik = k ik (2.2.17)
- Page 66 and 67: he=[3 2 2 4 1 1 ]; // конечн
- Page 68 and 69: 2.3. Обобщенные урав
- Page 70 and 71: -1 z1 1 -1 -1 z2 1 z′=Сt∙z∙C
- Page 72 and 73: По контурным токам,
- Page 74 and 75: ik+J ~ k = Y ~ kk·ek где J ~ kk
- Page 76 and 77: ники напряжения, а
- Page 78 and 79: 2.3.6. Алгоритм форми
- Page 80 and 81: (4) Перевод результа
- Page 82 and 83: формируется массив
- Page 84 and 85: function [C,A,g]=formc(ta,he) maxv1
- Page 86 and 87: чения источников д
- Page 88 and 89: После выполнения р
- Page 90 and 91: program kontur !******* 1. Объя
- Page 92 and 93: !нумерация узлов k=1;
- Page 94 and 95: call sgefa(C,vetvey,vetvey,ipvtC,in
- Page 96 and 97: мутирующие цепочки
- Page 98 and 99: 2.4.1. Цепь пересечен
- Page 100 and 101: Для каждой ветви эл
- Page 102 and 103: Матрица сопротивле
- Page 104 and 105: epk= -Cpkō ∙ Uō. (2.4.13) Ве
- Page 106 and 107: подпрограммами рас
- Page 108 and 109: Блок 5. Производитс
- Page 110 and 111: 2.4.7. Сравнение скор
- Page 112 and 113: 2.5. Расчет по частям
- Page 114 and 115: а) б) Рис. 2.5.2. Расчет
- Page 118 and 119: 3. Моделирование эл
- Page 120 and 121: 3.1. Резистивная схе
- Page 122 and 123: положен обратный п
- Page 124 and 125: 3.2. Элементарная сх
- Page 126 and 127: По известной топол
- Page 128 and 129: Матрица Cōs : p s ō 1 2 3
- Page 130 and 131: 3.4.2. Расчет поля при
- Page 132 and 133: Для схем с параллел
- Page 134 and 135: • окончательно узл
- Page 136 and 137: ……… kyzl=kyzl+1; B(kyzl,it)=U
- Page 138 and 139: U 0 4.2. Резистивные м
- Page 140 and 141: Метод структурных
- Page 142 and 143: (%i6) /*Отдельные комп
- Page 144 and 145: случае при непосре
- Page 146 and 147: Приложения. Порядо
- Page 148 and 149: При нажатии на ОК п
- Page 150 and 151: Модуль общих перем
- Page 152 and 153: Подпрограмма ввода
- Page 154 and 155: CALL GWGTXT (idtext(i),texts(i)) !
- Page 156 and 157: Подпрограмма расче
- Page 158 and 159: ! индуктивности и н
- Page 160 and 161: elseif (indexs==168) then Sym=char(
- Page 162 and 163: 16. Брамеллер А., Алл
- Page 164: Сохор Юрий Николае