results for au:Yao_J in:math

- Let $\mathbf{X}_n=(x_{ij})$ be a $k \times n$ data matrix with complex-valued, independent and standardized entries satisfying a Lindeberg-type moment condition. We consider simultaneously $R$ sample covariance matrices $\mathbf{B}_{nr}=\frac1n \mathbf{Q}_r \mathbf{X}_n \mathbf{X}_n^*\mathbf{Q}_r^\top,~1\le r\le R$, where the $\mathbf{Q}_{r}$'s are nonrandom real matrices with common dimensions $p\times k~(k\geq p)$. Assuming that both the dimension $p$ and the sample size $n$ grow to infinity, the limiting distributions of the eigenvalues of the matrices $\{\mathbf{B}_{nr}\}$ are identified, and as the main result of the paper, we establish a joint central limit theorem for linear spectral statistics of the $R$ matrices $\{\mathbf{B}_{nr}\}$. Next, this new CLT is applied to the problem of testing a high dimensional white noise in time series modelling. In experiments the derived test has a controlled size and is significantly faster than the classical permutation test, though it does have lower power. This application highlights the necessity of such joint CLT in the presence of several dependent sample covariance matrices. In contrast, all the existing works on CLT for linear spectral statistics of large sample covariance matrices deal with a single sample covariance matrix ($R=1$).
- Jan 10 2018 math.OC arXiv:1801.02752v1We study the equilibrium problem on general Riemannian manifolds. The results on existence of solutions and on the convex structure of the solution set are established. Our approach consists in relating the equilibrium problem to a suitable variational inequality problem on Riemannian manifolds, and is completely different from previous ones on this topic in the literature. As applications, the corresponding results for the mixed variational inequality and the Nash equilibrium are obtained. Moreover, we formulate and analyze the convergence of the proximal point algorithm for the equilibrium problem. In particular, correct proofs are provided for the results claimed in J. Math. Anal. Appl. 388, 61-77, 2012 (i.e., Theorems 3.5 and 4.9 there) regarding the existence of the mixed variational inequality and the domain of the resolvent for the equilibrium problem on Hadamard manifolds.
- Dec 29 2017 math.OC arXiv:1712.09743v1In this paper, we study second-order necessary and sufficient optimality conditions of Karush--Kuhn--Tucker-type for locally optimal solutions in the sense of Pareto to a class of multi-objective optimal control problems with mixed pointwise constraint. To deal with the problems, we first derive second-order optimality conditions for abstract multi-objective optimal control problems which satisfy the Robinson constraint qualification. We then apply the obtained results to our concrete problems. The proofs of obtained results are direct, self-contained without using scalarization techniques.
- A note on approximate Karush-Kuhn-Tucker conditions in locally Lipschitz multiobjective optimizationNov 27 2017 math.OC arXiv:1711.08551v1In the recent paper of Giorgi, Jiménez and Novo (J Optim Theory Appl 171:70--89, 2016), the authors introduced the so-called approximate Karush-Kuhn-Tucker (AKKT) condition for smooth multiobjective optimization problems and obtained some AKKT-type necessary optimality conditions and sufficient optimality conditions for weak efficient solutions of such a problem. In this note, we extend these optimality conditions to locally Lipschitz multiobjective optimization problems using Mordukhovich subdifferentials. Furthermore, we prove that, under some suitable additional conditions, an AKKT condition is also a KKT one.
- This paper studies the design and secrecy performance of linear multihop networks, in the presence of randomly distributed eavesdroppers in a large-scale two-dimensional space. Depending on whether there is feedback from the receiver to the transmitter, we study two transmission schemes: on-off transmission (OFT) and non-on-off transmission (NOFT). In the OFT scheme, transmission is suspended if the instantaneous received signal-to-noise ratio (SNR) falls below a given threshold, whereas there is no suspension of transmission in the NOFT scheme. We investigate the optimal design of the linear multiple network in terms of the optimal rate parameters of the wiretap code as well as the optimal number of hops. These design parameters are highly interrelated since more hops reduces the distance of per-hop communication which completely changes the optimal design of the wiretap coding rates. Despite the analytical difficulty, we are able to characterize the optimal designs and the resulting secure transmission throughput in mathematically tractable forms in the high SNR regime. Our numerical results demonstrate that our analytical results obtained in the high SNR regime are accurate at practical SNR values. Hence, these results provide useful guidelines for designing linear multihop networks with targeted physical layer security performance.
- In this paper, we study the secure transmission in multihop wireless networks with randomize-and-forward (RaF) relaying, in the presence of randomly distributed eavesdroppers. By considering adaptive encoder with on-off transmission (OFT) scheme, we investigate the optimal design of the wiretap code and routing strategies to maximize the secrecy rate while satisfying the secrecy outage probability (SOP) constraint. We derive the exact expressions for the optimal rate parameters of the wiretap code. Then the secure routing problem is solved by revising the classical Bellman-Ford algorithm. Simulation results are conducted to verify our analysis.
- Oct 12 2017 math.OC arXiv:1710.03989v2In the present paper, we discuss second-order optimality conditions in constrained vector optimization problems, where the objective functions and active constraint functions are locally Lipschitz at the referee point. By using the second-order upper generalized directional derivative and the second-order tangent set, we introduce a new type of second-order regularity condition in the sense of Abadie. Then we establish some second-order Karush-Kuhn-Tucker necessary optimality conditions in primal form for local (weak, Geoffrion properly) efficient solutions to the considered problem. Examples are given to illustrate the obtained results.
- Oct 02 2017 math.OC arXiv:1709.10227v1Generalized polyhedral convex optimization problems in locally convex Hausdorff topological vector spaces are studied systematically in this paper. We establish solution existence theorems, necessary and sufficient optimality conditions, weak and strong duality theorems. In particular, we show that the dual problem has the same structure as the primal problem, and the strong duality relation holds under three different sets of conditions.
- Aug 15 2017 math.PR arXiv:1708.03749v1This paper investigates the central limit theorem for linear spectral statistics of high dimensional sample covariance matrices of the form $\mathbf{B}_n=n^{-1}\sum_{j=1}^{n}\mathbf{Q}\mathbf{x}_j\mathbf{x}_j^{*}\mathbf{Q}^{*}$ where $\mathbf{Q}$ is a nonrandom matrix of dimension $p\times k$, and $\{\mathbf{x}_j\}$ is a sequence of independent $k$-dimensional random vector with independent entries, under the assumption that $p/n\to y>0$. A key novelty here is that the dimension $k\ge p$ can be arbitrary, possibly infinity. This new model of sample covariance matrices $\mathbf{B}_n$ covers most of the known models as its special cases. For example, standard sample covariance matrices are obtained with $k=p$ and $\mathbf{Q}=\mathbf{T}_n^{1/2}$ for some positive definite Hermitian matrix $\mathbf{T}_n$. Also with $k=\infty$ our model covers the case of repeated linear processes considered in recent high-dimensional time series literature. The CLT found in this paper substantially generalizes the seminal CLT in Bai and Silverstein (2004). Applications of this new CLT are proposed for testing the structure of a high-dimensional covariance matrix. The derived tests are then used to analyse a large fMRI data set regarding its temporary correlation structure.
- We present new connections among anomalous diffusion (AD), normal diffusion (ND) and the Central Limit Theorem. This is done by defining a point transformation to a new position variable, which we postulate to be Cartesian, motivated by considerations from super-symmetric quantum mechanics. Canonically quantizing in the new position and momentum variables according to Dirac gives rise to generalized negative semi-definite and self-adjoint Laplacian operators. These lead to new generalized Fourier transformations and associated probability distributions, which are form invariant under the corresponding transform. The new Laplacians also lead us to generalized diffusion equations, which imply a connection to the CLT. We show that the derived diffusion equations capture all of the Fractal and Non-Fractal Diffusion equations of O'Shaughnessy and Procaccia. However, we also obtain new equations that cannot (so far as we are able to tell) be expressed as examples of the O'Shaughnessy and Procaccia equations. These equations also possess asymptotics that are related to a CLT but with bi-modal distributions as limits. The results show, in part, that experimentally measuring the diffusion scaling law can determine the point transformation (for monomial point transformations). We also show that AD in the original, physical position is actually ND when viewed in terms of displacements in an appropriately transformed position variable. Finally, we show that there is a new, anomalous diffusion possible for bi-modal probability distributions that also display attractor behavior which is the consequence of an underlying CLT.
- Jul 14 2017 math.OC arXiv:1707.03955v1A parametric constrained convex optimal control problem, where the initial state is perturbed and the linear state equation contains a noise, is considered in this paper. Formulas for computing the subdifferential and the singular subdifferential of the optimal value function at a given parameter are obtained by means of some recent results on differential stability in mathematical programming. The computation procedures and illustrative examples are presented.
- May 22 2017 math.OC arXiv:1705.06892v1Generalized polyhedral convex sets, generalized polyhedral convex functions on locally convex Hausdorff topological vector spaces, and the related constructions such as sum of sets, sum of functions, directional derivative, infimal convolution, normal cone, conjugate function, subdifferential, are studied thoroughly in this paper. Among other things, we show how a generalized polyhedral convex set can be characterized via the finiteness of the number of its faces. In addition, it is proved that the infimal convolution of a generalized polyhedral convex function and a polyhedral convex function is a polyhedral convex function. The obtained results can be applied to scalar optimization problems described by generalized polyhedral convex sets and generalized polyhedral convex functions.
- We consider forecasting a single time series using high-dimensional predictors in the presence of a possible nonlinear forecast function. The sufficient forecasting (Fan et al., 2016) used sliced inverse regression to estimate lower-dimensional sufficient indices for nonparametric forecasting using factor models. However, Fan et al. (2016) is fundamentally limited to the inverse first-moment method, by assuming the restricted fixed number of factors, linearity condition for factors, and monotone effect of factors on the response. In this work, we study the inverse second-moment method using directional regression and the inverse third-moment method to extend the methodology and applicability of the sufficient forecasting. As the number of factors diverges with the dimension of predictors, the proposed method relaxes the distributional assumption of the predictor and enhances the capability of capturing the non-monotone effect of factors on the response. We not only provide a high-dimensional analysis of inverse moment methods such as exhaustiveness and rate of convergence, but also prove their model selection consistency. The power of our proposed methods is demonstrated in both simulation studies and an empirical study of forecasting monthly macroeconomic data from Q1 1959 to Q1 2016. During our theoretical development, we prove an invariance result for inverse moment methods, which make a separate contribution to the sufficient dimension reduction.
- Mar 14 2017 math.OC arXiv:1703.03966v1In this paper, we study constraint qualifications for the nonconvex inequality defined by a proper lower semicontinuous function. These constraint qualifications involve the generalized construction of normal cones and subdifferentials. Several conditions for these constraint qualifications are also provided therein. When restricted to the convex inequality, these constraint qualifications reduce to basic constraint qualification (BCQ) and strong BCQ studied in [SIAM J. Optim., 14(2004), 757-772] and [Math. Oper. Res., 30 (2005), 956-965].
- We propose a numerical recipe for risk evaluation defined by a backward stochastic differential equation. Using dual representation of the risk measure, we convert the risk valuation to a stochastic control problem where the control is a certain Radon-Nikodym derivative process. By exploring the maximum principle, we show that a piecewise-constant dual control provides a good approximation on a short interval. A dynamic programming algorithm extends the approximation to a finite time horizon. Finally, we illustrate the application of the procedure to risk management in conjunction with nested simulation.
- The theoretical analysis of detection and decoding of low-density parity-check (LDPC) codes transmitted over channels with two-dimensional (2D) interference and additive white Gaussian noise (AWGN) is provided in this paper. The detection and decoding system adopts the joint iterative detection and decoding scheme (JIDDS) in which the log-domain sum-product algorithm is adopted to decode the LDPC codes. The graph representations of the JIDDS are explained. Using the graph representations, we prove that the message-flow neighborhood of the detection and decoding system will be tree-like for a sufficiently long code length. We further confirm that the performance of the JIDDS will concentrate around the performance in which message-flow neighborhood is tree-like. Based on the tree-like message-flow neighborhood, we employ a modified density evolution algorithm to track the message densities during the iterations. A threshold is calculated using the density evolution algorithm which can be considered as the theoretical performance limit of the system. Simulation results demonstrate that the modified density evolution is effective in analyzing the performance of 2D interference systems.
- In additive white gaussian noise (AWGN) channel, chaos has been proved to be the optimal coherent communication waveform in the sense of using very simple matched filter to maximize the signal-to-noise ratio (SNR). Recently, Lyapunov exponent spectrum of the chaotic signals after being transmitted through a wireless channel has been shown to be unaltered, paving the way for wireless communication using chaos. In wireless communication systems, inter-symbol interference (ISI) caused by multipath propagation is one of the main obstacles to achieve high bit transmission rate and low bit error rate (BER). How to resist multipath effect is a fundamental problem in a chaos-based wireless communication system (CWCS). In this paper, implementation of a CWCS is presented. It is built to transmit chaotic signals generated by a hybrid dynamical system and then to filter the received signals by using the corresponding matched filter to decrease the noise effect and to detect the binary information. We find that the multipath effect can be effectively resisted by regrouping the return map of the received signal and by setting the corresponding threshold based on the available information. We show that the optimal threshold is a function of the channel parameters and of the transmitted information symbols. Practically, the channel parameters are time-variant, and the future information symbols are unavailable. In this case, a suboptimal threshold (SOT) is proposed, and the BER using the SOT is derived analytically. Simulation results show that the CWCS achieves a remarkable competitive performance even under inaccurate channel parameters.
- Jul 25 2016 math.OC arXiv:1607.06569v1In this paper, in terms of three types of generalized second-order derivatives of a nonsmooth function, we mainly study the corresponding second-order optimality conditions in a Hilbert space and prove the equivalence among these optimality conditions for paraconcave functions. As applications, we use these second-order optimality conditions to study strict local minimizers of order two and provide sufficient and/or necessary conditions for ensuring the local minimizer. This work extends and generalizes the study on second-order optimality conditions from the finite-dimensional space to the Hilbert space.
- Jul 05 2016 math.AP arXiv:1607.00581v1We investigate the existence and multiplicity of solutions to the following $p(x)$-Laplacian problem in $\mathbb{R}^{N}$ via critical point theory \beginequation* \left{ \beginarrayl -\bigtriangleup _p(x)u+V(x)\left\vert u\right\vert ^p(x)-2u=f(x,u),\text in \mathbbR^N, \{u∈W^1,p(⋅)(\mathbbR^N). \endarray \right. \endequation* We propose a new set of growth conditions which matches the variable exponent nature of the problem. Under this new set of assumptions, we manage to verify the Cerami compactness condition. Therefore, we succeed in proving the existence of multiple solutions to the above problem without the well-known Ambrosetti--Rabinowitz type growth condition. Meanwhile, we could also characterize the pointwise asymptotic behaviors of these solutions. In our main argument, the idea of localization, decomposition of the domain, regularity of weak solutions and comparison principle are crucial ingredients among others.
- Jul 05 2016 math.AP arXiv:1607.00584v1We investigate the following Dirichlet problem with variable exponents: \beginequation* \left{ \beginarrayl -\bigtriangleup _p(x)u=\lambda \alpha (x)\left\vert u\right\vert ^\alpha (x)-2u\left\vert v\right\vert ^\beta (x)+F_u(x,u,v),\text in \Omega , \\ -\bigtriangleup _q(x)v=\lambda \beta (x)\left\vert u\right\vert ^\alpha (x)\left\vert v\right\vert ^\beta (x)-2v+F_v(x,u,v),\text in \Omega , \{u=0=v,\text on ∂\Omega. \endarray \right. \endequation* We present here, in the system setting, a new set of growth conditions under which we manage to use a novel method to verify the Cerami compactness condition. By localization argument, decomposition technique and variational methods, we are able to show the existence of multiple solutions with constant sign for the problem without the well-known Ambrosetti--Rabinowitz type growth condition. More precisely, we manage to show that the problem admits four, six and infinitely many solutions respectively.
- May 04 2016 math.NT arXiv:1605.00813v1In this work we extend our study on a link between automaticity and certain algebraic power series over finite fields. Our starting point is a family of sequences in a finite field of characteristic $2$, recently introduced by the first author in connection with algebraic continued fractions. By including it in a large family of recurrent sequences in an arbitrary finite field, we prove its automaticity. Then we give a criterion on automatic sequences, generalizing a previous result and this allows us to present new families of automatic sequences in an arbitrary finite field.
- We prove that the famous diffusive Brusselator model can support more complicated spatial-temporal wave structure than the usual temporal-oscillation from a standard Hopf bifurcation. In our investigation, we discover that the diffusion term in the model is neither a usual parabolic stabilizer nor a destabilizer as in the Turing instability of uniform state, but rather plays the role of maintaining an equivariant Hopf bifurcation spectral mechanism. At the same time, we show that such a mechanism can occur around any nonzero wave number and this finding is also different from the former works where oscillations caused by diffusion can cause the growth of wave structure only at a particular wavelength. Our analysis also demonstrates that the complicated spatial-temporal oscillation is not solely driven by the inhomogeneity of the reactants.
- In this paper, we study the problem of secure routing in a multihop wireless ad-hoc network in the presence of randomly distributed eavesdroppers. Specifically, the locations of the eavesdroppers are modeled as a homogeneous Poisson point process (PPP) and the source-destination pair is assisted by intermediate relays using the decode-and-forward (DF) strategy. We analytically characterize the physical layer security performance of any chosen multihop path using the end-to-end secure connection probability (SCP) for both colluding and non-colluding eavesdroppers. To facilitate finding an efficient solution to secure routing, we derive accurate approximations of the SCP. Based on the SCP approximations, we study the secure routing problem which is defined as finding the multihop path having the highest SCP. A revised Bellman-Ford algorithm is adopted to find the optimal path in a distributed manner. Simulation results demonstrate that the proposed secure routing scheme achieves nearly the same performance as exhaustive search.
- Oct 01 2015 math.NT arXiv:1509.09075v1The aim of this note is to show the existence of a correspondance between certain algebraic continued fractions in fields of power series over a finite field and automatic sequences in the same finite field. this connection is illustrated by three families of examples and a counterexample.
- Sep 25 2015 math.OC arXiv:1509.07264v2We study some basic properties of the function $f_0:M\rightarrow\IR$ on Hadamard manifolds defined by $$ f_0(x):=\langle u_0,\exp_x_0^-1x\rangle\quad\mboxfor any $x\in M$. $$ A characterization for the function to be linear affine is given and a counterexample on Poincaré plane is provided, which in particular, shows that assertions (i) and (ii) claimed in \cite[Proposition 3.4]Papa2009 are not true, and that the function $f_0$ is indeed not quasi-convex. Furthermore, we discuss the convexity properties of the sub-level sets of the function on Riemannian manifolds with constant sectional curvatures.
- Aug 24 2015 math.OC arXiv:1508.05316v2We consider optimal control problems for diffusion processes, where the objective functional is defined by a time-consistent dynamic risk measure. We focus on coherent risk measures defined by $g$-evaluations. For such problems, we construct a family of time and space perturbed systems with piecewise-constant control functions. We obtain a regularized optimal value function by a special mollification procedure. This allows us to establish a bound on the difference between the optimal value functions of the original problem and of the problem with piecewise-constant controls.
- We consider forecasting a single time series when there is a large number of predictors and a possible nonlinear effect. The dimensionality was first reduced via a high-dimensional (approximate) factor model implemented by the principal component analysis. Using the extracted factors, we develop a novel forecasting method called the sufficient forecasting, which provides a set of sufficient predictive indices, inferred from high-dimensional predictors, to deliver additional predictive power. The projected principal component analysis will be employed to enhance the accuracy of inferred factors when a semi-parametric (approximate) factor model is assumed. Our method is also applicable to cross-sectional sufficient regression using extracted factors. The connection between the sufficient forecasting and the deep learning architecture is explicitly stated. The sufficient forecasting correctly estimates projection indices of the underlying factors even in the presence of a nonparametric forecasting function. The proposed method extends the sufficient dimension reduction to high-dimensional regimes by condensing the cross-sectional information through factor models. We derive asymptotic properties for the estimate of the central subspace spanned by these projection directions as well as the estimates of the sufficient predictive indices. We further show that the natural method of running multiple regression of target on estimated factors yields a linear estimate that actually falls into this central subspace. Our method and theory allow the number of predictors to be larger than the number of observations. We finally demonstrate that the sufficient forecasting improves upon the linear forecasting in both simulation studies and an empirical study of forecasting macroeconomic variables.
- Consider two $p$-variate populations, not necessarily Gaussian, with covariance matrices $\Sigma_1$ and $\Sigma_2$, respectively, and let $S_1$ and $S_2$ be the sample covariances matrices from samples of the populations with degrees of freedom $T$ and $n$, respectively. When the difference $\Delta$ between $\Sigma_1$ and $\Sigma_2$ is of small rank compared to $p,T$ and $n$, the Fisher matrix $F=S_2^{-1}S_1$ is called a \em spiked Fisher matrix. When $p,T$ and $n$ grow to infinity proportionally, we establish a phase transition for the extreme eigenvalues of $F$: when the eigenvalues of $\Delta$ (\em spikes) are above (or under) a critical value, the associated extreme eigenvalues of the Fisher matrix will converge to some point outside the support of the global limit (LSD) of other eigenvalues; otherwise, they will converge to the edge points of the LSD. Furthermore, we derive central limit theorems for these extreme eigenvalues of the spiked Fisher matrix. The limiting distributions are found to be Gaussian if and only if the corresponding population spike eigenvalues in $\Delta$ are \em simple. Numerical examples are provided to demonstrate the finite sample performance of the results. In addition to classical applications of a Fisher matrix in high-dimensional data analysis, we propose a new method for the detection of signals allowing an arbitrary covariance structure of the noise. Simulation experiments are conducted to illustrate the performance of this detector.
- The irrationality exponent of an irrational number $\xi$, which measures the approximation rate of $\xi$ by rationals, is in general extremely difficult to compute explicitly, unless we know the continued fraction expansion of $\xi$. Results obtained so far are rather fragmentary, and often treated case by case. In this work, we shall unify all the known results on the subject by showing that the irrationality exponents of large classes of automatic numbers and Mahler numbers (which are transcendental) are exactly equal to $2$. Our classes contain the Thue--Morse--Mahler numbers, the sum of the reciprocals of the Fermat numbers, the regular paperfolding numbers, which have been previously considered respectively by Bugeaud, Coons, and Guo, Wu and Wen, but also new classes such as the Stern numbers and so on. Among other ingredients, our proofs use results on Hankel determinants obtained recently by Han.
- Feb 16 2015 math.NA arXiv:1502.03828v1In this paper, we develop a multiscale finite element method for solving flows in fractured media. Our approach is based on Generalized Multiscale Finite Element Method (GMsFEM), where we represent the fracture effects on a coarse grid via multiscale basis functions. These multiscale basis functions are constructed in the offline stage via local spectral problems following GMsFEM. To represent the fractures on the fine grid, we consider two approaches (1) Discrete Fracture Model (DFM) (2) Embedded Fracture Model (EFM) and their combination. In DFM, the fractures are resolved via the fine grid, while in EFM the fracture and the fine grid block interaction is represented as a source term. In the proposed multiscale method, additional multiscale basis functions are used to represent the long fractures, while short-size fractures are collectively represented by a single basis functions. The procedure is automatically done via local spectral problems. In this regard, our approach shares common concepts with several approaches proposed in the literature as we discuss. Numerical results are presented where we demonstrate how one can adaptively add basis functions in the regions of interest based on error indicators. We also discuss the use of randomized snapshots (\citerandomized2014) which reduces the offline computational cost.
- Jan 28 2015 math.PR arXiv:1501.06641v1Let $(\varepsilon_{t})_{t>0}$ be a sequence of independent real random vectors of $p$-dimension and let $X_T=\sum_{t=s+1}^{s+T}\varepsilon_t\varepsilon^T_{t-s}/T$ be the lag-$s$ ($s$ is a fixed positive integer) auto-covariance matrix of $\varepsilon_t$. This paper investigates the limiting behavior of the singular values of $X_T$ under the so-called \em ultra-dimensional regime where $p\to\infty$ and $T\to\infty$ in a related way such that $p/T\to 0$. First, we show that the singular value distribution of $X_T$ after a suitable normalization converges to a nonrandom limit $G$ (quarter law) under the forth-moment condition. Second, we establish the convergence of its largest singular value to the right edge of $G$. Both results are derived using the moment method.
- Oct 06 2014 math.PR arXiv:1410.0752v2Let $(\varepsilon_{t})_{t>0}$ be a sequence of independent real random vectors of $p$-dimension and let $X_T= \sum_{t=s+1}^{s+T}\varepsilon_t\varepsilon^T_{t-s}/T$ be the lag-$s$ ($s$ is a fixed positive integer) auto-covariance matrix of $\varepsilon_t$. Since $X_T$ is not symmetric, we consider its singular values, which are the square roots of the eigenvalues of $X_TX^T_T$. Therefore, the purpose of this paper is to investigate the limiting behaviors of the eigenvalues of $X_TX^T_T$ in two aspects. First, we show that the empirical spectral distribution of its eigenvalues converges to a nonrandom limit $F$. Second, we establish the convergence of its largest eigenvalue to the right edge of $F$. Both results are derived using moment methods.
- We rigorously show that a class of systems of partial differential equations modeling wave bifurcations supports stationary equivariant bifurcation dynamics through deriving its full dynamics on the center manifold(s). A direct consequence of our result is that the oscillations of the dynamics are \textitnot due to rotation waves though the system exhibits Euclidean symmetries. The main difficulties of carrying out the program are: 1) the system under study contains multi bifurcation parameters and we do not know \textita priori how they come into play in the bifurcation dynamics. 2) the representation of the linear operator on the center space is a $2\times 2$ zero matrix, which makes the characteristic condition in the well-known normal form theorem trivial. We overcome the first difficulty by using projection method. We managed to overcome the second subtle difficulty by using a conjugate pair coordinate for the center space and applying duality and projection arguments. Due to the specific complex pair parametrization, we could naturally get a form of the center manifold reduction function, which makes the study of the current dynamics on the center manifold possible. The symmetry of the system plays an essential role in excluding the possibility of bifurcating rotation waves.
- May 12 2014 math.NT arXiv:1405.2171v1In 1986, some examples of algebraic, and nonquadratic, power series over a fi?nite prime ?field, having a continued fraction expansion with partial quotients all of degree one, were discovered by W. Mills and D. Robbins. In this note we show how these few examples are included in a very large family of continued fractions for certain algebraic power series over an arbitrary ?finite fi?eld of odd characteristic.
- Random Fisher matrices arise naturally in multivariate statistical analysis and understanding the properties of its eigenvalues is of primary importance for many hypothesis testing problems like testing the equality between two multivariate population covariance matrices, or testing the independence between sub-groups of a multivariate random vector. This paper is concerned with the properties of a large-dimensional Fisher matrix when the dimension of the population is proportionally large compared to the sample size. Most of existing works on Fisher matrices deal with a particular Fisher matrix where populations have i.i.d components so that the population covariance matrices are all identity. In this paper, we consider general Fisher matrices with arbitrary population covariance matrices. The first main result of the paper establishes the limiting distribution of the eigenvalues of a Fisher matrix while in a second main result, we provide a central limit theorem for a wide class of functionals of its eigenvalues. Some applications of these results are also proposed for testing hypotheses on high-dimensional covariance matrices.
- Feb 26 2014 math.PR arXiv:1402.6064v3In this paper, we derive a joint central limit theorem for random vector whose components are function of random sesquilinear forms. This result is a natural extension of the existing central limit theory on random quadratic forms. We also provide applications in random matrix theory related to large-dimensional spiked population models. For the first application, we find the joint distribution of grouped extreme sample eigenvalues correspond to the spikes. And for the second application, under the assumption that the population covariance matrix is diagonal with $k$ (fixed) simple spikes, we derive the asymptotic joint distribution of the extreme sample eigenvalue and its corresponding sample eigenvector projection.
- Jan 11 2014 math.AP arXiv:1401.2197v1Extending work of Texier and Zumbrun in the semilinear non-re ection symmetric case, we study O(2) transverse Hopf bifurcation, or \cellular instability," of viscous shock waves in a channel, for a class of quasilinear hyperbolicparabolic systems including the equations of thermoviscoelasticity. The main diffi?culties are to (i) obtain Fr?echet diff?erentiability of the time-T solution operator by appropriate hyperbolicparabolic energy estimates, and (ii) handle O(2) symmetry in the absence of either center manifold reduction (due to lack of spectral gap) or (due to nonstandard quasilinear hyperbolic-parabolic form) the requisite framework for treatment by spatial dynamics on the space of time-periodic functions, the two standard treatments for this problem. The latter issue is resolved by LyapunovSchmidt reduction of the time-T map, yielding a four-dimensional problem with O(2) plus approximate S1 symmetry, which we treat \by hand" using direct Implicit Function Theorem arguments. The former is treated by balancing information obtained in Lagrangian coordinates with that from an augmented system. Interestingly, this argument does not apply to gas dynamics or magnetohydrodynamics (MHD), due to the in?nite-dimensional family of Lagrangian symmetries corresponding to invariance under arbitrary volume-preserving diff?eomorphisms.
- The equivariant Hopf bifurcation dynamics of a class of system of partial differential equations is carefully studied. The connections between the current dynamics and fundamental concepts in hyperbolic conservation laws are explained. The unique approximation property of center manifold reduction function is used in the current work to determine certain parameter in the normal form. The current work generalizes the study of the second author ([J. Yao, $O(2)$-Hopf bifurcation for a model of cellular shock instability, Physica D, 269 (2014), 63-75.]) and supplies a class of examples of $O(2)$ Hopf bifurcation with two parameters arising from systems of partial differential equations.
- Sep 17 2013 math.PR arXiv:1309.3728v4Consider a $N\times n$ matrix $\Sigma_n=\frac{1}{\sqrt{n}}R_n^{1/2}X_n$, where $R_n$ is a nonnegative definite Hermitian matrix and $X_n$ is a random matrix with i.i.d. real or complex standardized entries. The fluctuations of the linear statistics of the eigenvalues \[\operatorname Tracef \bigl(\Sigma_n\Sigma_n^*\bigr)=\sum_i=1^Nf(\lambda_i),\qquad (\lambda_i)\ eigenvalues\ of\ \Sigma_n\Sigma_n^*,\]are shown to be Gaussian, in the regime where both dimensions of matrix $\Sigma_n$ go to infinity at the same pace and in the case where $f$ is of class $C^3$, that is, has three continuous derivatives. The main improvements with respect to Bai and Silverstein's CLT [Ann. Probab. 32 (2004) 553-605] are twofold: First, we consider general entries with finite fourth moment, but whose fourth cumulant is nonnull, that is, whose fourth moment may differ from the moment of a (real or complex) Gaussian random variable. As a consequence, extra terms proportional to $ \vert \mathcal{V}\vert ^2=\bigl|\mathbb{E}\bigl(X_{11}^n\bigr) ^2\bigr|^2$ and $\kappa=\mathbb{E}\bigl \vert X_{11}^n\bigr \vert ^4-\vert {\mathcal{V}}\vert ^2-2$ appear in the limiting variance and in the limiting bias, which not only depend on the spectrum of matrix $R_n$ but also on its eigenvectors. Second, we relax the analyticity assumption over $f$ by representing the linear statistics with the help of Helffer-Sjöstrand's formula. The CLT is expressed in terms of vanishing Lévy-Prohorov distance between the linear statistics' distribution and a Gaussian probability distribution, the mean and the variance of which depend upon $N$ and $n$ and may not converge.
- In this paper, we develop new statistical theory for probabilistic principal component analysis models in high dimensions. The focus is the estimation of the noise variance, which is an important and unresolved issue when the number of variables is large in comparison with the sample size. We first unveil the reasons of a widely observed downward bias of the maximum likelihood estimator of the variance when the data dimension is high. We then propose a bias-corrected estimator using random matrix theory and establish its asymptotic normality. The superiority of the new (bias-corrected) estimator over existing alternatives is first checked by Monte-Carlo experiments with various combinations of $(p, n)$ (dimension and sample size). In order to demonstrate further potential benefits from the results of the paper to general probability PCA analysis, we provide evidence of net improvements in two popular procedures (Ulfarsson and Solo, 2008; Bai and Ng, 2002) for determining the number of principal components when the respective variance estimator proposed by these authors is replaced by the bias-corrected estimator. The new estimator is also used to derive new asymptotics for the related goodness-of-fit statistic under the high-dimensional scheme.
- This paper proposes a CLT for linear spectral statistics of random matrix $S^{-1}T$ for a general non-negative definite and \bf non-random Hermitian matrix $T$.
- Apr 24 2013 math.PR arXiv:1304.6164v2In this note, we establish an asymptotic expansion for the centering parameter appearing in the central limit theorems for linear spectral statistic of large-dimensional sample covariance matrices when the population has a spiked covariance structure. As an application, we provide an asymptotic power function for the corrected likelihood ratio statistic for testing the presence of spike eigenvalues in the population covariance matrix. This result generalizes an existing formula from the literature where only one simple spike exists.
- In this paper, we propose corrections to the likelihood ratio test and John's test for sphericity in large-dimensions. New formulas for the limiting parameters in the CLT for linear spectral statistics of sample covariance matrices with general fourth moments are first established. Using these formulas, we derive the asymptotic distribution of the two proposed test statistics under the null. These asymptotics are valid for general population, i.e. not necessarily Gaussian, provided a finite fourth-moment. Extensive Monte-Carlo experiments are conducted to assess the quality of these tests with a comparison to several existing methods from the literature. Moreover, we also obtain their asymptotic power functions under the alternative of a spiked population model as a specific alternative.
- Oct 22 2012 math.AP arXiv:1210.5308v1In this paper, we investigate nonlinear Schr$\ddot{o}$dinger type equations in $R^N$ under the framework of variable exponent spaces. We propose new assumptions on the nonlinear term to yield bounded Palais-Smale sequences and then prove the special sequences we find converge to critical points respectively. The main arguments are based on the geometry supplied by Fountain Theorem. Consequently, we show that the equation has a sequence of solutions with high energies.
- For a multivariate linear model, Wilk's likelihood ratio test (LRT) constitutes one of the cornerstone tools. However, the computation of its quantiles under the null or the alternative requires complex analytic approximations and more importantly, these distributional approximations are feasible only for moderate dimension of the dependent variable, say $p\le 20$. On the other hand, assuming that the data dimension $p$ as well as the number $q$ of regression variables are fixed while the sample size $n$ grows, several asymptotic approximations are proposed in the literature for Wilk's $\bLa$ including the widely used chi-square approximation. In this paper, we consider necessary modifications to Wilk's test in a high-dimensional context, specifically assuming a high data dimension $p$ and a large sample size $n$. Based on recent random matrix theory, the correction we propose to Wilk's test is asymptotically Gaussian under the null and simulations demonstrate that the corrected LRT has very satisfactory size and power, surely in the large $p$ and large $n$ context, but also for moderately large data dimensions like $p=30$ or $p=50$. As a byproduct, we give a reason explaining why the standard chi-square approximation fails for high-dimensional data. We also introduce a new procedure for the classical multiple sample significance test in MANOVA which is valid for high-dimensional data.
- May 31 2012 math.AP arXiv:1205.6991v2We evaluate by direct calculation the Lopatinski determinant for ZND detonations in Majda's model for reacting flow, and show that on the nonstable (nonnegative real part) complex half-plane it has a single zero at the origin of multiplicity one, implying stability. Together with results of Zumbrun on the inviscid limit, this recovers the result of RoqueJoffre-Vila that viscous detonations of Majda's model also are stable for sufficiently small viscosity, for any fixed detonation strength, heat release, and rate of reaction.
- In this paper we study the nonlinear Neumann boundary value problem of the following equations -\textdiv(|∇u|^p_1(x)-2∇u)-\textdiv(|∇u|^p_2(x)-2∇u)+|u|^p_1(x)-2u+|u|^p_2(x)-2u=\lambda f(x,u) in a bounded smooth domain $\Omega\subset\mathbb{R}^{N}$ with Neumann boundary condition given by |∇u|^p_1(x)-2\frac∂u∂\nu+|∇u|^p_2(x)-2\frac∂u∂\nu=\mu g(x,u) on $\partial\Omega$. Under appropriate conditions on the source and boundary nonlinearities, we obtain a number of results on existence and multiplicity of solutions by variational methods in the framework of variable exponent Lebesgue and Sobolev spaces.
- May 10 2012 math.FA arXiv:1205.1854v1In this paper, we investigate the $(p_{1}(x), p_{2}(x))$-Laplace operator, the properties of the corresponding integral functional and weak solutions to the related differential equations. We show that the integral functional admits a derivative of type $(S_+)$ which induces a homeomorphism between duality space pairs. As applications of the above results, we gave some existence results of the $(p_{1}(x), p_{2}(x))$-Laplace equation -\textdiv(|∇u|^p_1(x)-2∇u)-\textdiv(|∇u|^p_2(x)-2∇u)=f(x,u) in a bounded smooth domain $\Omega\subset\mathbb{R}^{N}$ with Dirichlet boundary condition.
- This paper deals with the problem of estimating the covariance matrix of a series of independent multivariate observations, in the case where the dimension of each observation is of the same order as the number of observations. Although such a regime is of interest for many current statistical signal processing and wireless communication issues, traditional methods fail to produce consistent estimators and only recently results relying on large random matrix theory have been unveiled. In this paper, we develop the parametric framework proposed by Mestre, and consider a model where the covariance matrix to be estimated has a (known) finite number of eigenvalues, each of it with an unknown multiplicity. The main contributions of this work are essentially threefold with respect to existing results, and in particular to Mestre's work: To relax the (restrictive) separability assumption, to provide joint consistent estimates for the eigenvalues and their multiplicities, and to study the variance error by means of a Central Limit theorem.
- Nov 01 2011 math.AP arXiv:1110.6897v2Extending investigations of Antman & Malek-Madani, Schecter & Shearer, Slemrod, Barker & Lewicka & Zumbrun, and others, we investigate phase-transitional elasticity models of strain-gradient effect. We prove the existence of non-constant planar periodic standing waves in these models with strain-gradient effects by variational methods and phase-plane analysis, for deformations of arbitrary dimension and general, physical, viscosity and strain-gradient terms. Previous investigations considered one-dimensional phenomenological models with artificial viscosity/strain gradient effect, for which the existence reduces to a standard (scalar) nonlinear oscillator. For our variational analysis, we require that the mean vector of the unknowns over one period be in the elliptic region with respect to the corresponding pure inviscid elastic model. For our (1-D) phase-plane analysis, we have no such restriction, obtaining essentially complete information on the existence of non-constant periodic waves and bounding homoclinic/heteroclinic waves. Our variational framework has implications also for time-evolutionary stability, through the link between the action functional for the traveling-wave ODE and the relative mechanical energy for the time-evolutionary system. Finally, we show that spectral implies modulational nonlinear stability by using a change of variables introduced by Kotschote to transform our system to a strictly parabolic system to which general results of Johnson--Zumbrun apply. Previous such results were confined to one-dimensional deformations in models with artificial viscosity--strain-gradient coefficients.
- Estimating the number of spikes in a spiked model is an important problem in many areas such as signal processing. Most of the classical approaches assume a large sample size $n$ whereas the dimension $p$ of the observations is kept small. In this paper, we consider the case of high dimension, where $p$ is large compared to $n$. The approach is based on recent results of random matrix theory. We extend our previous results to a more difficult situation where some spikes are equal, and compare our algorithm to an existing benchmark method.
- In this note we develop an extension of the Marčenko-Pastur theorem to time series model with temporal correlations. The limiting spectral distribution (LSD) of the sample covariance matrix is characterised by an explicit equation for its Stieltjes transform depending on the spectral density of the time series. A numerical algorithm is then given to compute the density functions of these LSD's.
- This article provides a central limit theorem for a consistent estimator of population eigenvalues with large multiplicities based on sample covariance matrices. The focus is on limited sample size situations, whereby the number of available observations is known and comparable in magnitude to the observation dimension. An exact expression as well as an empirical, asymptotically accurate, approximation of the limiting variance is derived. Simulations are performed that corroborate the theoretical claims. A specific application to wireless sensor networks is developed.
- In a spiked population model, the population covariance matrix has all its eigenvalues equal to units except for a few fixed eigenvalues (spikes). Determining the number of spikes is a fundamental problem which appears in many scientific fields, including signal processing (linear mixture model) or economics (factor model). Several recent papers studied the asymptotic behavior of the eigenvalues of the sample covariance matrix (sample eigenvalues) when the dimension of the observations and the sample size both grow to infinity so that their ratio converges to a positive constant. Using these results, we propose a new estimator based on the difference between two consecutive sample eigenvalues.
- This paper considers nonlinear regular-singular stochastic optimal control of large insurance company. The company controls the reinsurance rate and dividend payout process to maximize the expected present value of the dividend pay-outs until the time of bankruptcy. However, if the optimal dividend barrier is too low to be acceptable, it will make the company result in bankruptcy soon. Moreover, although risk and return should be highly correlated, over-risking is not a good recipe for high return, the supervisors of the company have to impose their preferred risk level and additional charge on firm seeking services beyond or lower than the preferred risk level. These indeed are nonlinear regular-singular stochastic optimal problems under ruin probability constraints. This paper aims at solving this kind of the optimal problems, that is, deriving the optimal retention ratio,dividend payout level, optimal return function and optimal control strategy of the insurance company. As a by-product, the paper also sets a risk-based capital standard to ensure the capital requirement of can cover the total given risk, and the effect of the risk level on optimal retention ratio, dividend payout level and optimal control strategy are also presented.
- Based on a point of view that solvency and security are first, this paper considers regular-singular stochastic optimal control problem of a large insurance company facing positive transaction cost asked by reinsurer under solvency constraint. The company controls proportional reinsurance and dividend pay-out policy to maximize the expected present value of the dividend pay-outs until the time of bankruptcy. The paper aims at deriving the optimal retention ratio, dividend payout level, explicit value function of the insurance company via stochastic analysis and PDE methods. The results present the best equilibrium point between maximization of dividend pay-outs and minimization of risks. The paper also gets a risk-based capital standard to ensure the capital requirement of can cover the total given risk. We present numerical results to make analysis how the model parameters, such as, volatility, premium rate, and risk level, impact on risk-based capital standard, optimal retention ratio, optimal dividend payout level and the company's profit.
- Oct 27 2009 math.GT arXiv:0910.4949v2Let $\imath: M\to \RR^{p+2}$ be a smooth embedding from a connected, oriented, closed $p$-dimesional smooth manifold to $\RR^{p+2}$, then there is a spin structure $\imath^\sharp(\varsigma^{p+2})$ on $M$ canonically induced from the embedding. If an orientation-preserving diffeomorphism $\tau$ of $M$ extends over $\imath$ as an orientation-preserving topological homeomorphism of $\RR^{p+2}$, then $\tau$ preserves the induced spin structure. Let $\esg_\cat(\imath)$ be the subgroup of the $\cat$-mapping class group $\mcg_\cat(M)$ consisting of elements whose representatives extend over $\RR^{p+2}$ as orientation-preserving $\cat$-homeomorphisms, where $\cat=\topo$, $\pl$ or $\diff$. The invariance of $\imath^\sharp(\varsigma^{p+2})$ gives nontrivial lower bounds to $[\mcg_\cat(M):\esg_\cat(\imath)]$ in various special cases. We apply this to embedded surfaces in $\RR^4$ and embedded $p$-dimensional tori in $\RR^{p+2}$. In particular, in these cases the index lower bounds for $\esg_\topo(\imath)$ are achieved for unknotted embeddings.
- In this paper, we give an explanation to the failure of two likelihood ratio procedures for testing about covariance matrices from Gaussian populations when the dimension is large compared to the sample size. Next, using recent central limit theorems for linear spectral statistics of sample covariance matrices and of random F-matrices, we propose necessary corrections for these LR tests to cope with high-dimensional effects. The asymptotic distributions of these corrected tests under the null are given. Simulations demonstrate that the corrected LR tests yield a realized size close to nominal level for both moderate p (around 20) and high dimension, while the traditional LR tests with chi-square approximation fails. Another contribution from the paper is that for testing the equality between two covariance matrices, the proposed correction applies equally for non-Gaussian populations yielding a valid pseudo-likelihood ratio test.
- If there exists a diffeomorphism $f$ on a closed, orientable $n$-manifold $M$ such that the non-wandering set $\Omega(f)$ consists of finitely many orientable $(\pm)$ attractors derived from expanding maps, then $M$ must be a rational homology sphere; moreover all those attractors are of topological dimension $n-2$. Expanding maps are expanding on (co)homologies.
- In this paper we consider the realization of DE attractors by self-diffeomorphisms of manifolds. For any expanding self-map $\phi:M\to M$ of a connected, closed $p$-dimensional manifold $M$, one can always realize a $(p,q)$-type attractor derived from $\phi$ by a compactly-supported self-diffeomorphsm of $\RR^{p+q}$, as long as $q\geq p+1$. Thus lower codimensional realizations are more interesting, related to the knotting problem below the stable range. We show that for any expanding self-map $\phi$ of a standard smooth $p$-dimensional torus $T^p$, there is compactly-supported self-diffeomorphism of $\RR^{p+2}$ realizing an attractor derived from $\phi$. A key ingredient of the construction is to understand automorphisms of $T^p$ which extend over $\RR^{p+2}$ as a self-diffeomorphism via the standard unknotted embedding $\imath_p:T^p\hookrightarrow\RR^{p+2}$. We show that these automorphisms form a subgroup $E_{\imath_p}$ of $\Aut(T^p)$ of index at most $2^p-1$.
- Jun 17 2008 math.PR arXiv:0806.2503v1In a spiked population model, the population covariance matrix has all its eigenvalues equal to units except for a few fixed eigenvalues (spikes). This model is proposed by Johnstone to cope with empirical findings on various data sets. The question is to quantify the effect of the perturbation caused by the spike eigenvalues. A recent work by Baik and Silverstein establishes the almost sure limits of the extreme sample eigenvalues associated to the spike eigenvalues when the population and the sample sizes become large. This paper establishes the limiting distributions of these extreme sample eigenvalues. As another important result of the paper, we provide a central limit theorem on random sesquilinear forms.
- In the spiked population model introduced by Johnstone (2001),the population covariance matrix has all its eigenvalues equal to unit except for a few fixed eigenvalues (spikes). The question is to quantify the effect of the perturbation caused by the spike eigenvalues. Baik and Silverstein (2006) establishes the almost sure limits of the extreme sample eigenvalues associated to the spike eigenvalues when the population and the sample sizes become large. In a recent work (Bai and Yao, 2008), we have provided the limiting distributions for these extreme sample eigenvalues. In this paper, we extend this theory to a \em generalized spiked population model where the base population covariance matrix is arbitrary, instead of the identity matrix as in Johnstone's case. New mathematical tools are introduced for establishing the almost sure convergence of the sample eigenvalues generated by the spikes.
- In several application fields like daily pluviometry data modelling, or motion analysis from image sequences, observations contain two components of different nature. A first part is made with discrete values accounting for some symbolic information and a second part records a continuous (real-valued) measurement. We call such type of observations "mixed-state observations". This paper introduces spatial models suited for the analysis of these kinds of data. We consider multi-parameter auto-models whose local conditional distributions belong to a mixed state exponential family. Specific examples with exponential distributions are detailed, and we present some experimental results for modelling motion measurements from video sequences.
- This paper will be splited into two papers and submited later.
- For each composite number $n\ne 2^k$, there does not exist a single connected closed $(n+1)$-manifold such that any smooth, simply-connected, closed $n$-manifold can be topologically flat embedded into it. There is a single connected closed 5-manifold $W$ such that any simply-connected, 4-manifold $M$ can be topologically flat embedded into $W$ if $M$ is either closed and indefinite, or compact and with non-empty boundary.
- Mar 25 2005 math.PR arXiv:math/0503527v1Let Y be an Ornstein-Uhlenbeck diffusion governed by a stationary and ergodic Markov jump process X: dY_t=a(X_t)Y_t dt+\sigma(X_t) dW_t, Y_0=y_0. Ergodicity conditions for Y have been obtained. Here we investigate the tail propriety of the stationary distribution of this model. A characterization of either heavy or light tail case is established. The method is based on a renewal theorem for systems of equations with distributions on R.