亚洲男人的天堂2018av,欧美草比,久久久久久免费视频精选,国色天香在线看免费,久久久久亚洲av成人片仓井空

We consider a prototypical problem of Bayesian inference for a structured spiked model: a low-rank signal is corrupted by additive noise. While both information-theoretic and algorithmic limits are well understood when the noise is a Gaussian Wigner matrix, the more realistic case of structured noise still proves to be challenging. To capture the structure while maintaining mathematical tractability, a line of work has focused on rotationally invariant noise. However, existing studies either provide sub-optimal algorithms or are limited to special cases of noise ensembles. In this paper, using tools from statistical physics (replica method) and random matrix theory (generalized spherical integrals) we establish the first characterization of the information-theoretic limits for a noise matrix drawn from a general trace ensemble. Remarkably, our analysis unveils the asymptotic equivalence between the rotationally invariant model and a surrogate Gaussian one. Finally, we show how to saturate the predicted statistical limits using an efficient algorithm inspired by the theory of adaptive Thouless-Anderson-Palmer (TAP) equations.

相關內容

Ordinary differential equations (ODE) are a popular tool to model the spread of infectious diseases, yet they implicitly assume an exponential distribution to describe the flow from one infection state to another. However, scientific experience yields more plausible distributions where the likelihood of disease progression or recovery changes accordingly with the duration spent in a particular state of the disease. Furthermore, transmission dynamics depend heavily on the infectiousness of individuals. The corresponding nonlinear variation with the time individuals have already spent in an infectious state requires more realistic models. The previously mentioned items are particularly crucial when modeling dynamics at change points such as the implementation of nonpharmaceutical interventions. In order to capture these aspects and to enhance the accuracy of simulations, integro-differential equations (IDE) can be used. In this paper, we propose a generalized model based on integro-differential equations with eight infection states. The model allows for variable stay time distributions and generalizes the concept of ODE-based models as well as IDE-based age-of-infection models. In this, we include particular infection states for severe and critical cases to allow for surveillance of the clinical sector, avoiding bottlenecks and overloads in critical epidemic situations. We will extend a recently introduced nonstandard numerical scheme to solve a simpler IDE-based model. This scheme is adapted to our more advanced model and we prove important mathematical and biological properties for the numerical solution. Furthermore, we validate our approach numerically by demonstrating the convergence rate. Eventually, we also show that our novel model is intrinsically capable of better assessing disease dynamics upon the introduction of nonpharmaceutical interventions.

Zeitlin's model is a spatial discretization for the 2-D Euler equations on the flat 2-torus or the 2-sphere. Contrary to other discretizations, it preserves the underlying geometric structure, namely that the Euler equations describe Riemannian geodesics on a Lie group. Here we show how to extend Zeitlin's approach to the axisymmetric Euler equations on the 3-sphere. It is the first discretization of the 3-D Euler equations that fully preserves the geometric structure. Thus, this finite-dimensional model admits Riemannian curvature and Jacobi equations, which are discussed.

Data depth has emerged as an invaluable nonparametric measure for the ranking of multivariate samples. The main contribution of depth-based two-sample comparisons is the introduction of the Q statistic (Liu and Singh, 1993), a quality index. Unlike traditional methods, data depth does not require the assumption of normal distributions and adheres to four fundamental properties. Many existing two-sample homogeneity tests, which assess mean and/or scale changes in distributions often suffer from low statistical power or indeterminate asymptotic distributions. To overcome these challenges, we introduced a DEEPEAST (depth-explored same-attraction sample-to-sample central-outward ranking) technique for improving statistical power in two-sample tests via the same-attraction function. We proposed two novel and powerful depth-based test statistics: the sum test statistic and the product test statistic, which are rooted in Q statistics, share a "common attractor" and are applicable across all depth functions. We further proved the asymptotic distribution of these statistics for various depth functions. To assess the performance of power gain, we apply three depth functions: Mahalanobis depth (Liu and Singh, 1993), Spatial depth (Brown, 1958; Gower, 1974), and Projection depth (Liu, 1992). Through two-sample simulations, we have demonstrated that our sum and product statistics exhibit superior power performance, utilizing a strategic block permutation algorithm and compare favourably with popular methods in literature. Our tests are further validated through analysis on Raman spectral data, acquired from cellular and tissue samples, highlighting the effectiveness of the proposed tests highlighting the effective discrimination between health and cancerous samples.

We show that Pfaffians or contiguity relations of hypergeometric functions of several variables give a direct sampling algorithm from toric models in statistics, which is a Markov chain on a lattice generated by a matrix $A$. A correspondence among graphical models and $A$-hypergeometric system is discussed and we give a sum formula of special values of $A$-hypergeometric polynomials. Some hypergeometric series which are interesting in view of statistics are presented.

