In order to carry out tensor analysis in a neighborhood of a reference surface,the principal-direction orthogonal basis accompanying with Lame s coefficients or general curvilinear coordinate systems are widely used.A...In order to carry out tensor analysis in a neighborhood of a reference surface,the principal-direction orthogonal basis accompanying with Lame s coefficients or general curvilinear coordinate systems are widely used.A novel kind of field theory termed as the nonholonomic theory of the Principal-Direction Orthonormal Basis(PDOB)is presented systematically in the present paper,in which the formal Christoffel symbols are related directly to the principal and geodesic curvatures with respect to the principal directions of the surface.Furthermore,a systematic and simple way to determine the curvatures of the surface are presented with some examples.It provides a way to recognize qualitatively the bending property of a surface.展开更多
We present a simple demonstration on the orthonormality of Volkov solutions with emphasizing on the suffcient condition to the orthonormality. Properly aligning the external electromagnetic wave along the third-axis, ...We present a simple demonstration on the orthonormality of Volkov solutions with emphasizing on the suffcient condition to the orthonormality. Properly aligning the external electromagnetic wave along the third-axis, the Volkov solutions are eigenfunctions of the hermitian momentum ■1, ■2 and the light-cone hamiltonian operators with real eigenvalues, which can lead to a verification of the orthonormality in the context of quantum mechanics when the x3-integration of the external potential is not singularity as severe as δ(0). The hermiticity of the fermion field fourmomentum operators validates the application of the demonstration to the intense field quantum electrodynamic. The proof based on a direct calculation to the inner products of the solutions is recapitulated as well in a general manner without dependence on explicit representation of the Dirac matrices and spinors, which can be conducive to understand the suffcient condition and to the study of the polarized electron production where a convenient representation is selected elaborately to project out the spin-polarization.展开更多
A two-degree-of-freedom vibro-impact system having symmetrical rigid stops and subjected to periodic excitation is investigated in this paper. By introducing local maps between different stages of motion in the whole ...A two-degree-of-freedom vibro-impact system having symmetrical rigid stops and subjected to periodic excitation is investigated in this paper. By introducing local maps between different stages of motion in the whole impact process, the Poincare map of the system is constructed. Using the Poincare map and the Gram Schmidt orthonormalization, a method of calculating the spectrum of Lyapunov exponents of the above vibro-impact system is presented. Then the phase portraits of periodic and chaotic attractors for the system and the corresponding convergence diagrams of the spectrum of Lyapunov exponents are given out through the numerical simulations. To further identify the validity of the aforementioned computation method, the bifurcation diagram of the system with respect to the bifurcation parameter and the corresponding largest Lyapunov exponents are shown.展开更多
A parameter-free method based on orthonormal wavelet transforms is recommended for calculating the principal time scale of coherent structures in atmospheric boundary-layer measurements. First, the atmospheric turbule...A parameter-free method based on orthonormal wavelet transforms is recommended for calculating the principal time scale of coherent structures in atmospheric boundary-layer measurements. First, the atmospheric turbulent signal is decomposed into the small scale vortex that has approximate isotropy and the large scale vortex with the digital filter. Then, the large scale vortex is used to detect coherent structures with this method. The principal time scale and profile of coherent structures for velocity components (u, v, w above rice fields are obtained. In order to testify the validity of this method, the correlation of coherent structures and non-coherent structures are also calculated.展开更多
Motivated by representing multidimensional periodic nonlinear and non-stationary signals (functions), we study a class of orthonormal exponential basis for L2(Id) with I := [0, 1), whose exponential parts are pi...Motivated by representing multidimensional periodic nonlinear and non-stationary signals (functions), we study a class of orthonormal exponential basis for L2(Id) with I := [0, 1), whose exponential parts are piecewise linear spectral sequences with p-knots. It is widely applied in time-frequency analysis.展开更多
In this paper we show how to construct a scaling function and an orthonor-mal wavelet basis from a multiresolution approximation using an operator theoreticmethod.
In this paper,we propose a new linear multiuser receiver for synchronous code-division multiple-access (CDMA) systems,referred to as the orthogonal multiuser(OMU) receiver.Unlike the linear minimum mean-squared error(...In this paper,we propose a new linear multiuser receiver for synchronous code-division multiple-access (CDMA) systems,referred to as the orthogonal multiuser(OMU) receiver.Unlike the linear minimum mean-squared error(MMSE) receiver,the OMU receiver depends only on the signature vectors and does not require knowledge of the received amplitudes or the channel signal-to-noise ratio(SNR).Here we develop methods that construct an optimal set of vectors with a specified inner product structure,from a given set of vectors in a complex Hilbert space.The optimal vectors are chosen to minimize the sum of the squared norms of the errors between the constructed vectors and the given vectors.An algorithm has been developed using the principles of quantum parameters and some of its axioms and constraints.In place of the classical matched filter(MF) receiver we propose a modified receiver.This approach assumes that improving the accuracy will necessarily result in im- proved performance.The simulation results provided here suggests that in certain cases the OMU and POMU receivers can significantly increase the probability of correct detection with low error rate over the MF receiver.展开更多
The modal wave number tomography approach is used to obtain sound speed profile of water column in deep ocean. The approach consists of estimation of the local modal eigenvalues from complex pressure field and use of ...The modal wave number tomography approach is used to obtain sound speed profile of water column in deep ocean. The approach consists of estimation of the local modal eigenvalues from complex pressure field and use of these data as input to modal perturbative inversion method for obtaining the local sound speed profile. The empirical orthonormal function (EOF) is applied to reduce the parameter search space. The ocean environment used for numerical simulations includes the Munk profile as the unperturbed background speed profile and a weak Gaussian eddy as the sound speed profile perturbation. The results of numerical simulations show the method is capable of monitoring the oceanic interior structure.展开更多
Routh stability test is covered in almost all undergraduate control texts. It determines the stability or, a little beyond, the number of unstable roots of a polynomial in terms of the signs of certain entries of the ...Routh stability test is covered in almost all undergraduate control texts. It determines the stability or, a little beyond, the number of unstable roots of a polynomial in terms of the signs of certain entries of the Routh table constructed from the coefficients of the polynomial. The use of the Routh table, as far as the common textbooks show, is only limited to this function. We will show that the Routh table can actually be used to construct an orthonormal basis in the space of strictly proper rational functions with a common stable denominator. This orthonormal basis can then be used for many other purposes, including the computation of the H2 norm, the Hankel singular values and singular vectors, model reduction, H∞ optimization, etc. Keywords Routh stablity criterion - Orthonormal basis - Root-mean-squared value - Hankel operater - Nehari problem - Model reduction This work was supported by the Hong Kong Research Grants Council.展开更多
A continuous-time Model Predictive Controller was proposed using Kautz function in order to improve the performance of Load Frequency Control(LFC).A dynamic model of an interconnected power system was used for Model P...A continuous-time Model Predictive Controller was proposed using Kautz function in order to improve the performance of Load Frequency Control(LFC).A dynamic model of an interconnected power system was used for Model Predictive Controller(MPC)design.MPC predicts the future trajectory of the dynamic model by calculating the optimal closed loop feedback gain matrix.In this paper,the optimal closed loop feedback gain matrix was calculated using Kautz function.Being an Orthonormal Basis Function(OBF),Kautz function has an advantage of solving complex pole-based nonlinear system.Genetic Algorithm(GA)was applied to optimally tune the Kautz function-based MPC.A constraint based on phase plane analysis was implemented with the cost function in order to improve the robustness of the Kautz function-based MPC.The proposed method was simulated with three area interconnected power system and the efficiency of the proposed method was measured and exhibited by comparing with conventional Proportional and Integral(PI)controller and Linear Quadratic Regulation(LQR).展开更多
There are many forms of tensor theory which are quite different. Discussions are put forward and their relationships are found. The differences between them depend on whether there is metric on the space and the basis...There are many forms of tensor theory which are quite different. Discussions are put forward and their relationships are found. The differences between them depend on whether there is metric on the space and the basis is orthonormal.展开更多
Suppose M and N are two r×r and s×s dilation matrices,respectively.LetΓM andΓN represent the complete sets of representatives of distinct cosets of the quotient groups M-T Zr/Zr and N-T Zs/Zs,respectively....Suppose M and N are two r×r and s×s dilation matrices,respectively.LetΓM andΓN represent the complete sets of representatives of distinct cosets of the quotient groups M-T Zr/Zr and N-T Zs/Zs,respectively.Two methods for constructing nonseparable Ω-filter banks from M-filter banks and N-filter banks are presented,where Ω is a(r+s) ×(r+s) dilation matrix such that one of its complete sets of representatives of distinct cosets of the quotient groups Ω-T Zr+s/Zr+s areΓΩ={[γT h,ζ T q] T:γh∈ΓM,ζq∈ΓN}.Specially,Ω can be [MΘ0N],whereΘis a r×s integer matrix with M-1Θbeing also an integer matrix.Moreover,if the constructed filter bank satisfies Lawton's condition,which can be easy to verify,then it generates an orthonormal nonseparable Ω-wavelet basis for L2(Rr+s).Properties,including Lawton's condition,vanishing moments and regularity of the new Ω-filter banks or new Ω-wavelet basis are discussed then.Finally,a class of nonseparable Ω-wavelet basis for L2(Rr+1) are constructed and three other examples are given to illustrate the results.In particular,when M=N=2,all results obtained in this paper appeared in[1].展开更多
A new modification of the Homotopy Analysis Method (HAM) is presented for highly nonlinear ODEs on a semi-infinite domain. The main advantage of the modified HAM is that the number of terms in the series solution can ...A new modification of the Homotopy Analysis Method (HAM) is presented for highly nonlinear ODEs on a semi-infinite domain. The main advantage of the modified HAM is that the number of terms in the series solution can be greatly reduced;meanwhile the accuracy of the solution can be well retained. In this way, much less CPU is needed. Two typical examples are used to illustrate the efficiency of the proposed approach.展开更多
Let is the Walsh generalized system. In the paper constructed a weighted space , and series in the Walsh generalized system with monotonically decreasing coefficient such that for each function in the space one can fi...Let is the Walsh generalized system. In the paper constructed a weighted space , and series in the Walsh generalized system with monotonically decreasing coefficient such that for each function in the space one can find a subseries that converges to in the weighted and almost everywhere on [0,1].展开更多
Let An∈M2(ℤ)be integral matrices such that the infinite convolution of Dirac measures with equal weightsμ{A_(n),n≥1}δA_(1)^(-1)D*δA_(1)^(-1)A_(2)^(-2)D*…is a probability measure with compact support,where D={(0,...Let An∈M2(ℤ)be integral matrices such that the infinite convolution of Dirac measures with equal weightsμ{A_(n),n≥1}δA_(1)^(-1)D*δA_(1)^(-1)A_(2)^(-2)D*…is a probability measure with compact support,where D={(0,0)^(t),(1,0)^(t),(0,1)^(t)}is the Sierpinski digit.We prove that there exists a setΛ⊂ℝ2 such that the family{e2πi〈λ,x〉:λ∈Λ} is an orthonormal basis of L^(2)(μ{A_(n),n≥1})if and only if 1/3(1,-1)A_(n)∈Z^(2)for n≥2 under some metric conditions on A_(n).展开更多
In this paper,based on the imaginary time gradient flow model in the density functional theory,a scalar auxiliary variable(SAV)method is developed for the ground state calculation of a given electronic structure syste...In this paper,based on the imaginary time gradient flow model in the density functional theory,a scalar auxiliary variable(SAV)method is developed for the ground state calculation of a given electronic structure system.To handle the orthonormality constraint on those wave functions,two kinds of penalty terms are introduced in designing the modified energy functional in SAV,i.e.,one for the norm preserving of each wave function,another for the orthogonality between each pair of different wave functions.A numerical method consisting of a designed scheme and a linear finite element method is used for the discretization.Theoretically,the desired unconditional decay of modified energy can be obtained from our method,while computationally,both the original energy and modified energy decay behaviors can be observed successfully from a number of numerical experiments.More importantly,numerical results show that the orthonormality among those wave functions can be automatically preserved,without explicitly preserving orthogonalization operations.This implies the potential of our method in large-scale simulations in density functional theory.展开更多
The truncated version of the generalized minimal residual method (GMRES), the incomplete generalized minimal residual method (IGMRES), is studied. It is based on an incomplete orthogonalization of the Krylov vectors i...The truncated version of the generalized minimal residual method (GMRES), the incomplete generalized minimal residual method (IGMRES), is studied. It is based on an incomplete orthogonalization of the Krylov vectors in question, and gives an approximate or quasi_minimum residual solution over the Krylov subspace. A convergence analysis of this method is given, showing that in the non_restarted version IGMRES can behave like GMRES once the basis vectors of Krylov subspace generated by the incomplete orthogonalization are strongly linearly independent. Meanwhile, some relationships between the residual norms for IOM and IGMRES are established. Numerical experiments are reported to show convergence behavior of IGMRES and of its restarted version IGMRES( m ).展开更多
Large dimensional predictors are often introduced in regressions to attenuate the possible modeling bias. We consider the stable direction recovery in single-index models in which we solely assume the response Y is in...Large dimensional predictors are often introduced in regressions to attenuate the possible modeling bias. We consider the stable direction recovery in single-index models in which we solely assume the response Y is independent of the diverging dimensional predictors X when βτ 0 X is given, where β 0 is a p n × 1 vector, and p n →∞ as the sample size n →∞. We first explore sufficient conditions under which the least squares estimation β n0 recovers the direction β 0 consistently even when p n = o(√ n). To enhance the model interpretability by excluding irrelevant predictors in regressions, we suggest an e1-regularization algorithm with a quadratic constraint on magnitude of least squares residuals to search for a sparse estimation of β 0 . Not only can the solution β n of e1-regularization recover β 0 consistently, it also produces sufficiently sparse estimators which enable us to select "important" predictors to facilitate the model interpretation while maintaining the prediction accuracy. Further analysis by simulations and an application to the car price data suggest that our proposed estimation procedures have good finite-sample performance and are computationally efficient.展开更多
In this paper, we consider the problem of the existence of general non-separable variate orthonormal compactly supported wavelet basis when the symbol function has a special form. We prove that the general non-separab...In this paper, we consider the problem of the existence of general non-separable variate orthonormal compactly supported wavelet basis when the symbol function has a special form. We prove that the general non-separable variate orthonormal wavelet basis doesn't exist if the symbol function possesses a certain form. This helps us to explicate the difficulty of constructing the non-separable variate wavlet basis and to hint how to construct non-separable variate wavlet basis.展开更多
基金Project supported by the National Natural Science Foundation of China(11972120,11472082,12032016)。
文摘In order to carry out tensor analysis in a neighborhood of a reference surface,the principal-direction orthogonal basis accompanying with Lame s coefficients or general curvilinear coordinate systems are widely used.A novel kind of field theory termed as the nonholonomic theory of the Principal-Direction Orthonormal Basis(PDOB)is presented systematically in the present paper,in which the formal Christoffel symbols are related directly to the principal and geodesic curvatures with respect to the principal directions of the surface.Furthermore,a systematic and simple way to determine the curvatures of the surface are presented with some examples.It provides a way to recognize qualitatively the bending property of a surface.
基金Supported in part by the National Natural Science Foundation of China under Grants Nos.11475258,11205242,and 11675263
文摘We present a simple demonstration on the orthonormality of Volkov solutions with emphasizing on the suffcient condition to the orthonormality. Properly aligning the external electromagnetic wave along the third-axis, the Volkov solutions are eigenfunctions of the hermitian momentum ■1, ■2 and the light-cone hamiltonian operators with real eigenvalues, which can lead to a verification of the orthonormality in the context of quantum mechanics when the x3-integration of the external potential is not singularity as severe as δ(0). The hermiticity of the fermion field fourmomentum operators validates the application of the demonstration to the intense field quantum electrodynamic. The proof based on a direct calculation to the inner products of the solutions is recapitulated as well in a general manner without dependence on explicit representation of the Dirac matrices and spinors, which can be conducive to understand the suffcient condition and to the study of the polarized electron production where a convenient representation is selected elaborately to project out the spin-polarization.
基金supported by the National Natural Science Foundation of China (61071189)Innovation Scientists and Technicians Troop Construction of Henan Province of China (084100510012)the Natural Science Foundation for the Education Department of Henan Province of China (2008B510001)
文摘In this paper, a characterization of orthonormal wavelet families in Sobolev spaces H s (R) is established.
基金supported by the National Natural Science Foundation of China (Grant No. 10972059)the Natural Science Foundation of the Guangxi Zhuang Autonmous Region of China (Grant Nos. 0640002 and 2010GXNSFA013110)+1 种基金the Guangxi Youth Science Foundation of China (Grant No. 0832014)the Project of Excellent Innovating Team of Guangxi University
文摘A two-degree-of-freedom vibro-impact system having symmetrical rigid stops and subjected to periodic excitation is investigated in this paper. By introducing local maps between different stages of motion in the whole impact process, the Poincare map of the system is constructed. Using the Poincare map and the Gram Schmidt orthonormalization, a method of calculating the spectrum of Lyapunov exponents of the above vibro-impact system is presented. Then the phase portraits of periodic and chaotic attractors for the system and the corresponding convergence diagrams of the spectrum of Lyapunov exponents are given out through the numerical simulations. To further identify the validity of the aforementioned computation method, the bifurcation diagram of the system with respect to the bifurcation parameter and the corresponding largest Lyapunov exponents are shown.
基金Acknowledgments. This research is supported by the Knowledge Innovative Foundation of Chinese Academy of Science (No. KZCX2-204, No. KZ-CX-SW-01-01B), and the National Natural Science Foundation of China (No. 40035010). The authors thank Professors Huang
文摘A parameter-free method based on orthonormal wavelet transforms is recommended for calculating the principal time scale of coherent structures in atmospheric boundary-layer measurements. First, the atmospheric turbulent signal is decomposed into the small scale vortex that has approximate isotropy and the large scale vortex with the digital filter. Then, the large scale vortex is used to detect coherent structures with this method. The principal time scale and profile of coherent structures for velocity components (u, v, w above rice fields are obtained. In order to testify the validity of this method, the correlation of coherent structures and non-coherent structures are also calculated.
基金Supported in part by the President Fund of GUCASSupported in part by National Natural Foundation of China(Grant No.10631080)National Natural Foundation of Beijing (Grant No.1092004)
文摘Motivated by representing multidimensional periodic nonlinear and non-stationary signals (functions), we study a class of orthonormal exponential basis for L2(Id) with I := [0, 1), whose exponential parts are piecewise linear spectral sequences with p-knots. It is widely applied in time-frequency analysis.
基金The first author got support in part from the fund provided by the University of North Carolina at Charlotte.The second author got support from the National Natural Science Foundation of China and the Science Foundation of Zhejiang Province.
文摘In this paper we show how to construct a scaling function and an orthonor-mal wavelet basis from a multiresolution approximation using an operator theoreticmethod.
文摘In this paper,we propose a new linear multiuser receiver for synchronous code-division multiple-access (CDMA) systems,referred to as the orthogonal multiuser(OMU) receiver.Unlike the linear minimum mean-squared error(MMSE) receiver,the OMU receiver depends only on the signature vectors and does not require knowledge of the received amplitudes or the channel signal-to-noise ratio(SNR).Here we develop methods that construct an optimal set of vectors with a specified inner product structure,from a given set of vectors in a complex Hilbert space.The optimal vectors are chosen to minimize the sum of the squared norms of the errors between the constructed vectors and the given vectors.An algorithm has been developed using the principles of quantum parameters and some of its axioms and constraints.In place of the classical matched filter(MF) receiver we propose a modified receiver.This approach assumes that improving the accuracy will necessarily result in im- proved performance.The simulation results provided here suggests that in certain cases the OMU and POMU receivers can significantly increase the probability of correct detection with low error rate over the MF receiver.
文摘The modal wave number tomography approach is used to obtain sound speed profile of water column in deep ocean. The approach consists of estimation of the local modal eigenvalues from complex pressure field and use of these data as input to modal perturbative inversion method for obtaining the local sound speed profile. The empirical orthonormal function (EOF) is applied to reduce the parameter search space. The ocean environment used for numerical simulations includes the Munk profile as the unperturbed background speed profile and a weak Gaussian eddy as the sound speed profile perturbation. The results of numerical simulations show the method is capable of monitoring the oceanic interior structure.
文摘Routh stability test is covered in almost all undergraduate control texts. It determines the stability or, a little beyond, the number of unstable roots of a polynomial in terms of the signs of certain entries of the Routh table constructed from the coefficients of the polynomial. The use of the Routh table, as far as the common textbooks show, is only limited to this function. We will show that the Routh table can actually be used to construct an orthonormal basis in the space of strictly proper rational functions with a common stable denominator. This orthonormal basis can then be used for many other purposes, including the computation of the H2 norm, the Hankel singular values and singular vectors, model reduction, H∞ optimization, etc. Keywords Routh stablity criterion - Orthonormal basis - Root-mean-squared value - Hankel operater - Nehari problem - Model reduction This work was supported by the Hong Kong Research Grants Council.
文摘A continuous-time Model Predictive Controller was proposed using Kautz function in order to improve the performance of Load Frequency Control(LFC).A dynamic model of an interconnected power system was used for Model Predictive Controller(MPC)design.MPC predicts the future trajectory of the dynamic model by calculating the optimal closed loop feedback gain matrix.In this paper,the optimal closed loop feedback gain matrix was calculated using Kautz function.Being an Orthonormal Basis Function(OBF),Kautz function has an advantage of solving complex pole-based nonlinear system.Genetic Algorithm(GA)was applied to optimally tune the Kautz function-based MPC.A constraint based on phase plane analysis was implemented with the cost function in order to improve the robustness of the Kautz function-based MPC.The proposed method was simulated with three area interconnected power system and the efficiency of the proposed method was measured and exhibited by comparing with conventional Proportional and Integral(PI)controller and Linear Quadratic Regulation(LQR).
基金Sponsored by the National Natural Science Foundation of China(10372074)
文摘There are many forms of tensor theory which are quite different. Discussions are put forward and their relationships are found. The differences between them depend on whether there is metric on the space and the basis is orthonormal.
基金Supported by the National Natural Science Foundation of China(11071152)the Natural Science Foundation of Guangdong Province(10151503101000025)
文摘Suppose M and N are two r×r and s×s dilation matrices,respectively.LetΓM andΓN represent the complete sets of representatives of distinct cosets of the quotient groups M-T Zr/Zr and N-T Zs/Zs,respectively.Two methods for constructing nonseparable Ω-filter banks from M-filter banks and N-filter banks are presented,where Ω is a(r+s) ×(r+s) dilation matrix such that one of its complete sets of representatives of distinct cosets of the quotient groups Ω-T Zr+s/Zr+s areΓΩ={[γT h,ζ T q] T:γh∈ΓM,ζq∈ΓN}.Specially,Ω can be [MΘ0N],whereΘis a r×s integer matrix with M-1Θbeing also an integer matrix.Moreover,if the constructed filter bank satisfies Lawton's condition,which can be easy to verify,then it generates an orthonormal nonseparable Ω-wavelet basis for L2(Rr+s).Properties,including Lawton's condition,vanishing moments and regularity of the new Ω-filter banks or new Ω-wavelet basis are discussed then.Finally,a class of nonseparable Ω-wavelet basis for L2(Rr+1) are constructed and three other examples are given to illustrate the results.In particular,when M=N=2,all results obtained in this paper appeared in[1].
文摘A new modification of the Homotopy Analysis Method (HAM) is presented for highly nonlinear ODEs on a semi-infinite domain. The main advantage of the modified HAM is that the number of terms in the series solution can be greatly reduced;meanwhile the accuracy of the solution can be well retained. In this way, much less CPU is needed. Two typical examples are used to illustrate the efficiency of the proposed approach.
文摘Let is the Walsh generalized system. In the paper constructed a weighted space , and series in the Walsh generalized system with monotonically decreasing coefficient such that for each function in the space one can find a subseries that converges to in the weighted and almost everywhere on [0,1].
基金supported by the National Natural Science Foundation of China (Grant Nos. 12371087, 11971109,11971194, 11672074 and 12271185)supported by the program for Probability and Statistics:Theory and Application (Grant No. IRTL1704)+1 种基金the program for Innovative Research Team in Science and Technology in Fujian Province University (Grant No. IRTSTFJ)supported by Guangdong NSFC (Grant No. 2022A1515011124)
文摘Let An∈M2(ℤ)be integral matrices such that the infinite convolution of Dirac measures with equal weightsμ{A_(n),n≥1}δA_(1)^(-1)D*δA_(1)^(-1)A_(2)^(-2)D*…is a probability measure with compact support,where D={(0,0)^(t),(1,0)^(t),(0,1)^(t)}is the Sierpinski digit.We prove that there exists a setΛ⊂ℝ2 such that the family{e2πi〈λ,x〉:λ∈Λ} is an orthonormal basis of L^(2)(μ{A_(n),n≥1})if and only if 1/3(1,-1)A_(n)∈Z^(2)for n≥2 under some metric conditions on A_(n).
基金The first author would like to thank the support from the UM-Funded PhD Assistantship from University of MacaoThe second author was partially supported by Macao Young Scholar Program(AM201919)+5 种基金excellent youth project of Hunan Education Department(19B543)Hunan National Applied Mathematics Center of Hunan Provincial Science and Technology Department(2020ZYT003)The third author would like to thank financial support from National Natural Science Foundation of China(Grant Nos.11922120,11871489)FDCT of Macao SAR(Grant No.0082/2020/A2)University of Macao(Grant No.MYRG2020-00265-FST)Guangdong-Hong Kong-Macao Joint Laboratory for Data-Driven Fluid Mechanics and Engineering Applications(Grant No.2020B1212030001).
文摘In this paper,based on the imaginary time gradient flow model in the density functional theory,a scalar auxiliary variable(SAV)method is developed for the ground state calculation of a given electronic structure system.To handle the orthonormality constraint on those wave functions,two kinds of penalty terms are introduced in designing the modified energy functional in SAV,i.e.,one for the norm preserving of each wave function,another for the orthogonality between each pair of different wave functions.A numerical method consisting of a designed scheme and a linear finite element method is used for the discretization.Theoretically,the desired unconditional decay of modified energy can be obtained from our method,while computationally,both the original energy and modified energy decay behaviors can be observed successfully from a number of numerical experiments.More importantly,numerical results show that the orthonormality among those wave functions can be automatically preserved,without explicitly preserving orthogonalization operations.This implies the potential of our method in large-scale simulations in density functional theory.
文摘The truncated version of the generalized minimal residual method (GMRES), the incomplete generalized minimal residual method (IGMRES), is studied. It is based on an incomplete orthogonalization of the Krylov vectors in question, and gives an approximate or quasi_minimum residual solution over the Krylov subspace. A convergence analysis of this method is given, showing that in the non_restarted version IGMRES can behave like GMRES once the basis vectors of Krylov subspace generated by the incomplete orthogonalization are strongly linearly independent. Meanwhile, some relationships between the residual norms for IOM and IGMRES are established. Numerical experiments are reported to show convergence behavior of IGMRES and of its restarted version IGMRES( m ).
基金supported by National Natural Science Foundation of China (Grant No. 10701035)Chen Guang Project of Shanghai Education Development Foundation (Grant No. 2007CG33)+1 种基金supported by Research Grants Council of Hong KongFaculty Research Grant from Hong Kong Baptist University
文摘Large dimensional predictors are often introduced in regressions to attenuate the possible modeling bias. We consider the stable direction recovery in single-index models in which we solely assume the response Y is independent of the diverging dimensional predictors X when βτ 0 X is given, where β 0 is a p n × 1 vector, and p n →∞ as the sample size n →∞. We first explore sufficient conditions under which the least squares estimation β n0 recovers the direction β 0 consistently even when p n = o(√ n). To enhance the model interpretability by excluding irrelevant predictors in regressions, we suggest an e1-regularization algorithm with a quadratic constraint on magnitude of least squares residuals to search for a sparse estimation of β 0 . Not only can the solution β n of e1-regularization recover β 0 consistently, it also produces sufficiently sparse estimators which enable us to select "important" predictors to facilitate the model interpretation while maintaining the prediction accuracy. Further analysis by simulations and an application to the car price data suggest that our proposed estimation procedures have good finite-sample performance and are computationally efficient.
基金the National Natural Science Foundation of China (No.69982002) and theOpening Foundation of National Mobile Communications Res
文摘In this paper, we consider the problem of the existence of general non-separable variate orthonormal compactly supported wavelet basis when the symbol function has a special form. We prove that the general non-separable variate orthonormal wavelet basis doesn't exist if the symbol function possesses a certain form. This helps us to explicate the difficulty of constructing the non-separable variate wavlet basis and to hint how to construct non-separable variate wavlet basis.