Compactly Supported Tight Frames Associated with Refinable Functions

It is well known that in applied and computational mathematics, cardinal B-splines play an important role in geometric modeling (in computer-aided geometric design), statistical data representation (or modeling), solution of differential equations (in numerical analysis), and so forth. More recently, in the development of wavelet analysis, cardinal B-splines also serve as a canonical example of scaling functions that generate multiresolution analyses of L²(−∞,∞). However, although cardinal B-splines have compact support, their corresponding orthonormal wavelets (of Battle and Lemarie) have infinite duration. To preserve such properties as self-duality while requiring compact support, the notion of tight frames is probably the only replacement of that of orthonormal wavelets. In this paper, we study compactly supported tight frames Ψ={ψ¹,…,ψ^N} for L²(−∞,∞) that correspond to some refinable functions with compact support, give a precise existence criterion of Ψ in terms of an inequality condition on the Laurent polynomial symbols of the refinable functions, show that this condition is not always satisfied (implying the nonexistence of tight frames via the matrix extension approach), and give a constructive proof that when Ψ does exist, two functions with compact support are sufficient to constitute Ψ, while three guarantee symmetry/anti-symmetry, when the given refinable function is symmetric.     

Cesàro means of double Walsh-Fourier series

We investigate approximation properties of Cesàro (C; −α, −β)-means of double Walsh-Fourier series with α, β ∈ (0, 1). As an application, we obtain a sufficient condition for the convergence of the means σ _n,m ^/−α,−β (f; x, y) to f(x,y) in the L ^p -metric, p ∈ [1, ∞]. We also show that this sufficient condition cannot be improved.     

Augmented Lagrangian function, non-quadratic growth condition and exact penalization

In this paper, under the assumption that the perturbation function satisfies a growth condition, necessary and sufficient conditions for an exact penalty representation and a zero duality gap property between the primal problem and its augmented Lagrangian dual problem are established.     

Lower bounds for the cardinality of minimal blocking sets in projective spaces

In this paper new lower bounds for the cardinality of minimal m-blocking sets are determined. Let r₂(q) be the number such that q+r₂(q)+1 is the cardinality of the smallest non-trivial line-blocking set in a plane of order q. If B is a minimal m-blocking set in PG(n,q) that contains at most q^m+q^m−1+…+q+1+r₂(q)·(∑_i=2m−n′^m−1qⁱ) points for an integer n′ satisfying mn′2m, then the dimension of B is at most n′. If the dimension of B is n′, then the following holds. The cardinality of B equals q^m+q^m−1+…+q+1+r₂(q)(∑_i=2m−n′^m−1qⁱ). For n′=m the set B is an m-dimensional subspace and for n′=m+1 the set B is a cone with an (m−2)-dimensional vertex over a non-trivial line-blocking set of cardinality q+r₂(q)+1 in a plane skew to the vertex. This result is due to Heim (Mitt. Math. Semin. Giessen 226 (1996), 4–82). For n′>m+1 and q not a prime the number q is a square and for q16 the set B is a Baer cone. If q is odd and |B|<q^m+q^m−1+…+q+1+r₂(q)(q^m−1+q^m−2), it follows from this result that the subspace generated by B has dimension at most m+1. Furthermore we prove that in this case, if , then B is an m-dimensional subspace or a cone with an (m−2)-dimensional vertex over a non-trivial line-blocking set of cardinality q+r₂(q)+1 in a plane skew to the vertex. For q=p^3h, p7 and q not a square we show this assertion for |B|q^m+q^m−1+…+q+1+q^2/3·(q^m−1+…+1).     

Convergence of the Augmented Lagrangian Method for Nonlinear Optimization Problems over Second-Order Cones

The paper analyzes the rate of local convergence of the augmented Lagrangian method for nonlinear second-order cone optimization problems. Under the constraint nondegeneracy condition and the strong second order sufficient condition, we demonstrate that the sequence of iterate points generated by the augmented Lagrangian method locally converges to a local minimizer at a linear rate, whose ratio constant is proportional to 1/τ with penalty parameter τ not less than a threshold . Importantly and interestingly enough, the analysis does not require the strict complementarity condition. Supported by the National Natural Science Foundation of China under Project 10771026 and by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry.     

Second order optimality conditions for bilevel set optimization problems

In this paper, we introduce a new notion of augmenting function known as indicator augmenting function to establish a minmax type duality relation, existence of a path of solution converging to optimal value and a zero duality gap relation for a nonconvex primal problem and the corresponding Lagrangian dual problem. We also obtain necessary and sufficient conditions for an exact penalty representation in the framework of indicator augmented Lagrangian.     

Some characteristic functions of unimodal distributions

A sufficient condition of Pólya type is given for a function to be the Fourier transform of a unimodal distribution. This condition is used to show that exp(− ¦ x ¦³ − 3 ¦ x ¦) is the Fourier transform of a unimodal distribution.     

Pressure Stabilization Method for Steady Viscous Flows in a System of Pipes

