The emergency communication system based on rail is an unconventional emergency communication mode,it is a complement equipment for that conventional communication system can’t work while tunnel mine accident occurs....The emergency communication system based on rail is an unconventional emergency communication mode,it is a complement equipment for that conventional communication system can’t work while tunnel mine accident occurs.Medium of transmission channel is the widely existing rail in the tunnel.In this paper we analyzed the characteristics of the rail transmission channel,verified the feasibility that information is transmitted by vibration signal in rail,we proposed the realization plan of the system.Communication protocol and processing mechanism suitable for rail transmission are designed according to the characteristics of channel bandwidth and low data transmission.Information communication with low bit rate and low bit error is realized in the communication simulation model.In the simplified model,we realized to transmit recognition speech information,and the error rate of the key text information is low to accept.The most concerned problem of personnel location in the mine disaster rescue is proposed,the composite algorithm is based on the model of signal amplitude attenuation,key node information and data frame transmission delay.Location information of hitting point can be achieved within the simplified model of the experiment.Furthermore,we discuss the characteristics of vibration signals passing through different channels.展开更多
This paper presents a parallel composite local search algorithm based on multiple search neighborhoods to solve a special kind of timetable problem. The new algorithm can also effectively solve those problems that can...This paper presents a parallel composite local search algorithm based on multiple search neighborhoods to solve a special kind of timetable problem. The new algorithm can also effectively solve those problems that can be solved by general local search algorithms. Experimental results show that the new algorithm can generate better solutions than general local search algorithms.展开更多
A method to improve the low-velocity impact performance of composite laminate is proposed, and a multi-island genetic algorithm is used for the optimization of composite laminate stacking sequence under low-velocity i...A method to improve the low-velocity impact performance of composite laminate is proposed, and a multi-island genetic algorithm is used for the optimization of composite laminate stacking sequence under low-velocity impact loads based on a 2D dynamic impact finite element analysis. Low-velocity impact tests and compression-after impact(CAI) tests have been conducted to verify the effectiveness of optimization method. Experimental results show that the impact damage areas of the optimized laminate have been reduced by 42.1% compared to the baseline specimen, and the residual compression strength has been increased by 10.79%, from baseline specimen 156.97 MPa to optimized 173.91 MPa. The tests result shows that optimization method can effectively enhance the impact performances of the laminate.展开更多
Algorithm of fuzzy reasoning has been successful applied in fuzzy control,but its theoretical foundation of algorithms has not been thoroughly investigated. In this paper,structure of basic algorithms of fuzzy reasoni...Algorithm of fuzzy reasoning has been successful applied in fuzzy control,but its theoretical foundation of algorithms has not been thoroughly investigated. In this paper,structure of basic algorithms of fuzzy reasoning was studied, its rationality was discussed from the viewpoint of logic and mathematics, and three theorems were proved. These theorems shows that there always exists a mathe-~matical relation (that is, a bounded real function) between the premises and the conclusion for fuzzy reasoning, and in fact various algorithms of fuzzy reasoning are specific forms of this function. Thus these results show that algorithms of fuzzy reasoning are theoretically reliable.展开更多
A new optimization method for the optimization of stacking of composite glass fiber laminates is developed. The fiber orientation and angle of the layers of the cylindrical shells are sought considering the buckling l...A new optimization method for the optimization of stacking of composite glass fiber laminates is developed. The fiber orientation and angle of the layers of the cylindrical shells are sought considering the buckling load. The proposed optimization algorithm applies both finite element analysis and the mode-pursuing sampling (MPS)method. The algorithms suggest the optimal stacking sequence for achieving the maximal buckling load. The procedure is implemented by integrating ANSYS and MATLAB. The stacking sequence designing for the symmetric angle-ply three-layered and five-layered composite cylinder shells is presented to illustrate the optimization process, respectively. Compared with the genetic algorithms, the proposed optimization method is much faster and efficient for composite staking sequence plan.展开更多
With the continuous development of deep learning and artificial neural networks(ANNs), algorithmic composition has gradually become a hot research field. In order to solve the music-style problem in generating chord m...With the continuous development of deep learning and artificial neural networks(ANNs), algorithmic composition has gradually become a hot research field. In order to solve the music-style problem in generating chord music, a multi-style chord music generation(MSCMG) network is proposed based on the previous ANN for creation. A music-style extraction module and a style extractor are added by the network on the original basis;the music-style extraction module divides the entire music content into two parts, namely the music-style information Mstyleand the music content information Mcontent. The style extractor removes the music-style information entangled in the music content information. The similarity of music generated by different models is compared in this paper. It is also evaluated whether the model can learn music composition rules from the database. Through experiments, it is found that the model proposed in this paper can generate music works in the expected style. Compared with the long short term memory(LSTM) network, the MSCMG network has a certain improvement in the performance of music styles.展开更多
Global look-up table strategy proposed recently has been proven to be an efficient method to accelerate the interpolation, which is the most time-consuming part in the iterative sub-pixel digital image correlation (...Global look-up table strategy proposed recently has been proven to be an efficient method to accelerate the interpolation, which is the most time-consuming part in the iterative sub-pixel digital image correlation (DIC) algorithms. In this paper, a global look-up table strategy with cubic B-spline interpolation is developed for the DIC method based on the inverse compositional Gauss-Newton (IC-GN) algorithm. The performance of this strategy, including accuracy, precision, and computation efficiency, is evaluated through a theoretical and experimental study, using the one with widely employed bicubic interpolation as a benchmark. The global look-up table strategy with cubic B-spline interpolation improves significantly the accuracy of the IC-GN algorithm-based DIC method compared with the one using the bicubic interpolation, at a trivial price of computation efficiency.展开更多
Algorithmic composition is a very popular research field today. Bach's "two voice part invention" is the research object in this paper. The grammar and compositional rules of "invention" are in...Algorithmic composition is a very popular research field today. Bach's "two voice part invention" is the research object in this paper. The grammar and compositional rules of "invention" are introduced first. Then two soft computational methods,genetic algorithms and back propagation (BP) neural network technology,are combined to the experiment on assisting in composing "two voice part inventions". The system presented in this paper is quite effective and satisfactory.展开更多
This paper proposes the efficient model building in active appearance model(AAM) for the rotated face.Finding an exact region of the face is generally difficult due to different shapes and viewpoints.Unlike many paper...This paper proposes the efficient model building in active appearance model(AAM) for the rotated face.Finding an exact region of the face is generally difficult due to different shapes and viewpoints.Unlike many papers about the fitting method of AAM,this paper treats how images are chosen for fitting of the rotated face in modelling process.To solve this problem,databases of facial rotation and expression are selected and models are built using Procrustes method and principal component analysis(PCA).These models are applied in fitting methods like basic AAM fitting,inverse compositional alignment(ICA),project-out ICA,normalization ICA,robust normalization inverse compositional algorithm(RNIC)and efficient robust normalization algorithm(ERN).RNIC and ERN can fit the rotated face in images efficiently.The efficiency of model building is checked using sequence images made by ourselves.展开更多
This paper approaches melody harmonization with piano accompaniment as a machine learning task, in a probabilistic framework. An existing Sample score set of Chinese folk song with piano accompaniment is used to build...This paper approaches melody harmonization with piano accompaniment as a machine learning task, in a probabilistic framework. An existing Sample score set of Chinese folk song with piano accompaniment is used to build a model of piano accomimniment process. The mcdel can then be used to harmonize new mclody with piano accompaniment.展开更多
One important application of independent component analysis (ICA) is in image processing. A two dimensional (2-D) composite ICA algorithm framework for 2-D image independent component analysis (2-D ICA) is propo...One important application of independent component analysis (ICA) is in image processing. A two dimensional (2-D) composite ICA algorithm framework for 2-D image independent component analysis (2-D ICA) is proposed. The 2-D nature of the algorithm provides it an advantage of circumventing the roundabout transforming procedures between two dimensional (2-D) image deta and one-dimensional (l-D) signal. Moreover the combination of the Newton (fixed-point algorithm) and natural gradient algorithms in this composite algorithm increases its efficiency and robustness. The convincing results of a successful example in functional magnetic resonance imaging (fMRI) show the potential application of composite 2-D ICA in the brain activity detection.展开更多
In order to decrease the number of design variables and improve the efficiency of com- posite structure optimal design, a single-level composite structure optimization method based on a tapered model is presented. Com...In order to decrease the number of design variables and improve the efficiency of com- posite structure optimal design, a single-level composite structure optimization method based on a tapered model is presented. Compared with the conventional multi-level composite structure opti- mization method, this single-level method has many advantages. First, by using a distance variable and a ply group variable, the number of design variables is decreased evidently and independent with the density of sub-regions, which makes the single-level method very suitable for large-scale composite structures. Second, it is very convenient to optimize laminate thickness and stacking sequence in the same level, which probably improves the quality of optimal result. Third, ply con-tinuity can be guaranteed between sub-regions in the single-level method, which could reduce stress concentration and manufacturing difficulty. An example of a composite wing is used to demonstrate the advantages and competence of the single-level method proposed.展开更多
Large-scale service composition has become an important research topic in Service-Oriented Computing(SOC). Quality of Service(Qo S) has been mostly applied to represent nonfunctional properties of web services and...Large-scale service composition has become an important research topic in Service-Oriented Computing(SOC). Quality of Service(Qo S) has been mostly applied to represent nonfunctional properties of web services and to differentiate those with the same functionality. Many studies for measuring service composition in terms of Qo S have been completed. Among current popular optimization methods for service composition, the exhaustion method has some disadvantages such as requiring a large number of calculations and poor scalability. Similarly,the traditional evolutionary computation method has defects such as exhibiting slow convergence speed and falling easily into the local optimum. In order to solve these problems, an improved optimization algorithm, WS FOA(Web Service composition based on Fruit Fly Optimization Algorithm) for service composition, was proposed, on the basis of the modeling of service composition and the FOA. Simulated experiments demonstrated that the algorithm is effective, feasible, stable, and possesses good global searching ability.展开更多
Wireless Sensor Networks(WSNs) have many applications, such as climate monitoring systems, fire detection, smart homes, and smart cities. It is expected that WSNs will be integrated into the Internet of Things(IoT...Wireless Sensor Networks(WSNs) have many applications, such as climate monitoring systems, fire detection, smart homes, and smart cities. It is expected that WSNs will be integrated into the Internet of Things(IoT)and participate in various tasks. WSNs play an important role monitoring and reporting environment information and collecting surrounding context. In this paper we consider a WSN deployed for an application such as environment monitoring, and a mobile sink which acts as the gateway between the Internet and the WSN. Data gathering is a challenging problem in WSNs and in the IoT because the information has to be available quickly and effectively without delays and redundancies. In this paper we propose several distributed algorithms for composite event detection and reporting to a mobile sink. Once data is collected by the sink, it can be shared using the IoT infrastructure. We analyze the performance of our algorithms using WSNet simulator, which is specially designed for event-based WSNs. We measure various metrics such as average residual energy, percentage of composite events processed successfully at the sink, and the average number of hops to reach the sink.展开更多
Parallel computing techniques have been introduced into digital image correlation(DIC) in recent years and leads to a surge in computation speed. The graphics processing unit(GPU)-based parallel computing demonstrated...Parallel computing techniques have been introduced into digital image correlation(DIC) in recent years and leads to a surge in computation speed. The graphics processing unit(GPU)-based parallel computing demonstrated a surprising effect on accelerating the iterative subpixel DIC, compared with CPU-based parallel computing. In this paper, the performances of the two kinds of parallel computing techniques are compared for the previously proposed path-independent DIC method, in which the initial guess for the inverse compositional Gauss-Newton(IC-GN) algorithm at each point of interest(POI) is estimated through the fast Fourier transform-based cross-correlation(FFT-CC) algorithm. Based on the performance evaluation, a heterogeneous parallel computing(HPC) model is proposed with hybrid mode of parallelisms in order to combine the computing power of GPU and multicore CPU. A scheme of trial computation test is developed to optimize the configuration of the HPC model on a specific computer. The proposed HPC model shows excellent performance on a middle-end desktop computer for real-time subpixel DIC with high resolution of more than 10000 POIs per frame.展开更多
In this paper, a meta-structure of piano accompaniment figure (meta-structure for short) is proposed to harmonize a melodic piece of music so as to construct a multi-voice music. Here we approach melody harmonizatio...In this paper, a meta-structure of piano accompaniment figure (meta-structure for short) is proposed to harmonize a melodic piece of music so as to construct a multi-voice music. Here we approach melody harmonization with piano accompaniment as a machine learning task in a probabilistic framework. A series of piano accompaniment figures are collected from the massive existing sample scores and converted into a set of meta-structure. After the procedure of samples training, a model is formulated to generate a proper piano accompaniment figure for a harmonizing unit in the context. This model is flexible in harmonizing a melody with piano accompaniment. The experimental results are evaluated and discussed.展开更多
基金The authors would like to thank National Natural Science Foundation of China for the grant of the project(41574137)Furthermore,they would like to specially thank Prof.Guo Yong for his contributions and his support in this paper.
文摘The emergency communication system based on rail is an unconventional emergency communication mode,it is a complement equipment for that conventional communication system can’t work while tunnel mine accident occurs.Medium of transmission channel is the widely existing rail in the tunnel.In this paper we analyzed the characteristics of the rail transmission channel,verified the feasibility that information is transmitted by vibration signal in rail,we proposed the realization plan of the system.Communication protocol and processing mechanism suitable for rail transmission are designed according to the characteristics of channel bandwidth and low data transmission.Information communication with low bit rate and low bit error is realized in the communication simulation model.In the simplified model,we realized to transmit recognition speech information,and the error rate of the key text information is low to accept.The most concerned problem of personnel location in the mine disaster rescue is proposed,the composite algorithm is based on the model of signal amplitude attenuation,key node information and data frame transmission delay.Location information of hitting point can be achieved within the simplified model of the experiment.Furthermore,we discuss the characteristics of vibration signals passing through different channels.
文摘This paper presents a parallel composite local search algorithm based on multiple search neighborhoods to solve a special kind of timetable problem. The new algorithm can also effectively solve those problems that can be solved by general local search algorithms. Experimental results show that the new algorithm can generate better solutions than general local search algorithms.
基金Funded by the National Natural Science Foundation of China(No.51275393)the Fundamental Research Funds for the Central Universities(No.xjj2017160)
文摘A method to improve the low-velocity impact performance of composite laminate is proposed, and a multi-island genetic algorithm is used for the optimization of composite laminate stacking sequence under low-velocity impact loads based on a 2D dynamic impact finite element analysis. Low-velocity impact tests and compression-after impact(CAI) tests have been conducted to verify the effectiveness of optimization method. Experimental results show that the impact damage areas of the optimized laminate have been reduced by 42.1% compared to the baseline specimen, and the residual compression strength has been increased by 10.79%, from baseline specimen 156.97 MPa to optimized 173.91 MPa. The tests result shows that optimization method can effectively enhance the impact performances of the laminate.
文摘Algorithm of fuzzy reasoning has been successful applied in fuzzy control,but its theoretical foundation of algorithms has not been thoroughly investigated. In this paper,structure of basic algorithms of fuzzy reasoning was studied, its rationality was discussed from the viewpoint of logic and mathematics, and three theorems were proved. These theorems shows that there always exists a mathe-~matical relation (that is, a bounded real function) between the premises and the conclusion for fuzzy reasoning, and in fact various algorithms of fuzzy reasoning are specific forms of this function. Thus these results show that algorithms of fuzzy reasoning are theoretically reliable.
基金Innovation Team Development Program of Ministry of Education of China (No. IRT0763)National Natural Science Foundation of China (No. 50205028).
文摘A new optimization method for the optimization of stacking of composite glass fiber laminates is developed. The fiber orientation and angle of the layers of the cylindrical shells are sought considering the buckling load. The proposed optimization algorithm applies both finite element analysis and the mode-pursuing sampling (MPS)method. The algorithms suggest the optimal stacking sequence for achieving the maximal buckling load. The procedure is implemented by integrating ANSYS and MATLAB. The stacking sequence designing for the symmetric angle-ply three-layered and five-layered composite cylinder shells is presented to illustrate the optimization process, respectively. Compared with the genetic algorithms, the proposed optimization method is much faster and efficient for composite staking sequence plan.
基金National Natural Science Foundation of China (No.61801106)。
文摘With the continuous development of deep learning and artificial neural networks(ANNs), algorithmic composition has gradually become a hot research field. In order to solve the music-style problem in generating chord music, a multi-style chord music generation(MSCMG) network is proposed based on the previous ANN for creation. A music-style extraction module and a style extractor are added by the network on the original basis;the music-style extraction module divides the entire music content into two parts, namely the music-style information Mstyleand the music content information Mcontent. The style extractor removes the music-style information entangled in the music content information. The similarity of music generated by different models is compared in this paper. It is also evaluated whether the model can learn music composition rules from the database. Through experiments, it is found that the model proposed in this paper can generate music works in the expected style. Compared with the long short term memory(LSTM) network, the MSCMG network has a certain improvement in the performance of music styles.
基金financially supported by the National Natural Science Foundation of China(11202081,11272124,and 11472109)the State Key Lab of Subtropical Building Science,South China University of Technology(2014ZC17)
文摘Global look-up table strategy proposed recently has been proven to be an efficient method to accelerate the interpolation, which is the most time-consuming part in the iterative sub-pixel digital image correlation (DIC) algorithms. In this paper, a global look-up table strategy with cubic B-spline interpolation is developed for the DIC method based on the inverse compositional Gauss-Newton (IC-GN) algorithm. The performance of this strategy, including accuracy, precision, and computation efficiency, is evaluated through a theoretical and experimental study, using the one with widely employed bicubic interpolation as a benchmark. The global look-up table strategy with cubic B-spline interpolation improves significantly the accuracy of the IC-GN algorithm-based DIC method compared with the one using the bicubic interpolation, at a trivial price of computation efficiency.
基金National Natural Science Foundation of China (No.60975076)
文摘Algorithmic composition is a very popular research field today. Bach's "two voice part invention" is the research object in this paper. The grammar and compositional rules of "invention" are introduced first. Then two soft computational methods,genetic algorithms and back propagation (BP) neural network technology,are combined to the experiment on assisting in composing "two voice part inventions". The system presented in this paper is quite effective and satisfactory.
基金Next-Generation Information Computing Development Program through the National Research Foundation of Korea(NRF)funded by the Ministry of Education,Science and Technology(No.2012M3C4A7032182)The MSIP(Ministry of Science,ICT&Future Planning),Korea,under the ITRC(Information Technology Research Center)support program(NIPA-2013-H0301-13-2006)supervised by the NIPA(National IT Industry Promotion Agency)
文摘This paper proposes the efficient model building in active appearance model(AAM) for the rotated face.Finding an exact region of the face is generally difficult due to different shapes and viewpoints.Unlike many papers about the fitting method of AAM,this paper treats how images are chosen for fitting of the rotated face in modelling process.To solve this problem,databases of facial rotation and expression are selected and models are built using Procrustes method and principal component analysis(PCA).These models are applied in fitting methods like basic AAM fitting,inverse compositional alignment(ICA),project-out ICA,normalization ICA,robust normalization inverse compositional algorithm(RNIC)and efficient robust normalization algorithm(ERN).RNIC and ERN can fit the rotated face in images efficiently.The efficiency of model building is checked using sequence images made by ourselves.
基金National Natural Science Foundation of China (No.60975076)
文摘This paper approaches melody harmonization with piano accompaniment as a machine learning task, in a probabilistic framework. An existing Sample score set of Chinese folk song with piano accompaniment is used to build a model of piano accomimniment process. The mcdel can then be used to harmonize new mclody with piano accompaniment.
基金Supported by the 973 Project (No.2003CB716106), NSFC (No.90208003, 30200059), TRAPOYT, Doctor Training Fund of MOE, PRC, Key Research Project of Science and Technology of MOE, Fok Ying Tong Education Foundation (No.91041)
文摘One important application of independent component analysis (ICA) is in image processing. A two dimensional (2-D) composite ICA algorithm framework for 2-D image independent component analysis (2-D ICA) is proposed. The 2-D nature of the algorithm provides it an advantage of circumventing the roundabout transforming procedures between two dimensional (2-D) image deta and one-dimensional (l-D) signal. Moreover the combination of the Newton (fixed-point algorithm) and natural gradient algorithms in this composite algorithm increases its efficiency and robustness. The convincing results of a successful example in functional magnetic resonance imaging (fMRI) show the potential application of composite 2-D ICA in the brain activity detection.
基金supported by National Natural Science Foundation of China(No.1110216/A020312)Foundation Sciences of Northwestern Polytechnical University(No.JC20120210)
文摘In order to decrease the number of design variables and improve the efficiency of com- posite structure optimal design, a single-level composite structure optimization method based on a tapered model is presented. Compared with the conventional multi-level composite structure opti- mization method, this single-level method has many advantages. First, by using a distance variable and a ply group variable, the number of design variables is decreased evidently and independent with the density of sub-regions, which makes the single-level method very suitable for large-scale composite structures. Second, it is very convenient to optimize laminate thickness and stacking sequence in the same level, which probably improves the quality of optimal result. Third, ply con-tinuity can be guaranteed between sub-regions in the single-level method, which could reduce stress concentration and manufacturing difficulty. An example of a composite wing is used to demonstrate the advantages and competence of the single-level method proposed.
基金supported by the National Natural Science Foundation of China (Nos. 61402006 and 61202227)the Natural Science Foundation of Anhui Province of China (No. 1408085MF132)+2 种基金the Science and Technology Planning Project of Anhui Province of China (No. 1301032162)the College Students Scientific Research Training Program (No. KYXL2014060)the 211 Project of Anhui University (No. 02303301)
文摘Large-scale service composition has become an important research topic in Service-Oriented Computing(SOC). Quality of Service(Qo S) has been mostly applied to represent nonfunctional properties of web services and to differentiate those with the same functionality. Many studies for measuring service composition in terms of Qo S have been completed. Among current popular optimization methods for service composition, the exhaustion method has some disadvantages such as requiring a large number of calculations and poor scalability. Similarly,the traditional evolutionary computation method has defects such as exhibiting slow convergence speed and falling easily into the local optimum. In order to solve these problems, an improved optimization algorithm, WS FOA(Web Service composition based on Fruit Fly Optimization Algorithm) for service composition, was proposed, on the basis of the modeling of service composition and the FOA. Simulated experiments demonstrated that the algorithm is effective, feasible, stable, and possesses good global searching ability.
文摘Wireless Sensor Networks(WSNs) have many applications, such as climate monitoring systems, fire detection, smart homes, and smart cities. It is expected that WSNs will be integrated into the Internet of Things(IoT)and participate in various tasks. WSNs play an important role monitoring and reporting environment information and collecting surrounding context. In this paper we consider a WSN deployed for an application such as environment monitoring, and a mobile sink which acts as the gateway between the Internet and the WSN. Data gathering is a challenging problem in WSNs and in the IoT because the information has to be available quickly and effectively without delays and redundancies. In this paper we propose several distributed algorithms for composite event detection and reporting to a mobile sink. Once data is collected by the sink, it can be shared using the IoT infrastructure. We analyze the performance of our algorithms using WSNet simulator, which is specially designed for event-based WSNs. We measure various metrics such as average residual energy, percentage of composite events processed successfully at the sink, and the average number of hops to reach the sink.
基金supported by the National Natural Science Foundation of China(Grant Nos.11972013 and 12272145)the Ministry of Science and Technology of China(Grant No.2018YFF01014200).
基金supported by the National Natural Science Foundation of China(Grant Nos.11772131,11772132,11772134&11472109)the Natural Science Foundation of Guangdong Province,China(Grant Nos.2015A030308017,2015A030311046&2015B010131009)+2 种基金the Opening fund of State Key Laboratory of Nonlinear Mechanics(LNM)CASthe State Key Lab of Subtropical Building Science,South China University of Technology(Grant Nos.2014ZC17&2017ZD096)
文摘Parallel computing techniques have been introduced into digital image correlation(DIC) in recent years and leads to a surge in computation speed. The graphics processing unit(GPU)-based parallel computing demonstrated a surprising effect on accelerating the iterative subpixel DIC, compared with CPU-based parallel computing. In this paper, the performances of the two kinds of parallel computing techniques are compared for the previously proposed path-independent DIC method, in which the initial guess for the inverse compositional Gauss-Newton(IC-GN) algorithm at each point of interest(POI) is estimated through the fast Fourier transform-based cross-correlation(FFT-CC) algorithm. Based on the performance evaluation, a heterogeneous parallel computing(HPC) model is proposed with hybrid mode of parallelisms in order to combine the computing power of GPU and multicore CPU. A scheme of trial computation test is developed to optimize the configuration of the HPC model on a specific computer. The proposed HPC model shows excellent performance on a middle-end desktop computer for real-time subpixel DIC with high resolution of more than 10000 POIs per frame.
基金Supported by the National Natural Science Foundation of China under Grant No.60975076
文摘In this paper, a meta-structure of piano accompaniment figure (meta-structure for short) is proposed to harmonize a melodic piece of music so as to construct a multi-voice music. Here we approach melody harmonization with piano accompaniment as a machine learning task in a probabilistic framework. A series of piano accompaniment figures are collected from the massive existing sample scores and converted into a set of meta-structure. After the procedure of samples training, a model is formulated to generate a proper piano accompaniment figure for a harmonizing unit in the context. This model is flexible in harmonizing a melody with piano accompaniment. The experimental results are evaluated and discussed.