Root-Clustering Criteria (I): The Composite-Matrix Approach

In a series of articles we investigate different approachesto root clustering criteria. In particular, given A C^nxn andan algebraic region A C, we are looking for a criterion (asingle set of inequalities) ensuring the inclusion (A) A. Here,in the first article, we collect previous ideas on transformableregions and composite matrices, state new results, and correctand refine others.     

A Simple Proof for a Spectral Factorization Theorem

It is shown using simple methods that the transform Z(z), z C, of the coefficient sequence of the Wold decomposition ofany full-rank wide-sense stationary purely non-deterministicstochastic process satisfies (i) Z(z) H² (D) and (ii) Z^–1(z) H(D). Further it is shown that all spectral factors satisfying(i) and (ii) are equal up to right multiplication by orthogonalmatrices, and that among these the normalized (Z(0) =I) spectralfactors are equal to the transform of the Wold decomposition.An elementary proof of Youla's Theorem is then given togetherwith a simple proof that the rows of a Cholesky factor of abanded block Toeplitz matrix converge to the coefficients ofa stable matrix polynomial. Computer and Automation Institute of the Hungarian Academyof Sciences, Budapes, Hungary.     

On the Compactness of the Nonlinear Evolution Operator in a Banach Space

Let X be a real Banach space and let A(t): X 2^x be dissipativefor all t(0, T). Assume that {A(t)} generates an evolution operatorU(t, s) of type (D, , f). Necessary and sufficient conditionsare given for the compactness of U(t, s) for 0 s < t <T.     

Metric Entropy of Convex Hulls in Hilbert Spaces

We show in this note the following statement which is an improvementover a result of R. M. Dudley and which is also of independentinterest. Let X be a set of a Hilbert space with the propertythat there are constants , >0, and for each n N, the setX can be covered by at most n balls of radius n^–. Then,for each n N, the convex hull of X can be covered by 2ⁿ ballsof radius . The estimate is best possible for all n N, apart from the value c=c(, , X).In other words, let N(, X), >0, be the minimal number ofballs of radius covering the set X. Then the above result isequivalent to saying that if N(, X)=O(^–1/) as 0, thenfor the convex hull conv (X) of X, N(, conv (X)) =O(exp(^–2/(12))). Moreover, we give an interplay between several coveringparameters based on coverings by balls (entropy numbers) andcoverings by cylindrical sets (Kolmogorov numbers). 1991 MathematicsSubject Classification 41A46.     

A necessary and sufficient algebraic condition for the controllability of a thermoelastic plate equation

In this paper we give a necessary and sufficient algebraic conditionfor the approximate controllability of the following thermoelasticplate equation with Dirichlet boundary conditions w_tt + ²w + = a₁(x)u₁ + ... + a_m(x)u_m, t 0, x , _t – ß – w_t = d₁(x)u₁ + ... + d_m(x)u_m,t 0, x , = w = w = 0, t 0, x , where 0, ß > 0, is a sufficiently regular boundeddomain in R^N, a_i, d_i, L² (; R), the control functions u_i L²(0, t₁; R); i = 1, 2, ..., m. This condition is easy to checkand is given by Rank [P_jBA_jP_jBA²_jP_jB ... A^3_j–1_jP_jB] = 3_j,BU=b₁U₁+...+b_mU_m,b_i=[0, a_i, d_i], A_j=[0, –²_j, 0, 1, 0, –_j, 0, _j, –ß_j]P_j, j1, where _j, S are the eigenvalues of – with Dirichlet boundarycondition and P_j, S are the projections on the correspondingeigenspace.     

Sparse Systems of Functions Closed on Large Sets in RN

Let E(Z) = {e^inx}_nZ denote the trigonometrical exponential system.It is well known that E(Z) forms an orthogonal basis in thespace L²(0, 2). In 1964, H. Landau discovered that the trigonometricalsystem has the following property: certain small perturbationsof E(Z) yield exponential systems which are complete in L² onany finite union of 2-periodic translations of any interval(, 2–), 0 < < .     

Approximation of Ground State Eigenvalues and Eigenfunctions of Dirichlet Laplacians

Let be a bounded connected open set in R^N, N 2, and let –0be the Dirichlet Laplacian defined in L²(). Let > 0 be thesmallest eigenvalue of –, and let > 0 be its correspondingeigenfunction, normalized by ||||₂ = 1. For sufficiently small>0 we let R() be a connected open subset of satisfying Let – 0 be the Dirichlet Laplacian on R(), and let >0and >0 be its ground state eigenvalue and ground state eigenfunction,respectively, normalized by ||||₂=1. For functions f definedon , we let Sf denote the restriction of f to R(). For functionsg defined on R(), we let Tg be the extension of g to satisfying 1991 Mathematics SubjectClassification 47F05.     

