Optimal control of fractional systems: a diffusive formulation

Optimal control of fractional linear systems on a finite horizon can be classically formulated using the adjoint system. But the adjoint of a causal fractional integral or derivative operator happens to be an anti-causal operator: hence, the adjoint equations are not easy to solve in the first place. Using an equivalent diffusive realization helps transform the original problem into a coupled system of PDEs, for which the adjoint system can be more easily derived and properly studied

Wild monodromy and automorphisms of curves

Let $R$ be a complete discrete valuation ring of mixed characteristic $(0,p)$ with field of fractions $K$ containing the $p$-th roots of unity. This paper is concerned with semi-stable models of $p$-cyclic covers of the projective line C \la \PK. We start by providing a new construction of a semi-stable model of $C$ in the case of an equidistant branch locus. If the cover is given by the Kummer equation $Z^p=f(X_0)$ we define what we called the monodromy polynomial ${\mathcal L}(Y)$ of $f(X_0)$; a polynomial with coefficients in $K$. Its zeros are key to obtaining a semi-stable model of $C$. As a corollary we obtain an upper bound for the minimal extension $K'/K$ over which a stable model of the curve $C$ exists. Consider the polynomial ${\cal L}(Y)\prod(Y^p-f(y_i))$ where the $y_i$ range over the zeros of ${\cal L}(Y)$. We show that the splitting field of this polynomial always contains $K'$, and that in some instances the two fields are equal.Comment: The final version of this article will be published in the Duke Mathematical Journal, published by Duke University pres

Diffusive approximation of a time-fractional Burger's equation in nonlinear acoustics

A fractional time derivative is introduced into the Burger's equation to model losses of nonlinear waves. This term amounts to a time convolution product, which greatly penalizes the numerical modeling. A diffusive representation of the fractional derivative is adopted here, replacing this nonlocal operator by a continuum of memory variables that satisfy local-in-time ordinary differential equations. Then a quadrature formula yields a system of local partial differential equations, well-suited to numerical integration. The determination of the quadrature coefficients is crucial to ensure both the well-posedness of the system and the computational efficiency of the diffusive approximation. For this purpose, optimization with constraint is shown to be a very efficient strategy. Strang splitting is used to solve successively the hyperbolic part by a shock-capturing scheme, and the diffusive part exactly. Numerical experiments are proposed to assess the efficiency of the numerical modeling, and to illustrate the effect of the fractional attenuation on the wave propagation.Comment: submitted to Siam SIA

Diffusive representations for fractional Laplacian: systems theory framework and numerical issues

Bridging the gap between an abstract definition of pseudo-differential operators, such as (-\Delta)^{\gamma} for - 1/2 < \gamma < 1/2, and a concrete way to represent them has proved difficult; deriving stable numerical schemes for such operators is not an easy task either. Thus, the framework of diffusive representations, as already developed for causal fractional integrals and derivatives, is being applied to fractional Laplacian: it can be seen as an extension of the Wiener-­Hopf factorization and splitting techniques to irrational transfer functions

Systems control theory applied to natural and synthetic musical sounds

Systems control theory is a far developped field which helps to study stability, estimation and control of dynamical systems. The physical behaviour of musical instruments, once described by dynamical systems, can then be controlled and numerically simulated for many purposes. The aim of this paper is twofold: first, to provide the theoretical background on linear system theory, both in continuous and discrete time, mainly in the case of a finite number of degrees of freedom ; second, to give illustrative examples on wind instruments, such as the vocal tract represented as a waveguide, and a sliding flute

Stable Realization of a Delay System Modeling a Convergent Acoustic Cone

This paper deals with the physical modeling and the digital time simulation of acoustic pipes. We will study the simplified case of a single convergent cone. It is modeled by a linear system made of delays and a transfer function which represents the wave reflection at the entry of the cone. According to [1], the input/output relation of this system is causal and stable whereas the reflection function is unstable. In the continuous time-domain, a first state space representation of this delay system is done. Then, we use a change of state to separate the unobservable subspace and its orthogonal complement, which is observable. Whereas the unobservable part is unstable, it is proved that the observable part is stable, using the D-Subdivision method. Thus, isolating this latter observable subspace, to build the minimal realization, defines a stable system. Finally, a discrete-time version of this system is derived and is proved to be stable using the Jury criterion. The main contribution of this work is neither the minimal realization of the system nor the proofs of stability, but it is rather the solving of an old problem of acoustics which has heen achieved using standard tools of automatic control

Digital waveguide modeling for wind instruments: building a state-space representation based on the Webster-Lokshin model

This paper deals with digital waveguide modeling of wind instruments. It presents the application of state-space representations for the refined acoustic model of Webster-Lokshin. This acoustic model describes the propagation of longitudinal waves in axisymmetric acoustic pipes with a varying cross-section, visco-thermal losses at the walls, and without assuming planar or spherical waves. Moreover, three types of discontinuities of the shape can be taken into account (radius, slope and curvature). The purpose of this work is to build low-cost digital simulations in the time domain based on the Webster-Lokshin model. First, decomposing a resonator into independent elementary parts and isolating delay operators lead to a Kelly-Lochbaum network of input/output systems and delays. Second, for a systematic assembling of elements, their state-space representations are derived in discrete time. Then, standard tools of automatic control are used to reduce the complexity of digital simulations in the time domain. The method is applied to a real trombone, and results of simulations are presented and compared with measurements. This method seems to be a promising approach in term of modularity, complexity of calculation and accuracy, for any acoustic resonators based on tubes

On the singularities of fractional differential systems, using a mathematical limiting process based on physical grounds

Fractional systems are associated with irrational transfer functions for which nonunique analytic continuations are available (from some right-half Laplace plane to a maximal domain). They involve continuous sets of singularities, namely cuts, which link fixed branching points with an arbitrary path. In this paper, an academic example of the 1D heat equation and a realistic model of an acoustic pipe on bounded domains are considered. Both involve a transfer function with a unique analytic continuation and singularities of pole type. The set of singularities degenerates into uniquely defined cuts when the length of the physical domain becomes infinite. From a mathematical point of view, both the convergence in Hardy spaces of some right-half complex plane and the pointwise convergence are studied and proved

Puzzles in pipes with negative curvature: from the Webster PDE to stable numerical simulation in real time

Minimal realizations of a class of delay-differential systems are derived for the digital simulation of waveguides, modelled by the Webster horn equation. Studying their stability is an interesting issue, since negative curvatures could lead to unstable systems. Spectral properties of Toeplitz matrix play a key role in this work

State-space representation for digital waveguide networks of lossy flared acoustic pipes

This paper deals with digital waveguide modeling of wind instruments. It presents the application of state-space representations to the acoustic model of Webster-Lokshin. This acoustic model describes the propagation of longitudinal waves in axisymmetric acoustic pipes with a varying cross-section, visco-thermal losses at the walls, and without assuming planar or spherical waves. Moreover, three types of discontinuities of the shape can be taken into account (radius, slope and curvature), which can lead to a good fit of the original shape of pipe. The purpose of this work is to build low-cost digital simulations in the time domain, based on the Webster-Lokshin model. First, decomposing a resonator into independent elementary parts and isolating delay operators lead to a network of input/output systems and delays, of Kelly-Lochbaum network type. Second, for a systematic assembling of elements, their state-space representations are derived in discrete time. Then, standard tools of automatic control are used to reduce the complexity of digital simulations in time domain. In order to validate the method, simulations are presented and compared with measurements