In this paper,we focus on inferring graph Laplacian matrix from the spatiotemporal signal which is defined as“time-vertex signal”.To realize this,we first represent the signals on a joint graph which is the Cartesia...In this paper,we focus on inferring graph Laplacian matrix from the spatiotemporal signal which is defined as“time-vertex signal”.To realize this,we first represent the signals on a joint graph which is the Cartesian product graph of the time-and vertex-graphs.By assuming the signals follow a Gaussian prior distribution on the joint graph,a meaningful representation that promotes the smoothness property of the joint graph signal is derived.Furthermore,by decoupling the joint graph,the graph learning framework is formulated as a joint optimization problem which includes signal denoising,timeand vertex-graphs learning together.Specifically,two algorithms are proposed to solve the optimization problem,where the discrete second-order difference operator with reversed sign(DSODO)in the time domain is used as the time-graph Laplacian operator to recover the signal and infer a vertex-graph in the first algorithm,and the time-graph,as well as the vertex-graph,is estimated by the other algorithm.Experiments on both synthetic and real-world datasets demonstrate that the proposed algorithms can effectively infer meaningful time-and vertex-graphs from noisy and incomplete data.展开更多
The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding freque...The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...展开更多
Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions...Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.展开更多
In the present paper, the random interfacial waves in N-layer density-stratified fluids moving at different steady uniform speeds are researched by using an expansion technique, and the second-order asymptotic solutio...In the present paper, the random interfacial waves in N-layer density-stratified fluids moving at different steady uniform speeds are researched by using an expansion technique, and the second-order asymptotic solutions of the random displacements of the density interfaces and the associated velocity potentials in N-layer fluid are presented based on the small amplitude wave theory. The obtained results indicate that the wave-wave second-order nonlinear interactions of the wave components and the second-order nonlinear interactions between the waves and currents are described. As expected, the solutions include those derived by Chen （2006） as a special case where the steady uniform currents of the N-layer fluids are taken as zero, and the solutions also reduce to those obtained by Song （2005） for second-order solutions for random interfacial waves with steady uniform currents if N = 2.展开更多
In the graph signal processing(GSP)framework,distributed algorithms are highly desirable in processing signals defined on large-scale networks.However,in most existing distributed algorithms,all nodes homogeneously pe...In the graph signal processing(GSP)framework,distributed algorithms are highly desirable in processing signals defined on large-scale networks.However,in most existing distributed algorithms,all nodes homogeneously perform the local computation,which calls for heavy computational and communication costs.Moreover,in many real-world networks,such as those with straggling nodes,the homogeneous manner may result in serious delay or even failure.To this end,we propose active network decomposition algorithms to select non-straggling nodes(normal nodes)that perform the main computation and communication across the network.To accommodate the decomposition in different kinds of networks,two different approaches are developed,one is centralized decomposition that leverages the adjacency of the network and the other is distributed decomposition that employs the indicator message transmission between neighboring nodes,which constitutes the main contribution of this paper.By incorporating the active decomposition scheme,a distributed Newton method is employed to solve the least squares problem in GSP,where the Hessian inverse is approximately evaluated by patching a series of inverses of local Hessian matrices each of which is governed by one normal node.The proposed algorithm inherits the fast convergence of the second-order algorithms while maintains low computational and communication cost.Numerical examples demonstrate the effectiveness of the proposed algorithm.展开更多
This paper investigates the optimal control problem of spacecraft reorientation subject to attitude forbidden constraints,angular velocity saturation and actuator saturation simultaneously.A second-order cone programm...This paper investigates the optimal control problem of spacecraft reorientation subject to attitude forbidden constraints,angular velocity saturation and actuator saturation simultaneously.A second-order cone programming(SOCP)technology is developed to solve the strong nonlinear and non-convex control problem in real time.Specifically,the nonlinear attitude kinematic and dynamic are transformed and relaxed to a standard affine system,and linearization and L1 penalty technique are adopted to convexify non-convex inequality constraints.With the proposed quadratic performance index of angular velocity,the optimal control solution is obtained with high accuracy using the successive SOCP algorithm.Finally,the effectiveness of the algorithm is validated by numerical simulation.展开更多
The optimal planning and design of an integrated energy system(IES)is of great significance to facilitate distributed renewable energy(DRE)technology and improve the overall energy efficiency of the energy system.With...The optimal planning and design of an integrated energy system(IES)is of great significance to facilitate distributed renewable energy(DRE)technology and improve the overall energy efficiency of the energy system.With the increased penetration of distributed generation(DG),the power supply and load sides of an IES present more increased levels of uncertainties.Demand response(DR)and the energy storage system(ESS)serve as important means to shift energy supply and use across time to counter the indeterminate variations.However,the current IES planning methods are unable to effectively deal with the uncertainties of DREs and loads,and to optimize the operations of DG-DR-ESS due to the enormous possible combinations.In this paper,a new method for the optimal planning and design of an integrated energy system has been introduced and verified.The new method consists of three integrated elements.First,the method of the probability scenario has been used to model the uncertainties of the DREs and loads so as to better characterize the impact of uncertainty on the planning and design of the IES.Secondly,the optimal operation of the IES under different probability scenarios is ensured using the second-order cone optimization for quick solutions due to the simplicity of this sub-problem,serving as the bottom-level optimization.Thirdly,the optimal planning and design of IES through optimal sizing of the power generating components and ESS are performed using a special meta-model based global optimization method due to the complex,black-box,and computation intensive nature of this top-level optimization in a nested,bi-level global optimization problem.The combined approach takes full account of the interrelated operations of DG-DR-ESS under different design configurations to support a better optimal planning and design of the IES.The simulation has been carried out on an IES system modified from the IEEE 33-node distribution system.The simulation results show that the proposed method and model are effective.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.61966007)Key Laboratory of Cognitive Radio and Information Processing,Ministry of Education(No.CRKL180106,No.CRKL180201)+1 种基金Guangxi Key Laboratory of Wireless Wideband Communication and Signal Processing,Guilin University of Electronic Technology(No.GXKL06180107,No.GXKL06190117)Guangxi Colleges and Universities Key Laboratory of Satellite Navigation and Position Sensing.
文摘In this paper,we focus on inferring graph Laplacian matrix from the spatiotemporal signal which is defined as“time-vertex signal”.To realize this,we first represent the signals on a joint graph which is the Cartesian product graph of the time-and vertex-graphs.By assuming the signals follow a Gaussian prior distribution on the joint graph,a meaningful representation that promotes the smoothness property of the joint graph signal is derived.Furthermore,by decoupling the joint graph,the graph learning framework is formulated as a joint optimization problem which includes signal denoising,timeand vertex-graphs learning together.Specifically,two algorithms are proposed to solve the optimization problem,where the discrete second-order difference operator with reversed sign(DSODO)in the time domain is used as the time-graph Laplacian operator to recover the signal and infer a vertex-graph in the first algorithm,and the time-graph,as well as the vertex-graph,is estimated by the other algorithm.Experiments on both synthetic and real-world datasets demonstrate that the proposed algorithms can effectively infer meaningful time-and vertex-graphs from noisy and incomplete data.
基金Supported by the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20060056036)
文摘The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...
基金supported by the National Natural Science Foundation of China(11001037,11102037 and 11290143)the Fundamental Research Funds for the Central Universities
文摘Isoparametric quadrilateral elements are widely used in the finite element method, but the accuracy of the isoparametric quadrilateral elements will drop obviously deteriorate due to mesh distortions. Spline functions have some properties of simplicity and conformality. Two 8-node quadrilateral elements have been developed using the trian- gular area coordinates and the B-net method, which can ex- actly model the quadratic field for both convex and concave quadrangles. Some appropriate examples are employed to evaluate the performance of the proposed elements. The nu- merical results show that the two spline elements can obtain solutions which are highly accurate and insensitive to mesh distortions.
基金supported by the Natural Science Foundation of Inner Mongolia,China (Grant No 200711020116)Open Fund of the Key Laboratory of Ocean Circulation and Waves,Chinese Academy of Sciences (Grant No KLOCAW0805)+1 种基金the Key Program of the Scientific Research Plan of Inner Mongolia University of Technology,China (Grant No ZD200608)the National Science Fund for Distinguished Young Scholars of China (Grant No 40425015)
文摘In the present paper, the random interfacial waves in N-layer density-stratified fluids moving at different steady uniform speeds are researched by using an expansion technique, and the second-order asymptotic solutions of the random displacements of the density interfaces and the associated velocity potentials in N-layer fluid are presented based on the small amplitude wave theory. The obtained results indicate that the wave-wave second-order nonlinear interactions of the wave components and the second-order nonlinear interactions between the waves and currents are described. As expected, the solutions include those derived by Chen （2006） as a special case where the steady uniform currents of the N-layer fluids are taken as zero, and the solutions also reduce to those obtained by Song （2005） for second-order solutions for random interfacial waves with steady uniform currents if N = 2.
基金supported by National Natural Science Foundation of China(Grant No.61761011)Natural Science Foundation of Guangxi(Grant No.2020GXNSFBA297078).
文摘In the graph signal processing(GSP)framework,distributed algorithms are highly desirable in processing signals defined on large-scale networks.However,in most existing distributed algorithms,all nodes homogeneously perform the local computation,which calls for heavy computational and communication costs.Moreover,in many real-world networks,such as those with straggling nodes,the homogeneous manner may result in serious delay or even failure.To this end,we propose active network decomposition algorithms to select non-straggling nodes(normal nodes)that perform the main computation and communication across the network.To accommodate the decomposition in different kinds of networks,two different approaches are developed,one is centralized decomposition that leverages the adjacency of the network and the other is distributed decomposition that employs the indicator message transmission between neighboring nodes,which constitutes the main contribution of this paper.By incorporating the active decomposition scheme,a distributed Newton method is employed to solve the least squares problem in GSP,where the Hessian inverse is approximately evaluated by patching a series of inverses of local Hessian matrices each of which is governed by one normal node.The proposed algorithm inherits the fast convergence of the second-order algorithms while maintains low computational and communication cost.Numerical examples demonstrate the effectiveness of the proposed algorithm.
基金This work was supported by the National Natural Science Foundation of China(Nos.61960206011,61633003)the Beijing Natural Science Foundation(No.JQ19017)。
文摘This paper investigates the optimal control problem of spacecraft reorientation subject to attitude forbidden constraints,angular velocity saturation and actuator saturation simultaneously.A second-order cone programming(SOCP)technology is developed to solve the strong nonlinear and non-convex control problem in real time.Specifically,the nonlinear attitude kinematic and dynamic are transformed and relaxed to a standard affine system,and linearization and L1 penalty technique are adopted to convexify non-convex inequality constraints.With the proposed quadratic performance index of angular velocity,the optimal control solution is obtained with high accuracy using the successive SOCP algorithm.Finally,the effectiveness of the algorithm is validated by numerical simulation.
基金This work was supported by the National Natural Science Foundation of China(51607170)the Key Front Science Project of Chinese Academy of Sciences(QYZDB-SSW-JSC024)the International Collaboration Programs of the Chinese Academy of Sciences and the Foreign Expert Affairs.
文摘The optimal planning and design of an integrated energy system(IES)is of great significance to facilitate distributed renewable energy(DRE)technology and improve the overall energy efficiency of the energy system.With the increased penetration of distributed generation(DG),the power supply and load sides of an IES present more increased levels of uncertainties.Demand response(DR)and the energy storage system(ESS)serve as important means to shift energy supply and use across time to counter the indeterminate variations.However,the current IES planning methods are unable to effectively deal with the uncertainties of DREs and loads,and to optimize the operations of DG-DR-ESS due to the enormous possible combinations.In this paper,a new method for the optimal planning and design of an integrated energy system has been introduced and verified.The new method consists of three integrated elements.First,the method of the probability scenario has been used to model the uncertainties of the DREs and loads so as to better characterize the impact of uncertainty on the planning and design of the IES.Secondly,the optimal operation of the IES under different probability scenarios is ensured using the second-order cone optimization for quick solutions due to the simplicity of this sub-problem,serving as the bottom-level optimization.Thirdly,the optimal planning and design of IES through optimal sizing of the power generating components and ESS are performed using a special meta-model based global optimization method due to the complex,black-box,and computation intensive nature of this top-level optimization in a nested,bi-level global optimization problem.The combined approach takes full account of the interrelated operations of DG-DR-ESS under different design configurations to support a better optimal planning and design of the IES.The simulation has been carried out on an IES system modified from the IEEE 33-node distribution system.The simulation results show that the proposed method and model are effective.