This paper deals with a singular perturbation of the stationary Stokes and Navier-Stokes systems. The term ε²Δp is added to the continuity equation, where ε is a small parameter. For a domain with cylindrical outlets to infinity and exponentially decaying data, existence and uniqueness of solutions under flux conditions at infinity are established for the linear problem and also for the nonlinear problem in the case of small data. Asymptotically exact estimates are proved for ε tending to zero. For sufficiently regular data, these estimates imply the convergence in H _loc ^5/2−δ for the velocity parts and in H _loc ^3/2−δ for the pressure parts, respectively. Bibliography: 17 titles.Dedicated to V. A. Solonnikov on the occasion of his 70th birthday__________Published in Zapiski Nauchnykh Seminarov POMI, Vol. 306, 2003, pp. 107–133.     

Application of Hilbert＇s Projective Metric on Symmetric Cones

Let Ω be a symmetric cone. In this note, we introduce Hilbert＇s projective metric on Ω in terms of Jordan algebras and we apply it to prove that, given a linear invertible transformation g such that g（Ω） = Ω and a real number p, ｜p｜ 〉 1, there exists a unique element x ∈ Ω satisfying g（x） = x^p.     

Asymptotic properties of Cramér-Smirnov statistics—A new approach

Let Y₁,…, Y_n be independent identically distributed random variables with distribution function F(x, θ), θ = (θ′₁, θ′₂), where θ_i (i = 1, 2) is a vector of p_i components, p = p₁ + p₂ and for θI, an open interval in ^p, F(x, θ) is continuous. In the present paper the author shows that the asymptotic distribution of modified Cramér-Smirnov statistic under H_n: θ₁ = θ₁₀ + n^−1/2γ, θ₂ unspecified, where γ is a given vector independent of n, is the distribution of a sum of weighted noncentral χ₁² variables whose weights are eigenvalues of a covariance function of a Gaussian process and noncentrality parameters are Fourier coefficients of the mean function of the Gaussian process. Further, the author exploits the special form of the covariance function by using perturbation theory to obtain the noncentrality parameters and the weights. The technique is applicable to other goodness-of-fit statistics such as U² [G. S. Watson, Biometrika 48 (1961), 109–114].     

𝒪-Operators of Loday Algebras and Analogues of the Classical Yang–Baxter Equation

We introduce notions of 𝒪-operators of the Loday algebras including the dendriform algebras and quadri-algebras as a natural generalization of Rota–Baxter operators. The invertible 𝒪-operators give a sufficient and necessary condition on the existence of the 2 ⁿ⁺¹ operations on an algebra with the 2 ⁿ operations in an associative cluster. The analogues of the classical Yang–Baxter equation in these algebras can be understood as the 𝒪-operators associated to certain dual bimodules. As a byproduct, the constraint conditions (invariances) of nondegenerate bilinear forms on these algebras are given.     

Perturbation theory in large order

For many quantum mechanical models, the behavior of perturbation theory in large order is strikingly simple. For example, in the quantum anharmonic oscillator, which is defined by−″ + (x²/4+εx⁴/4−E)y=0, y(±∞)=0,the perturbation coefficients A_n in the expansion for the ground-state energysimplify dramatically as n → ∞:.We use the methods of applied mathematics to investigate the nature of perturbation theory in quantum mechanics and we show that its large-order behavior is determined by the semiclassical content of the theory. In quantum field theory the perturbation coefficients are computed by summing Feynman graphs. We present a statistical procedure in a simple λ⁴ model for summing the set of all graphs as the number of vertices → ∞. Finally, we discuss the connection between the large-order behavior of perturbation theory in quantum electrodynamics and the value of α, the charge on the electron.     

A Two-Parameter Family of Orthogonal Polynomials with Respect to a Jacobi-Type Weight on the Unit Circle

A two-parameter family of polynomials is introduced by a recursion formula. The polynomials are orthogonal on the unit circle with respect to the weight ω_{α, β}(θ) = |(1 − z)^α(1 + z)^β|², α, β > − , z = e^iθ. Explicit representation, norm estimates, shift identities, and explicit connection to Jacobi polynomials on the real interval [−1, 1] is presented.     

Minimax estimation of the mean of spherically symmetric distributions under general quadratic loss

For X one observation on a p-dimensional (p ≥ 4) spherically symmetric (s.s.) distribution about θ, minimax estimators whose risks dominate the risk of X (the best invariant procedure) are found with respect to general quadratic loss, L(δ, θ) = (δ − θ)′ D(δ − θ) where D is a known p × p positive definite matrix. For C a p × p known positive definite matrix, conditions are given under which estimators of the form δ_a,r,C,D(X) = (I − (ar(|X|²)) D^−1/2CD^1/2 |X|⁻²)X are minimax with smaller risk than X. For the problem of estimating the mean when n observations X₁, X₂, …, X_n are taken on a p-dimensional s.s. distribution about θ, any spherically symmetric translation invariant estimator, δ(X₁, X₂, …, X_n), with have a s.s. distribution about θ. Among the estimators which have these properties are best invariant estimators, sample means and maximum likelihood estimators. Moreover, under certain conditions, improved robust estimators can be found.     

