文摘
In this paper, we investigate interpolatory model order reduction for large-scale bilinear descriptor systems. Recently, it was shown in Wyatt et al. (2013) for linear descriptor systems that directly extending the standard rational interpolation conditions used in d="mmlsi3" class="mathmlsrc">class="formulatext stixSupport mathImg" data-mathURL="/science?_ob=MathURL&_method=retrieve&_eid=1-s2.0-S0167691116301141&_mathId=si3.gif&_user=111111111&_pii=S0167691116301141&_rdoc=1&_issn=01676911&md5=bd020f5421a488634c8ea377c98b15b4" title="Click to view the MathML source">H2class="mathContainer hidden">class="mathCode"> optimal model reduction to descriptor systems, in general, yields an unbounded error in the d="mmlsi3" class="mathmlsrc">class="formulatext stixSupport mathImg" data-mathURL="/science?_ob=MathURL&_method=retrieve&_eid=1-s2.0-S0167691116301141&_mathId=si3.gif&_user=111111111&_pii=S0167691116301141&_rdoc=1&_issn=01676911&md5=bd020f5421a488634c8ea377c98b15b4" title="Click to view the MathML source">H2class="mathContainer hidden">class="mathCode">-norm. This is due to the possible mismatch of the polynomial part of the original and reduced-order systems. This conclusion holds for nonlinear systems as well. In this paper, we deal with bilinear descriptor systems and aim to pay attention to the polynomial part of the bilinear descriptor system along with interpolation. To this end, we have shown in Goyal et al. (2015) how to determine the polynomial part of each subsystem of the bilinear descriptor system explicitly, by assuming special structures of the system matrices. Considering the same structured bilinear descriptor systems, in this paper we first show how to achieve multipoint interpolation of the underlying Volterra series of bilinear descriptor systems while retaining the polynomial part of each subsystem of the bilinear system. Then, we extend the interpolation based first-order necessary conditions for d="mmlsi3" class="mathmlsrc">class="formulatext stixSupport mathImg" data-mathURL="/science?_ob=MathURL&_method=retrieve&_eid=1-s2.0-S0167691116301141&_mathId=si3.gif&_user=111111111&_pii=S0167691116301141&_rdoc=1&_issn=01676911&md5=bd020f5421a488634c8ea377c98b15b4" title="Click to view the MathML source">H2class="mathContainer hidden">class="mathCode"> optimality to bilinear descriptor systems and propose an iterative scheme to obtain an d="mmlsi3" class="mathmlsrc">class="formulatext stixSupport mathImg" data-mathURL="/science?_ob=MathURL&_method=retrieve&_eid=1-s2.0-S0167691116301141&_mathId=si3.gif&_user=111111111&_pii=S0167691116301141&_rdoc=1&_issn=01676911&md5=bd020f5421a488634c8ea377c98b15b4" title="Click to view the MathML source">H2class="mathContainer hidden">class="mathCode"> optimal reduced-order system. By means of various numerical examples, we demonstrate the efficiency of the proposed model order reduction technique and compare it with reduced bilinear systems obtained by using linear IRKA, the Loewner method for bilinear systems and POD-based approximations.