Refine
Year of publication
Language
- English (1103)
- German (10)
- Multiple languages (1)
Keywords
- optimal control (27)
- stability (14)
- integer programming (11)
- Stochastic programming (9)
- finite elements (9)
- mixed integer programming (9)
- Hamiltonian matrix (8)
- finite element method (8)
- model reduction (8)
- state constraints (8)
Dissipative Hamiltonian (DH) systems are an important concept in energy based modeling of dynamical
systems. One of the major advantages of the DH formulation is that the system encodes system
properties in an algebraic way in the system. Making use of the structure,
it is easy to see that DH systems are stable. In this paper
the question is discussed when a linear constant coefficient DH system is on the boundary
of the region of asymptotic stability, i.e., when it has purely imaginary eigenvalues,
or how much it has to be perturbed to be on this boundary. For unstructured systems this distance to instability (stability radius) is well-understood. In this paper,
explicit formulas for this distance under structure preserving perturbations are determined.
It is also shown (via numerical examples) that under structured perturbations the asymptotical
stability of a DH system is much more robust than for unstructured perturbations, since the
distance can be much larger.
Canonical forms for matrix triples $(A,G,\hat G)$, where
$A$ is arbitrary rectangular and $G$, $\hat G$ are either real symmetric
or skew symmetric, or complex Hermitian or skew Hermitian, are derived.
These forms generalize classical product Schur forms as well as
singular value decompositions.
An new proof for the complex case is given, where there is no need to
distinguish whether $G$ and $\hat G$ are Hermitian or skew Hermitian.
This proof is independent from the results in Bolschakov/Reichstein 1995, where
a similar canonical form has been obtained for the complex case,
and it allows generalization to the real case. Here,
the three cases, i.e., that
$G$ and $\hat G$ are both symmetric, both skew symmetric or one each,
are treated separately.
The long standing problem is discussed of how to deflate the part associated with the eigenvalue infinity in a structured matrix pencil using structure preserving unitary transformations. We derive such a deflation procedure and apply this new technique to symmetric, Hermitian or alternating pencils and in a modified form to (anti)-palindromic pencils. We present a detailed error and perturbation analysis of this and other deflation procedures and demonstrate the properties of the new algorithm with several numerical examples.
The paper provides a structural analysis of the feasible set defined by linear probabilistic constraints. Emphasis is laid on single (individual) probabilistic constraints. A classical convexity result by Van de Panne/Popp and Kataoka is extended to a broader class of distributions and to more general functions of the decision vector. The range of probability levels for which convexity can be expected is exactly identified. Apart from convexity, also nontriviality and compactness of the
feasible set are precisely characterized at the same time. The relation between feasible sets with negative and with nonnegative right-hand side is revealed. Finally, an existence result is formulated for the more difficult case of joint probabilistic constraints.
Actin is a major structural protein of the eukaryotic cytoskeleton and enables cell motility.
Here, we present a model of the actin filament (F-actin) that incorporates the global structure
of the recently published model by Oda et al. but also conserves internal stereochemistry. A
comparison is made using molecular dynamics simulation of the model with other recent F-
actin models. A number of structural determents such as the protomer propeller angle, the
number of hydrogen bonds and the structural variation among the protomers are analyzed.
The MD comparison is found to reflect the evolution in quality of actin models over the last
six years. In addition, simulations of the model are carried out in states with both ADP or
ATP bound and local hydrogen-bonding differences characterized. The results point to the
significance of a direct interaction of Gln137 with ATP for activation of ATPase activity after
the G-to-F-actin transition.
Diffusion Weighted Imaging has become and will certainly continue to be an important tool in medical research and diagnostics. Data obtained with Diffusion Weighted Imaging are characterized by a high noise level. Thus, estimation of quantities like anisotropy indices or the main diffusion direction may be significantly compromised by noise in clinical or neuroscience applications. Here, we present a new package dti for R, which provides functions for the analysis of diffusion weighted data within the diffusion tensor model. This includes smoothing by a recently proposed structural adaptive smoothing procedure based on the Propagation-Separation approach in the context of the widely used Diffusion Tensor Model. We extend the procedure and show, how a correction for Rician bias can be incorporated. We use a heteroscedastic nonlinear regression model to estimate the diffusion tensor. The smoothing procedure naturally adapts to different structures of different size and thus avoids oversmoothing edges and fine structures. We illustrate the usage and capabilities of the package through some examples.
Functional Magnetic Resonance Imaging inherently involves noisy measurements and a severe multiple test
problem. Smoothing is usually used to reduce the effective number of multiple
comparisons and to locally integrate the signal and hence increase the
signal-to-noise ratio. Here, we provide a new structural adaptive segmentation
algorithm (AS)
that naturally combines the signal detection with noise reduction in one procedure.
Moreover, the new method
is closely related to a recently proposed structural adaptive smoothing
algorithm and preserves shape and spatial extent of activation areas without
blurring the borders.
In this paper, we consider the characterization of strong stationary solutions to
equilibrium problems with equilibrium constraints (EPECs). Assuming that the underlying
generalized equation satisfies strong regularity in the sense of Robinson, an explicit
multiplier-based stationarity condition can be derived. This is applied then
to an equilibrium model arising from ISO-regulated electricity spot markets.
Research on flows over time has been conducted mainly in two separate and mainly independent approaches, namely \emph{discrete} and \emph{continuous} models, depending on whether a discrete or continuous representation of time is used. Recently, Borel flows have been introduced to build a bridge between these two models.
In this paper, we consider the maximum Borel flow problem formulated in a network where capacities on arcs are given as Borel measures and storage might be allowed at the nodes of the network. This problem is formulated as a linear program in a space of measures. We define a dual problem and prove a strong duality result. We show that strong duality is closely related to a MaxFlow-MinCut Theorem.
We present a novel algorithm for automatic parameterization of tube-like surfaces of arbitrary genus such as the surfaces of knots, trees, blood vessels, neurons, or any tubular graph with a globally consistent stripe texture. We use the principal curvature frame field of the underlying tube-like surface to guide the creation of a global, topologically consistent stripe parameterization of the surface. Our algorithm extends the QuadCover algorithm and is based, first, on the use of so-called projective vector fields instead of frame fields, and second, on different types of branch points. That does not only simplify the mathematical theory, but also reduces computation time by the decomposition of the underlying stiffness matrices.
In the first part of this article, we have shown how time-dependent optimal control for partial
differential equations can be realized in a modern high-level modeling and simulation package. In this second part we extend our approach to (state) constrained problems. "Pure" state constraints in a function space
setting lead to non-regular Lagrange multipliers (if they exist), i.e. the Lagrange multipliers are in general Borel
measures. This will be overcome by different regularization techniques.
To implement inequality constraints, active set methods and interior point methods (or barrier methods) are widely in use. We show how these techniques can be realized in the modeling and simulation package Comsol
Multiphysics.
In contrast to the first part, only the one-shot-approach based on space-time elements is considered. We implemented a projection method based on active sets as well as a barrier method and compare these methods
by a specialized PDE optimization program, and a program that optimizes the discrete version of the given problem.
We show how time-dependent optimal control for partial differential equations can be realized in a modern high-level modeling and simulation package. We summarize the general formulation for distributed and boundary control for initial-boundary value problems for parabolic PDEs and derive the optimality system including the adjoint equation. The main difficulty therein is that the latter has to be integrated backwards in time. This implies that complicated implementation effort is necessary to couple state and adjoint equations to compute an optimal solution. Furthermore a large amount of computational effort or storage is required to provide the needed information (i.e the trajectories) of the state and adjoint variables. We show how this can be realized in the modeling and simulation package COMSOL MULTIPHYSICS, taking advantage of built-in discretization, solver and post-processing technologies and thus minimizing the implementation effort. We present two strategies: The treatment of the coupled optimality system in the space-time cylinder, and the iterative approach by sequentially solving state and adjoint system and updating the controls. Numerical examples show the elegance of the implementation and the efficiency of the two strategies.
Boolean modeling frameworks have long since proved their worth for capturing and analyzing essential characteristics of complex systems.
Hybrid approaches aim at exploiting the advantages of Boolean formalisms while refining expressiveness. In this paper, we present a formalism that augments Boolean models with stochastic aspects. More specifically, biological reactions effecting a system in a given state are associated
with probabilities, resulting in dynamical behavior represented as a Markov chain. Using this approach, we model and analyze the cytokinin
response network of Arabidopsis thaliana with a focus on clarifying the character of an important feedback mechanism.
We present a mixed-integer multistage stochastic programming model for the short term unit commitment of a hydro-thermal power system under uncertainty in load, inflow to reservoirs, and prices for fuel and delivery contracts. The model is implemented for uncertain load and tested on realistic data from a German power utility. Load scenario trees are generated by a procedure consisting of two steps: (i) Simulation of load scenarios using an explicit respresentation of the load distribution and (ii) construction of a tree out of these scenarios. The dimension of the corresponding mixed-integer programs ranges up to 200,000 binary and 350,000 continuous variables. The model is solved by a Lagrangian-based decomposition strategy exploiting the loose coupling structure. Solving the Lagrangian dual by a proximal bundle method leads to a successive decomposition into single unit subproblems, which are solved by specific algorithms. Finally, Lagrangian heuristics are used to construct nearly optimal first stage decisions.
We consider the preemptive and non-preemptive problems of scheduling jobs with precedence constraints on parallel machines with the
objective to minimize the sum of~(weighted) completion times. We investigate an online model in which the scheduler learns about a
job when all its predecessors have completed. For scheduling on a single machine, we show matching lower and upper bounds of~$\Theta(n)$ and~$\Theta(\sqrt{n})$ for jobs with general and equal weights, respectively. We also derive corresponding results on parallel machines.
Our result for arbitrary job weights holds even in the more general stochastic online scheduling model where, in addition to the limited information about the job set, processing times are uncertain. For a
large class of processing time distributions, we derive also an improved performance guarantee if weights are equal.
We consider a non-preemptive, stochastic parallel machine
scheduling model with the goal to minimize the weighted completion
times of jobs. In contrast to the classical stochastic model where jobs
with their processing time distributions are known beforehand, we assume
that jobs appear one by one, and every job must be assigned
to a machine online. We propose a simple online scheduling policy for
that model, and prove a performance guarantee that matches the currently
best known performance guarantee for stochastic parallel machine
scheduling. For the more general model with job release dates we derive
an analogous result, and for NBUE distributed processing times we
even improve upon the previously best known performance guarantee for
stochastic parallel machine scheduling. Moreover, we derive some lower
bounds on approximation.
We consider empirical approximations of two-stage stochastic mixed-integer linear programs and derive central
limit theorems for the objectives and optimal values. The limit theorems are based on empirical process theory
and the functional delta method. We also show how these limit theorems can be used to derive confidence intervals
for optimal values via a certain modification of the bootstrapping method.
We analyze an interactive model of credit ratings where external shocks, initially
affecting only a small number of firms, spread by a contagious chain reaction to the
entire economy. Counterparty relationships along with discrete adjustments of credit
ratings generate a transition mechanism that allows the financial distress of one firm
to spill over to its business partners. Such a contagious infectious of financial distress
constitutes a source of intrinsic risk for large portfolios of credit sensitive securities that
cannot be “diversified away.” We provide a complete characterization of the fluctuations
of credit ratings in large economies when adjustments follow a threshold rule. We also
analyze the effects of downgrading cascades on aggregate losses of credit portfolios. We
show that the loss distribution has a power-law tail if the interaction between different
companies is strong enough.
Stochastic Optimization of Electricity Portfolios: Scenario Tree Modeling and Risk Management
(2008)
We present recent developments in the field of stochastic programming with regard to application in power management. In particular we discuss issues of scenario tree modeling, i.e., appropriate discrete approximations of the underlying stochastic parameters. Moreover, we suggest risk avoidance strategies via the incorporation of
so-called polyhedral risk functionals into stochastic programs. This approach, motivated through tractability of the resulting problems, is a constructive framework providing particular flexibility with respect to the dynamic aspects of risk.
A strategy for controlling the stepsize in the numerical integration of stochastic
differential equations (SDEs) is presented. It is based on estimating the p-th mean of
local errors. The strategy leads to deterministic stepsize sequences that are identical
for all paths. For the family of Euler schemes for SDEs with small noise we derive
computable estimates for the dominating term of the p-th mean of local errors
and show that the strategy becomes efficient for reasonable stepsizes. Numerical
experience is reported for test examples including scalar SDEs and a stochastic
circuit model.
The package fmri is provided for analysis of single run functional
Magnetic Resonance Imaging data. It implements structural adaptive smoothing
methods with signal detection for adaptive noise reduction which avoids blurring
of edges of activation areas. fmri provides fmri analysis from time series
modeling to signal detection and publication-ready images.
New types of stationary solutions of a one-dimensional driven sixth-order Cahn-Hilliard type equation that arises as a model for epitaxially growing nano-structures such as quantum dots, are derived by an extension of the method of matched asymptotic expansions that retains exponentially small terms. This method yields analytical expressions for far-field behavior as well as the widths of the humps of these spatially non-monotone solutions in the limit of small driving force strength which is the deposition rate in case of epitaxial growth. These solutions extend the family of the
monotone kink and antikink solutions. The hump spacing is related to solutions of the Lambert $W$ function.
Using phase space analysis for the corresponding fifth-order dynamical
system, we use a numerical technique that enables the efficient and accurate tracking of the solution branches, where the asymptotic solutions are used as initial input.
Additionally, our approach is first demonstrated for the related but simpler driven fourth-order Cahn-Hilliard equation, also known as the convective Cahn-Hilliard equation.
We give sufficient conditions for a non-zero sum discounted stochastic game with
compact and convex action spaces and with norm-continuous transition probabilities,
but with possibly unbounded state space, to have a Nash equilibrium in homogeneous
Markov strategies that depends in a Lipschitz continuous manner on the current state. If
the underlying state space is compact this yields the existence of a stationary equilibrium.
Stochastic games with weakly interacting players provide a probabilistic framework within
which to study strategic behavior in models of non-market interactions.
We discuss a parallel library of efficient algorithms for model reduction of largescale
systems with state-space dimension up to O(104). We survey the numerical
algorithms underlying the implementation of the chosen model reduction methods.
The approach considered here is based on state-space truncation of the system
matrices and includes absolute and relative error methods for both stable and unstable
systems. In contrast to serial implementations of these methods, we employ
Newton-type iterative algorithms for the solution of the major computational tasks.
Experimental results report the numerical accuracy and the parallel performance of
our approach on a cluster of Intel Pentium II processors.
A state-constrained optimal control problem arising in the context of sublimation crystal growth is considered. The presence of pointwise state-constraints and nonlocal radiation interface conditions
constitutes the major issue of this problem. A regularity result of the state is presented that allows to
derive the optimality condition.
We consider a control- and state-constrained optimal control problem
governed by a semilinear
elliptic equation with nonlocal interface conditions.
These conditions occur during the
modeling of diffuse-gray conductive-radiative heat transfer.
The nonlocal radiation interface condition and the pointwise state-constraints
represent the particular features of this problem. To deal with the
state-constraints, continuity of the state is shown which allows to
derive first-order necessary conditions. Afterwards, we establish second-order
sufficient conditions that account for strongly active sets and
ensure local optimality in an $L^2$-neighborhood.
In optimal control problems with nonlinear time-dependent 3D PDEs, full 4D discretizations are usually prohibitive due to the storage requirement. For this reason gradient and Newton type methods working on the reduced functional are often employed. The computation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. The state enters into the adjoint equation, again requiring the storage of a full 4D data set. We propose a lossy compression algorithm using an inexact but cheap predictor for the state data, with additional entropy coding of prediction errors. As the data is used inside a discretized, iterative algorithm, lossy compression
maintaining a certain error bound turns out to be sufficient.
We consider first order optimality conditions for state constrained optimal control problems. In particular we study the case where the state equation has not enough regularity to admit existence of a Slater point in function space. We overcome this difficulty by a special transformation. Under a density condition we show existence of Lagrange multipliers, which have a representation via measures and additional regularity properties.
We discuss the eigenvalue problem for
general and structured matrix polynomials which may
be singular and may have eigenvalues at infinity.
We derive staircase
condensed forms that allow deflation of the infinite eigenvalue and
singular structure of the matrix polynomial.
The remaining reduced order staircase form leads to
new types of linearizations which determine the finite eigenvalues and
and corresponding eigenvectors. The new linearizations
also simplify the construction of structure preserving linearizations.
Whenever the invariant stationary density of metastable dynamical systems decomposes into almost invariant partial densities, its computation as eigenvector of some transition probability matrix is an ill-conditioned problem. In order to avoid this computational difficulty, we suggest to apply an aggregation/disaggregation method which only addresses wellconditioned sub-problems and thus results in a stable algorithm. In contrast to existing methods, the aggregation step is done via a sampling algorithm which covers only small patches of the sampling space. Finally, the theoretical analysis is illustrated by two biomolecular examples.
The aim of this paper is to study the behaviour of a weak solution to Navier-Stokes equations for isothermal fluids with a nonlinear stress tensor for time going to infinity. In an analogous way as in [18], we construct a suitable function which approximates the density for time going to infinity. Using properties of this function, we can prove the strong convergence of the density to its limit state. The behaviour of the velocity field and kinetic energy is mentioned as well.
This work studies the stability and the stochastic properties of neural activity evoked by external
stimulation. The underlying model describes the spatiotemporal dynamics of neural populations
involving both synaptic delay and axonal transmission delay. We show, that the linear model
recasts to a set of affne delay differential equations in spatial Fourier space. Besides a stability
study for general kernels and general external stimulation, the power spectrum of evoked activity
is derived analytically in case of external Gaussian noise. Further applications to specific kernels
reveal critical
uctuations at Hopf- and Turing bifurcations and allow the numerical detection of
1/f fluctuations near the stability threshold.
We study linear dissipative Hamiltonian (DH) systems with real constant coefficients that arise in energy based modeling of dynamical
systems. In this paper we analyze when such a system is on the boundary
of the region of asymptotic stability, i.e., when it has purely imaginary eigenvalues,
or how much the dissipation term has to be perturbed to be on this boundary. For unstructured systems the explicit construction of the \emph{real distance to instability (real stability radius)} has been a challenging problem. In this paper, we analyze this real distance under different structured perturbations to the dissipation term that preserve the DH structure and we derive explicit formulas for this distance in terms of low rank perturbations. We also show (via numerical examples) that under real structured perturbations to the dissipation the asymptotical
stability of a DH system is much more robust than for unstructured perturbations.
Classical stability properties of solutions
that are well-known for ordinary differential
equations (ODEs) are generalized to differential-algebraic equations (DAEs).
A new test equation is derived for the analysis of numerical methods applied
to DAEs with respect to the stability of the numerical approximations.
Morevover, a stabilization technique is developed to improve the stability of classical DAE integration methods. The stability regions for these stabilized discretization methods are determined and it is shown that they much better reproduce the stability properties known for the ODE case
than in the unstabilized form.
Movies that depict the stability regions for several methods are included for interactive use.
We analyse stability aspects of linear multistage stochastic programs with polyhedral risk measures in the objective. In particular, we consider sensitivity of the optimal value with respect perturbations of the underlying stochastic input process. An existing stability result for multistage stochastic programs with expectation objective is carried forward to the case of polyhedral risk-averse objectives. Beside Lr-distances these results also involve filtration distances of the perturbations of the stochastic process. We discuss additional requirements for the
polyhedral risk measures such that the problem dependent filtration distances can be bounded by problem independent ones. Stability and such bounds are the basis for scenario tree approximation techniques used in practical problem solving.
Quantitative stability of linear multistage stochastic programs is studied. It
is shown that the infima of such programs behave (locally) Lipschitz continuous
with respect to the sum of an Lr-distance and of a distance measure for the filtrations
of the original and approximate stochastic (input) processes. Various issues
of the result are discussed and an illustrative example is given. Consequences for
the reduction of scenario trees are also discussed.
Stability of Linear Stochastic Difference Equations in Strategically Controlled Random Environments
(2004)
We consider the stochastic sequence fYtgt2N defined recursively by the linear relation
Yt+1 = AtYt+Bt in a random environment. The environment is described by the stochastic
process f(At;Bt)gt2N and is under the simultaneous control of several agents playing a
discounted stochastic game. We formulate sufficient conditions on the game which ensure
the existence of Nash equilibrium in Markov strategies which has the additional property
that, in equilibrium, the process fYtgt2N converges in distribution to a stationary regime.
The Lang-Kobayashi model is a system of delay differential equations (DDEs) describing the dynamics of a semiconductor laser under delayed optical feedback. In this paper, we study the stability of so called external cavity modes (ECMs), which are harmonic oscillations corresponding to stationary lasing states. We focus on experimentally relevant situations, when the delay is large compared to the internal time scales of the laser. In this case, both the number of ECMs and the number of critical eigenvalues grows to infinity. Applying a newly developed asymptotic description for the spectrum of linearized DDEs with long delay, we are able to overcome this difficulty and to give a complete description of the stability properties of all ECMs. In particular, we distinguish between different types of weak and strong instabilities and calculate bifurcation diagrams that indicate the regions with different stability properties and the transitions between them.
An analysis of convex stochastic programs is provided if the underlying probability distribution is subjected to (small) perturbations. It is shown, in particular, that epsilon-approximate solution sets of convex stochastic programs behave Lipschitz continuous with respect to certain distances of probability distributions that are generated by the relevant integrands. It is shown that these results apply to linear two-stage stochastic programs with random recourse. Consequences are discussed on associating Fortet-Mourier metrics to two-stage models and on the asymptotic behavior of empirical estimates of such models, respectively.
We consider convex optimization problems with $k$th order stochastic dominance constraints for $k\ge 2$. We discuss distances of random variables that are relevant for the dominance relation and establish quantitative stability results for optimal values and solution sets in terms of a suitably selected probability metrics.Moreover, we provide conditions ensuring that the optimal value function is Hadamard directionally differentiable. Finally, we discuss some implications of the results for empirical (Monte Carlo,
sample average) approximations of dominance constrained optimization models.
Stability and Sensitivity of Optimization Problems with First Order Stochastic Dominance Constraints
(2007)
We analyze the stability and sensitivity of stochastic optimization problems with stochastic dominance constraints of first order. We consider general perturbations of the underlying probability measures in the space of regular measures equipped with a suitable discrepancy distance. We show that the graph of the feasible set mapping is closed under rather general assumptions. We obtain conditions for the continuity of the optimal value and upper-semicontinuity of the optimal solutions, as well as quantitative stability estimates of Lipschitz type.
Furthermore, we analyze the sensitivity of the optimal value and obtain upper and lower bounds for the directional
derivatives of the optimal value. The estimates are formulated in terms of the dual utility functions associated with the
dominance constraints.
By extending the stability analysis of [17] for multistage stochastic programs we show that their solution sets behave stable with respect to the sum of an Lr-distance and a filtration distance. Based on such stability results we suggest a scenario tree generation method for the (multivariate) stochastic input process. It starts with a fan of individual scenarios and consists of a recursive deletion and branching procedure which is controlled by bounding the approximation error. Some numerical experience for generating scenario trees in electricity portfolio management is reported.
A stability analysis is presented for neural field equations in the presence
of axonal delays and for a general class of connectivity kernels and synap-
tic properties. Sufficient conditions are given for the stability of equilibrium
solutions. It is shown that the delays play a crucial role in non-stationary
bifurcations of equilibria, whereas the stationary bifurcations depend only on
the kernel. Bounds are determined for the frequencies of bifurcating periodic
solutions. A perturbative scheme is used to calculate the types of bifurca-
tions leading to spatial patterns, oscillatory solutions, and traveling waves.
For high transmission speeds a simple method is derived that allows the de-
termination of the bifurcation type by visual inspection of the Fourier trans-
forms of the connectivity kernel and its first moment. Results are numerically
illustrated on a class of neurologically plausible second order systems with
combinations of Gaussian excitatory and inhibitory connections.
In this paper, we discuss stability properties of positive descriptor systems in the continuous-time as well as in the discrete-time case. We present different characterisations of positivity and establish generalised stability criteria for the case of positive descriptor systems. We show that if the spectral projector onto the right finite deflating subspace of the matrix pair $(E,A)$ is non-negative, then all stability criteria for standard positive systems take a comparably simple form in the positive descriptor case. Furthermore, we provide sufficient conditions that guarantee entry-wise non-negativity along with positive semi-definiteness of solutions of generalised projected Lyapunov equations. As an application of the framework established throughout this paper, we exemplarily generalise two criteria for the stability of two switched standard positive systems under arbitrary switching to the descriptor case.
In this paper we discuss the stability and model order reduction of coupled linear
time-invariant systems. Sufficient conditions for a closed-loop system to be asymptotically stable are
given. We present a model reduction approach for coupled systems based on reducing the order of the
subsystems and coupling the reduced-order subsystems by the same interconnection matrices as for
the original model. Such an approach allows to obtain error bounds for the reduced-order closed-loop
system in terms of the errors in the reduced-order subsystems. Model reduction of coupled systems
with unstable subsystems is also considered. Numerical examples are given.
The paper introduces an identification problem arising in modern regional hyperthermia, a cancer
therapy aiming at heating the tumor by microwave radiation. The task is to identify the highly
individual perfusion, which affects the resulting temperature distribution, from MR measurements.
The identification problem is formulated as an optimization problem. Existence of a solution and
optimality conditions are analyzed. Different regularizations and problem variants are considered. For
the numerical solution, a standard SQP method is used. Sufficient conditions for the convergence of
the method are derived. Finally, numerical examples on artificial as well as clinical data are presented.
This work studies dynamical properties of spatially extended neu-
ronal ensembles. We first derive an evolution equation from tem-
poral properties and statistical distributions of synapses and somata.
The obtained integro-differential equation considers both synaptic and
axonal propagation delay, while spatial synaptic connectivities ex-
hibit gamma-distributed distributions. This familiy of connectivity
kernels also covers the cases of divergent, finite, and negligible self-
connections. The work derives conditions for both stationary and
nonstationary instabilities for gamma-distributed kernels.It turns out
that the stability conditions can be formulated in terms of the mean spatial interaction ranges and the mean spatial interaction times. In
addition, a numerical study examines the evoked spatiotemporal re-
sponse activity caused by short local stimuli and reveals maximum
response activity after the mean interaction time at a distance from
stimulus offset location equal to the mean interaction range. These
findings propose new insights to neuronal mechanisms of experimen-
tally observed evoked brain activity.
Chimera states are particular trajectories
in systems of phase oscillators with non-local coupling
that display a spatio-temporal pattern of coherent and incoherent motion.
We present here a detailed analysis
of the spectral properties for such trajectories.
First, we study numerically their Lyapunov spectrum
and its behavior for an increasing number of oscillators.
The spectra demonstrate the hyperchaotic nature of the chimera states
and show a correspondence of the Lyapunov dimension
with the number of incoherent oscillators.
Then, we pass to the thermodynamic limit equation
and present an analytic approach
to the spectrum of a corresponding linearized evolution operator.
We show that in this setting, the chimera state is neutrally stable
and that the continuous spectrum coincides with the limit
of the hyperchaotic Lyapunov spectrum obtained for the finite size systems.
Lyapunov and exponential dichotomy spectral theory is extended
from ordinary differential equations (ODEs) to nonautonomous
differential-algebraic equations (DAEs). By using orthogonal
changes of variables, the original DAE system is transformed into
appropriate condensed forms, for which concepts such as Lyapunov
exponents, Bohl exponents, exponential dichotomy and spectral
intervals of various kinds can be analyzed via the resulting
underlying ODE. Some essential differences between the spectral
theory for ODEs and that for DAEs are pointed out. Numerical
methods for computing the spectral intervals associated with
Lyapunov and Sacker-Sell (exponential dichotomy) spectra are
derived by modifying and extending those methods proposed for ODEs. Perturbation theory and error analysis are discussed, as
well. Finally, some numerical examples are presented to illustrate
the theoretical results and the properties of the numerical
methods.