Download - Wolfram Research

252 3. Reference Manual
Show equation system.
In[4]:= tran // DisplayForm
Out[4]//DisplayForm=
I$AC$D1t I$V0t ⩵⩵ 0, I$AC$D1t
V$2t C1 V$2 ′ t ⩵⩵ 0,
R1
V$1t ⩵⩵ V0, I$AC$D1t ⩵⩵ 1. 1 38.6635 V$1tV$2t IS$D1,
V$1t, V$2t, I$V0t, I$AC$D1t,
DesignPoint V0 2. Sin1000000 t, R1 100.,
C1 1. 10 7 , BV$D1 ∞, GLEAK$D1 0, GLEAKSW$D1 0,
IBV$D1 0.001, IS$D1 1. 10 14 , ISW$D1 0, PERIM$D1 0,
RS$D1 0, CD$D1 0, CJO$D1 0, CJSW$D1 0
Calculate operating point. In[5]:= op = First[NDAESolve[tran, {t, 0.}]]
Out[5]= V$2 0., V$1 0., I$V0 0., I$AC$D1 0.
Linearize equations about
the operating point.
Extract equation system.
Display design point
stored in the DAEObject.
Note that new parameters
have been generated.
In[6]:= ac = ACEquations[tran, op, {V0 −> 1}]
Out[6]= DAEAC, 4 4
In[7]:= GetEquations[ac]
Out[7]=
I$AC$D1 I$V0 ⩵⩵ 0, I$AC$D1
1 C1 s V$2 ⩵⩵ 0,
R1
V$1 ⩵⩵ V0, I$AC$D1 38.6635 38.6635 V$1$opV$2$op IS$D1 V$1
38.6635 38.6635 V$1$opV$2$op IS$D1 V$2 ⩵⩵ 0
In[8]:= GetDesignPoint[ac]
Out[8]=
R1 100., C1 1. 10 7 , BV$D1 ∞, GLEAK$D1 0,
GLEAKSW$D1 0, IBV$D1 0.001, IS$D1 1. 10 14 ,
ISW$D1 0, PERIM$D1 0, RS$D1 0, CD$D1 0, CJO$D1 0,
CJSW$D1 0, V0 1, V$1$op 0., V$2$op 0., I$V0$op 0.,
I$AC$D1$op 0., V0$op 2. Sin1000000 t$op
3.5.3 DCEquations
DCEquations[dae]
transforms a dynamic Transient DAEObject into a static
DC DAEObject
Command structure of DCEquations.
Given a Transient DAEObject, DCEquations transforms the dynamic equations into a static system
and returns a DC DAEObject. All derivatives with respect to the IndependentVariable (as stored
in the DAEObject) are replaced by zero. The initial time is inserted into both the equations and the
design point. Depending on the value of the ModeValues option, the values of independent sources
stored in the DesginPoint option of the DAEObject are replaced by the corresponding DC values
stored in the ModeValues option of the DAEObject.
Note that, given a Netlist or Circuit object, you can set up DC equations directly using
CircuitEquations (Section 3.5.1).