Decentralized optimization has become a standard paradigm for solving large-scale decision-making problems and training large machine learning models without centralizing data. However, this paradigm introduces new privacy and security risks, with malicious agents potentially able to infer private data or impair the model accuracy. Over the past decade, significant advancements have been made in developing secure decentralized optimization and learning frameworks and algorithms. This survey provides a comprehensive tutorial on these advancements. We begin with the fundamentals of decentralized optimization and learning, highlighting centralized aggregation and distributed consensus as key modules exposed to security risks in federated and distributed optimization, respectively. Next, we focus on privacy-preserving algorithms, detailing three cryptographic tools and their integration into decentralized optimization and learning systems. Additionally, we examine resilient algorithms, exploring the design and analysis of resilient aggregation and consensus protocols that support these systems. We conclude the survey by discussing current trends and potential future directions.

In this paper, we present a class of high-order and efficient compact difference schemes for nonlinear convection diffusion equations, which can preserve both bounds and mass. For the one-dimensional problem, we first introduce a high-order compact Strang splitting scheme (denoted as HOC-Splitting), which is fourth-order accurate in space and second-order accurate in time. Then, by incorporating the Lagrange multiplier approach with the HOC-Splitting scheme, we construct two additional bound-preserving or/and mass-conservative HOC-Splitting schemes that do not require excessive computational cost and can automatically ensure the uniform bounds of the numerical solution. These schemes combined with an alternating direction implicit (ADI) method are generalized to the two-dimensional problem, which further enhance the computational efficiency for large-scale modeling and simulation. Besides, we present an optimal-order error estimate for the bound-preserving ADI scheme in the discrete $L_2$ norm. Finally, ample numerical examples are presented to verify the theoretical results and demonstrate the accuracy, efficiency, and effectiveness in preserving bounds or/and mass of the proposed schemes.

We prove sharp wavenumber-explicit error bounds for first- or second-type-N\'ed\'elec-element (a.k.a. edge-element) conforming discretisations, of arbitrary (fixed) order, of the variable-coefficient time-harmonic Maxwell equations posed in a bounded domain with perfect electric conductor (PEC) boundary conditions. The PDE coefficients are allowed to be piecewise regular and complex-valued; this set-up therefore includes scattering from a PEC obstacle and/or variable real-valued coefficients, with the radiation condition approximated by a perfectly matched layer (PML). In the analysis of the $h$-version of the finite-element method, with fixed polynomial degree $p$, applied to the time-harmonic Maxwell equations, the $\textit{asymptotic regime}$ is when the meshwidth, $h$, is small enough (in a wavenumber-dependent way) that the Galerkin solution is quasioptimal independently of the wavenumber, while the $\textit{preasymptotic regime}$ is the complement of the asymptotic regime. The results of this paper are the first preasymptotic error bounds for the time-harmonic Maxwell equations using first-type N\'ed\'elec elements or higher-than-lowest-order second-type N\'ed\'elec elements. Furthermore, they are the first wavenumber-explicit results, even in the asymptotic regime, for Maxwell scattering problems with a non-empty scatterer.

This work proposes a novel numerical scheme for solving the high-dimensional Hamilton-Jacobi-Bellman equation with a functional hierarchical tensor ansatz. We consider the setting of stochastic control, whereby one applies control to a particle under Brownian motion. In particular, the existence of diffusion presents a new challenge to conventional tensor network methods for deterministic optimal control. To overcome the difficulty, we use a general regression-based formulation where the loss term is the Bellman consistency error combined with a Sobolev-type penalization term. We propose two novel sketching-based subroutines for obtaining the tensor-network approximation to the action-value functions and the value functions, which greatly accelerate the convergence for the subsequent regression phase. We apply the proposed approach successfully to two challenging control problems with Ginzburg-Landau potential in 1D and 2D with 64 variables.

Given a finite set of matrices with integer entries, the matrix mortality problem asks if there exists a product of these matrices equal to the zero matrix. We consider a special case of this problem where all entries of the matrices are nonnegative. This case is equivalent to the NFA mortality problem, which, given an NFA, asks for a word $w$ such that the image of every state under $w$ is the empty set. The size of the alphabet of the NFA is then equal to the number of matrices in the set. We study the length of shortest such words depending on the size of the alphabet. We show that for an NFA with $n$ states this length can be at least $2^n - 1$ for an alphabet of size $n$, $2^{(n - 4)/2}$ for an alphabet of size $3$ and $2^{(n - 2)/3}$ for an alphabet of size $2$. We also discuss further open problems related to mortality of NFAs and DFAs.

The strong convergence of an explicit full-discrete scheme is investigated for the stochastic Burgers-Huxley equation driven by additive space-time white noise, which possesses both Burgers-type and cubic nonlinearities. To discretize the continuous problem in space, we utilize a spectral Galerkin method. Subsequently, we introduce a nonlinear-tamed exponential integrator scheme, resulting in a fully discrete scheme. Within the framework of semigroup theory, this study provides precise estimations of the Sobolev regularity, $L^\infty$ regularity in space, and H\"older continuity in time for the mild solution, as well as for its semi-discrete and full-discrete approximations. Building upon these results, we establish moment boundedness for the numerical solution and obtain strong convergence rates in both spatial and temporal dimensions. A numerical example is presented to validate the theoretical findings.

北京阿比特科技有限公司