Disassembly sequence planning (DSP) plays a significant role in maintenance planning of the aircraft. It is used during the design stage for the analysis of maintainability of the aircraft. To solve product disassem...Disassembly sequence planning (DSP) plays a significant role in maintenance planning of the aircraft. It is used during the design stage for the analysis of maintainability of the aircraft. To solve product disassembly sequence planning problems efficiently, a product disassembly hybrid graph model, which describes the connection, non-connection and precedence relationships between the product parts, is established based on the characteristic of disassembly. Farther, the optimization model is provided to optimize disassembly sequence. And the solution methodology based on the genetic/simulated annealing algorithm with binaxy-tree algorithm is given. Finally, an example is analyzed in detail, and the result shows that the model is correct and efficient.展开更多
The similarity measure is crucial to the performance of spectral clustering. The Gaussian kernel function based on the Euclidean distance is usual y adopted as the similarity measure. However, the Euclidean distance m...The similarity measure is crucial to the performance of spectral clustering. The Gaussian kernel function based on the Euclidean distance is usual y adopted as the similarity measure. However, the Euclidean distance measure cannot ful y reveal the complex distribution data, and the result of spectral clustering is very sensitive to the scaling parameter. To solve these problems, a new manifold distance measure and a novel simulated anneal-ing spectral clustering (SASC) algorithm based on the manifold distance measure are proposed. The simulated annealing based on genetic algorithm (SAGA), characterized by its rapid convergence to the global optimum, is used to cluster the sample points in the spectral mapping space. The proposed algorithm can not only reflect local and global consistency better, but also reduce the sensitivity of spectral clustering to the kernel parameter, which improves the algorithm’s clustering performance. To efficiently apply the algorithm to image segmentation, the Nystrom method is used to reduce the computation complexity. Experimental results show that compared with traditional clustering algorithms and those popular spectral clustering algorithms, the proposed algorithm can achieve better clustering performances on several synthetic datasets, texture images and real images.展开更多
In rough communication, because each agent has a different language and cannot provide precise communication to each other, the concept translated among multi-agents will loss some information and this results in a le...In rough communication, because each agent has a different language and cannot provide precise communication to each other, the concept translated among multi-agents will loss some information and this results in a less or rougher concept. With different translation sequences, the problem of information loss is varied. To get the translation sequence, in which the jth agent taking part in rough communication gets maximum information, a simulated annealing algorithm is used. Analysis and simulation of this algorithm demonstrate its effectiveness.展开更多
Evolutionary computation techniques have mostly been used to solve various optimization problems, and it is well known that graph isomorphism problem (GIP) is a nondeterministic polynomial problem. A simulated annea...Evolutionary computation techniques have mostly been used to solve various optimization problems, and it is well known that graph isomorphism problem (GIP) is a nondeterministic polynomial problem. A simulated annealing (SA) algorithm for detecting graph isomorphism is proposed, and the proposed SA algorithm is well suited to deal with random graphs with large size. To verify the validity of the proposed SA algorithm, simulations are performed on three pairs of small graphs and four pairs of large random graphs with edge densities 0.5, 0.1, and 0.01, respectively. The simulation results show that the proposed SA algorithm can detect graph isomorphism with a high probability.展开更多
At present,simultaneous localization and mapping(SLAM) for an autonomous underwater vehicle(AUV)is a research hotspot.Aiming at the problem of non-linear model and non-Gaussian noise in AUV motion,an improved method o...At present,simultaneous localization and mapping(SLAM) for an autonomous underwater vehicle(AUV)is a research hotspot.Aiming at the problem of non-linear model and non-Gaussian noise in AUV motion,an improved method of variance reduction fast simultaneous localization and mapping(FastSLAM) with simulated annealing is proposed to solve the problems of particle degradation,particle depletion and particle loss in traditional FastSLAM,which lead to the reduction of AUV location estimation accuracy.The adaptive exponential fading factor is generated by the anneal function of simulated annealing algorithm to improve the effective particle number and replace resampling.By increasing the weight of small particles and decreasing the weight of large particles,the variance of particle weight can be reduced,the number of effective particles can be increased,and the accuracy of AUV location and feature location estimation can be improved to some extent by retaining more information carried by particles.The experimental results based on trial data show that the proposed simulated annealing variance reduction FastSLAM method avoids particle degradation,maintains the diversity of particles,weakened the degeneracy and improves the accuracy and stability of AUV navigation and localization system.展开更多
Taking the ratio of heat transfer area to net power and heat recovery efficiency into account, a multi-objective mathematical model was developed for organic Rankine cycle (ORC). Working fluids considered were R123,...Taking the ratio of heat transfer area to net power and heat recovery efficiency into account, a multi-objective mathematical model was developed for organic Rankine cycle (ORC). Working fluids considered were R123, R134a, R141b, R227ea and R245fa. Under the given conditions, the parameters including evaporating and condensing pressures, working fluid and cooling water velocities were optimized by simulated annealing algorithm. The results show that the optimal evaporating pressure increases with the heat source temperature increasing. Compared with other working fluids, R123 is the best choice for the temperature range of 100--180℃ and R141 b shows better performance when the temperature is higher than 180 ℃. Economic characteristic of system decreases rapidly with the decrease of heat source temperature. ORC system is uneconomical for the heat source temperature lower than 100℃.展开更多
Accurate forecasting of wind velocity can improve the economic dispatch and safe operation of the power system. Support vector machine (SVM) has been proved to be an efficient approach for forecasting. According to th...Accurate forecasting of wind velocity can improve the economic dispatch and safe operation of the power system. Support vector machine (SVM) has been proved to be an efficient approach for forecasting. According to the analysis with support vector machine method, the drawback of determining the parameters only by experts' experience should be improved. After a detailed description of the methodology of SVM and simulated annealing, an improved algorithm was proposed for the automatic optimization of parameters using SVM method. An example has proved that the proposed method can efficiently select the parameters of the SVM method. And by optimizing the parameters, the forecasting accuracy of the max wind velocity increases by 34.45%, which indicates that the new SASVM model improves the forecasting accuracy.展开更多
In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position inf...In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.展开更多
A hybrid feature selection and classification strategy was proposed based on the simulated annealing genetic algonthrn and multiple instance learning (MIL). The band selection method was proposed from subspace decom...A hybrid feature selection and classification strategy was proposed based on the simulated annealing genetic algonthrn and multiple instance learning (MIL). The band selection method was proposed from subspace decomposition, which combines the simulated annealing algorithm with the genetic algorithm in choosing different cross-over and mutation probabilities, as well as mutation individuals. Then MIL was combined with image segmentation, clustering and support vector machine algorithms to classify hyperspectral image. The experimental results show that this proposed method can get high classification accuracy of 93.13% at small training samples and the weaknesses of the conventional methods are overcome.展开更多
Under the condition of the designated collection ratio and the interfused ratio of mullock, to ensure the least energy consumption, the parameters of collecting head (the feed speed, the axes height of collecting hea...Under the condition of the designated collection ratio and the interfused ratio of mullock, to ensure the least energy consumption, the parameters of collecting head (the feed speed, the axes height of collecting head, and the rotate speed) are chosen as the optimized parameters. According to the force on the cutting pick, the collecting size of the cobalt crust and bedrock and the optimized energy consumption of the collecting head, the optimized design model of collecting head is built. Taking two hundred groups seabed microtopography for grand in the range of depth displacement from 4.5 to 5.5 era, then making use of the improved simulated annealing genetic algorithm (SAGA), the corresponding optimized result can be obtained. At the same time, in order to speed up the controlling of collecting head, the optimization results are analyzed using the regression analysis method, and the conclusion of the second parameter of the seabed microtopography is drawn.展开更多
Most conventional ceramic based aluminum metal matrix composites(MMCs) are either heavy,costly or combination of both. In order to reduce cost and weight,while at the same time maintaining quality,cow horn particles(C...Most conventional ceramic based aluminum metal matrix composites(MMCs) are either heavy,costly or combination of both. In order to reduce cost and weight,while at the same time maintaining quality,cow horn particles(CHp) was used with aluminum alloy A356 to produce MMC for brake drum application and other engineering uses. The aim of this research is to model the age hardening process of the produced composite using response surface methodology(RSM) and artificial neural network(ANN),and to use the developed ANN as fitness function for a simulated annealing optimization algorithm(SA-NN system) for optimization of age hardening process parameters. The results show that ANN modeled the age hardening data excellently and better than RSM with a correlation coefficient of experimental response with ANN predictions being 0.9921 as against 0.9583 for the RSM. The SA-NN system optimized process parameters were in very close agreement with the experimental values with the maximum relative error of 1.2%,minimum of 0.35% and average of 0.71%.展开更多
According to time-sharing valuation principle (TSVP) of power supply, the relationships of current density and current efficiency at different acidities are obtained based on the processed data of electrolytic deposit...According to time-sharing valuation principle (TSVP) of power supply, the relationships of current density and current efficiency at different acidities are obtained based on the processed data of electrolytic deposition process of zinc (EDPZ) with the least square method (LSM). Thus an optimal model of time-sharing power supply system for EDPZ is established, which has been optimized by use of an improved efficient simulated annealing algorithm (SAA). Practical results show that industrial and mining enterprises can obtain enormous economic benefits every year.展开更多
The vehicle routing problem (VRP) can be described as the problem of designing the optimal delivery or collection routes from one or several depots to a number of geographically scattered customers, subject to load co...The vehicle routing problem (VRP) can be described as the problem of designing the optimal delivery or collection routes from one or several depots to a number of geographically scattered customers, subject to load constraints. The routing decision involves determining which of the demand s will be satisfied by each vehicle and what route each vehicle will follow in s erving its assigned demand in order to minimize total delivery cost. In this pap er, a methodology for the design of VRP by integrating optimization and simulate d annealing (SA) approach is presented hierarchically. To express the problem of vehicle routing, a new mathematical formulation is first conducted. The objecti ve function involves both the delivery cost and the vehicle acquisition cost wit h load constraints. A heuristic is then proposed to solve this problem by using SA procedure in conjunction with any solution procedure of travelling salesman p roblem (TSP). The initial configuration is arranged as one vehicle route ser ving one customer. The SA searching procedure is then developed to combine custo mer to any one of the vehicle routes existed in the system if the capacity and c ost are attractive. An important concept of this proposed heuristic is that it attempts to minimize total number of vehicle required in the system on the b asis of the fixed cost and the variable cost view points. In addition, this appr oach can be easily adapted to accommodate many additional problem complexities.展开更多
A two-pass annealing/quenching internal spinning process with small-end rotations is proposed to form a curved generatrix conical thin-walled shell.That is,annealing at 360°C for 2 h followed by the 1st pass spin...A two-pass annealing/quenching internal spinning process with small-end rotations is proposed to form a curved generatrix conical thin-walled shell.That is,annealing at 360°C for 2 h followed by the 1st pass spinning,and finally quenching in ice water after holding for 1 h at 498°C followed by the 2nd pass spinning.ABAQUS finite element software is used to simulate the internal spinning process of the products formed under different forming parameters.The distribution laws of spinning force,the stress and strain under different forming processes were compared and analyzed.The mechanical properties and microstructure of the products are subsequently analyzed.The results show that the strain and the residual stress in the skin area of the formed products under two-pass spinning process more uniform,and the hardness and the mechanical performance are improved.The microstructure of the products formed with the 0.15 mm thickness reduction at the 2nd pass is excellent.And the second phase grain size distributed uniformly in the range of 36μm.Whereas,the second phase particles are broken seriously and the size distribution inhomogeneity is increased when the thickness reduction in the skin area is greater than 0.20 mm at the 2nd pass spinning process.展开更多
This paper deals with identification of subset autoregressive time series model. A simulated annealing based on approach for deter min ing the optimal subset of regression terms included in the subset autogressive m...This paper deals with identification of subset autoregressive time series model. A simulated annealing based on approach for deter min ing the optimal subset of regression terms included in the subset autogressive model is introduced. Numerical examples are given to show the effectiveness of the proposed algorithm.展开更多
基金supported by the National High Technology Research and Development Program of China(2006AA04Z427).
文摘Disassembly sequence planning (DSP) plays a significant role in maintenance planning of the aircraft. It is used during the design stage for the analysis of maintainability of the aircraft. To solve product disassembly sequence planning problems efficiently, a product disassembly hybrid graph model, which describes the connection, non-connection and precedence relationships between the product parts, is established based on the characteristic of disassembly. Farther, the optimization model is provided to optimize disassembly sequence. And the solution methodology based on the genetic/simulated annealing algorithm with binaxy-tree algorithm is given. Finally, an example is analyzed in detail, and the result shows that the model is correct and efficient.
基金supported by the National Natural Science Foundationof China(61272119)
文摘The similarity measure is crucial to the performance of spectral clustering. The Gaussian kernel function based on the Euclidean distance is usual y adopted as the similarity measure. However, the Euclidean distance measure cannot ful y reveal the complex distribution data, and the result of spectral clustering is very sensitive to the scaling parameter. To solve these problems, a new manifold distance measure and a novel simulated anneal-ing spectral clustering (SASC) algorithm based on the manifold distance measure are proposed. The simulated annealing based on genetic algorithm (SAGA), characterized by its rapid convergence to the global optimum, is used to cluster the sample points in the spectral mapping space. The proposed algorithm can not only reflect local and global consistency better, but also reduce the sensitivity of spectral clustering to the kernel parameter, which improves the algorithm’s clustering performance. To efficiently apply the algorithm to image segmentation, the Nystrom method is used to reduce the computation complexity. Experimental results show that compared with traditional clustering algorithms and those popular spectral clustering algorithms, the proposed algorithm can achieve better clustering performances on several synthetic datasets, texture images and real images.
基金the Natural Science Foundation of Shandong Province (Y2006A12)the Scientific ResearchDevelopment Project of Shandong Provincial Education Department(J06P01)the Doctoral Foundation of University of Jinan(B0633).
文摘In rough communication, because each agent has a different language and cannot provide precise communication to each other, the concept translated among multi-agents will loss some information and this results in a less or rougher concept. With different translation sequences, the problem of information loss is varied. To get the translation sequence, in which the jth agent taking part in rough communication gets maximum information, a simulated annealing algorithm is used. Analysis and simulation of this algorithm demonstrate its effectiveness.
基金the National Natural Science Foundation of China (60373089, 60674106, and 60533010)the National High Technology Research and Development "863" Program (2006AA01Z104)
文摘Evolutionary computation techniques have mostly been used to solve various optimization problems, and it is well known that graph isomorphism problem (GIP) is a nondeterministic polynomial problem. A simulated annealing (SA) algorithm for detecting graph isomorphism is proposed, and the proposed SA algorithm is well suited to deal with random graphs with large size. To verify the validity of the proposed SA algorithm, simulations are performed on three pairs of small graphs and four pairs of large random graphs with edge densities 0.5, 0.1, and 0.01, respectively. The simulation results show that the proposed SA algorithm can detect graph isomorphism with a high probability.
基金supported by the National Science Fund of China under Grants 61603034China Postdoctoral Science Foundation under Grant 2019M653870XB+1 种基金Beijing Municipal Natural Science Foundation (3182027)Fundamental Research Funds for the Central Universities,China,FRF-GF-17-B44,and XJS191315
文摘At present,simultaneous localization and mapping(SLAM) for an autonomous underwater vehicle(AUV)is a research hotspot.Aiming at the problem of non-linear model and non-Gaussian noise in AUV motion,an improved method of variance reduction fast simultaneous localization and mapping(FastSLAM) with simulated annealing is proposed to solve the problems of particle degradation,particle depletion and particle loss in traditional FastSLAM,which lead to the reduction of AUV location estimation accuracy.The adaptive exponential fading factor is generated by the anneal function of simulated annealing algorithm to improve the effective particle number and replace resampling.By increasing the weight of small particles and decreasing the weight of large particles,the variance of particle weight can be reduced,the number of effective particles can be increased,and the accuracy of AUV location and feature location estimation can be improved to some extent by retaining more information carried by particles.The experimental results based on trial data show that the proposed simulated annealing variance reduction FastSLAM method avoids particle degradation,maintains the diversity of particles,weakened the degeneracy and improves the accuracy and stability of AUV navigation and localization system.
基金Project(2009GK2009) supported by Science and Technology Department Funds of Hunan Province,ChinaProject(08C26224302178) supported by Innovation Fund for Technology Based Firms of China
文摘Taking the ratio of heat transfer area to net power and heat recovery efficiency into account, a multi-objective mathematical model was developed for organic Rankine cycle (ORC). Working fluids considered were R123, R134a, R141b, R227ea and R245fa. Under the given conditions, the parameters including evaporating and condensing pressures, working fluid and cooling water velocities were optimized by simulated annealing algorithm. The results show that the optimal evaporating pressure increases with the heat source temperature increasing. Compared with other working fluids, R123 is the best choice for the temperature range of 100--180℃ and R141 b shows better performance when the temperature is higher than 180 ℃. Economic characteristic of system decreases rapidly with the decrease of heat source temperature. ORC system is uneconomical for the heat source temperature lower than 100℃.
基金Project(71071052) supported by the National Natural Science Foundation of ChinaProject(JB2011097) supported by the Fundamental Research Funds for the Central Universities of China
文摘Accurate forecasting of wind velocity can improve the economic dispatch and safe operation of the power system. Support vector machine (SVM) has been proved to be an efficient approach for forecasting. According to the analysis with support vector machine method, the drawback of determining the parameters only by experts' experience should be improved. After a detailed description of the methodology of SVM and simulated annealing, an improved algorithm was proposed for the automatic optimization of parameters using SVM method. An example has proved that the proposed method can efficiently select the parameters of the SVM method. And by optimizing the parameters, the forecasting accuracy of the max wind velocity increases by 34.45%, which indicates that the new SASVM model improves the forecasting accuracy.
基金Projects(40804027,41074085) supported by the National Natural Science Foundation of ChinaProject(09JJ3048) supported by the Natural Science Foundation of Hunan Province,ChinaProject(200805331082) supported by the Research Fund for the Doctoral Program of Higher Education,China
文摘In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.
文摘A hybrid feature selection and classification strategy was proposed based on the simulated annealing genetic algonthrn and multiple instance learning (MIL). The band selection method was proposed from subspace decomposition, which combines the simulated annealing algorithm with the genetic algorithm in choosing different cross-over and mutation probabilities, as well as mutation individuals. Then MIL was combined with image segmentation, clustering and support vector machine algorithms to classify hyperspectral image. The experimental results show that this proposed method can get high classification accuracy of 93.13% at small training samples and the weaknesses of the conventional methods are overcome.
基金Project(50875265) supported by the National Natural Science Foundation of ChinaProject(20080440992) supported by the Postdoctoral Science Foundation of ChinaProject(2009SK3159) supported by the Technology Support Plan of Hunan Province,China
文摘Under the condition of the designated collection ratio and the interfused ratio of mullock, to ensure the least energy consumption, the parameters of collecting head (the feed speed, the axes height of collecting head, and the rotate speed) are chosen as the optimized parameters. According to the force on the cutting pick, the collecting size of the cobalt crust and bedrock and the optimized energy consumption of the collecting head, the optimized design model of collecting head is built. Taking two hundred groups seabed microtopography for grand in the range of depth displacement from 4.5 to 5.5 era, then making use of the improved simulated annealing genetic algorithm (SAGA), the corresponding optimized result can be obtained. At the same time, in order to speed up the controlling of collecting head, the optimization results are analyzed using the regression analysis method, and the conclusion of the second parameter of the seabed microtopography is drawn.
文摘Most conventional ceramic based aluminum metal matrix composites(MMCs) are either heavy,costly or combination of both. In order to reduce cost and weight,while at the same time maintaining quality,cow horn particles(CHp) was used with aluminum alloy A356 to produce MMC for brake drum application and other engineering uses. The aim of this research is to model the age hardening process of the produced composite using response surface methodology(RSM) and artificial neural network(ANN),and to use the developed ANN as fitness function for a simulated annealing optimization algorithm(SA-NN system) for optimization of age hardening process parameters. The results show that ANN modeled the age hardening data excellently and better than RSM with a correlation coefficient of experimental response with ANN predictions being 0.9921 as against 0.9583 for the RSM. The SA-NN system optimized process parameters were in very close agreement with the experimental values with the maximum relative error of 1.2%,minimum of 0.35% and average of 0.71%.
文摘According to time-sharing valuation principle (TSVP) of power supply, the relationships of current density and current efficiency at different acidities are obtained based on the processed data of electrolytic deposition process of zinc (EDPZ) with the least square method (LSM). Thus an optimal model of time-sharing power supply system for EDPZ is established, which has been optimized by use of an improved efficient simulated annealing algorithm (SAA). Practical results show that industrial and mining enterprises can obtain enormous economic benefits every year.
文摘The vehicle routing problem (VRP) can be described as the problem of designing the optimal delivery or collection routes from one or several depots to a number of geographically scattered customers, subject to load constraints. The routing decision involves determining which of the demand s will be satisfied by each vehicle and what route each vehicle will follow in s erving its assigned demand in order to minimize total delivery cost. In this pap er, a methodology for the design of VRP by integrating optimization and simulate d annealing (SA) approach is presented hierarchically. To express the problem of vehicle routing, a new mathematical formulation is first conducted. The objecti ve function involves both the delivery cost and the vehicle acquisition cost wit h load constraints. A heuristic is then proposed to solve this problem by using SA procedure in conjunction with any solution procedure of travelling salesman p roblem (TSP). The initial configuration is arranged as one vehicle route ser ving one customer. The SA searching procedure is then developed to combine custo mer to any one of the vehicle routes existed in the system if the capacity and c ost are attractive. An important concept of this proposed heuristic is that it attempts to minimize total number of vehicle required in the system on the b asis of the fixed cost and the variable cost view points. In addition, this appr oach can be easily adapted to accommodate many additional problem complexities.
基金Project(51775479)supported by the National Natural Science Foundation of ChinaProject(E2017203046)supported by the Natural Science Foundation of Hebei Province,China
文摘A two-pass annealing/quenching internal spinning process with small-end rotations is proposed to form a curved generatrix conical thin-walled shell.That is,annealing at 360°C for 2 h followed by the 1st pass spinning,and finally quenching in ice water after holding for 1 h at 498°C followed by the 2nd pass spinning.ABAQUS finite element software is used to simulate the internal spinning process of the products formed under different forming parameters.The distribution laws of spinning force,the stress and strain under different forming processes were compared and analyzed.The mechanical properties and microstructure of the products are subsequently analyzed.The results show that the strain and the residual stress in the skin area of the formed products under two-pass spinning process more uniform,and the hardness and the mechanical performance are improved.The microstructure of the products formed with the 0.15 mm thickness reduction at the 2nd pass is excellent.And the second phase grain size distributed uniformly in the range of 36μm.Whereas,the second phase particles are broken seriously and the size distribution inhomogeneity is increased when the thickness reduction in the skin area is greater than 0.20 mm at the 2nd pass spinning process.
文摘This paper deals with identification of subset autoregressive time series model. A simulated annealing based on approach for deter min ing the optimal subset of regression terms included in the subset autogressive model is introduced. Numerical examples are given to show the effectiveness of the proposed algorithm.