Minimax estimation of means of multivariate normal mixtures

Assume X = (X₁, …, X_p)′ is a normal mixture distribution with density w.r.t. Lebesgue measure, , where Σ is a known positive definite matrix and F is any known c.d.f. on (0, ∞). Estimation of the mean vector under an arbitrary known quadratic loss function _Q(θ, a) = (a − θ)′ Q(a − θ), Q a positive definite matrix, is considered. An unbiased estimator of risk is obatined for an arbitrary estimator, and a sufficient condition for estimators to be minimax is then achieved. The result is applied to modifying all the Stein estimators for the means of independent normal random variables to be minimax estimators for the problem considered here. In particular the results apply to the Stein class of limited translation estimators.     

On the convergence of a smoothed penalty algorithm for semi-infinite programming

For semi-infinite programming (SIP), we consider a class of smoothed penalty functions, which approximate the exact $l_\rho (0<\rho \le 1)$ penalty functions. On base of the smoothed penalty function, we present a feasible penalty algorithm for solving SIP. Without any boundedness condition or coercive condition, we establish the global convergence theorem. Then we present a perturbation theorem for this algorithm and obtain a necessary and sufficient condition for the convergence to the optimal value of SIP. Under Mangasarian–Fromovitz constrained qualification condition, we further discuss the convergence properties for the algorithm based upon a subclass of smooth approximations to the exact $l_\rho $ penalty function. Finally, numerical results are given.     

On the approximate behavior of the posterior distribution for an extreme multivariate observation

The behavior of the posterior for a large observation is considered. Two basic situations are discussed; location vectors and natural parameters.Let X = (X₁, X₂, …, X_n) be an observation from a multivariate exponential distribution with that natural parameter Θ = (Θ₁, Θ₂, …, Θ_n). Let θ_x^* be the posterior mode. Sufficient conditions are presented for the distribution of Θ − θ_x^* given X = x to converge to a multivariate normal with mean vector 0 as |x| tends to infinity. These same conditions imply that E(Θ | X = x) − θ_x^* converges to the zero vector as |x| tends to infinity.The posterior for an observation X = (X₁, X₂, …, X_n is considered for a location vector Θ = (Θ₁, Θ₂, …, Θ_n) as x gets large along a path, γ, in Rⁿ. Sufficient conditions are given for the distribution of γ(t) − Θ given X = γ(t) to converge in law as t → ∞. Slightly stronger conditions ensure that γ(t) − E(Θ | X = γ(t)) converges to the mean of the limiting distribution.These basic results about the posterior mean are extended to cover other estimators. Loss functions which are convex functions of absolute error are considered. Let δ be a Bayes estimator for a loss function of this type. Generally, if the distribution of Θ − E(Θ | X = γ(t)) given X = γ(t) converges in law to a symmetric distribution as t → ∞, it is shown that δ(γ(t)) − E(Θ | X = γ(t)) → 0 as t → ∞.     

The diffraction problem for a half plane with different face impedances revisited

The problem of the diffraction of an electromagnetic wave by a half plane with different face impedances is dealt with, following a rigorous approach based on the [L₂⁺( ]₂ theory of systems of Wiener-Hopf equations with piecewise continuous presymbols. The corresponding operator is defined in spaces of physically admissible solutions, the Sobolev spaces H_α⁺( )×H_α−1⁺( ) for , and its Fredholm characteristics are determined. For it is shown that the operators are invertible and their inverses are calculated. In the final section the inverse of a related operator presented by Meister and Speck is also obtained.     

Jacobson’s lemma and Cline’s formula for generalized inverses in a ring with involution

Abstract

In a ring with involution, we prove that a Drazin invertible element is pseudo core invertible if and only if its spectral idempotent is {1, 4}-invertible. As its applications, we obtain necessary and sufficient conditions for 1???ba (resp., ba) being pseudo core invertible while 1???ab (resp., ab) has pseudo core inverse, and the pseudo core inverse of 1???ba (resp., ba) is given in terms of 1???ab (resp., ab). Inspired by the above idea, Jacobson’s lemma for Moore-Penrose inverse is considered.     

Robust M-estimators of location vectors

Let X₁,…,X_n be i.i.d. random vectors in R^m, let θεR^m be an unknown location parameter, and assume that the restriction of the distribution of X₁−θ to a sphere of radius d belongs to a specified neighborhood of distributions spherically symmetric about 0. Under regularity conditions on and d, the parameter θ in this model is identifiable, and consistent M-estimators of θ (i.e., solutions of Σ_i=1ⁿψ(|X_i− |)(X_i− )=0) are obtained by using “re-descenders,” i.e., ψ's wh satisfy ψ(x)=0 for x≥c. An iterative method for solving for is shown to produce consistent and asymptotically normal estimates of θ under all distributions in . The following asymptotic robustness problem is considered: finding the ψ which is best among the re-descenders according to Huber's minimax variance criterion.