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

case 'R' then //если 'R' то резистор
Z(i,i)=newshema(k+1);
case 'C' then //если 'С' то емкость
kc=kc+1; nomc(kc)=i; //формирование списка емкостей
Z(i,i)=dt/newshema(k+1);
case 'L' then //если 'L' то индуктивность
kL=kL+1; nomL(kL)=i; //формирование списка индуктивностей
Z(i,i)=newshema(k+1)/dt;
else //если не распознана ветвь
'Abort!!! There is no model '+sim(i)
abort
end;
// формирование векторов источников ветвей
e(i)=newshema(k+2); J(i)=newshema(k+3); k=k+5;
end;
/ / формирование топологических матриц
yz1=max (ta); yz2=max (he); yzlov=max (yz1,yz2) //определение числа узлов схемы
[C,A,g]= formc (ta',he');
/ / в ычисление матриц сопротивлений и проводимостей схемы
Zr=C'*Z*C;
Yk=Zr(yzlov:arc,yzlov:arc)^(- 1);
zok=Zr(1:yzlov-1,yzlov:arc); zko=Zr(yzlov:arc,1:yzlov-1);
zo =Zr(1:yzlov-1,1:yzlov-1); zor=zo-zok*Yk*zko;
for n=1:NN+1, / /Расчет в цикле по времени
er=C'*e; Jr=A'*J;
ek=er(yzlov:arc); eo=er(1:yzlov-1);
Jk=Jr(yzlov:arc); Jo=Jr(1:yzlov-1);
ik=Yk*ek-Yk*zko*Jo-Jk; Uo=zor*Jo+zok*Yk*ek-eo;
iv=C(1:arc,yzlov:arc)*ik; Uv=A(1:arc,1:yzlov-1)*Uo;
/ / начальные условия для следующего шага
if kL 0 then
for j=1:kL
J(nomL(j))=- iv(nomL(j)); //начальные условия для индуктивностей
end;
end;
if kc0 then
for j=1:kc
e(nomc(j))=- Uv(nomc(j)); //начальные условия для емкостей
end;
end;
/ / в ывод в массивы для построения графиков
ivL(n)=iv(4); uvc(n)=Uv(2);
end; / / о кончание цикла по времени
plot2d (t, [uvc,ivL] ) // построение графиков
xgrid(2)
Ниже приводится функция для расчета топологических матриц. Входными
данными для функции является список начальных узлов ветвей ta
и список конечных узлов ветвей he. В результате вычисляются ортогональные
матрицы: контурная С и узловая А, кроме того формируется
информация о топологическом графе схемы g, позволяющая отобразить
граф схемы по введенному списку начальных и конечных узлов схемы.
81