http://www.pspicelib.narod.ru

MatchSpice      www.pspicelib.narod.ru         www.pspice.comtv.ru
Трансляция уравнений и выражений в Mathematica

 В мире есть три конкурирующих математических САПР высокого уровня,  доступных  нашему пользователю.
Это Maple, Mathematica, MahCAD. Что лучше - решает пользователь.
Глупо не пользоваться.
MSpice - pspice совместимый пакет для аналитических расчётов электронных схем, работающий под Maple.
Глупо не пользоваться.
Такая система как Mathematica имеет беспрецендентную поддержку. В интернете можно найти кучу готовых и отлаженных примеров.
Глупо не пользоваться.
 Тогда туда нужно уметь передавать математические выражения и уравнения  из MSpice
Такя функция появилась в MSpice.
 Глупо не пользоваться.
 

Рис. 3. Схема усилителя
 [Maple Bitmap]

>    restart:
#read(`F:/PMAPLE/ESolver.m`);with(MSpice):
read(`F:/MSPICE/MSpice.m`);with(MSP):read(`F:/MSPICE/ESolver.m`);with(SolveNet):
Приборы:=Oдинаковые.AC.en: ESolve(OE,`02 V-усилитель на 1-м транзисторе.CKT`);

`Cистема Кирхгофа-Laplace solve { ODE , ISc }`

-(V3-V4)*beta/Rbe+I1-(V2-V1)/`Rн`

-V4/R3+I2+(V3-V4)*beta/Rbe-(V4-V3)/Rbe-(V4-V3)*s*Cbe

-V3/R2+(V4-V3)/Rbe+(V4-V3)*s*Cbe+IIN-(V3-V1)/R1

{V3, V2, V4}

`` 

MSpice v7.72:   www.pspicelib.narod.ru

Заданы узлы: {V1} Источники: [I1, I2, IIN, VB1]

Решения V_NET: solve(ODE,ISc):  assign(%);

J_NET: [Je, JT, JVB1, JR1, JR2, Jk, JR3, JRн, JCbe, JRbe, J1, J2, JIN]

``

При задании опции OE в диррективе ESolve выводятся уравнения цепи по Кирхгофу, но не решаются.
Такой режим нужен, когда вы используете какие то нетрадиционные способы поиска решения. Например хотите решить эту задачу в Mathematica.
Их можно посмотреть.

>    ODE:=ODE; ISc:=ISc;
 

ODE := {-V3/R2+(V4-V3)/Rbe+(V4-V3)*s*Cbe+IIN-(V3-VB1)/R1, -V4/R3+I2+(V3-V4)*beta/Rbe-(V4-V3)/Rbe-(V4-V3)*s*Cbe, -(V3-V4)*beta/Rbe+I1-(V2-VB1)/`Rн`}

ISc := {V3, V2, V4}

Транслируем в формат Mathematica. Затем через буфер копируем в окно Mthematica
И уже там получаем решение

>    convMma(ODE,ISc);

`Уравнениe в формате Mathematica`

`` 

"Solve[{-V3/R2+(V4-V3)/Rbe+(V4-V3)*s*Cbe+IIN-(V3-VB1)/R1 == 0, -V4/R3+I2+(V3-V4)*beta/Rbe-(V4-V3)/Rbe-(V4-V3)*s*Cbe == 0, -(V3-V4)*beta/Rbe+I1-(V2-VB1)/Rн == 0}, {V3, V2, V4}]"

Решение уравнений в Mathematica

[Maple Bitmap]

>   

MSpiceo    www.pspicelib.narod.ru
Преобразование выражений и уравнений в формат Mathematica

В мире есть три конкурирующих математических САПР высокого уровня,  доступных  нашему пользователю.
Это Maple, Mathematica, MahCAD. Что лучше - решает пользователь.
Глупо не пользоваться.
MSpice - pspice совместимый пакет для аналитических расчётов электронных схем, работающий под Maple.
Глупо не пользоваться.

Такая система как Mathematica имеет беспрецендентную поддержку. В интернете можно найти кучу готовых и отлаженных примеров.
Глупо не пользоваться.
Тогда туда нужно уметь передавать математические выражения и уравнения  из MSpice
Такя функция появилась в MSpice.
Глупо не пользоваться.  
 


[Maple Metafile]

>    restart:with(inttrans):
#read(`F:/PMAPLE/ESolver.m`);with(MSpice):
read(`F:/MSPICE/MSpice.m`);with(MSP):read(`F:/MSPICE/ESolver.m`);with(SolveNet):
ESolve(NQ,`M2/OP-PSpiceFiles/SCHEMATIC1/SCHEMATIC1.net`);

`` 

MSpice v7.72:   www.pspicelib.narod.ru

Заданы узлы: {VarDC} Источники: [нет источников]

Решения V_NET: [VOUT2, VOUT, VOUT3, VN682972, VN6830131, VN683176, V6, V2, V5, V1, V3, V4]

J_NET: [JSciense]

``

>    VOUT3:=VOUT3;  

VOUT3 := 2*sin(exp(t)+cos(t))*ln(exp(t)^2+exp(t)*cos(t))/ln(10)*cos(exp(t)+cos(t))*(exp(t)-sin(t))+sin(exp(t)+cos(t))^2*(2*exp(t)^2+exp(t)*cos(t)-exp(t)*sin(t))/(exp(t)^2+exp(t)*cos(t))/ln(10)

Преобразование выражения в формат Mathematica.

>    convHma(VOUT3);

`Выражение в формате Mathematica`

`` 

"2*Sin[Exp[t]+Cos[t]]*Ln[Exp[t]^2+Exp[t]*Cos[t]]/Ln[10]*Cos[Exp[t]+Cos[t]]*(Exp[t]-Sin[t])+Sin[Exp[t]+Cos[t]]^2*(2*Exp[t]^2+Exp[t]*Cos[t]-Exp[t]*Sin[t])/(Exp[t]^2+Exp[t]*Cos[t])/Ln[10]"

2*sin(exp(t)+cos(t))*Ln(exp(t)^2+exp(t)*cos(t))/Ln(10)*cos(exp(t)+cos(t))*(exp(t)-sin(t))+sin(exp(t)+cos(t))^2*(2*exp(t)^2+exp(t)*cos(t)-exp(t)*sin(t))/(exp(t)^2+exp(t)*cos(t))/Ln(10)

>   

Передача формулы  в Mathematica

[Maple Bitmap]

Hosted by uCoz