Our passivity correction algorithms are rigorous, and we guarantee our models to be passive. It's also worth pointing out that ensuring a system's passivity at discrete frequency points doesn't ensure that it will be passive when simulated in a variable time step simulator, since the simulation algorithms have to interpolate between those discrete frequency points and extrapolate beyond them. To be rigorous you have to ensure that the equations underlying the model are passive, and that the simulation algorithm will map a passive system to a stable simulation. -- Steve TDA Systems MeasureXtractor is EDN Innovativion of the Year finalist! http://www.reed-electronics.com/ednmag ------------------------------------------- Steven D. Corey, Ph.D. Time Domain Analysis Systems, Inc. "The Interconnect Analysis Company." http://www.tdasystems.com email: steven.corey@xxxxxxxxxxxxxx phone: (503) 246-2272 fax: (503) 246-2282 ------------------------------------------- Scott McMorrow wrote: > Michael and Ray, > First of all, I am not as smart as I sometimes sound. I had a good > friend of mine, Rob Hinz, develop the code for passivity checking and > passivity correction based upon some of my simple thoughts on the > subject . We use what we call the "Wack a Mole" method for correction > of non-passive s-parameters. You can imagine from that name what we > might do. > > You will want to compute the minimum real eigenvalue of Identity matrix > - Smatrix * Smatrix' on the complex data. > > In matlab this is done as follows > > mineigv = min(real(eig(eye(size(S))-S*S'))) > > > The matrix is passive the value is positive real > THe matrix is conditionally passive if the value is zero > The matrix is non-passive if the value is negative. > > > The original reference is: > > R. A. Rohrer, "Lumped network passivity criteria," IEEE Trans. Circuit > Theory, vol. > CT-15, no. 1, pp24-29, March 1968 > > > At this time, no one has developed a mathematiclly rigorous method for > passivity correction that I am aware of. > > > regards, > > scott > > > > > > > Mirmak, Michael wrote: > > >>To echo Raymond (forgive the questions from the matrix-impaired): >> >>Scott, >> >>If I have an 8-port S-parameter measurement (S11, S28, etc.), from a >>Touchstone or CITI file for example, would the given equation apply to >>the complex data or just to eh magnitudes? >> >>More specifically, taking the equation literally, my eigenvalues are >>calculated as: >> >>Identity matrix (8x8) - Smatrix * Smatrix' >> >>Where Smatrix is an 8x8 matrix of complex S-parameter values, Smatrix' >>is its transpose, and the identity matrix is just that, an 8x8 matrix >>filled with zeros on the off-diagonals and ones on the diagonals. Is it >>still correct to say that a complex result from that calculation, if >>having a negative real part, indicates non-passivity? >> >>More specifically still, I get an 8x8 matrix result of "eigenvalues" > >>from your equation. However, if I request eigenvalues from Matlab using > >>the same 8x8 Smatrix, I get a result which is 8x1 only. Why the >>difference? >> >>Thanks! >> >>- MM >> >> >> >>-----Original Message----- >>From: si-list-bounce@xxxxxxxxxxxxx [mailto:si-list-bounce@xxxxxxxxxxxxx] >>On Behalf Of Raymond Anderson >>Sent: Friday, February 20, 2004 12:26 PM >>To: si-list@xxxxxxxxxxxxx >>Subject: [SI-LIST] Re: Convergence problem after inserting detailed >>connector model >> >> >>Scott McMorrow wrote: >> >> >> >> >>>One way to check this is to run the model in the frequency domain and=20 >>>compute the eigenvalues of the S-parameters at each frequency point=20 >>>with I - S*S'. If there are any negative eigenvalues then the matrix=20 >>>is not passive. Dr Liu has had some experience with this, since his=20 >>>software will simulate with standard SPICE models and with=20 >>>s-parameters. He has found that not only can s-parameters be=20 >>>non-passive, but also that lumped element models with coupling can be=20 >>>non-passive. The culprit is usually mutual L. A simple test for this=20 >>>in a multi-section connector model is to comment out all K elements,=20 >>>re-run the simulation and see if it is stable. If the simulation is=20 >>>stable then it is most likely non-passive. More formally it should be=20 >>>run through an eigenvalue check, which Dr Liu reports to me that he has >>> >>> >> >> >> >> >>>done on a number of occasions. In that case, the model is bad and >>> >>> >> >>needs=20 >> >> >> >>>to go back to the developer. >>>=20 >>> >>> >>> >> >>Does anyone have any pointers to papers, books, or other info on how to >>compute the eigenvalues of the s-parameters as described above. i.e., >>are >>there any tools available or published matlab or whatever code available >>to >>automate the process? >> >>-Ray Anderson >>Sun Microsystems Inc. >> >>------------------------------------------------------------------ >>To unsubscribe from si-list: >>si-list-request@xxxxxxxxxxxxx with 'unsubscribe' in the Subject field >> >>or to administer your membership from a web page, go to: >>//www.freelists.org/webpage/si-list >> >>For help: >>si-list-request@xxxxxxxxxxxxx with 'help' in the Subject field >> >>List technical documents are available at: >> http://www.si-list.org >> >>List archives are viewable at: =20 >> //www.freelists.org/archives/si-list >>or at our remote archives: >> http://groups.yahoo.com/group/si-list/messages >>Old (prior to June 6, 2001) list archives are viewable at: >> http://www.qsl.net/wb6tpu >>=20 >> >>------------------------------------------------------------------ >>To unsubscribe from si-list: >>si-list-request@xxxxxxxxxxxxx with 'unsubscribe' in the Subject field >> >>or to administer your membership from a web page, go to: >>//www.freelists.org/webpage/si-list >> >>For help: >>si-list-request@xxxxxxxxxxxxx with 'help' in the Subject field >> >>List technical documents are available at: >> http://www.si-list.org >> >>List archives are viewable at: >> //www.freelists.org/archives/si-list >>or at our remote archives: >> http://groups.yahoo.com/group/si-list/messages >>Old (prior to June 6, 2001) list archives are viewable at: >> http://www.qsl.net/wb6tpu >> >> >> >> > > -- ------------------------------------------------------------------ To unsubscribe from si-list: si-list-request@xxxxxxxxxxxxx with 'unsubscribe' in the Subject field or to administer your membership from a web page, go to: //www.freelists.org/webpage/si-list For help: si-list-request@xxxxxxxxxxxxx with 'help' in the Subject field List technical documents are available at: http://www.si-list.org List archives are viewable at: //www.freelists.org/archives/si-list or at our remote archives: http://groups.yahoo.com/group/si-list/messages Old (prior to June 6, 2001) list archives are viewable at: http://www.qsl.net/wb6tpu