We employ the block negative dislocation model to invert the distribution of fault coupling and slip rate deficit on the different segments of the Tanlu(Tancheng-Lujiang) fault zone, according to the GPS horizontal ve...We employ the block negative dislocation model to invert the distribution of fault coupling and slip rate deficit on the different segments of the Tanlu(Tancheng-Lujiang) fault zone, according to the GPS horizontal velocity field from 1991 to 2007(the first phase) and 2013 to 2018(the second phase). By comparing the deformation characteristics results, we discuss the relationship between the deformation characteristics with the M earthquake in Japan. The results showed that the fault coupling rate of the northern section of Tancheng in the second phase reduced compared with that in the first phase. However, the results of the two phases showed that the northern section of Juxian still has a high coupling rate, a deep blocking depth, and a dextral compressive deficit, which is the enrapture section of the 1668 Tancheng earthquake. At the same time, the area strain results show that the strain rate of the central and eastern regions of the second phase is obviously enhanced compared with that of the first phase. The occurrence of the great earthquake in Japan has played a specific role in alleviating the strain accumulation in the middle and south sections of the Tanlu fault zone. The results of the maximum shear strain show that the shear strain in the middle section of the Tanlu fault zone in the second phase is weaker than that in the first phase, and the maximum shear strain in the southern section is stronger than that in the first phase. The fault coupling coefficient of the south Sihong to Jiashan section is high, and it is also the unruptured section of historical earthquakes. At the same time, small earthquakes in this area are not active and accumulate stress easily, so the future earthquake risk deserves attention.展开更多
<div style="text-align:justify;"> In order to speed up the global optimization-based mesh smoothing, an enhanced steepest descent method is presented in the paper. Numerical experiment results show tha...<div style="text-align:justify;"> In order to speed up the global optimization-based mesh smoothing, an enhanced steepest descent method is presented in the paper. Numerical experiment results show that the method performs better than the steepest descent method in the global smoothing. We also presented a physically-based interpretation to explain why the method works better than the steepest descent method. </div>展开更多
The classical variational inequality problem with a Lipschitzian and strongly monotone operator on a nonempty closed convex subset in a real Hilbert space is studied. A new three-step relaxed hybrid steepest-descent m...The classical variational inequality problem with a Lipschitzian and strongly monotone operator on a nonempty closed convex subset in a real Hilbert space is studied. A new three-step relaxed hybrid steepest-descent method for this class of variational inequalities is introduced. Strong convergence of this method is established under suitable assumptions imposed on the algorithm parameters.展开更多
Vertical records are critically important when determining the rupture model of an earthquake, especially a thrust earthquake. Due to the relatively low fitness level of near-field vertical displacements, the precisio...Vertical records are critically important when determining the rupture model of an earthquake, especially a thrust earthquake. Due to the relatively low fitness level of near-field vertical displacements, the precision of previous rupture models is relatively low, and the seismic hazard evaluated thereafter should be further updated. In this study, we applied three-component displacement records from GPS stations in and around the source region of the 2013 MW6.6 Lushan earthquake to re-investigate the rupture model.To improve the resolution of the rupture model, records from both continuous and campaign GPS stations were gathered, and secular deformations of the GPS movements were removed from the records of the campaign stations to ensure their reliability. The rupture model was derived by the steepest descent method(SDM), which is based on a layered velocity structure. The peak slip value was about 0.75 m, with a seismic moment release of 9.89 × 1018 N·m, which was equivalent to an MW6.6 event. The inferred fault geometry coincided well with the aftershock distribution of the Lushan earthquake. Unlike previous rupture models, a secondary slip asperity existed at a shallow depth and even touched the ground surface. Based on the distribution of the co-seismic ruptures of the Lushan and Wenchuan earthquakes, post-seismic relaxation of the Wenchuan earthquake, and tectonic loading process, we proposed that the seismic hazard is quite high and still needs special attention in the seismic gap between the two earthquakes.展开更多
In this paper, we consider the generalized variational inequality GVI(F, g, C), where F and g are mappings from a Hilbert space into itself and C is the fixed point set of a nonexpansive mapping. We propose two iter...In this paper, we consider the generalized variational inequality GVI(F, g, C), where F and g are mappings from a Hilbert space into itself and C is the fixed point set of a nonexpansive mapping. We propose two iterative algorithms to find approximate solutions of the GVI(F,g, C). Strong convergence results are established and applications to constrained generalized pseudo-inverse are included.展开更多
By extending the classical analysis techniques due to Samokish, Faddeev and Faddee- va, and Longsine and McCormick among others, we prove the convergence of the precon- ditioned steepest descent with implicit deflati...By extending the classical analysis techniques due to Samokish, Faddeev and Faddee- va, and Longsine and McCormick among others, we prove the convergence of the precon- ditioned steepest descent with implicit deflation (PSD-id) method for solving Hermitian- definite generalized eigenvalue problems. Furthermore, we derive a nonasymptotie estimate of the rate of convergence of the PSD-id method. We show that with a proper choice of the shift, the indefinite shift-and-invert preconditioner is a locally accelerated preconditioner, and is asymptotically optimal which leads to superlinear convergence Numerical examples are presented to verify the theoretical results on the convergence behavior of the PSD- id method for solving ill-conditioned Hermitian-definite generalized eigenvalue problems arising from electronic structure calculations. While rigorous and full-scale convergence proofs of preconditioned block steepest descent methods in practical use still largely eludes us, we believe the theoretical results presented in this paper shed light on an improved understanding of the convergence behavior of these block methods.展开更多
Spectral conjugate gradient method is an algorithm obtained by combination of spectral gradient method and conjugate gradient method,which is characterized with global convergence and simplicity of spectral gradient m...Spectral conjugate gradient method is an algorithm obtained by combination of spectral gradient method and conjugate gradient method,which is characterized with global convergence and simplicity of spectral gradient method,and small storage of conjugate gradient method.Besides,the spectral conjugate gradient method was proved that the search direction at each iteration is a descent direction of objective function even without relying on any line search method.Spectral conjugate gradient method is applied to full waveform inversion for numerical tests on Marmousi model.The authors give a comparison on numerical results obtained by steepest descent method,conjugate gradient method and spectral conjugate gradient method,which shows that the spectral conjugate gradient method is superior to the other two methods.展开更多
There are many phenomena that generate polygonal tessellations on surfaces of 3D objects. One interesting example is the jackfruit, a multiple fruit found in the tropics. A recent study found the best-fit spherical Vo...There are many phenomena that generate polygonal tessellations on surfaces of 3D objects. One interesting example is the jackfruit, a multiple fruit found in the tropics. A recent study found the best-fit spherical Voronoi diagram from a photo of jackfruit skin, but the optimization was relative to the radius of the sphere and the height of the spikes. In this study, we propose a method for adjusting the position of the center of the sphere in addition to these parameters. Experiments were conducted using both ideal and real data. However, convergence with real data has not been confirmed due to relaxation of the convergence condition.展开更多
In this paper,we address interesting soliton resolution,asymptotic stability of N-soliton solutions and the Painleve asymptotics for the Kundu-Eckhaus(KE)equation with nonzero boundary conditions iq_(t)+q_(xx)-2(l|q|^...In this paper,we address interesting soliton resolution,asymptotic stability of N-soliton solutions and the Painleve asymptotics for the Kundu-Eckhaus(KE)equation with nonzero boundary conditions iq_(t)+q_(xx)-2(l|q|^(2)-1)q+4β^(2)(lql^(4)-1)q+4iβ(lql^(2))_(x)q=0,q(x,0)=q_(0)(x)-±1,x→±∞.The key to proving these results is to establish the formulation of a Riemann-Hilbert(RH)problem associated with the above Cauchy problem and find its connection with the RH problem of the defocusing NLS equation.With the■-steepest descent method and the results of the defocusing NLS equation,we find complete leading order approximation formulas for the defocusing KE equation on the whole(x,t)half-plane including soliton resolution and asymptotic stability of N-soliton solutions in a solitonic region,Zakharov-Shabat asymptotics in a solitonless region and the Painlevéasymptotics in two transition regions.展开更多
We introduce the fractional-order global optimal backpropagation machine,which is trained by an improved fractionalorder steepest descent method(FSDM).This is a fractional-order backpropagation neural network(FBPNN),a...We introduce the fractional-order global optimal backpropagation machine,which is trained by an improved fractionalorder steepest descent method(FSDM).This is a fractional-order backpropagation neural network(FBPNN),a state-of-the-art fractional-order branch of the family of backpropagation neural networks(BPNNs),different from the majority of the previous classic first-order BPNNs which are trained by the traditional first-order steepest descent method.The reverse incremental search of the proposed FBPNN is in the negative directions of the approximate fractional-order partial derivatives of the square error.First,the theoretical concept of an FBPNN trained by an improved FSDM is described mathematically.Then,the mathematical proof of fractional-order global optimal convergence,an assumption of the structure,and fractional-order multi-scale global optimization of the FBPNN are analyzed in detail.Finally,we perform three(types of)experiments to compare the performances of an FBPNN and a classic first-order BPNN,i.e.,example function approximation,fractional-order multi-scale global optimization,and comparison of global search and error fitting abilities with real data.The higher optimal search ability of an FBPNN to determine the global optimal solution is the major advantage that makes the FBPNN superior to a classic first-order BPNN.展开更多
基金supported by the National Natural Science Foundation of China (Grand number 41802224)the Youth Program of Seismological Science and Technology Spark Program of China Earthquake Administration (Grand No. XH23019YC)the Joint Open Fund of National Geophysical Observation and Research Station in Mengcheng, Anhui Province (Grand No. MENGO-202114)。
文摘We employ the block negative dislocation model to invert the distribution of fault coupling and slip rate deficit on the different segments of the Tanlu(Tancheng-Lujiang) fault zone, according to the GPS horizontal velocity field from 1991 to 2007(the first phase) and 2013 to 2018(the second phase). By comparing the deformation characteristics results, we discuss the relationship between the deformation characteristics with the M earthquake in Japan. The results showed that the fault coupling rate of the northern section of Tancheng in the second phase reduced compared with that in the first phase. However, the results of the two phases showed that the northern section of Juxian still has a high coupling rate, a deep blocking depth, and a dextral compressive deficit, which is the enrapture section of the 1668 Tancheng earthquake. At the same time, the area strain results show that the strain rate of the central and eastern regions of the second phase is obviously enhanced compared with that of the first phase. The occurrence of the great earthquake in Japan has played a specific role in alleviating the strain accumulation in the middle and south sections of the Tanlu fault zone. The results of the maximum shear strain show that the shear strain in the middle section of the Tanlu fault zone in the second phase is weaker than that in the first phase, and the maximum shear strain in the southern section is stronger than that in the first phase. The fault coupling coefficient of the south Sihong to Jiashan section is high, and it is also the unruptured section of historical earthquakes. At the same time, small earthquakes in this area are not active and accumulate stress easily, so the future earthquake risk deserves attention.
文摘<div style="text-align:justify;"> In order to speed up the global optimization-based mesh smoothing, an enhanced steepest descent method is presented in the paper. Numerical experiment results show that the method performs better than the steepest descent method in the global smoothing. We also presented a physically-based interpretation to explain why the method works better than the steepest descent method. </div>
基金Project supported by the Key Science Foundation of Education Department of Sichuan Province of China (No.2003A081)Sichuan Province Leading Academic Discipline Project (No.SZD0406)
文摘The classical variational inequality problem with a Lipschitzian and strongly monotone operator on a nonempty closed convex subset in a real Hilbert space is studied. A new three-step relaxed hybrid steepest-descent method for this class of variational inequalities is introduced. Strong convergence of this method is established under suitable assumptions imposed on the algorithm parameters.
基金supported by the grant from the National Sichuan-Yunnan earthquake prediction experimental field of CEA (grant No. 2016CESE0204)
文摘Vertical records are critically important when determining the rupture model of an earthquake, especially a thrust earthquake. Due to the relatively low fitness level of near-field vertical displacements, the precision of previous rupture models is relatively low, and the seismic hazard evaluated thereafter should be further updated. In this study, we applied three-component displacement records from GPS stations in and around the source region of the 2013 MW6.6 Lushan earthquake to re-investigate the rupture model.To improve the resolution of the rupture model, records from both continuous and campaign GPS stations were gathered, and secular deformations of the GPS movements were removed from the records of the campaign stations to ensure their reliability. The rupture model was derived by the steepest descent method(SDM), which is based on a layered velocity structure. The peak slip value was about 0.75 m, with a seismic moment release of 9.89 × 1018 N·m, which was equivalent to an MW6.6 event. The inferred fault geometry coincided well with the aftershock distribution of the Lushan earthquake. Unlike previous rupture models, a secondary slip asperity existed at a shallow depth and even touched the ground surface. Based on the distribution of the co-seismic ruptures of the Lushan and Wenchuan earthquakes, post-seismic relaxation of the Wenchuan earthquake, and tectonic loading process, we proposed that the seismic hazard is quite high and still needs special attention in the seismic gap between the two earthquakes.
文摘In this paper, we consider the generalized variational inequality GVI(F, g, C), where F and g are mappings from a Hilbert space into itself and C is the fixed point set of a nonexpansive mapping. We propose two iterative algorithms to find approximate solutions of the GVI(F,g, C). Strong convergence results are established and applications to constrained generalized pseudo-inverse are included.
文摘By extending the classical analysis techniques due to Samokish, Faddeev and Faddee- va, and Longsine and McCormick among others, we prove the convergence of the precon- ditioned steepest descent with implicit deflation (PSD-id) method for solving Hermitian- definite generalized eigenvalue problems. Furthermore, we derive a nonasymptotie estimate of the rate of convergence of the PSD-id method. We show that with a proper choice of the shift, the indefinite shift-and-invert preconditioner is a locally accelerated preconditioner, and is asymptotically optimal which leads to superlinear convergence Numerical examples are presented to verify the theoretical results on the convergence behavior of the PSD- id method for solving ill-conditioned Hermitian-definite generalized eigenvalue problems arising from electronic structure calculations. While rigorous and full-scale convergence proofs of preconditioned block steepest descent methods in practical use still largely eludes us, we believe the theoretical results presented in this paper shed light on an improved understanding of the convergence behavior of these block methods.
文摘Spectral conjugate gradient method is an algorithm obtained by combination of spectral gradient method and conjugate gradient method,which is characterized with global convergence and simplicity of spectral gradient method,and small storage of conjugate gradient method.Besides,the spectral conjugate gradient method was proved that the search direction at each iteration is a descent direction of objective function even without relying on any line search method.Spectral conjugate gradient method is applied to full waveform inversion for numerical tests on Marmousi model.The authors give a comparison on numerical results obtained by steepest descent method,conjugate gradient method and spectral conjugate gradient method,which shows that the spectral conjugate gradient method is superior to the other two methods.
基金Partly Supported by the Grant-in-Aid for Basic Research of MEXT(No.24360039)
文摘There are many phenomena that generate polygonal tessellations on surfaces of 3D objects. One interesting example is the jackfruit, a multiple fruit found in the tropics. A recent study found the best-fit spherical Voronoi diagram from a photo of jackfruit skin, but the optimization was relative to the radius of the sphere and the height of the spikes. In this study, we propose a method for adjusting the position of the center of the sphere in addition to these parameters. Experiments were conducted using both ideal and real data. However, convergence with real data has not been confirmed due to relaxation of the convergence condition.
基金supported by the National Science Foundation of China(Grant No.12271104,51879045)。
文摘In this paper,we address interesting soliton resolution,asymptotic stability of N-soliton solutions and the Painleve asymptotics for the Kundu-Eckhaus(KE)equation with nonzero boundary conditions iq_(t)+q_(xx)-2(l|q|^(2)-1)q+4β^(2)(lql^(4)-1)q+4iβ(lql^(2))_(x)q=0,q(x,0)=q_(0)(x)-±1,x→±∞.The key to proving these results is to establish the formulation of a Riemann-Hilbert(RH)problem associated with the above Cauchy problem and find its connection with the RH problem of the defocusing NLS equation.With the■-steepest descent method and the results of the defocusing NLS equation,we find complete leading order approximation formulas for the defocusing KE equation on the whole(x,t)half-plane including soliton resolution and asymptotic stability of N-soliton solutions in a solitonic region,Zakharov-Shabat asymptotics in a solitonless region and the Painlevéasymptotics in two transition regions.
基金Project supported by the National Key Research and Development Program of China(No.2018YFC0830300)the National Natural Science Foundation of China(No.61571312)。
文摘We introduce the fractional-order global optimal backpropagation machine,which is trained by an improved fractionalorder steepest descent method(FSDM).This is a fractional-order backpropagation neural network(FBPNN),a state-of-the-art fractional-order branch of the family of backpropagation neural networks(BPNNs),different from the majority of the previous classic first-order BPNNs which are trained by the traditional first-order steepest descent method.The reverse incremental search of the proposed FBPNN is in the negative directions of the approximate fractional-order partial derivatives of the square error.First,the theoretical concept of an FBPNN trained by an improved FSDM is described mathematically.Then,the mathematical proof of fractional-order global optimal convergence,an assumption of the structure,and fractional-order multi-scale global optimization of the FBPNN are analyzed in detail.Finally,we perform three(types of)experiments to compare the performances of an FBPNN and a classic first-order BPNN,i.e.,example function approximation,fractional-order multi-scale global optimization,and comparison of global search and error fitting abilities with real data.The higher optimal search ability of an FBPNN to determine the global optimal solution is the major advantage that makes the FBPNN superior to a classic first-order BPNN.