BOUNDARY CONCENTRATION IN RADIAL SOLUTIONS TO A SYSTEM OF SEMILINEAR ELLIPTIC EQUATIONS

We study concentration phenomena for the system in the unit ball B₁ of ³ with Dirichlet boundaryconditions. Here , , > 0 and p > 1. We prove the existenceof positive radial solutions (, ) such that concentrates ata distance (/2)|log | away from the boundary B₁ as the parameter tends to 0. The approach is based on a combination of Lyapunov–Schmidtreduction procedure together with a variational method.     

Approximation of exact controllability problem involving parabolic differential equations

In this paper we examine computation of optimal control u^* ofthe exact controllability problem (referred to as the constraintproblem) governed by the following type of linear parabolicdifferential equations: (y/t) + Ay = u in Q y = 0 on y(0) = y₀ on where A is the second-order elliptic differential operator, is a bounded domain in ^k with smooth boundary , Q = (0, T)x , = (0, T) x and T > 0. This is achieved by approximatingu^* through a sequence {u_n} of controls corresponding to unconstrainedproblems involving a penalty function arising from the controllabilityconstraint.     

Spectra of graph neighborhoods and scattering

Let (G)_>0 be a family of ‘-thin’ Riemannian manifoldsmodeled on a finite metric graph G, for example, the -neighborhoodof an embedding of G in some Euclidean space with straight edges.We study the asymptotic behavior of the spectrum of the Laplace–Beltramioperator on G, as 0, for various boundary conditions. We obtaincomplete asymptotic expansions for the kth eigenvalue and theeigenfunctions, uniformly for kC^–1, in terms of scatteringdata on a non-compact limit space. We then use this to determinethe quantum graph which is to be regarded as the limit object,in a spectral sense, of the family (G). Our method is a directconstruction of approximate eigenfunctions from the scatteringand graph data, and the use of a priori estimates to show thatall eigenfunctions are obtained in this way.     

Turning Point Problems and Resonance

Consider the boundary value problem: ²yⁿ + (xp(x) + ²f(x, ))y'+ g(x, )y = 0, y(a) = A, y(b) = B, where a < 0 < b, p(x)< p(x) < 0, and p, f, and g are analytic. We investigatethe solution of this problem for small positive values of theparameter . If-g(0, 0)/p(0) c where c N = {0, 1, 2, 3,...},then so-called resonance does not occur, and y = o(ⁿ) on closedsubintervals of (a, b), for any n N, with expected boundarylayer behaviour at the end-points. If -g(0, 0)/p(0) = c, c N, then further transformations of dependent and independentvariables may still expose resonance or non-resonance. The setof necessary conditions that is developed is compared to otherauthors' criteria, most notably, Olver's sufficiency condition,and the necessary conditions of Cook & Eckhaus, Lakin, andMatkowsky. Finally, it is proved that these conditions are necessaryfor resonance.     

The Unmixing Problem

High-gain adaptive stabilization of minimum-phase linear systemsraises the so-called unmixing problem: finding a finite setA of matrices, with the property that for any M GL(n), (AM)C₊ for at least one A . First, some preliminary results will be given, leading up to a non-constructive proofof the existence of such finiteA's. For the cases n = 2 and3, a fairly satisfying solution is presented. In higher dimensions,the situation is more complicated, and negative results aregiven. Based on the Euler angles for O(n), an explicit solutionis given for arbitrary n.     

COMPRESSIONS, CONVEX GEOMETRY AND THE FREIMAN-BILU THEOREM

We note a link between combinatorial results of Bollobásand Leader concerning sumsets in the grid, the Brunn–Minkowskitheorem and a result of Freiman and Bilu concerning the structureof sets A with small doubling. Our main result is the following. If > 0 and if A is a finitenon-empty subset of a torsion-free abelian group with |A + A| K|A|, then A may be covered by e^KO(1) progressions of dimension log ₂ K + and size at most |A|.     

The Critical Conditions in a Reactive Slab with Slowly Varying Surface Temperature

Liouville's non-linear partial differential equation is consideredfor an infinite rectangular strip domain with a slowly varyingboundary condition. The equation describes a layer of chemicallyreactive material under conditions where the resistance to surfaceheat transfer is negligible and the ambient temperature variesslowly along the surface. Symmetrical heating by a zero orderexothermic reaction is assumed. If is a small dimensionlesstemperature difference between regions where the surface temperatureis effectively constant, a perturbation series solution in may be determined provided the Frank-Kamenetskii parameter satisfies _c(). It is shown that a plausible value for thecritical parameter is _c() = _c(0) e^–e,where _c(0) = 0.878.The corresponding critical temperature distribution is shownto have a dependence on different from that for subcriticalcases.     

