23.01.2015 Views

Download - Wolfram Research

Download - Wolfram Research

Download - Wolfram Research

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

3.8 Pole/Zero Analysis 319<br />

Solve the GEP A B.<br />

Specify initial guesses for<br />

the eigenvectors.<br />

In[6]:= LREigenpair[A, B, −1., InitialGuess −> {{0, 1}, {1, 0}}]<br />

Out[6]=<br />

1.94319, 0.288369, 0.957519,<br />

0.957519, 0.288369, 5.05207 10 10 , 8.70217 10 10 , 4<br />

Use the default projection<br />

vectors.<br />

In[7]:= LREigenpair[A, B, 1., MaxResidual −> 1.*^−15]<br />

Out[7]=<br />

0.514618, 0.992822, 0.1196,<br />

0.1196, 0.992822, 1.24127 10 16 , 1.77722 10 16 , 4<br />

Note that the following setting for ProjectionVectors increases the number of iterations.<br />

Use right eigenvectors as<br />

projection vectors.<br />

In[8]:= LREigenpair[A, B, 1., MaxResidual −> 1.*^−15,<br />

ProjectionVectors −><br />

{RightEigenvector, RightEigenvector}]<br />

Out[8]=<br />

0.514618, 0.992822, 0.1196,<br />

0.1196, 0.992822, 1.24127 10 16 , 4.47545 10 16 , 5<br />

The setting Tolerance −> 1.0 allows the value of the solution Λ to lie within and times the<br />

value of the initial guess theta0.<br />

Search for an eigenvalue<br />

in the interval <br />

In[9]:= LREigenpair[A, B, −100., Tolerance −> 1.]<br />

LREigenpair::tolx:<br />

Eigenvalue lies outside the specified tolerance region.<br />

Out[9]=<br />

1.94319, 0.288369, 0.957519,<br />

0.957519, 0.288369, 2.70894 10 11 , 4.66522 10 11 , 6<br />

With the default setting Tolerance −> Infinity, any solution of the GEP is accepted without a<br />

warning.<br />

Do not specify a tolerance<br />

region.<br />

Read in a PSpice netlist<br />

and small-signal data.<br />

Set up a system of<br />

symbolic AC equations.<br />

In[10]:= LREigenpair[A, B, −100., Tolerance −> Infinity]<br />

Out[10]=<br />

1.94319, 0.288369, 0.957519,<br />

0.957519, 0.288369, 2.70894 10 11 , 4.66522 10 11 , 6<br />

In[11]:= buffer = ReadNetlist[<br />

"AnalogInsydes/DemoFiles/Buffer.cir",<br />

"AnalogInsydes/DemoFiles/Buffer.out",<br />

Simulator −> "PSpice"]<br />

Out[11]= Circuit <br />

In[12]:= mnabuffersym = CircuitEquations[buffer,<br />

AnalysisMode −> AC, ElementValues −> Symbolic]<br />

Out[12]= DAEAC, 18 18

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!