摘要

Since the conditional nonlinear optimal perturbations (CNOPs) are formulated mathematically as a constrained minimization problem, the existing numerical methods for computing the CNOPs were designed to solve the original constrained optimization problems by directly adapting some (constrained) optimization algorithms. Although such an approach is natural and convenient, it often results in inefficient and expensive numerical methods for computing the CNOPs. In this article we propose an alternative approach which is based on the idea of first transforming the original CNOP problem into a special unconstrained optimization problem (known as a nonlinear least-squares problem) using a penalty strategy and then utilizing a more efficient Gauss-Newton iterative method, to solve the transformed (unconstrained) nonlinear least-squares problem. Compared with the existing numerical methods, the proposed Gauss-Newton iterative method has some desired advantages: it is easy to implement and portable because it does not rely on any existing constrained optimization package; it is also fast, and highly efficient. These advantages and potential merits of the proposed method are demonstrated by several sets of evaluation experiments based on the viscous Burgers equation, the T21L3 quasi-geostrophic model and the shallow-water equation model, respectively.