Convergence of Univariate Quasi-Interpolation Using Multiquadrics

Quasi-interpolants to a function f: RR on an infinite regularmesh of spacing h can be defined by where :RR is a function with fast decay for large argument. In the approach employing the radial-basis-function : RR, thefunction is a finite linear combination of basis functions(|•–jh|) (jZ). Choosing Hardy's multiquadrics (r)=(r²+c²)^?,we show that sufficiently fast-decaying exist that render quasi-interpolationexact for linear polynomials f. Then, approximating f C²(R),we obtain uniform convergence of s to f as (h, c)0, and convergenceof s' to f' as (h, c²/h)0. However, when c stays bounded awayfrom 0 as h0, there are f C(R) for which s does not convergeto f as h0. We also show that, for all which vanish at infinity but arenot integrable over R, there are no finite linear combinations of the given basis functions allowing the construction of admissiblequasi-interpolants. This includes the case of the inverse multiquadncs(r)=(r²+c²)^–?.     

Rate of convergence of numerical approximations to homoclinic bifurcation points

For x=f (x, ), x Rⁿ, R, having a hyperbolic or semihyperbolicequilibrium p(), we study the numerical approximation of parametervalues ^* at which there is an orbit homoclinic to p(). We approximate^* by solving a finite-interval boundary value problem on J=[T_–,T₊], T_–<0<T₊, with boundary conditions that sayx(T_–) and x(T₊) are in approximations to appropriate invariantmanifolds of p(). A phase condition is also necessary to makethe solution unique. Using a lemma of Xiao-Biao Lin, we improve,for certain phase conditions, existing estimates on the rateof convergence of the computed homoclinic bifurcation parametervalue , to the true value ^*. The estimates we obtain agree withthe rates of convergence observed in numerical experiments.Unfortunately, the phase condition most commonly used in numericalwork is not covered by our results.     

Finite-Element Approximation of a Plasma Equilibrium Problem

The plasma problem studied is: given R⁺ find (, d, u) R ?R ? H¹() such that Let ₁ < ₂ be the first two eigenvalues of the associatedlinear eigenvalue problem: find $$\left(\lambda ,\phi \right)\in\mathrm{R;}\times {\hbox{ H }}_{0}^{1}\left(\Omega \right)$$such that For ₀(0,₂) it is well known that there exists a unique solution(₀, d₀, u₀) to the above problem. We show that the standard continuous piecewise linear Galerkinfinite-element approximatinon $$\left({\lambda }_{0},{\hbox{d }}_{0}^{k},{u}_{0}^{h}\right)$$, for ₀(0,₂), converges atthe optimal rate in the H¹, L², and L norms as h, the mesh length,tends to 0. In addition, we show that dist (, ^h)Ch² ln 1/h,where $${\Gamma }^{\left(h\right)}=\left\{x\in \Omega :{u}_{0}^{\left(h\right)}\left(x\right)=0\right\}$$.Finally we consider a more practical approximation involvingnumerical integration.     

Local error estimates for radial basis function interpolation of scattered data

Introducing a suitable variational formulation for the localerror of scattered data interpolation by radial basis functions(r), the error can be bounded by a term depending on the Fouriertransform of the interpolated function f and a certain ‘Krigingfunction’, which allows a formulation as an integral involvingthe Fourier transform of . The explicit construction of locallywell-behaving admissible coefficient vectors makes the Krigingfunction bounded by some power of the local density h of datapoints. This leads to error estimates for interpolation of functionsf whose Fourier transform f is ‘dominated’ by thenonnegative Fourier transform of (x) = (||x||) in the sense . Approximation orders are arbitrarily high for interpolationwith Hardy multiquadrics, inverse multiquadrics and Gaussiankernels. This was also proven in recent papers by Madych andNelson, using a reproducing kernel Hilbert space approach andrequiring the same hypothesis as above on f, which limits thepractical applicability of the results. This work uses a differentand simpler analytic technique and allows to handle the casesof interpolation with (r) = r^s for s R, s > 1, s 2N, and(r) = r^s log r for s 2N, which are shown to have accuracy O(h^s/2)     

The Set of {omega}-Limit Points of A Map of the Circle is Closed

Let f be a continuous self-map of the unit circle, S¹. The -limitpoints (x) of a point x are the set of all limit points of thesequence of iterates of f acting on x. We shall show that theset of all -limit points _xS¹(x) a closed set in S¹.