Any of the list members who are interested in passivity enforcement as it pertains to modeling may be interested in the Matlab code just released by Dr. Bjorn Gustavsen of SINTEF Energy Research in Trondheim, Norway. He writes: > Hi Raymond, > > You may want to announce on your list-server that a package for > passivity enforcement (QPpassive.zip) has now been added to the site > http://www.energy.sintef.no/produkt/VECTFIT/home.asp > > Best regards, > Bjorn This code utilizes quadratic programming techniques. From Bjorn's web page: ------------------------------------------------------------------------ This package contains a post-processing procedure for enforcing passivity of a model given on pole-residue form. This is done by introducing the passivity criterion as a constraint equation in the least squares problem. The routine requires the user to provide frequency samples for the least squares problem and the passivity constraint. A large number of options is available. Download: QPpassive.zip* * <http://www.energy.sintef.no/produkt/VECTFIT/QPpassive.zip> *File* *Purpose* QPpassive.pdf User manual passivity.pdf Document describing the theory (IEEE paper) QPpassive.m Passivity enforcement routine rot.m Utility routine, rotates eigenvectors to minimize their imaginary part interchange.m Utility routine, interchanges eigenpairs to obtain eigenvalues as smooth functions of frequency ex1.m, ex2.m, ex3a.m, ex3b.m, ex4a.m, ex4b.m Example cases described in QPpassive.pdf ex1.mat, ex2.mat, ex3.mat, ex4.mat Workspaces containing variables for examples (ex1.m,...) Note: QPpassive.m uses the routine quadprog.m of the Matlab Optimization Toolbox. ---------------------------------------------------------------------------------------------- -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 FAQ wiki page is located at: http://si-list.org/wiki/wiki.pl?Si-List_FAQ 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