In this paper, two different methods are used to study the cyclic structure solution and the optimal approximation of the quaternion Stein equation AXB - X = F . Firstly, the matrix equation equivalent to the ta...In this paper, two different methods are used to study the cyclic structure solution and the optimal approximation of the quaternion Stein equation AXB - X = F . Firstly, the matrix equation equivalent to the target structure matrix is constructed by using the complex decomposition of the quaternion matrix, to obtain the necessary and sufficient conditions for the existence of the cyclic solution of the equation and the expression of the general solution. Secondly, the Stein equation is converted into the Sylvester equation by adding the necessary parameters, and the condition for the existence of a cyclic solution and the expression of the equation’s solution are then obtained by using the real decomposition of the quaternion matrix and the Kronecker product of the matrix. At the same time, under the condition that the solution set is non-empty, the optimal approximation solution to the given quaternion circulant matrix is obtained by using the property of Frobenius norm property. Numerical examples are given to verify the correctness of the theoretical results and the feasibility of the proposed method. .展开更多
This paper offers an extensive overview of the utilization of sequential approximate optimization approaches in the context of numerically simulated large-scale continuum structures.These structures,commonly encounter...This paper offers an extensive overview of the utilization of sequential approximate optimization approaches in the context of numerically simulated large-scale continuum structures.These structures,commonly encountered in engineering applications,often involve complex objective and constraint functions that cannot be readily expressed as explicit functions of the design variables.As a result,sequential approximation techniques have emerged as the preferred strategy for addressing a wide array of topology optimization challenges.Over the past several decades,topology optimization methods have been advanced remarkably and successfully applied to solve engineering problems incorporating diverse physical backgrounds.In comparison to the large-scale equation solution,sensitivity analysis,graphics post-processing,etc.,the progress of the sequential approximation functions and their corresponding optimizersmake sluggish progress.Researchers,particularly novices,pay special attention to their difficulties with a particular problem.Thus,this paper provides an overview of sequential approximation functions,related literature on topology optimization methods,and their applications.Starting from optimality criteria and sequential linear programming,the other sequential approximate optimizations are introduced by employing Taylor expansion and intervening variables.In addition,recent advancements have led to the emergence of approaches such as Augmented Lagrange,sequential approximate integer,and non-gradient approximation are also introduced.By highlighting real-world applications and case studies,the paper not only demonstrates the practical relevance of these methods but also underscores the need for continued exploration in this area.Furthermore,to provide a comprehensive overview,this paper offers several novel developments that aim to illuminate potential directions for future research.展开更多
In the research on spatial hearing and virtual auditory space,it is important to effectively model the head-related transfer functions(HRTFs).Based on the analysis of the HRTFs’spectrum and some perspectives of psych...In the research on spatial hearing and virtual auditory space,it is important to effectively model the head-related transfer functions(HRTFs).Based on the analysis of the HRTFs’spectrum and some perspectives of psychoacoustics,this paper applied multiple demes’parallel and real-valued coding genetic algorithm(GA)to approxi-mate the HRTFs’zero-pole model.Using the logarithmic magnitude’s error criterion for the human auditory sense,the results show that the performance of the GA is on the average 39%better than that of the traditional Prony method,and 46%better than that of the Yule-Walker algo-rithm.展开更多
The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward c...The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations.展开更多
Dykstra’s alternating projection algorithm was proposed to treat the problem of finding the projection of a given point onto the intersection of some closed convex sets. In this paper, we first apply Dykstra’s alter...Dykstra’s alternating projection algorithm was proposed to treat the problem of finding the projection of a given point onto the intersection of some closed convex sets. In this paper, we first apply Dykstra’s alternating projection algorithm to compute the optimal approximate symmetric positive semidefinite solution of the matrix equations AXB = E, CXD = F. If we choose the initial iterative matrix X<sub>0</sub> = 0, the least Frobenius norm symmetric positive semidefinite solution of these matrix equations is obtained. A numerical example shows that the new algorithm is feasible and effective.展开更多
This paper estimates an off-policy integral reinforcement learning(IRL) algorithm to obtain the optimal tracking control of unknown chaotic systems. Off-policy IRL can learn the solution of the HJB equation from the...This paper estimates an off-policy integral reinforcement learning(IRL) algorithm to obtain the optimal tracking control of unknown chaotic systems. Off-policy IRL can learn the solution of the HJB equation from the system data generated by an arbitrary control. Moreover, off-policy IRL can be regarded as a direct learning method, which avoids the identification of system dynamics. In this paper, the performance index function is first given based on the system tracking error and control error. For solving the Hamilton–Jacobi–Bellman(HJB) equation, an off-policy IRL algorithm is proposed.It is proven that the iterative control makes the tracking error system asymptotically stable, and the iterative performance index function is convergent. Simulation study demonstrates the effectiveness of the developed tracking control method.展开更多
Uncertainty in determining optimum conjunctive water use lies not only on variability of hydrological cycle and climate but also on lack of adequate data and perfect knowledge about groundwater-surface water system in...Uncertainty in determining optimum conjunctive water use lies not only on variability of hydrological cycle and climate but also on lack of adequate data and perfect knowledge about groundwater-surface water system interactions, errors in historic data and inherent variability of system parameters both in space and time. Simulation-optimization models are used for conjunctive water use management under uncertain conditions. However, direct application of such approach whereby all realizations are considered at every-iteration of the optimization process leads to a highly computational time-consuming optimization problem as the number of realizations increases. Hence, this study proposes a novel approach—a Retrospective Optimization Approximation (ROA) approach. In this approach, a simulation model was used to determine aquifer system responses (draw-downs) which were assembled as response matrices and incorporated in the optimization model (procedure) as coefficients in the constraints. The sample optimization sub-problems generated, were solved and analyzed through ROA-Active-Set procedure implemented under MATLAB code. The ROA-Active Set procedure solves and evaluates a sequence of sample path optimization sub-problems in an increasing number of realizations. The methodology was applied to a real-world conjunctive water use management problem found in Great Letaba River basin, South Africa. In the River basin, surface water source contributes 87% of the existing un-optimized total conjunctive water use withdrawal rate (6512.04 m<sup>3</sup>/day) and the remaining 13% is contributed by groundwater source. Through ROA approach, results indicate that the optimum percentages contribution of the surface and subsurface sources to the total water demand are 58% and 42% respectively. This implies that the existing percentage contribution can be increased or reduced by ±29% that is groundwater source can be increased by 29% while the surface water source contribution can be reduced by 29%. This reveals that the existing conjunctive water use practice is unsustainable wherein surface water system is overstressed while groundwater system is under-utilized. Through k-means sampling technique ROA-Active Set procedure was able to attain a converged maximum expected total optimum conjunctive water use withdrawal rate of 4.35 × 10<sup>4</sup> m<sup>3</sup>/day within a relatively few numbers of iterations (6 to 8 iterations) in about 2.30 Hrs. In conclusion, results demonstrated that ROA approach is capable of managing real-world regional aquifers sustainable conjunctive water use practice under hydro-geological uncertainty conditions.展开更多
In this article, the generalized reflexive solution of matrix equations (AX = B, XC = D) is considered. With special properties of generalized reflexive matrices, the necessary and sufficient conditions for the solv...In this article, the generalized reflexive solution of matrix equations (AX = B, XC = D) is considered. With special properties of generalized reflexive matrices, the necessary and sufficient conditions for the solvability and the general expression of the solution are obtained. Moreover, the related optimal approximation problem to a given matrix over the solution set is solved.展开更多
Appealing to the Clifford analysis and matching pursuits, we study the adaptive decompositions of functions of several variables of finite energy under the dictionaries consisting of shifted Cauchy kernels. This is a ...Appealing to the Clifford analysis and matching pursuits, we study the adaptive decompositions of functions of several variables of finite energy under the dictionaries consisting of shifted Cauchy kernels. This is a realization of matching pursuits among shifted Cauchy kernels in higher-dimensional spaces. It offers a method to process signals in arbitrary dimensions.展开更多
In this paper we first consider the existence and the general form of solution to the following generalized inverse eigenvalue problem(GIEP): given a set of n-dimension complex vectors {x j}m j=1 and a set of co...In this paper we first consider the existence and the general form of solution to the following generalized inverse eigenvalue problem(GIEP): given a set of n-dimension complex vectors {x j}m j=1 and a set of complex numbers {λ j}m j=1, find two n×n centrohermitian matrices A,B such that {x j}m j=1 and {λ j}m j=1 are the generalized eigenvectors and generalized eigenvalues of Ax=λBx, respectively. We then discuss the optimal approximation problem for the GIEP. More concretely, given two arbitrary matrices, , ∈C n×n, we find two matrices A and B such that the matrix (A*,B*) is closest to (,) in the Frobenius norm, where the matrix (A*,B*) is the solution to the GIEP. We show that the expression of the solution of the optimal approximation is unique and derive the expression for it.展开更多
The matrix equations (AX, XBH)=(C, DH) have been widely used in structural design, parameter identification, linear optimal control, and so on. But few researches studied the reflexive solutions. A new approach for th...The matrix equations (AX, XBH)=(C, DH) have been widely used in structural design, parameter identification, linear optimal control, and so on. But few researches studied the reflexive solutions. A new approach for the reflexive solutions to the matrix equations was proposed. By applying the canonical correlation decomposition (CCD) of matrix pairs, the necessary and sufficient conditions for the existence and the general expression for the reflexive solutions of the matrix equations (AX, XBH)=(C, DH) were established. In addition, by using the methods of space decomposition, the expression of the optimal approximation solution to a given matrix was derived.展开更多
A novel numerical method is presented to update mass and stiffness matrices simultaneously with measured vibration data by means of the combined acceleration and displacement output feedback.By the method,the required...A novel numerical method is presented to update mass and stiffness matrices simultaneously with measured vibration data by means of the combined acceleration and displacement output feedback.By the method,the required displacement and acceleration output feedback gain matrices are determined,and thus the optimal approximation mass matrix and stiffness matrix which satisfy the required orthogonality relation and eigenvalue equation are found.The proposed method is computationally efficient and the updated mass and stiffness matrices are also symmetric and have the compact expressions.The numerical example shows that the proposed method is reliable and attractive.展开更多
The least-square solutions of inverse problem for anti-symmetric and skew-symmetric matrices are studied. In addition, the problem of using anti-symmetric and skew-symmetric matrices to construct the optimal approxima...The least-square solutions of inverse problem for anti-symmetric and skew-symmetric matrices are studied. In addition, the problem of using anti-symmetric and skew-symmetric matrices to construct the optimal approximation to a given matrix is discussed, the necessary and sufficient conditions for the problem are derived, and the expression of the solution is provided. A numerical example is given to show the effectiveness of the proposed method.展开更多
A Lagrangian relaxation(LR) approach was presented which is with machine capacity relaxation and operation precedence relaxation for solving a flexible job shop(FJS) scheduling problem from the steelmaking-refining-co...A Lagrangian relaxation(LR) approach was presented which is with machine capacity relaxation and operation precedence relaxation for solving a flexible job shop(FJS) scheduling problem from the steelmaking-refining-continuous casting process. Unlike the full optimization of LR problems in traditional LR approaches, the machine capacity relaxation is optimized asymptotically, while the precedence relaxation is optimized approximately due to the NP-hard nature of its LR problem. Because the standard subgradient algorithm(SSA) cannot solve the Lagrangian dual(LD) problem within the partial optimization of LR problem, an effective deflected-conditional approximate subgradient level algorithm(DCASLA) was developed, named as Lagrangian relaxation level approach. The efficiency of the DCASLA is enhanced by a deflected-conditional epsilon-subgradient to weaken the possible zigzagging phenomena. Computational results and comparisons show that the proposed methods improve significantly the efficiency of the LR approach and the DCASLA adopting capacity relaxation strategy performs best among eight methods in terms of solution quality and running time.展开更多
The Quantum Approximate Optimization Algorithm(QAOA)is an algorithmic framework for finding approximate solutions to combinatorial optimization problems.It consists of interleaved unitary transformations induced by tw...The Quantum Approximate Optimization Algorithm(QAOA)is an algorithmic framework for finding approximate solutions to combinatorial optimization problems.It consists of interleaved unitary transformations induced by two operators labelled the mixing and problem Hamiltonians.To fit this framework,one needs to transform the original problem into a suitable form and embed it into these two Hamiltonians.In this paper,for the well-known NP-hard Traveling Salesman Problem(TSP),we encode its constraints into the mixing Hamiltonian rather than the conventional approach of adding penalty terms to the problem Hamiltonian.Moreover,we map edges(routes)connecting each pair of cities to qubits,which decreases the search space significantly in comparison to other approaches.As a result,our method can achieve a higher probability for the shortest round-trip route with only half the number of qubits consumed compared to IBM Q’s approach.We argue the formalization approach presented in this paper would lead to a generalized framework for finding,in the context of QAOA,high-quality approximate solutions to NP optimization problems.展开更多
Control systems designed by the principle of matching gives rise to problems of evaluating the peak output. This paper proposes a practical method for computing the peak output of linear time-invariant and non-anticip...Control systems designed by the principle of matching gives rise to problems of evaluating the peak output. This paper proposes a practical method for computing the peak output of linear time-invariant and non-anticipative systems for a class of possible sets that are characterized with many bounding conditions on the two- and/or the infinity-norms of the inputs and their derivatives. The original infinite-dimensional convex optimization problem is approximated as a large-scale convex programme defined in a Euclidean space, which are associated with sparse matrices and thus can be solved efficiently in practice. The numerical results show that the method performs satisfactorily, and that using a possible set with many bounding conditions can help to reduce the design conservatism and thereby yield a better match.展开更多
The multiple knapsack problem denoted by MKP (B,S,m,n) can be defined as fol- lows.A set B of n items and a set Sof m knapsacks are given such thateach item j has a profit pjand weightwj,and each knapsack i has a ca...The multiple knapsack problem denoted by MKP (B,S,m,n) can be defined as fol- lows.A set B of n items and a set Sof m knapsacks are given such thateach item j has a profit pjand weightwj,and each knapsack i has a capacity Ci.The goal is to find a subset of items of maximum profit such that they have a feasible packing in the knapsacks.MKP(B,S,m,n) is strongly NP- Complete and no polynomial- time approximation algorithm can have an approxima- tion ratio better than0 .5 .In the last ten years,semi- definite programming has been empolyed to solve some combinatorial problems successfully.This paper firstly presents a semi- definite re- laxation algorithm (MKPS) for MKP (B,S,m,n) .It is proved that MKPS have a approxima- tion ratio better than 0 .5 for a subclass of MKP (B,S,m,n) with n≤ 1 0 0 ,m≤ 5 and maxnj=1{ wj} minmi=1{ Ci} ≤ 2 3 .展开更多
In this paper,the generalized inverse eigenvalue problem for the(P,Q)-conjugate matrices and the associated approximation problem are discussed by using generalized singular value decomposition(GSVD).Moreover,the ...In this paper,the generalized inverse eigenvalue problem for the(P,Q)-conjugate matrices and the associated approximation problem are discussed by using generalized singular value decomposition(GSVD).Moreover,the least residual problem of the above generalized inverse eigenvalue problem is studied by using the canonical correlation decomposition(CCD).The solutions to these problems are derived.Some numerical examples are given to illustrate the main results.展开更多
In this paper,the Hermitian reflexive(Anti-Hermitian reflexive)least-squares so-lutions of matrix equations(AX = B,XC = D)are considered.With special properties of partitioned matrices and Hermitian reflexive(Ant...In this paper,the Hermitian reflexive(Anti-Hermitian reflexive)least-squares so-lutions of matrix equations(AX = B,XC = D)are considered.With special properties of partitioned matrices and Hermitian reflexive(Anti-Hermitian reflexive)matrices,the general expression of the solution is obtained.Moreover,the related optimal approximation problem to a given matrix over the solution set is considered.展开更多
Let SE denote the least-squares symmetric solution set of the matrix equation A×B = C, where A, B and C are given matrices of suitable size. To find the optimal approximate solution in the set SE to a given matri...Let SE denote the least-squares symmetric solution set of the matrix equation A×B = C, where A, B and C are given matrices of suitable size. To find the optimal approximate solution in the set SE to a given matrix, we give a new feasible method based on the projection theorem, the generalized SVD and the canonical correction decomposition.展开更多
文摘In this paper, two different methods are used to study the cyclic structure solution and the optimal approximation of the quaternion Stein equation AXB - X = F . Firstly, the matrix equation equivalent to the target structure matrix is constructed by using the complex decomposition of the quaternion matrix, to obtain the necessary and sufficient conditions for the existence of the cyclic solution of the equation and the expression of the general solution. Secondly, the Stein equation is converted into the Sylvester equation by adding the necessary parameters, and the condition for the existence of a cyclic solution and the expression of the equation’s solution are then obtained by using the real decomposition of the quaternion matrix and the Kronecker product of the matrix. At the same time, under the condition that the solution set is non-empty, the optimal approximation solution to the given quaternion circulant matrix is obtained by using the property of Frobenius norm property. Numerical examples are given to verify the correctness of the theoretical results and the feasibility of the proposed method. .
基金financially supported by the National Key R&D Program (2022YFB4201302)Guang Dong Basic and Applied Basic Research Foundation (2022A1515240057)the Huaneng Technology Funds (HNKJ20-H88).
文摘This paper offers an extensive overview of the utilization of sequential approximate optimization approaches in the context of numerically simulated large-scale continuum structures.These structures,commonly encountered in engineering applications,often involve complex objective and constraint functions that cannot be readily expressed as explicit functions of the design variables.As a result,sequential approximation techniques have emerged as the preferred strategy for addressing a wide array of topology optimization challenges.Over the past several decades,topology optimization methods have been advanced remarkably and successfully applied to solve engineering problems incorporating diverse physical backgrounds.In comparison to the large-scale equation solution,sensitivity analysis,graphics post-processing,etc.,the progress of the sequential approximation functions and their corresponding optimizersmake sluggish progress.Researchers,particularly novices,pay special attention to their difficulties with a particular problem.Thus,this paper provides an overview of sequential approximation functions,related literature on topology optimization methods,and their applications.Starting from optimality criteria and sequential linear programming,the other sequential approximate optimizations are introduced by employing Taylor expansion and intervening variables.In addition,recent advancements have led to the emergence of approaches such as Augmented Lagrange,sequential approximate integer,and non-gradient approximation are also introduced.By highlighting real-world applications and case studies,the paper not only demonstrates the practical relevance of these methods but also underscores the need for continued exploration in this area.Furthermore,to provide a comprehensive overview,this paper offers several novel developments that aim to illuminate potential directions for future research.
基金supported by the National Basic Research of China(No.2002CB312102)。
文摘In the research on spatial hearing and virtual auditory space,it is important to effectively model the head-related transfer functions(HRTFs).Based on the analysis of the HRTFs’spectrum and some perspectives of psychoacoustics,this paper applied multiple demes’parallel and real-valued coding genetic algorithm(GA)to approxi-mate the HRTFs’zero-pole model.Using the logarithmic magnitude’s error criterion for the human auditory sense,the results show that the performance of the GA is on the average 39%better than that of the traditional Prony method,and 46%better than that of the Yule-Walker algo-rithm.
基金This paper is supported by the National Natural Science Foundation of China (No. 40474047)State Key Laboratory of Geological Processes and Mineral Resources (No. GPMR200654)the Focused Subject Program of Beijing (No. XK104910598).
文摘The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations.
文摘Dykstra’s alternating projection algorithm was proposed to treat the problem of finding the projection of a given point onto the intersection of some closed convex sets. In this paper, we first apply Dykstra’s alternating projection algorithm to compute the optimal approximate symmetric positive semidefinite solution of the matrix equations AXB = E, CXD = F. If we choose the initial iterative matrix X<sub>0</sub> = 0, the least Frobenius norm symmetric positive semidefinite solution of these matrix equations is obtained. A numerical example shows that the new algorithm is feasible and effective.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61304079 and 61374105)the Beijing Natural Science Foundation,China(Grant Nos.4132078 and 4143065)+2 种基金the China Postdoctoral Science Foundation(Grant No.2013M530527)the Fundamental Research Funds for the Central Universities,China(Grant No.FRF-TP-14-119A2)the Open Research Project from State Key Laboratory of Management and Control for Complex Systems,China(Grant No.20150104)
文摘This paper estimates an off-policy integral reinforcement learning(IRL) algorithm to obtain the optimal tracking control of unknown chaotic systems. Off-policy IRL can learn the solution of the HJB equation from the system data generated by an arbitrary control. Moreover, off-policy IRL can be regarded as a direct learning method, which avoids the identification of system dynamics. In this paper, the performance index function is first given based on the system tracking error and control error. For solving the Hamilton–Jacobi–Bellman(HJB) equation, an off-policy IRL algorithm is proposed.It is proven that the iterative control makes the tracking error system asymptotically stable, and the iterative performance index function is convergent. Simulation study demonstrates the effectiveness of the developed tracking control method.
文摘Uncertainty in determining optimum conjunctive water use lies not only on variability of hydrological cycle and climate but also on lack of adequate data and perfect knowledge about groundwater-surface water system interactions, errors in historic data and inherent variability of system parameters both in space and time. Simulation-optimization models are used for conjunctive water use management under uncertain conditions. However, direct application of such approach whereby all realizations are considered at every-iteration of the optimization process leads to a highly computational time-consuming optimization problem as the number of realizations increases. Hence, this study proposes a novel approach—a Retrospective Optimization Approximation (ROA) approach. In this approach, a simulation model was used to determine aquifer system responses (draw-downs) which were assembled as response matrices and incorporated in the optimization model (procedure) as coefficients in the constraints. The sample optimization sub-problems generated, were solved and analyzed through ROA-Active-Set procedure implemented under MATLAB code. The ROA-Active Set procedure solves and evaluates a sequence of sample path optimization sub-problems in an increasing number of realizations. The methodology was applied to a real-world conjunctive water use management problem found in Great Letaba River basin, South Africa. In the River basin, surface water source contributes 87% of the existing un-optimized total conjunctive water use withdrawal rate (6512.04 m<sup>3</sup>/day) and the remaining 13% is contributed by groundwater source. Through ROA approach, results indicate that the optimum percentages contribution of the surface and subsurface sources to the total water demand are 58% and 42% respectively. This implies that the existing percentage contribution can be increased or reduced by ±29% that is groundwater source can be increased by 29% while the surface water source contribution can be reduced by 29%. This reveals that the existing conjunctive water use practice is unsustainable wherein surface water system is overstressed while groundwater system is under-utilized. Through k-means sampling technique ROA-Active Set procedure was able to attain a converged maximum expected total optimum conjunctive water use withdrawal rate of 4.35 × 10<sup>4</sup> m<sup>3</sup>/day within a relatively few numbers of iterations (6 to 8 iterations) in about 2.30 Hrs. In conclusion, results demonstrated that ROA approach is capable of managing real-world regional aquifers sustainable conjunctive water use practice under hydro-geological uncertainty conditions.
基金supported by National Natural Science Foundation of China (10571047)and by Scientific Research Fund of Hunan Provincial Education Department of China Grant(06C235)+1 种基金by Central South University of Forestry and Technology (06Y017)by Specialized Research Fund for the Doctoral Program of Higher Education (20060532014)
文摘In this article, the generalized reflexive solution of matrix equations (AX = B, XC = D) is considered. With special properties of generalized reflexive matrices, the necessary and sufficient conditions for the solvability and the general expression of the solution are obtained. Moreover, the related optimal approximation problem to a given matrix over the solution set is solved.
基金supported by Macao FDCT(098/2012/A3)Research Grant of the University of Macao(UL017/08-Y4/MAT/QT01/FST)+1 种基金National Natural Science Funds for Young Scholars(10901166)Sun Yat-sen University Operating Costs of Basic ResearchProjects to Cultivate Young Teachers(11lgpy99)
文摘Appealing to the Clifford analysis and matching pursuits, we study the adaptive decompositions of functions of several variables of finite energy under the dictionaries consisting of shifted Cauchy kernels. This is a realization of matching pursuits among shifted Cauchy kernels in higher-dimensional spaces. It offers a method to process signals in arbitrary dimensions.
文摘In this paper we first consider the existence and the general form of solution to the following generalized inverse eigenvalue problem(GIEP): given a set of n-dimension complex vectors {x j}m j=1 and a set of complex numbers {λ j}m j=1, find two n×n centrohermitian matrices A,B such that {x j}m j=1 and {λ j}m j=1 are the generalized eigenvectors and generalized eigenvalues of Ax=λBx, respectively. We then discuss the optimal approximation problem for the GIEP. More concretely, given two arbitrary matrices, , ∈C n×n, we find two matrices A and B such that the matrix (A*,B*) is closest to (,) in the Frobenius norm, where the matrix (A*,B*) is the solution to the GIEP. We show that the expression of the solution of the optimal approximation is unique and derive the expression for it.
基金National Natural Science Foundation of China ( No. 60875007)
文摘The matrix equations (AX, XBH)=(C, DH) have been widely used in structural design, parameter identification, linear optimal control, and so on. But few researches studied the reflexive solutions. A new approach for the reflexive solutions to the matrix equations was proposed. By applying the canonical correlation decomposition (CCD) of matrix pairs, the necessary and sufficient conditions for the existence and the general expression for the reflexive solutions of the matrix equations (AX, XBH)=(C, DH) were established. In addition, by using the methods of space decomposition, the expression of the optimal approximation solution to a given matrix was derived.
文摘A novel numerical method is presented to update mass and stiffness matrices simultaneously with measured vibration data by means of the combined acceleration and displacement output feedback.By the method,the required displacement and acceleration output feedback gain matrices are determined,and thus the optimal approximation mass matrix and stiffness matrix which satisfy the required orthogonality relation and eigenvalue equation are found.The proposed method is computationally efficient and the updated mass and stiffness matrices are also symmetric and have the compact expressions.The numerical example shows that the proposed method is reliable and attractive.
文摘The least-square solutions of inverse problem for anti-symmetric and skew-symmetric matrices are studied. In addition, the problem of using anti-symmetric and skew-symmetric matrices to construct the optimal approximation to a given matrix is discussed, the necessary and sufficient conditions for the problem are derived, and the expression of the solution is provided. A numerical example is given to show the effectiveness of the proposed method.
基金Projects(51435009,51575212,61573249,61371200)supported by the National Natural Science Foundation of ChinaProjects(2015T80798,2014M552040,2014M561250,2015M571328)supported by Postdoctoral Science Foundation of ChinaProject(L2015372)supported by Liaoning Province Education Administration,China
文摘A Lagrangian relaxation(LR) approach was presented which is with machine capacity relaxation and operation precedence relaxation for solving a flexible job shop(FJS) scheduling problem from the steelmaking-refining-continuous casting process. Unlike the full optimization of LR problems in traditional LR approaches, the machine capacity relaxation is optimized asymptotically, while the precedence relaxation is optimized approximately due to the NP-hard nature of its LR problem. Because the standard subgradient algorithm(SSA) cannot solve the Lagrangian dual(LD) problem within the partial optimization of LR problem, an effective deflected-conditional approximate subgradient level algorithm(DCASLA) was developed, named as Lagrangian relaxation level approach. The efficiency of the DCASLA is enhanced by a deflected-conditional epsilon-subgradient to weaken the possible zigzagging phenomena. Computational results and comparisons show that the proposed methods improve significantly the efficiency of the LR approach and the DCASLA adopting capacity relaxation strategy performs best among eight methods in terms of solution quality and running time.
基金This work is supported by the Natural Science Foundation,China(Grant No.61802002)Natural Science Foundation of Anhui Province,China(Grant No.1708085MF162).
文摘The Quantum Approximate Optimization Algorithm(QAOA)is an algorithmic framework for finding approximate solutions to combinatorial optimization problems.It consists of interleaved unitary transformations induced by two operators labelled the mixing and problem Hamiltonians.To fit this framework,one needs to transform the original problem into a suitable form and embed it into these two Hamiltonians.In this paper,for the well-known NP-hard Traveling Salesman Problem(TSP),we encode its constraints into the mixing Hamiltonian rather than the conventional approach of adding penalty terms to the problem Hamiltonian.Moreover,we map edges(routes)connecting each pair of cities to qubits,which decreases the search space significantly in comparison to other approaches.As a result,our method can achieve a higher probability for the shortest round-trip route with only half the number of qubits consumed compared to IBM Q’s approach.We argue the formalization approach presented in this paper would lead to a generalized framework for finding,in the context of QAOA,high-quality approximate solutions to NP optimization problems.
文摘Control systems designed by the principle of matching gives rise to problems of evaluating the peak output. This paper proposes a practical method for computing the peak output of linear time-invariant and non-anticipative systems for a class of possible sets that are characterized with many bounding conditions on the two- and/or the infinity-norms of the inputs and their derivatives. The original infinite-dimensional convex optimization problem is approximated as a large-scale convex programme defined in a Euclidean space, which are associated with sparse matrices and thus can be solved efficiently in practice. The numerical results show that the method performs satisfactorily, and that using a possible set with many bounding conditions can help to reduce the design conservatism and thereby yield a better match.
基金Supported by the National Natural Science Foundation of China(1 9971 0 78)
文摘The multiple knapsack problem denoted by MKP (B,S,m,n) can be defined as fol- lows.A set B of n items and a set Sof m knapsacks are given such thateach item j has a profit pjand weightwj,and each knapsack i has a capacity Ci.The goal is to find a subset of items of maximum profit such that they have a feasible packing in the knapsacks.MKP(B,S,m,n) is strongly NP- Complete and no polynomial- time approximation algorithm can have an approxima- tion ratio better than0 .5 .In the last ten years,semi- definite programming has been empolyed to solve some combinatorial problems successfully.This paper firstly presents a semi- definite re- laxation algorithm (MKPS) for MKP (B,S,m,n) .It is proved that MKPS have a approxima- tion ratio better than 0 .5 for a subclass of MKP (B,S,m,n) with n≤ 1 0 0 ,m≤ 5 and maxnj=1{ wj} minmi=1{ Ci} ≤ 2 3 .
基金Supported by the Key Discipline Construction Project of Tianshui Normal University
文摘In this paper,the generalized inverse eigenvalue problem for the(P,Q)-conjugate matrices and the associated approximation problem are discussed by using generalized singular value decomposition(GSVD).Moreover,the least residual problem of the above generalized inverse eigenvalue problem is studied by using the canonical correlation decomposition(CCD).The solutions to these problems are derived.Some numerical examples are given to illustrate the main results.
文摘In this paper,the Hermitian reflexive(Anti-Hermitian reflexive)least-squares so-lutions of matrix equations(AX = B,XC = D)are considered.With special properties of partitioned matrices and Hermitian reflexive(Anti-Hermitian reflexive)matrices,the general expression of the solution is obtained.Moreover,the related optimal approximation problem to a given matrix over the solution set is considered.
基金The work of this author was supported in part by Natural Science Foundation of Hunan Province (No. 03JJY6028).
文摘Let SE denote the least-squares symmetric solution set of the matrix equation A×B = C, where A, B and C are given matrices of suitable size. To find the optimal approximate solution in the set SE to a given matrix, we give a new feasible method based on the projection theorem, the generalized SVD and the canonical correction decomposition.