Time series forecasting plays an important role in various fields, such as energy, finance, transport, and weather. Temporal convolutional networks (TCNs) based on dilated causal convolution have been widely used in t...Time series forecasting plays an important role in various fields, such as energy, finance, transport, and weather. Temporal convolutional networks (TCNs) based on dilated causal convolution have been widely used in time series forecasting. However, two problems weaken the performance of TCNs. One is that in dilated casual convolution, causal convolution leads to the receptive fields of outputs being concentrated in the earlier part of the input sequence, whereas the recent input information will be severely lost. The other is that the distribution shift problem in time series has not been adequately solved. To address the first problem, we propose a subsequence-based dilated convolution method (SDC). By using multiple convolutional filters to convolve elements of neighboring subsequences, the method extracts temporal features from a growing receptive field via a growing subsequence rather than a single element. Ultimately, the receptive field of each output element can cover the whole input sequence. To address the second problem, we propose a difference and compensation method (DCM). The method reduces the discrepancies between and within the input sequences by difference operations and then compensates the outputs for the information lost due to difference operations. Based on SDC and DCM, we further construct a temporal subsequence-based convolutional network with difference (TSCND) for time series forecasting. The experimental results show that TSCND can reduce prediction mean squared error by 7.3% and save runtime, compared with state-of-the-art models and vanilla TCN.展开更多
A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the c...A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.展开更多
Prestack reverse time migration (RTM) is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of t...Prestack reverse time migration (RTM) is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference (FD) method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.展开更多
Micro-satellite cluster enables a whole new class of missions for communications, remote sensing, and scientific research for both civilian and military purposes. Synchronizing the time of the satellites in a cluster ...Micro-satellite cluster enables a whole new class of missions for communications, remote sensing, and scientific research for both civilian and military purposes. Synchronizing the time of the satellites in a cluster is important for both cluster sensing capabilities and its autonomous operating. However, the existing time synchronization methods are not suitable for microsatellite cluster, because it requires too many human interventions and occupies too much ground control resource. Although, data post-process may realize the equivalent time synchronization, it requires processing time and powerful computing ability on the ground, which cannot be implemented by cluster itself. In order to autonomously establish and maintain the time benchmark in a cluster, we propose a compact time difference compensation system(TDCS), which is a kind of time control loop that dynamically adjusts the satellite reference frequency according to the time difference. Consequently, the time synchronization in the cluster can be autonomously achieved on-orbit by synchronizing the clock of other satellites to a chosen one's. The experimental result shows that the standard deviation of time synchronization is about 102 ps when the carrier to noise ratio(CNR) is 95 d BHz, and the standard deviation of corresponding frequency difference is approximately0.36 Hz.展开更多
The advancement of small satellites is promoting the development of distributed satellite systems,and for the latter,it is essential to coordinate the spatial and temporal relations between mutually visible satellites...The advancement of small satellites is promoting the development of distributed satellite systems,and for the latter,it is essential to coordinate the spatial and temporal relations between mutually visible satellites.By now,dual one-way ranging(DOWR)and two-way time transfer(TWTT)are generally integrated in the same software and hardware system to meet the limitations of small satellites in terms of size,weight and power(SWaP)consumption.However,studies show that pseudo-noise regenerative ranging(PNRR)performs better than DOWR if some advanced implementation technologies are employed.Besides,PNRR has no requirement on time synchronization.To apply PNRR to small satellites,and meanwhile,meet the demand for time difference measurement,we propose the round-way time difference measurement,which can be combined with PNRR to form a new integrated system without exceeding the limits of SWaP.The new integrated system can provide distributed small satellite systems with on-orbit high-accuracy and high-precision distance measurement and time difference measurement in real time.Experimental results show that the precision of ranging is about 1.94 cm,and that of time difference measurement is about 78.4 ps,at the signal to noise ratio of 80 dBHz.展开更多
Magneto-rhelological(MR) dampers are devices that employ rheological fluids to modify their mechanical properties. Their mechanical simplicity, high dynamic range, lower power requirements, large force capacity, rob...Magneto-rhelological(MR) dampers are devices that employ rheological fluids to modify their mechanical properties. Their mechanical simplicity, high dynamic range, lower power requirements, large force capacity, robustness and safe maimer of operation in cases of failure have made them attractive devices for semi-active real-time control in civil, aerospace and automotive applications. Time response characteristic is one of the most important technical performances of MR dampers, and response time directly affects the control frequency, application range and the actual effect of MR dampers. In this study, one kind of finite difference solution for predicting the response time ofmagneto-rheological dampers from "off-state" to "on-state" is put forward. A laminar flow model is used to describe the flow in the MR valve, and a bi-viscous fluid flow model is utilized to describe the relationship of shear stress and shear rate of MR fluid. An explicit difference format is used to diseretize the Novior-Stoks equation, and stability condition of this algorithm is built by Von-Neumann stability criterion. The pressure gradient along the flow duct is solved by a dichotomy algorithm with iteration, and the changing curve of the damping force versus time of MR damper is obtained as well. According to the abovementioned numerical algorithm, the damping forces versus time curves from "off-state" to "on-state" of a cylindrical piston type MR damper are computed. Moreover, the MR damper is installed in a material test system(MTS), the magnetic field in the wire circles of the MR damper is "triggered" when the MR damper is imposed to do a constant speed motion, and the damping force curves are recorded. The comparison between numerical results and experimental results indicates that this finite difference algorithm can be used to estimate the response time delay of MR dampers.展开更多
The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integ...The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.展开更多
The transmission coefficients of electromagnetic (EM) waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain (SO-FDTD) method is used in the ...The transmission coefficients of electromagnetic (EM) waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain (SO-FDTD) method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.展开更多
Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order tim...Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.展开更多
In order to improve the accuracy and engineering feasibility of four-Satellite localization system, the frequency difference measurement is introduced to the four-Satellite TDOA (Time Difference of Arrival) localizati...In order to improve the accuracy and engineering feasibility of four-Satellite localization system, the frequency difference measurement is introduced to the four-Satellite TDOA (Time Difference of Arrival) localization algorithm. The TDOA/FDOA (Frequency Difference of Arrival) localization algorithm is used to optimize the GDOP (geometric dilution of precision) of four-Satellite localization. The simulation results show that the absolute position measurement accuracy has little influence on TDOA/FDOA localization accuracy as compared with TDOA localization. Under the same conditions, TDOA/FDOA localization has better accuracy and its GDOP shows more uniform distribution in diamond configuration case. The localization accuracy of four-Satellite TDOA/FDOA is better than the localization accuracy of four-Satellite TDOA.展开更多
The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems wit...The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.展开更多
A structure-dependent explicit method with enhanced stability properties is proposed in this study. In general, the method offers unconditional stability for structural systems except those with a particular instantan...A structure-dependent explicit method with enhanced stability properties is proposed in this study. In general, the method offers unconditional stability for structural systems except those with a particular instantaneous stiffness hardening behavior. In addition, it is second-order accurate and displays no overshooting in high frequency responses. Numerical experiments reveal that the proposed method saves a substantial amount of computational effort in solving inertial problems where only the low frequency responses are of interest, when compared to a general second-order accurate integration method.展开更多
A new absorbing boundary condition (ABC) for frequency dependent finite difference time domain algorithm for the arbitrary dispersive media is presented. The concepts of the digital systems are introduced to the (F...A new absorbing boundary condition (ABC) for frequency dependent finite difference time domain algorithm for the arbitrary dispersive media is presented. The concepts of the digital systems are introduced to the (FD) 2TD method. On the basis of digital filter designing and vector algebra, the absorbing boundary condition under arbitrary angle of incidence are derived. The transient electromagnetic problems in two dimensions and three dimensions are calculated and the validity of the ABC is verified.展开更多
This paper introduces a real-time high precision measurement of phase difference based on Field Programmable Gate Array(FPGA) technology,which has been successfully applied to laser grating interference measurement ...This paper introduces a real-time high precision measurement of phase difference based on Field Programmable Gate Array(FPGA) technology,which has been successfully applied to laser grating interference measurement and real-time feedback of plasma electron density in HL-2A tokamak.It can track the changes of electron density while setting the starting point of the density curve to zero.In a laboratory test,the measuring accuracy of phase difference is less than 0.1°,the time resolution is 80 ns,and the feedback delay is 180 μs.展开更多
The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of c...The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.展开更多
A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-di...A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones, respectively, thus leading both the time and spatial derivatives to the second-order accuracy. The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift. Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach, respectively, and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor. From comparison of numerical results with theoretical or experimental ones, it is found that they are in reasonable agreement.展开更多
Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer(CPML) is given for solving the truncation problem of the open...Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer(CPML) is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain(ECT-CFDTD) method is used to simulate the wave propagation inside a perfect electric conductor(PEC) waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.展开更多
In this paper, we approximate the solution to time-fractional telegraph equation by two kinds of difference methods: the Grünwald formula and Caputo fractional difference.
Analyzed the relation between time delay difference and time delay estimation errors, based on the principles of three-point locating, a reformed threshold method for time delay difference estimation of impulse signal...Analyzed the relation between time delay difference and time delay estimation errors, based on the principles of three-point locating, a reformed threshold method for time delay difference estimation of impulse signals, called as amendment estimation for short, is developed by introducing channel equalization technique to its conventional version, named as direct estimation in this paper, to improve the estimation stability. After inherent relationship between time delay and phase shift of signals is analyzed, an integer period error compensation method utilized the diversities of both contribution share and contribution mode of concerned estimates is proposed under the condition of high precision phase lag estimation. Finally, a cooperative multi-threshold estimation method composed of amendment and direct estimations to process impulse signals with three thresholds is established. In sea trials data tests of passive locating, this method improves the estimation precision of time delay difference efficiently. The experiments verify the theoretical predictions.展开更多
基金supported by the National Key Research and Development Program of China(No.2018YFB2101300)the National Natural Science Foundation of China(Grant No.61871186)the Dean’s Fund of Engineering Research Center of Software/Hardware Co-Design Technology and Application,Ministry of Education(East China Normal University).
文摘Time series forecasting plays an important role in various fields, such as energy, finance, transport, and weather. Temporal convolutional networks (TCNs) based on dilated causal convolution have been widely used in time series forecasting. However, two problems weaken the performance of TCNs. One is that in dilated casual convolution, causal convolution leads to the receptive fields of outputs being concentrated in the earlier part of the input sequence, whereas the recent input information will be severely lost. The other is that the distribution shift problem in time series has not been adequately solved. To address the first problem, we propose a subsequence-based dilated convolution method (SDC). By using multiple convolutional filters to convolve elements of neighboring subsequences, the method extracts temporal features from a growing receptive field via a growing subsequence rather than a single element. Ultimately, the receptive field of each output element can cover the whole input sequence. To address the second problem, we propose a difference and compensation method (DCM). The method reduces the discrepancies between and within the input sequences by difference operations and then compensates the outputs for the information lost due to difference operations. Based on SDC and DCM, we further construct a temporal subsequence-based convolutional network with difference (TSCND) for time series forecasting. The experimental results show that TSCND can reduce prediction mean squared error by 7.3% and save runtime, compared with state-of-the-art models and vanilla TCN.
基金the National Natural Science Foundation of China(No.11701103)the Young Top-notch Talent Program of Guangdong Province of China(No.2017GC010379)+4 种基金the Natural Science Foundation of Guangdong Province of China(No.2022A1515012147)the Project of Science and Technology of Guangzhou of China(No.202102020704)the Opening Project of Guangdong Province Key Laboratory of Computational Science at the Sun Yat-sen University of China(2021023)the Science and Technology Development Fund,Macao SAR(File No.0005/2019/A)the University of Macao of China(File Nos.MYRG2020-00035-FST,MYRG2018-00047-FST).
文摘A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.
基金This research was supported by the National Nature Science Foundation of China (No. 41074100) and the Program for NewCentury Excellent Talents in the University of the Ministry of Education of China (No. NCET- 10-0812).
文摘Prestack reverse time migration (RTM) is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference (FD) method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.
基金supported by the National Natural Science Foundation of China(61401389)the Joint Fund of the Ministry of Education of China(6141A02033310)
文摘Micro-satellite cluster enables a whole new class of missions for communications, remote sensing, and scientific research for both civilian and military purposes. Synchronizing the time of the satellites in a cluster is important for both cluster sensing capabilities and its autonomous operating. However, the existing time synchronization methods are not suitable for microsatellite cluster, because it requires too many human interventions and occupies too much ground control resource. Although, data post-process may realize the equivalent time synchronization, it requires processing time and powerful computing ability on the ground, which cannot be implemented by cluster itself. In order to autonomously establish and maintain the time benchmark in a cluster, we propose a compact time difference compensation system(TDCS), which is a kind of time control loop that dynamically adjusts the satellite reference frequency according to the time difference. Consequently, the time synchronization in the cluster can be autonomously achieved on-orbit by synchronizing the clock of other satellites to a chosen one's. The experimental result shows that the standard deviation of time synchronization is about 102 ps when the carrier to noise ratio(CNR) is 95 d BHz, and the standard deviation of corresponding frequency difference is approximately0.36 Hz.
基金supported by the National Natural Science Foundation of China(61401389).
文摘The advancement of small satellites is promoting the development of distributed satellite systems,and for the latter,it is essential to coordinate the spatial and temporal relations between mutually visible satellites.By now,dual one-way ranging(DOWR)and two-way time transfer(TWTT)are generally integrated in the same software and hardware system to meet the limitations of small satellites in terms of size,weight and power(SWaP)consumption.However,studies show that pseudo-noise regenerative ranging(PNRR)performs better than DOWR if some advanced implementation technologies are employed.Besides,PNRR has no requirement on time synchronization.To apply PNRR to small satellites,and meanwhile,meet the demand for time difference measurement,we propose the round-way time difference measurement,which can be combined with PNRR to form a new integrated system without exceeding the limits of SWaP.The new integrated system can provide distributed small satellite systems with on-orbit high-accuracy and high-precision distance measurement and time difference measurement in real time.Experimental results show that the precision of ranging is about 1.94 cm,and that of time difference measurement is about 78.4 ps,at the signal to noise ratio of 80 dBHz.
基金supported by National Natural Science Foundation of China (Grant No. 50675106)
文摘Magneto-rhelological(MR) dampers are devices that employ rheological fluids to modify their mechanical properties. Their mechanical simplicity, high dynamic range, lower power requirements, large force capacity, robustness and safe maimer of operation in cases of failure have made them attractive devices for semi-active real-time control in civil, aerospace and automotive applications. Time response characteristic is one of the most important technical performances of MR dampers, and response time directly affects the control frequency, application range and the actual effect of MR dampers. In this study, one kind of finite difference solution for predicting the response time ofmagneto-rheological dampers from "off-state" to "on-state" is put forward. A laminar flow model is used to describe the flow in the MR valve, and a bi-viscous fluid flow model is utilized to describe the relationship of shear stress and shear rate of MR fluid. An explicit difference format is used to diseretize the Novior-Stoks equation, and stability condition of this algorithm is built by Von-Neumann stability criterion. The pressure gradient along the flow duct is solved by a dichotomy algorithm with iteration, and the changing curve of the damping force versus time of MR damper is obtained as well. According to the abovementioned numerical algorithm, the damping forces versus time curves from "off-state" to "on-state" of a cylindrical piston type MR damper are computed. Moreover, the MR damper is installed in a material test system(MTS), the magnetic field in the wire circles of the MR damper is "triggered" when the MR damper is imposed to do a constant speed motion, and the damping force curves are recorded. The comparison between numerical results and experimental results indicates that this finite difference algorithm can be used to estimate the response time delay of MR dampers.
文摘The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.
基金Project supported partly by the Open Research Program in State Key Laboratory of Millimeter Waves of China(Grant No.K200802)partly by the National Natural Science Foundation of China(Grant No.60971122)
文摘The transmission coefficients of electromagnetic (EM) waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain (SO-FDTD) method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.
基金Supported by the Discipline Construction and Teaching Research Fund of LUTcte(20140089)
文摘Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.
文摘In order to improve the accuracy and engineering feasibility of four-Satellite localization system, the frequency difference measurement is introduced to the four-Satellite TDOA (Time Difference of Arrival) localization algorithm. The TDOA/FDOA (Frequency Difference of Arrival) localization algorithm is used to optimize the GDOP (geometric dilution of precision) of four-Satellite localization. The simulation results show that the absolute position measurement accuracy has little influence on TDOA/FDOA localization accuracy as compared with TDOA localization. Under the same conditions, TDOA/FDOA localization has better accuracy and its GDOP shows more uniform distribution in diamond configuration case. The localization accuracy of four-Satellite TDOA/FDOA is better than the localization accuracy of four-Satellite TDOA.
基金supported by the National Natural Science Foundation of China(Grant Nos.52271278 and 52111530137)the Natural Science Found of Jiangsu Province(Grant No.BK20221389)the Newton Advanced Fellowships(Grant No.NAF\R1\180304)by the Royal Society.
文摘The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.
基金The Science Council,Chinese Taipei Under Grant No.NSC-99-2221-E-027-029
文摘A structure-dependent explicit method with enhanced stability properties is proposed in this study. In general, the method offers unconditional stability for structural systems except those with a particular instantaneous stiffness hardening behavior. In addition, it is second-order accurate and displays no overshooting in high frequency responses. Numerical experiments reveal that the proposed method saves a substantial amount of computational effort in solving inertial problems where only the low frequency responses are of interest, when compared to a general second-order accurate integration method.
文摘A new absorbing boundary condition (ABC) for frequency dependent finite difference time domain algorithm for the arbitrary dispersive media is presented. The concepts of the digital systems are introduced to the (FD) 2TD method. On the basis of digital filter designing and vector algebra, the absorbing boundary condition under arbitrary angle of incidence are derived. The transient electromagnetic problems in two dimensions and three dimensions are calculated and the validity of the ABC is verified.
基金supported by National Natural Science Foundation of China(Nos.11375195,11075048)the National Magnetic Confinement Fusion Science Program of China(No.2013GB104003)
文摘This paper introduces a real-time high precision measurement of phase difference based on Field Programmable Gate Array(FPGA) technology,which has been successfully applied to laser grating interference measurement and real-time feedback of plasma electron density in HL-2A tokamak.It can track the changes of electron density while setting the starting point of the density curve to zero.In a laboratory test,the measuring accuracy of phase difference is less than 0.1°,the time resolution is 80 ns,and the feedback delay is 180 μs.
文摘The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.
基金This work wasjointlysupported by the National Natural Science Foundation of China(Grant No.40106008) the National Natural Science Fundfor Distinguished Young Scholars(Grant No.40225014)
文摘A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones, respectively, thus leading both the time and spatial derivatives to the second-order accuracy. The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift. Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach, respectively, and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor. From comparison of numerical results with theoretical or experimental ones, it is found that they are in reasonable agreement.
基金supported by the National Natural Science Foundation of China(Grant No.61231003)
文摘Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer(CPML) is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain(ECT-CFDTD) method is used to simulate the wave propagation inside a perfect electric conductor(PEC) waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.
文摘In this paper, we approximate the solution to time-fractional telegraph equation by two kinds of difference methods: the Grünwald formula and Caputo fractional difference.
文摘Analyzed the relation between time delay difference and time delay estimation errors, based on the principles of three-point locating, a reformed threshold method for time delay difference estimation of impulse signals, called as amendment estimation for short, is developed by introducing channel equalization technique to its conventional version, named as direct estimation in this paper, to improve the estimation stability. After inherent relationship between time delay and phase shift of signals is analyzed, an integer period error compensation method utilized the diversities of both contribution share and contribution mode of concerned estimates is proposed under the condition of high precision phase lag estimation. Finally, a cooperative multi-threshold estimation method composed of amendment and direct estimations to process impulse signals with three thresholds is established. In sea trials data tests of passive locating, this method improves the estimation precision of time delay difference efficiently. The experiments verify the theoretical predictions.