Acoustic fields with impedance boundary conditions have high engineering applications, such as noise control and evaluation of sound insulation materials, and can be approximated by three-dimensional Helmholtz boundar...Acoustic fields with impedance boundary conditions have high engineering applications, such as noise control and evaluation of sound insulation materials, and can be approximated by three-dimensional Helmholtz boundary value problems. Finite difference method is widely applied to solving these problems due to its ease of use. However, when the wave number is large, the pollution effects are still a major difficulty in obtaining accurate numerical solutions. We develop a fast algorithm for solving three-dimensional Helmholtz boundary problems with large wave numbers. The boundary of computational domain is discrete based on high-order compact difference scheme. Using the properties of the tensor product and the discrete Fourier sine transform method, the original problem is solved by splitting it into independent small tridiagonal subsystems. Numerical examples with impedance boundary conditions are used to verify the feasibility and accuracy of the proposed algorithm. Results demonstrate that the algorithm has a fourth- order convergence in and -norms, and costs less CPU calculation time and random access memory.展开更多
This paper is concerned with the Diophantine properties of the sequence {ξθn}, where 1 ≤ξ 〈 θ and θ is a rational or an algebraic integer. We establish a combinatorial proposition which can be used to study suc...This paper is concerned with the Diophantine properties of the sequence {ξθn}, where 1 ≤ξ 〈 θ and θ is a rational or an algebraic integer. We establish a combinatorial proposition which can be used to study such two cases in the same manner. It is shown that the decay rate of the Fourier transforms of self-similar measures μλ with λ = θ-1 as the uniform contractive ratio is logarithmic. This generalizes some results of Kershner and Bufetov-Solomyak, who consider the case of Bernoulli convolutions. As an application, we prove that μλ ahaost every x is normal to any base b ≥ 2, which implies that there exist infinitely many absolute normal numbers on the corresponding self-similar set. This can be seen as a complementary result of the well-known Cassels-Schmidt theorem.展开更多
In the history of mathematics different methods have been used to detect if a number is prime or not. In this paper a new one will be shown. It will be demonstrated that if the following equation is zero for a certain...In the history of mathematics different methods have been used to detect if a number is prime or not. In this paper a new one will be shown. It will be demonstrated that if the following equation is zero for a certain number p, this number p would be prime. And being m an integer number higher than (the lowest, the most efficient the operation). . If the result is an integer, this result will tell us how many permutations of two divisors, the input number has. As you can check, no recurrent division by odd or prime numbers is done, to check if the number is prime or has divisors. To get to this point, we will do the following. First, we will create a domain with all the composite numbers. This is easy, as you can just multiply one by one all the integers (greater or equal than 2) in that domain. So, you will get all the composite numbers (not getting any prime) in that domain. Then, we will use the Fourier transform to change from this original domain (called discrete time domain in this regards) to the frequency domain. There, we can check, using Parseval’s theorem, if a certain number is there or not. The use of Parseval’s theorem leads to the above integral. If the number p that we want to check is not in the domain, the result of the integral is zero and the number is a prime. If instead, the result is an integer, this integer will tell us how many permutations of two divisors the number p has. And, in consequence information how many factors, the number p has. So, for any number p lower than 2m?- 1, you can check if it is prime or not, just making the numerical definite integration. We will apply this integral in a computer program to check the efficiency of the operation. We will check, if no further developments are done, the numerical integration is inefficient computing-wise compared with brute-force checking. To be added, is the question regarding the level of accuracy needed (number of decimals and number of steps in the numerical integration) to have a reliable result for large numbers. This will be commented on the paper, but a separate study will be needed to have detailed conclusions. Of course, the best would be that in the future, an analytical result (or at least an approximation) for the summation or for the integration is achieved.展开更多
A new interpolation algorithm for Head-Related Transfer Function (HRTF) is proposed to realize 3D sound reproduction via headphones in arbitrary spatial direction. HRTFs are modeled as a weighted sum of spherical harm...A new interpolation algorithm for Head-Related Transfer Function (HRTF) is proposed to realize 3D sound reproduction via headphones in arbitrary spatial direction. HRTFs are modeled as a weighted sum of spherical harmonics on a spherical surface. Truncated Singular Value Decomposition (SVD) is adopted to calculate the weights of the model. The truncation number is chosen according to Frobenius norm ratio and the partial condition number. Compared with other interpolated methods, our proposed approach not only is continuous but exploits global information of available directions. The HRTF from any desired direction can be obtained more accurately and robustly. Reconstructed and interpolated results demonstrate that our proposed algorithm acquired better performance.展开更多
The study of large-scale atmospheric turbulence and transport processes is of vital importance in the general circulation of the atmosphere. The governing equations of the power and cross-spectra for the atmospheric m...The study of large-scale atmospheric turbulence and transport processes is of vital importance in the general circulation of the atmosphere. The governing equations of the power and cross-spectra for the atmospheric motion and transports in the domain of wave number frequency space have been derived. The contributions of the nonlinear interactions of the atmospheric waves in velocity and temperature fields to the conversion of kinetic and potential energies and to the meridional transports of angular momentum and sensible heat in the atmosphere have been discussed.展开更多
文摘Acoustic fields with impedance boundary conditions have high engineering applications, such as noise control and evaluation of sound insulation materials, and can be approximated by three-dimensional Helmholtz boundary value problems. Finite difference method is widely applied to solving these problems due to its ease of use. However, when the wave number is large, the pollution effects are still a major difficulty in obtaining accurate numerical solutions. We develop a fast algorithm for solving three-dimensional Helmholtz boundary problems with large wave numbers. The boundary of computational domain is discrete based on high-order compact difference scheme. Using the properties of the tensor product and the discrete Fourier sine transform method, the original problem is solved by splitting it into independent small tridiagonal subsystems. Numerical examples with impedance boundary conditions are used to verify the feasibility and accuracy of the proposed algorithm. Results demonstrate that the algorithm has a fourth- order convergence in and -norms, and costs less CPU calculation time and random access memory.
基金Project(1834201)supported by the National Natural Science Foundation of ChinaProject(2020YJ0076)supported by the Sichuan Science and Technology Program,China+1 种基金Project(2682020CX35)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(2020M673280)supported by the Postdoctoral Science Foundation,China。
文摘This paper is concerned with the Diophantine properties of the sequence {ξθn}, where 1 ≤ξ 〈 θ and θ is a rational or an algebraic integer. We establish a combinatorial proposition which can be used to study such two cases in the same manner. It is shown that the decay rate of the Fourier transforms of self-similar measures μλ with λ = θ-1 as the uniform contractive ratio is logarithmic. This generalizes some results of Kershner and Bufetov-Solomyak, who consider the case of Bernoulli convolutions. As an application, we prove that μλ ahaost every x is normal to any base b ≥ 2, which implies that there exist infinitely many absolute normal numbers on the corresponding self-similar set. This can be seen as a complementary result of the well-known Cassels-Schmidt theorem.
文摘In the history of mathematics different methods have been used to detect if a number is prime or not. In this paper a new one will be shown. It will be demonstrated that if the following equation is zero for a certain number p, this number p would be prime. And being m an integer number higher than (the lowest, the most efficient the operation). . If the result is an integer, this result will tell us how many permutations of two divisors, the input number has. As you can check, no recurrent division by odd or prime numbers is done, to check if the number is prime or has divisors. To get to this point, we will do the following. First, we will create a domain with all the composite numbers. This is easy, as you can just multiply one by one all the integers (greater or equal than 2) in that domain. So, you will get all the composite numbers (not getting any prime) in that domain. Then, we will use the Fourier transform to change from this original domain (called discrete time domain in this regards) to the frequency domain. There, we can check, using Parseval’s theorem, if a certain number is there or not. The use of Parseval’s theorem leads to the above integral. If the number p that we want to check is not in the domain, the result of the integral is zero and the number is a prime. If instead, the result is an integer, this integer will tell us how many permutations of two divisors the number p has. And, in consequence information how many factors, the number p has. So, for any number p lower than 2m?- 1, you can check if it is prime or not, just making the numerical definite integration. We will apply this integral in a computer program to check the efficiency of the operation. We will check, if no further developments are done, the numerical integration is inefficient computing-wise compared with brute-force checking. To be added, is the question regarding the level of accuracy needed (number of decimals and number of steps in the numerical integration) to have a reliable result for large numbers. This will be commented on the paper, but a separate study will be needed to have detailed conclusions. Of course, the best would be that in the future, an analytical result (or at least an approximation) for the summation or for the integration is achieved.
基金Supported by Shanghai Natural Science Foundation, Shanghai Leading Academic Discipline Project, and STCSM of China (No. 08ZR1408300, S30108, and 08DZ2231100)
文摘A new interpolation algorithm for Head-Related Transfer Function (HRTF) is proposed to realize 3D sound reproduction via headphones in arbitrary spatial direction. HRTFs are modeled as a weighted sum of spherical harmonics on a spherical surface. Truncated Singular Value Decomposition (SVD) is adopted to calculate the weights of the model. The truncation number is chosen according to Frobenius norm ratio and the partial condition number. Compared with other interpolated methods, our proposed approach not only is continuous but exploits global information of available directions. The HRTF from any desired direction can be obtained more accurately and robustly. Reconstructed and interpolated results demonstrate that our proposed algorithm acquired better performance.
文摘The study of large-scale atmospheric turbulence and transport processes is of vital importance in the general circulation of the atmosphere. The governing equations of the power and cross-spectra for the atmospheric motion and transports in the domain of wave number frequency space have been derived. The contributions of the nonlinear interactions of the atmospheric waves in velocity and temperature fields to the conversion of kinetic and potential energies and to the meridional transports of angular momentum and sensible heat in the atmosphere have been discussed.