A neural network with a feed forward topology and Bayesian regularization training algorithm is used to predict the austenite formation temperatures (At1 and A13) by considering the percentage of alloying elements i...A neural network with a feed forward topology and Bayesian regularization training algorithm is used to predict the austenite formation temperatures (At1 and A13) by considering the percentage of alloying elements in chemical composition of steel. The data base used here involves a large variety of different steel types such as struc- tural steels, stainless steels, rail steels, spring steels, high temperature creep resisting steels and tool steels. Scatter diagrams and mean relative error (MRE) statistical criteria are used to compare the performance of developed neural network with the results of Andrew% empirical equations and a feed forward neural network with "gradient descent with momentum" training algorithm. The results showed that Bayesian regularization neural network has the best performance. Also, due to the satisfactory results of the developed neural network, it was used to investigate the effect of the chemical composition on Ac1 and At3 temperatures. Results are in accordance with materials science theories.展开更多
For geophysical inversion problems,deterministic inversion methods can easily fall into local optimal solutions,while stochastic optimization methods can theoretically converge to global optimal solutions.These proble...For geophysical inversion problems,deterministic inversion methods can easily fall into local optimal solutions,while stochastic optimization methods can theoretically converge to global optimal solutions.These problems have always been a concern for researchers.Among many stochastic optimization methods,particle swarm optimization(PSO)has been applied to solve geophysical inversion problems due to its simple principle and the fact that only a few parameters require adjustment.To overcome the nonuniqueness of inversion,model constraints can be added to PSO optimization.However,using fixed regularization parameters in PSO iteration is equivalent to keeping the default model constraint at a certain level,yielding an inversion result that is considerably affected by the model constraint.This study proposes a hybrid method that combines the regularized least squares method(RLSM)with the PSO method.The RLSM is used to improve the global optimal particle and accelerate convergence,while the adaptive regularization strategy is used to update the regularization parameters to avoid the influence of model constraints on the inversion results.Further,the inversion results of the RLSM and hybrid algorithm are compared and analyzed by considering the audio magnetotelluric synthesis and field data as examples.Experiments show that the proposed hybrid method is superior to the RLSM.Furthermore,compared with the standard PSO algorithm,the hybrid algorithm needs a broader model space but a smaller particle swarm and fewer iteration steps,thus reducing the prior conditions and the computational cost used in the inversion.展开更多
Bioluminescence tomography(BLT)is an important noninvasive optical molecular imaging modality in preclinical research.To improve the image quality,reconstruction algorithms have to deal with the inherent ill-posedness...Bioluminescence tomography(BLT)is an important noninvasive optical molecular imaging modality in preclinical research.To improve the image quality,reconstruction algorithms have to deal with the inherent ill-posedness of BLT inverse problem.The sparse characteristic of bioluminescent sources in spatial distribution has been widely explored in BLT and many L1-regularized methods have been investigated due to the sparsity-inducing properties of L1 norm.In this paper,we present a reconstruction method based on L_(1/2) regularization to enhance sparsity of BLT solution and solve the nonconvex L_(1/2) norm problem by converting it to a series of weighted L1 homotopy minimization problems with iteratively updated weights.To assess the performance of the proposed reconstruction algorithm,simulations on a heterogeneous mouse model are designed to compare it with three representative sparse reconstruction algorithms,including the weighted interior-point,L1 homotopy,and the Stagewise Orthogonal Matching Pursuit algorithm.Simulation results show that the proposed method yield stable reconstruction results under different noise levels.Quantitative comparison results demonstrate that the proposed algorithm outperforms the competitor algorithms in location accuracy,multiple-source resolving and image quality.展开更多
Tomographic synthetic aperture radar(TomoSAR)imaging exploits the antenna array measurements taken at different elevation aperture to recover the reflectivity function along the elevation direction.In these years,for ...Tomographic synthetic aperture radar(TomoSAR)imaging exploits the antenna array measurements taken at different elevation aperture to recover the reflectivity function along the elevation direction.In these years,for the sparse elevation distribution,compressive sensing(CS)is a developed favorable technique for the high-resolution elevation reconstruction in TomoSAR by solving an L_(1) regularization problem.However,because the elevation distribution in the forested area is nonsparse,if we want to use CS in the recovery,some basis,such as wavelet,should be exploited in the sparse L_(1/2) representation of the elevation reflectivity function.This paper presents a novel wavelet-based L_(2) regularization CS-TomoSAR imaging method of the forested area.In the proposed method,we first construct a wavelet basis,which can sparsely represent the elevation reflectivity function of the forested area,and then reconstruct the elevation distribution by using the L_(1/2) regularization technique.Compared to the wavelet-based L_(1) regularization TomoSAR imaging,the proposed method can improve the elevation recovered quality efficiently.展开更多
Bayesian empirical likelihood is a semiparametric method that combines parametric priors and nonparametric likelihoods, that is, replacing the parametric likelihood function in Bayes theorem with a nonparametric empir...Bayesian empirical likelihood is a semiparametric method that combines parametric priors and nonparametric likelihoods, that is, replacing the parametric likelihood function in Bayes theorem with a nonparametric empirical likelihood function, which can be used without assuming the distribution of the data. It can effectively avoid the problems caused by the wrong setting of the model. In the variable selection based on Bayesian empirical likelihood, the penalty term is introduced into the model in the form of parameter prior. In this paper, we propose a novel variable selection method, L<sub>1/2</sub> regularization based on Bayesian empirical likelihood. The L<sub>1/2</sub> penalty is introduced into the model through a scale mixture of uniform representation of generalized Gaussian prior, and the posterior distribution is then sampled using MCMC method. Simulations demonstrate that the proposed method can have better predictive ability when the error violates the zero-mean normality assumption of the standard parameter model, and can perform variable selection.展开更多
In view of the composition analysis and identification of ancient glass products, L1 regularization, K-Means cluster analysis, elbow rule and other methods were comprehensively used to build logical regression, cluste...In view of the composition analysis and identification of ancient glass products, L1 regularization, K-Means cluster analysis, elbow rule and other methods were comprehensively used to build logical regression, cluster analysis, hyper-parameter test and other models, and SPSS, Python and other tools were used to obtain the classification rules of glass products under different fluxes, sub classification under different chemical compositions, hyper-parameter K value test and rationality analysis. Research can provide theoretical support for the protection and restoration of ancient glass relics.展开更多
Seismic data regularization is an important preprocessing step in seismic signal processing. Traditional seismic acquisition methods follow the Shannon–Nyquist sampling theorem, whereas compressive sensing(CS) prov...Seismic data regularization is an important preprocessing step in seismic signal processing. Traditional seismic acquisition methods follow the Shannon–Nyquist sampling theorem, whereas compressive sensing(CS) provides a fundamentally new paradigm to overcome limitations in data acquisition. Besides the sparse representation of seismic signal in some transform domain and the 1-norm reconstruction algorithm, the seismic data regularization quality of CS-based techniques strongly depends on random undersampling schemes. For 2D seismic data, discrete uniform-based methods have been investigated, where some seismic traces are randomly sampled with an equal probability. However, in theory and practice, some seismic traces with different probability are required to be sampled for satisfying the assumptions in CS. Therefore, designing new undersampling schemes is imperative. We propose a Bernoulli-based random undersampling scheme and its jittered version to determine the regular traces that are randomly sampled with different probability, while both schemes comply with the Bernoulli process distribution. We performed experiments using the Fourier and curvelet transforms and the spectral projected gradient reconstruction algorithm for 1-norm(SPGL1), and ten different random seeds. According to the signal-to-noise ratio(SNR) between the original and reconstructed seismic data, the detailed experimental results from 2D numerical and physical simulation data show that the proposed novel schemes perform overall better than the discrete uniform schemes.展开更多
By defining fuzzy valued simple functions and giving L1(μ) approximations of fuzzy valued integrably bounded functions by such simple functions, the paper analyses by L1(μ)-norm the approximation capability of four-...By defining fuzzy valued simple functions and giving L1(μ) approximations of fuzzy valued integrably bounded functions by such simple functions, the paper analyses by L1(μ)-norm the approximation capability of four-layer feedforward regular fuzzy neural networks to the fuzzy valued integrably bounded function F : Rn → FcO(R). That is, if the transfer functionσ: R→R is non-polynomial and integrable function on each finite interval, F may be innorm approximated by fuzzy valued functions defined as to anydegree of accuracy. Finally some real examples demonstrate the conclusions.展开更多
We propose an ?~1 regularized method for numerical differentiation using empirical eigenfunctions. Compared with traditional methods for numerical differentiation, the output of our method can be considered directly ...We propose an ?~1 regularized method for numerical differentiation using empirical eigenfunctions. Compared with traditional methods for numerical differentiation, the output of our method can be considered directly as the derivative of the underlying function. Moreover,our method could produce sparse representations with respect to empirical eigenfunctions.Numerical results show that our method is quite effective.展开更多
L(d, 1)-labeling is a kind of graph coloring problem from frequency assignment in radio networks, in which adjacent nodes must receive colors that are at least d apart while nodes at distance two from each other must ...L(d, 1)-labeling is a kind of graph coloring problem from frequency assignment in radio networks, in which adjacent nodes must receive colors that are at least d apart while nodes at distance two from each other must receive different colors. We focus on L(d, 1)-labeling of regular tilings for d≥3 since the cases d=0, 1 or 2 have been researched by Calamoneri and Petreschi. For all three kinds of regular tilings, we give their L (d, 1)-labeling numbers for any integer d≥3. Therefore, combined with the results given by Calamoneri and Petreschi, the L(d, 1)-labeling numbers of regular tilings for any nonnegative integer d may be determined completely.展开更多
The paper discusses the core parameters of the 3 D and 4 D variational merging based on L1 norm regularization,namely optimization characteristic correlation length of background error covariance matrix and regulariza...The paper discusses the core parameters of the 3 D and 4 D variational merging based on L1 norm regularization,namely optimization characteristic correlation length of background error covariance matrix and regularization parameter. Classical 3 D/4 D variational merging is based on the theory that error follows Gaussian distribution. It involves the solution of the objective functional gradient in minimization iteration,which requires the data to have continuity and differentiability. Classic 3 D/4 D-dimensional variational merging method was extended,and L1 norm was used as the constraint coupling to the classical variational merged model. Experiment was carried out by using linear advection-diffusion equation as four-dimensional prediction model,and parameter optimization of this method is discussed. Considering the strong temporal and spatial variation of water vapor,this method is further applied to the precipitable water vapor( PWV) merging by calculating reanalysis data and GNSS retrieval.Parameters were adjusted gradually to analyze the influence of background field on the merging result,and the experiment results show that the mathematical algorithm adopted in this paper is feasible.展开更多
The properties of generalized flip Markov chains on connected regular digraphs are discussed.The 1-Flipper operation on Markov chains for undirected graphs is generalized to that for multi-digraphs.The generalized 1-F...The properties of generalized flip Markov chains on connected regular digraphs are discussed.The 1-Flipper operation on Markov chains for undirected graphs is generalized to that for multi-digraphs.The generalized 1-Flipper operation preserves the regularity and weak connectivity of multi-digraphs.The generalized 1-Flipper operation is proved to be symmetric.Moreover,it is presented that a series of random generalized 1-Flipper operations eventually lead to a uniform probability distribution over all connected d-regular multi-digraphs without loops.展开更多
Let us consider the following elliptic systems of second order-D_α(A_i~α(x, u, Du))=B_4(x, u, Du), i=1, …, N, x∈Q(?)R^n, n≥3 (1) and supposeⅰ) |A_i~α(x, u, Du)|≤L(1+|Du|);ⅱ) (1+|p|)^(-1)A_i~α(x, u, p)are H(?...Let us consider the following elliptic systems of second order-D_α(A_i~α(x, u, Du))=B_4(x, u, Du), i=1, …, N, x∈Q(?)R^n, n≥3 (1) and supposeⅰ) |A_i~α(x, u, Du)|≤L(1+|Du|);ⅱ) (1+|p|)^(-1)A_i~α(x, u, p)are H(?)lder-continuous functions with some exponent δ on (?)×R^N uniformly with respect to p, i.e.ⅲ) A_i~α(x, u, p) are differentiable function in p with bounded and continuous derivativesⅳ)ⅴ) for all u∈H_(loc)~1(Ω, R^N)∩L^(n(γ-1)/(2-γ))(Ω, R^N), B(x, u, Du)is ineasurable and |B(x, u, p)|≤a(|p|~γ+|u|~τ)+b(x), where 1+2/n<γ<2, τ≤max((n+2)/(n-2), (γ-1)/(2-γ)-ε), (?)ε>0, b(x)∈L2n/(n+2), n^2/(n+2)+e(Ω), (?)ε>0.Remarks. Only bounded open set Q will be considered in this paper; for all p≥1, λ≥0, which is clled a Morrey Space.Let assumptions ⅰ)-ⅳ) hold, Giaquinta and Modica have proved the regularity of both the H^1 weak solutions of (1) under controllable growth condition |B|≤α(|p|~γ+|u|^((n+2)/(n-2))+b, 0<γ≤1+2/n and the H^1∩L~∞ weak solutions of (1) under natural growth condition |B|≤α|p|~2+b with a smallness condition 2aM<λ(|u|≤M), which implys that the H^1∩L~∞ weak solutions have the same regularty in the case of 1+2/n<γ<2. In the case of γ=2, many counterexamples (see [2] showed that u must be in H^1L~∞, while in the case of 1+2/n<γ<2, we consider the H^1∩L^n(γ-1)/(2-γ) weak solutions of (1), weaken the instability conditions upon them (from L~∞ to L^n(γ-1)/(2-γ) and obtain the same regularity results. Finally we show that the exponent n(γ-1)/(2-γ) can not be docreased anymore for the sake of the regularity results.Delinition 1. We call u∈H^1∩L^n(γ-1)/(2-γ)(Q, R^N) be a weak solution of (1), providod that where We use the convention that repeated indices are summed. i, j go from 1 to N ann α, β from 1 to n.展开更多
On a compact Riemann surface with finite punctures P_(1),…P_(k),we define toric curves as multivalued,totallyunramified holomorphic maps to P^(n)with monodromy in a maximal torus of PSU(n+1).Toric solutions to SU(n+1...On a compact Riemann surface with finite punctures P_(1),…P_(k),we define toric curves as multivalued,totallyunramified holomorphic maps to P^(n)with monodromy in a maximal torus of PSU(n+1).Toric solutions to SU(n+1)Todasystems on X\{P_(1);…;P_(k)}are recognized by the associated toric curves in.We introduce character n-ensembles as-tuples of meromorphic one-forms with simple poles and purely imaginary periods,generating toric curves on minus finitelymany points.On X,we establish a correspondence between character-ensembles and toric solutions to the SU(n+1)system with finitely many cone singularities.Our approach not only broadens seminal solutions with two conesingularities on the Riemann sphere,as classified by Jost-Wang(Int.Math.Res.Not.,2002,(6):277-290)andLin-Wei-Ye(Invent.Math.,2012,190(1):169-207),but also advances beyond the limits of Lin-Yang-Zhong’s existencetheorems(J.Differential Geom.,2020,114(2):337-391)by introducing a new solution class.展开更多
In this paper, we built upon the estimating primaries by sparse inversion (EPSI) method. We use the 3D curvelet transform and modify the EPSI method to the sparse inversion of the biconvex optimization and Ll-norm r...In this paper, we built upon the estimating primaries by sparse inversion (EPSI) method. We use the 3D curvelet transform and modify the EPSI method to the sparse inversion of the biconvex optimization and Ll-norm regularization, and use alternating optimization to directly estimate the primary reflection coefficients and source wavelet. The 3D curvelet transform is used as a sparseness constraint when inverting the primary reflection coefficients, which results in avoiding the prediction subtraction process in the surface-related multiples elimination (SRME) method. The proposed method not only reduces the damage to the effective waves but also improves the elimination of multiples. It is also a wave equation- based method for elimination of surface multiple reflections, which effectively removes surface multiples under complex submarine conditions.展开更多
基金Financial support of mechanical engineering center of excellence at Roudbar Azad University
文摘A neural network with a feed forward topology and Bayesian regularization training algorithm is used to predict the austenite formation temperatures (At1 and A13) by considering the percentage of alloying elements in chemical composition of steel. The data base used here involves a large variety of different steel types such as struc- tural steels, stainless steels, rail steels, spring steels, high temperature creep resisting steels and tool steels. Scatter diagrams and mean relative error (MRE) statistical criteria are used to compare the performance of developed neural network with the results of Andrew% empirical equations and a feed forward neural network with "gradient descent with momentum" training algorithm. The results showed that Bayesian regularization neural network has the best performance. Also, due to the satisfactory results of the developed neural network, it was used to investigate the effect of the chemical composition on Ac1 and At3 temperatures. Results are in accordance with materials science theories.
基金supported by the National Natural Science Foundation of China(NSFC)[grant number 41374133]
文摘For geophysical inversion problems,deterministic inversion methods can easily fall into local optimal solutions,while stochastic optimization methods can theoretically converge to global optimal solutions.These problems have always been a concern for researchers.Among many stochastic optimization methods,particle swarm optimization(PSO)has been applied to solve geophysical inversion problems due to its simple principle and the fact that only a few parameters require adjustment.To overcome the nonuniqueness of inversion,model constraints can be added to PSO optimization.However,using fixed regularization parameters in PSO iteration is equivalent to keeping the default model constraint at a certain level,yielding an inversion result that is considerably affected by the model constraint.This study proposes a hybrid method that combines the regularized least squares method(RLSM)with the PSO method.The RLSM is used to improve the global optimal particle and accelerate convergence,while the adaptive regularization strategy is used to update the regularization parameters to avoid the influence of model constraints on the inversion results.Further,the inversion results of the RLSM and hybrid algorithm are compared and analyzed by considering the audio magnetotelluric synthesis and field data as examples.Experiments show that the proposed hybrid method is superior to the RLSM.Furthermore,compared with the standard PSO algorithm,the hybrid algorithm needs a broader model space but a smaller particle swarm and fewer iteration steps,thus reducing the prior conditions and the computational cost used in the inversion.
基金supported by the National Natural Science Foundation of China(No.61401264,11574192)the Natural Science Research Plan Program in Shaanxi Province of China(No.2015JM6322)the Fundamental Research Funds for the Central Universities(No.GK201603025).
文摘Bioluminescence tomography(BLT)is an important noninvasive optical molecular imaging modality in preclinical research.To improve the image quality,reconstruction algorithms have to deal with the inherent ill-posedness of BLT inverse problem.The sparse characteristic of bioluminescent sources in spatial distribution has been widely explored in BLT and many L1-regularized methods have been investigated due to the sparsity-inducing properties of L1 norm.In this paper,we present a reconstruction method based on L_(1/2) regularization to enhance sparsity of BLT solution and solve the nonconvex L_(1/2) norm problem by converting it to a series of weighted L1 homotopy minimization problems with iteratively updated weights.To assess the performance of the proposed reconstruction algorithm,simulations on a heterogeneous mouse model are designed to compare it with three representative sparse reconstruction algorithms,including the weighted interior-point,L1 homotopy,and the Stagewise Orthogonal Matching Pursuit algorithm.Simulation results show that the proposed method yield stable reconstruction results under different noise levels.Quantitative comparison results demonstrate that the proposed algorithm outperforms the competitor algorithms in location accuracy,multiple-source resolving and image quality.
基金This work was supported by the Fundamental Research Funds for the Central Universities(NE2020004)the National Natural Science Foundation of China(61901213)+3 种基金the Natural Science Foundation of Jiangsu Province(BK20190397)the Aeronautical Science Foundation of China(201920052001)the Young Science and Technology Talent Support Project of Jiangsu Science and Technology Associationthe Foundation of Graduate Innovation Center in Nanjing University of Aeronautics and Astronautics(kfjj20200419).
文摘Tomographic synthetic aperture radar(TomoSAR)imaging exploits the antenna array measurements taken at different elevation aperture to recover the reflectivity function along the elevation direction.In these years,for the sparse elevation distribution,compressive sensing(CS)is a developed favorable technique for the high-resolution elevation reconstruction in TomoSAR by solving an L_(1) regularization problem.However,because the elevation distribution in the forested area is nonsparse,if we want to use CS in the recovery,some basis,such as wavelet,should be exploited in the sparse L_(1/2) representation of the elevation reflectivity function.This paper presents a novel wavelet-based L_(2) regularization CS-TomoSAR imaging method of the forested area.In the proposed method,we first construct a wavelet basis,which can sparsely represent the elevation reflectivity function of the forested area,and then reconstruct the elevation distribution by using the L_(1/2) regularization technique.Compared to the wavelet-based L_(1) regularization TomoSAR imaging,the proposed method can improve the elevation recovered quality efficiently.
文摘Bayesian empirical likelihood is a semiparametric method that combines parametric priors and nonparametric likelihoods, that is, replacing the parametric likelihood function in Bayes theorem with a nonparametric empirical likelihood function, which can be used without assuming the distribution of the data. It can effectively avoid the problems caused by the wrong setting of the model. In the variable selection based on Bayesian empirical likelihood, the penalty term is introduced into the model in the form of parameter prior. In this paper, we propose a novel variable selection method, L<sub>1/2</sub> regularization based on Bayesian empirical likelihood. The L<sub>1/2</sub> penalty is introduced into the model through a scale mixture of uniform representation of generalized Gaussian prior, and the posterior distribution is then sampled using MCMC method. Simulations demonstrate that the proposed method can have better predictive ability when the error violates the zero-mean normality assumption of the standard parameter model, and can perform variable selection.
文摘In view of the composition analysis and identification of ancient glass products, L1 regularization, K-Means cluster analysis, elbow rule and other methods were comprehensively used to build logical regression, cluster analysis, hyper-parameter test and other models, and SPSS, Python and other tools were used to obtain the classification rules of glass products under different fluxes, sub classification under different chemical compositions, hyper-parameter K value test and rationality analysis. Research can provide theoretical support for the protection and restoration of ancient glass relics.
基金financially supported by The 2011 Prospective Research Project of SINOPEC(P11096)
文摘Seismic data regularization is an important preprocessing step in seismic signal processing. Traditional seismic acquisition methods follow the Shannon–Nyquist sampling theorem, whereas compressive sensing(CS) provides a fundamentally new paradigm to overcome limitations in data acquisition. Besides the sparse representation of seismic signal in some transform domain and the 1-norm reconstruction algorithm, the seismic data regularization quality of CS-based techniques strongly depends on random undersampling schemes. For 2D seismic data, discrete uniform-based methods have been investigated, where some seismic traces are randomly sampled with an equal probability. However, in theory and practice, some seismic traces with different probability are required to be sampled for satisfying the assumptions in CS. Therefore, designing new undersampling schemes is imperative. We propose a Bernoulli-based random undersampling scheme and its jittered version to determine the regular traces that are randomly sampled with different probability, while both schemes comply with the Bernoulli process distribution. We performed experiments using the Fourier and curvelet transforms and the spectral projected gradient reconstruction algorithm for 1-norm(SPGL1), and ten different random seeds. According to the signal-to-noise ratio(SNR) between the original and reconstructed seismic data, the detailed experimental results from 2D numerical and physical simulation data show that the proposed novel schemes perform overall better than the discrete uniform schemes.
基金Supported by the National Natural Science Foundation of China(No:69872039)
文摘By defining fuzzy valued simple functions and giving L1(μ) approximations of fuzzy valued integrably bounded functions by such simple functions, the paper analyses by L1(μ)-norm the approximation capability of four-layer feedforward regular fuzzy neural networks to the fuzzy valued integrably bounded function F : Rn → FcO(R). That is, if the transfer functionσ: R→R is non-polynomial and integrable function on each finite interval, F may be innorm approximated by fuzzy valued functions defined as to anydegree of accuracy. Finally some real examples demonstrate the conclusions.
基金Supported by the National Nature Science Foundation of China(Grant Nos.11301052,11301045,11271060,11601064,11671068)the Fundamental Research Funds for the Central Universities(Grant No.DUT16LK33)the Fundamental Research of Civil Aircraft(Grant No.MJ-F-2012-04)
文摘We propose an ?~1 regularized method for numerical differentiation using empirical eigenfunctions. Compared with traditional methods for numerical differentiation, the output of our method can be considered directly as the derivative of the underlying function. Moreover,our method could produce sparse representations with respect to empirical eigenfunctions.Numerical results show that our method is quite effective.
文摘L(d, 1)-labeling is a kind of graph coloring problem from frequency assignment in radio networks, in which adjacent nodes must receive colors that are at least d apart while nodes at distance two from each other must receive different colors. We focus on L(d, 1)-labeling of regular tilings for d≥3 since the cases d=0, 1 or 2 have been researched by Calamoneri and Petreschi. For all three kinds of regular tilings, we give their L (d, 1)-labeling numbers for any integer d≥3. Therefore, combined with the results given by Calamoneri and Petreschi, the L(d, 1)-labeling numbers of regular tilings for any nonnegative integer d may be determined completely.
基金Supported by Open Foundation Project of Shenyang Institute of Atmospheric Environment,China Meteorological Administration(2016SYIAE14)Natural Science Foundation of Anhui Province,China(1708085QD89)National Natural Science Foundation of China(41805080)
文摘The paper discusses the core parameters of the 3 D and 4 D variational merging based on L1 norm regularization,namely optimization characteristic correlation length of background error covariance matrix and regularization parameter. Classical 3 D/4 D variational merging is based on the theory that error follows Gaussian distribution. It involves the solution of the objective functional gradient in minimization iteration,which requires the data to have continuity and differentiability. Classic 3 D/4 D-dimensional variational merging method was extended,and L1 norm was used as the constraint coupling to the classical variational merged model. Experiment was carried out by using linear advection-diffusion equation as four-dimensional prediction model,and parameter optimization of this method is discussed. Considering the strong temporal and spatial variation of water vapor,this method is further applied to the precipitable water vapor( PWV) merging by calculating reanalysis data and GNSS retrieval.Parameters were adjusted gradually to analyze the influence of background field on the merging result,and the experiment results show that the mathematical algorithm adopted in this paper is feasible.
基金National Natural Science Foundation of China(No.11671258)。
文摘The properties of generalized flip Markov chains on connected regular digraphs are discussed.The 1-Flipper operation on Markov chains for undirected graphs is generalized to that for multi-digraphs.The generalized 1-Flipper operation preserves the regularity and weak connectivity of multi-digraphs.The generalized 1-Flipper operation is proved to be symmetric.Moreover,it is presented that a series of random generalized 1-Flipper operations eventually lead to a uniform probability distribution over all connected d-regular multi-digraphs without loops.
基金This work is supported in part by the Foundation of Zhongshan University, Advanced Research Center.
文摘Let us consider the following elliptic systems of second order-D_α(A_i~α(x, u, Du))=B_4(x, u, Du), i=1, …, N, x∈Q(?)R^n, n≥3 (1) and supposeⅰ) |A_i~α(x, u, Du)|≤L(1+|Du|);ⅱ) (1+|p|)^(-1)A_i~α(x, u, p)are H(?)lder-continuous functions with some exponent δ on (?)×R^N uniformly with respect to p, i.e.ⅲ) A_i~α(x, u, p) are differentiable function in p with bounded and continuous derivativesⅳ)ⅴ) for all u∈H_(loc)~1(Ω, R^N)∩L^(n(γ-1)/(2-γ))(Ω, R^N), B(x, u, Du)is ineasurable and |B(x, u, p)|≤a(|p|~γ+|u|~τ)+b(x), where 1+2/n<γ<2, τ≤max((n+2)/(n-2), (γ-1)/(2-γ)-ε), (?)ε>0, b(x)∈L2n/(n+2), n^2/(n+2)+e(Ω), (?)ε>0.Remarks. Only bounded open set Q will be considered in this paper; for all p≥1, λ≥0, which is clled a Morrey Space.Let assumptions ⅰ)-ⅳ) hold, Giaquinta and Modica have proved the regularity of both the H^1 weak solutions of (1) under controllable growth condition |B|≤α(|p|~γ+|u|^((n+2)/(n-2))+b, 0<γ≤1+2/n and the H^1∩L~∞ weak solutions of (1) under natural growth condition |B|≤α|p|~2+b with a smallness condition 2aM<λ(|u|≤M), which implys that the H^1∩L~∞ weak solutions have the same regularty in the case of 1+2/n<γ<2. In the case of γ=2, many counterexamples (see [2] showed that u must be in H^1L~∞, while in the case of 1+2/n<γ<2, we consider the H^1∩L^n(γ-1)/(2-γ) weak solutions of (1), weaken the instability conditions upon them (from L~∞ to L^n(γ-1)/(2-γ) and obtain the same regularity results. Finally we show that the exponent n(γ-1)/(2-γ) can not be docreased anymore for the sake of the regularity results.Delinition 1. We call u∈H^1∩L^n(γ-1)/(2-γ)(Q, R^N) be a weak solution of (1), providod that where We use the convention that repeated indices are summed. i, j go from 1 to N ann α, β from 1 to n.
基金supported by the National Natural Science Foundation of China(11931009,12271495,11971450,and 12071449)Anhui Initiative in Quantum Information Technologies(AHY150200)the Project of Stable Support for Youth Team in Basic Research Field,Chinese Academy of Sciences(YSBR-001).
文摘On a compact Riemann surface with finite punctures P_(1),…P_(k),we define toric curves as multivalued,totallyunramified holomorphic maps to P^(n)with monodromy in a maximal torus of PSU(n+1).Toric solutions to SU(n+1)Todasystems on X\{P_(1);…;P_(k)}are recognized by the associated toric curves in.We introduce character n-ensembles as-tuples of meromorphic one-forms with simple poles and purely imaginary periods,generating toric curves on minus finitelymany points.On X,we establish a correspondence between character-ensembles and toric solutions to the SU(n+1)system with finitely many cone singularities.Our approach not only broadens seminal solutions with two conesingularities on the Riemann sphere,as classified by Jost-Wang(Int.Math.Res.Not.,2002,(6):277-290)andLin-Wei-Ye(Invent.Math.,2012,190(1):169-207),but also advances beyond the limits of Lin-Yang-Zhong’s existencetheorems(J.Differential Geom.,2020,114(2):337-391)by introducing a new solution class.
基金supported by the National Science and Technology Major Project (No.2011ZX05023-005-008)
文摘In this paper, we built upon the estimating primaries by sparse inversion (EPSI) method. We use the 3D curvelet transform and modify the EPSI method to the sparse inversion of the biconvex optimization and Ll-norm regularization, and use alternating optimization to directly estimate the primary reflection coefficients and source wavelet. The 3D curvelet transform is used as a sparseness constraint when inverting the primary reflection coefficients, which results in avoiding the prediction subtraction process in the surface-related multiples elimination (SRME) method. The proposed method not only reduces the damage to the effective waves but also improves the elimination of multiples. It is also a wave equation- based method for elimination of surface multiple reflections, which effectively removes surface multiples under complex submarine conditions.