This paper investigates the anomaly-resistant decentralized state estimation(SE) problem for a class of wide-area power systems which are divided into several non-overlapping areas connected through transmission lines...This paper investigates the anomaly-resistant decentralized state estimation(SE) problem for a class of wide-area power systems which are divided into several non-overlapping areas connected through transmission lines. Two classes of measurements(i.e., local measurements and edge measurements) are obtained, respectively, from the individual area and the transmission lines. A decentralized state estimator, whose performance is resistant against measurement with anomalies, is designed based on the minimum error entropy with fiducial points(MEEF) criterion. Specifically, 1) An augmented model, which incorporates the local prediction and local measurement, is developed by resorting to the unscented transformation approach and the statistical linearization approach;2) Using the augmented model, an MEEF-based cost function is designed that reflects the local prediction errors of the state and the measurement;and 3) The local estimate is first obtained by minimizing the MEEF-based cost function through a fixed-point iteration and then updated by using the edge measuring information. Finally, simulation experiments with three scenarios are carried out on the IEEE 14-bus system to illustrate the validity of the proposed anomaly-resistant decentralized SE scheme.展开更多
Kalman filter is commonly used in data filtering and parameters estimation of nonlinear system,such as projectile's trajectory estimation and control.While there is a drawback that the prior error covariance matri...Kalman filter is commonly used in data filtering and parameters estimation of nonlinear system,such as projectile's trajectory estimation and control.While there is a drawback that the prior error covariance matrix and filter parameters are difficult to be determined,which may result in filtering divergence.As to the problem that the accuracy of state estimation for nonlinear ballistic model strongly depends on its mathematical model,we improve the weighted least squares method(WLSM)with minimum model error principle.Invariant embedding method is adopted to solve the cost function including the model error.With the knowledge of measurement data and measurement error covariance matrix,we use gradient descent algorithm to determine the weighting matrix of model error.The uncertainty and linearization error of model are recursively estimated by the proposed method,thus achieving an online filtering estimation of the observations.Simulation results indicate that the proposed recursive estimation algorithm is insensitive to initial conditions and of good robustness.展开更多
Performance of the Adaptive Coding and Modulation(ACM) strongly depends on the retrieved Channel State Information(CSI),which can be obtained using the channel estimation techniques relying on pilot symbol transmissio...Performance of the Adaptive Coding and Modulation(ACM) strongly depends on the retrieved Channel State Information(CSI),which can be obtained using the channel estimation techniques relying on pilot symbol transmission.Earlier analysis of methods of pilot-aided channel estimation for ACM systems were relatively little.In this paper,we investigate the performance of CSI prediction using the Minimum Mean Square Error(MMSE)channel estimator for an ACM system.To solve the two problems of MMSE:high computational operations and oversimplified assumption,we then propose the Low-Complexity schemes(LC-MMSE and Recursion LC-MMSE(R-LC-MMSE)).Computational complexity and Mean Square Error(MSE) are presented to evaluate the efficiency of the proposed algorithm.Both analysis and numerical results show that LC-MMSE performs close to the wellknown MMSE estimator with much lower complexity and R-LC-MMSE improves the application of MMSE estimation to specific circumstances.展开更多
Gauss-Markov model is frequently used in data analysis; the analysis and estimation of its parameters is always a hot issue. Based on the information theory and from the viewpoint of optimal information on description...Gauss-Markov model is frequently used in data analysis; the analysis and estimation of its parameters is always a hot issue. Based on the information theory and from the viewpoint of optimal information on description—minimum description length, this paper discusses a case: where there is multi-collinearity in the coefficient matrix, principal component estimation is used to estimate and select the original parameters, so as to reduce its multi-collinearity and improve its credibility. From the viewpoint of minimum description length, this paper discusses the approach of selecting principal components and uses this approach to solve a practical problem.展开更多
Based on the influence of circuit element tolerances to the k-fault diagnosis, a method of fault diagnosis is presented which is called minimum tolerance estimation algorithm and has clear physical meaning. Using this...Based on the influence of circuit element tolerances to the k-fault diagnosis, a method of fault diagnosis is presented which is called minimum tolerance estimation algorithm and has clear physical meaning. Using this’method, an effective estimation of the equivalent fault sources can be obtained with less computing time. It is especially worthwhile to point out that an adaptive sub-optimum algorithm, which comes from the above method, requires even less computing-labor and is particularly suitable to more complicated circuits as well as real-time fault location.展开更多
The novel compensating method directly demodulates the signals without the carrier recovery processes, in which the carrier with original modulation frequency is used as the local coherent carrier. In this way, the ph...The novel compensating method directly demodulates the signals without the carrier recovery processes, in which the carrier with original modulation frequency is used as the local coherent carrier. In this way, the phase offsets due to frequency shift are linear. Based on this premise, the compensation processes are: firstly, the phase offsets between the baseband neighbor-symbols after clock recovery is unbiasedly estimated among the reference symbols; then, the receiving signals symbols are adjusted by the phase estimation value; finally, the phase offsets after adjusting are compensated by the least mean squares (LMS) algorithm. In order to express the compensation processes and ability clearly, the quadrature phase shift keying (QPSK) modulation signals are regarded as examples for Matlab simulation. BER simulations are carried out using the Monte-Carlo method. The learning curves are obtained to study the algorithm's convergence ability. The constellation figures are also simulated to observe the compensation results directly.展开更多
The solution of the grey model(GM(1,1)model)generally involves equal-precision observations,and the(co)variance matrix is established from the prior information.However,the data are generally available with unequal-pr...The solution of the grey model(GM(1,1)model)generally involves equal-precision observations,and the(co)variance matrix is established from the prior information.However,the data are generally available with unequal-precision measurements in reality.To deal with the errors of all observations for GM(1,1)model with errors-in-variables(EIV)structure,we exploit the total least-squares(TLS)algorithm to estimate the parameters of GM(1,1)model in this paper.Ignoring that the effect of the improper prior stochastic model and the homologous observations may degrade the accuracy of parameter estimation,we further present a nonlinear total least-squares variance component estimation approach for GM(1,1)model,which resorts to the minimum norm quadratic unbiased estimation(MINQUE).The practical and simulative experiments indicate that the presented approach has significant merits in improving the predictive accuracy in comparison with control methods.展开更多
In multi-user multiple input multiple output (MU-MIMO) systems, the outdated channel state information at the transmit- ter caused by channel time variation has been shown to greatly reduce the achievable ergodic su...In multi-user multiple input multiple output (MU-MIMO) systems, the outdated channel state information at the transmit- ter caused by channel time variation has been shown to greatly reduce the achievable ergodic sum capacity. A simple yet effec- tive solution to this problem is presented by designing a channel extrapolator relying on Karhunen-Loeve (KL) expansion of time- varying channels. In this scheme, channel estimation is done at the base station (BS) rather than at the user terminal (UT), which thereby dispenses the channel parameters feedback from the UT to the BS. Moreover, the inherent channel correlation and the parsimonious parameterization properties of the KL expan- sion are respectively exploited to reduce the channel mismatch error and the computational complexity. Simulations show that the presented scheme outperforms conventional schemes in terms of both channel estimation mean square error (MSE) and ergodic capacity.展开更多
Higher transmission rate is one of the technological features of promi-nently used wireless communication namely Multiple Input Multiple Output-Orthogonal Frequency Division Multiplexing(MIMO–OFDM).One among an effec...Higher transmission rate is one of the technological features of promi-nently used wireless communication namely Multiple Input Multiple Output-Orthogonal Frequency Division Multiplexing(MIMO–OFDM).One among an effective solution for channel estimation in wireless communication system,spe-cifically in different environments is Deep Learning(DL)method.This research greatly utilizes channel estimator on the basis of Convolutional Neural Network Auto Encoder(CNNAE)classifier for MIMO-OFDM systems.A CNNAE classi-fier is one among Deep Learning(DL)algorithm,in which video signal is fed as input by allotting significant learnable weights and biases in various aspects/objects for video signal and capable of differentiating from one another.Improved performances are achieved by using CNNAE based channel estimation,in which extension is done for channel selection as well as achieve enhanced performances numerically,when compared with conventional estimators in quite a lot of scenar-ios.Considering reduction in number of parameters involved and re-usability of weights,CNNAE based channel estimation is quite suitable and properlyfits to the video signal.CNNAE classifier weights updation are done with minimized Sig-nal to Noise Ratio(SNR),Bit Error Rate(BER)and Mean Square Error(MSE).展开更多
The channel estimation technique is investigated in OFDM communication systems with multi-antenna Amplify-and-Forward(AF) relay.The Space-Time Block Code(STBC) is applied at the transmitter of the relay to obtain dive...The channel estimation technique is investigated in OFDM communication systems with multi-antenna Amplify-and-Forward(AF) relay.The Space-Time Block Code(STBC) is applied at the transmitter of the relay to obtain diversity gain.According to the transmission characteristics of OFDM symbols on multiple antennas,a pilot-aided Linear Minimum Mean-Square-Error(LMMSE) channel estimation algorithm with low complexity is designed.Simulation results show that,the proposed LMMSE estimator outperforms least-square estimator and approaches the optimal estimator without error in the performance of Symbol Error Ratio(SER) under several modulation modes,and has a good estimation effect in the realistic relay communication scenario.展开更多
In present paper, the properties of the Bayes Shrinkage estimator is studied for the measure of dispersion of an inverse Gaussian model under the Minimax estimation criteria.
Our purpose is twofold: to present a prototypical example of the conditioning technique to obtain the best estimator of a parameter and to show that th</span><span style="font-family:Verdana;">is...Our purpose is twofold: to present a prototypical example of the conditioning technique to obtain the best estimator of a parameter and to show that th</span><span style="font-family:Verdana;">is technique resides in the structure of an inner product space. Th</span><span style="font-family:Verdana;">e technique uses conditioning </span></span><span style="font-family:Verdana;">of</span><span style="font-family:Verdana;"> an unbiased estimator </span><span style="font-family:Verdana;">on</span><span style="font-family:Verdana;"> a sufficient statistic. This procedure is founded upon the conditional variance formula, which leads to an inner product space and a geometric interpretation. The example clearly illustrates the dependence on the sampling methodology. These advantages show the power and centrality of this process.展开更多
In optical techniques,noise signal is a classical problem in medical image processing.Recently,there has been considerable interest in using the wavelet transform with Bayesian estimation as a powerful tool for recove...In optical techniques,noise signal is a classical problem in medical image processing.Recently,there has been considerable interest in using the wavelet transform with Bayesian estimation as a powerful tool for recovering image from noisy data.In wavelet domain,if Bayesian estimator is used for denoising problem,the solution requires a prior knowledge about the distribution of wavelet coeffcients.Indeed,wavelet coeffcients might be better modeled by super Gaussian density.The super Gaussian density can be generated by Gaussian scale mixture(GSM).So,we present new minimum mean square error(MMSE)estimator for spherically-contoured GSM with Maxwell distribution in additive white Gaussian noise(AWGN).We compare our proposed method to current state-of-the-art method applied on standard test image and we quantify achieved performance improvement.展开更多
The paper aims to discuss three interesting issues of statistical inferences for a common risk ratio (RR) in sparse meta-analysis data. Firstly, the conventional log-risk ratio estimator encounters a number of problem...The paper aims to discuss three interesting issues of statistical inferences for a common risk ratio (RR) in sparse meta-analysis data. Firstly, the conventional log-risk ratio estimator encounters a number of problems when the number of events in the experimental or control group is zero in sparse data of a 2 × 2 table. The adjusted log-risk ratio estimator with the continuity correction points based upon the minimum Bayes risk with respect to the uniform prior density over (0, 1) and the Euclidean loss function is proposed. Secondly, the interest is to find the optimal weights of the pooled estimate that minimize the mean square error (MSE) of subject to the constraint on where , , . Finally, the performance of this minimum MSE weighted estimator adjusted with various values of points is investigated to compare with other popular estimators, such as the Mantel-Haenszel (MH) estimator and the weighted least squares (WLS) estimator (also equivalently known as the inverse-variance weighted estimator) in senses of point estimation and hypothesis testing via simulation studies. The results of estimation illustrate that regardless of the true values of RR, the MH estimator achieves the best performance with the smallest MSE when the study size is rather large and the sample sizes within each study are small. The MSE of WLS estimator and the proposed-weight estimator adjusted by , or , or are close together and they are the best when the sample sizes are moderate to large (and) while the study size is rather small.展开更多
Joint power control has advantages of multi-user detection and power control; and it can combat the multi-access interference and the near-far problem. A novel adaptive joint power control algorithm with channel estim...Joint power control has advantages of multi-user detection and power control; and it can combat the multi-access interference and the near-far problem. A novel adaptive joint power control algorithm with channel estimation in a CDMA cellular system was designed. Simulation results show that the algorithm can control the power not only quickly but also precisely with a time change. The method is useful for increasing system capacity.展开更多
基金supported in part by the National Natural Science Foundation of China(61933007, U21A2019, 62273005, 62273088, 62303301)the Program of Shanghai Academic/Technology Research Leader of China (20XD1420100)+2 种基金the Hainan Province Science and Technology Special Fund of China(ZDYF2022SHFZ105)the Natural Science Foundation of Anhui Province of China (2108085MA07)the Alexander von Humboldt Foundation of Germany。
文摘This paper investigates the anomaly-resistant decentralized state estimation(SE) problem for a class of wide-area power systems which are divided into several non-overlapping areas connected through transmission lines. Two classes of measurements(i.e., local measurements and edge measurements) are obtained, respectively, from the individual area and the transmission lines. A decentralized state estimator, whose performance is resistant against measurement with anomalies, is designed based on the minimum error entropy with fiducial points(MEEF) criterion. Specifically, 1) An augmented model, which incorporates the local prediction and local measurement, is developed by resorting to the unscented transformation approach and the statistical linearization approach;2) Using the augmented model, an MEEF-based cost function is designed that reflects the local prediction errors of the state and the measurement;and 3) The local estimate is first obtained by minimizing the MEEF-based cost function through a fixed-point iteration and then updated by using the edge measuring information. Finally, simulation experiments with three scenarios are carried out on the IEEE 14-bus system to illustrate the validity of the proposed anomaly-resistant decentralized SE scheme.
基金This work is supported by Postgraduate Research&Practice Innovation Program of Jiangsu Province(KYCX18_0467)Jiangsu Province,China.During the revision of this paper,the author is supported by China Scholarship Council(No.201906840021)China to continue some research related to data processing.
文摘Kalman filter is commonly used in data filtering and parameters estimation of nonlinear system,such as projectile's trajectory estimation and control.While there is a drawback that the prior error covariance matrix and filter parameters are difficult to be determined,which may result in filtering divergence.As to the problem that the accuracy of state estimation for nonlinear ballistic model strongly depends on its mathematical model,we improve the weighted least squares method(WLSM)with minimum model error principle.Invariant embedding method is adopted to solve the cost function including the model error.With the knowledge of measurement data and measurement error covariance matrix,we use gradient descent algorithm to determine the weighting matrix of model error.The uncertainty and linearization error of model are recursively estimated by the proposed method,thus achieving an online filtering estimation of the observations.Simulation results indicate that the proposed recursive estimation algorithm is insensitive to initial conditions and of good robustness.
基金supported by the 2011 China Aerospace Science and Technology Foundationthe Certain Ministry Foundation under Grant No.20212HK03010
文摘Performance of the Adaptive Coding and Modulation(ACM) strongly depends on the retrieved Channel State Information(CSI),which can be obtained using the channel estimation techniques relying on pilot symbol transmission.Earlier analysis of methods of pilot-aided channel estimation for ACM systems were relatively little.In this paper,we investigate the performance of CSI prediction using the Minimum Mean Square Error(MMSE)channel estimator for an ACM system.To solve the two problems of MMSE:high computational operations and oversimplified assumption,we then propose the Low-Complexity schemes(LC-MMSE and Recursion LC-MMSE(R-LC-MMSE)).Computational complexity and Mean Square Error(MSE) are presented to evaluate the efficiency of the proposed algorithm.Both analysis and numerical results show that LC-MMSE performs close to the wellknown MMSE estimator with much lower complexity and R-LC-MMSE improves the application of MMSE estimation to specific circumstances.
基金Project(40074001) supported by National Natural Science Foundation of China Project (SD2003 -10) supported by the Open ResearchFund Programof the Key Laboratory of Geomatics and Digital Technilogy ,Shandong Province
文摘Gauss-Markov model is frequently used in data analysis; the analysis and estimation of its parameters is always a hot issue. Based on the information theory and from the viewpoint of optimal information on description—minimum description length, this paper discusses a case: where there is multi-collinearity in the coefficient matrix, principal component estimation is used to estimate and select the original parameters, so as to reduce its multi-collinearity and improve its credibility. From the viewpoint of minimum description length, this paper discusses the approach of selecting principal components and uses this approach to solve a practical problem.
基金Supported by the National Natural Science Foundation of Chilla
文摘Based on the influence of circuit element tolerances to the k-fault diagnosis, a method of fault diagnosis is presented which is called minimum tolerance estimation algorithm and has clear physical meaning. Using this’method, an effective estimation of the equivalent fault sources can be obtained with less computing time. It is especially worthwhile to point out that an adaptive sub-optimum algorithm, which comes from the above method, requires even less computing-labor and is particularly suitable to more complicated circuits as well as real-time fault location.
基金supported by the National Natural Science Foundation of China(60532030)
文摘The novel compensating method directly demodulates the signals without the carrier recovery processes, in which the carrier with original modulation frequency is used as the local coherent carrier. In this way, the phase offsets due to frequency shift are linear. Based on this premise, the compensation processes are: firstly, the phase offsets between the baseband neighbor-symbols after clock recovery is unbiasedly estimated among the reference symbols; then, the receiving signals symbols are adjusted by the phase estimation value; finally, the phase offsets after adjusting are compensated by the least mean squares (LMS) algorithm. In order to express the compensation processes and ability clearly, the quadrature phase shift keying (QPSK) modulation signals are regarded as examples for Matlab simulation. BER simulations are carried out using the Monte-Carlo method. The learning curves are obtained to study the algorithm's convergence ability. The constellation figures are also simulated to observe the compensation results directly.
基金supported by the National Natural Science Foundation of China(No.41874001 and No.41664001)Support Program for Outstanding Youth Talents in Jiangxi Province(No.20162BCB23050)National Key Research and Development Program(No.2016YFB0501405)。
文摘The solution of the grey model(GM(1,1)model)generally involves equal-precision observations,and the(co)variance matrix is established from the prior information.However,the data are generally available with unequal-precision measurements in reality.To deal with the errors of all observations for GM(1,1)model with errors-in-variables(EIV)structure,we exploit the total least-squares(TLS)algorithm to estimate the parameters of GM(1,1)model in this paper.Ignoring that the effect of the improper prior stochastic model and the homologous observations may degrade the accuracy of parameter estimation,we further present a nonlinear total least-squares variance component estimation approach for GM(1,1)model,which resorts to the minimum norm quadratic unbiased estimation(MINQUE).The practical and simulative experiments indicate that the presented approach has significant merits in improving the predictive accuracy in comparison with control methods.
基金supported by the National Natural Science Foundation of China (6096200161071088)+2 种基金the Natural Science Foundation of Fujian Province of China (2012J05119)the Fundamental Research Funds for the Central Universities (11QZR02)the Research Fund of Guangxi Key Lab of Wireless Wideband Communication & Signal Processing (21104)
文摘In multi-user multiple input multiple output (MU-MIMO) systems, the outdated channel state information at the transmit- ter caused by channel time variation has been shown to greatly reduce the achievable ergodic sum capacity. A simple yet effec- tive solution to this problem is presented by designing a channel extrapolator relying on Karhunen-Loeve (KL) expansion of time- varying channels. In this scheme, channel estimation is done at the base station (BS) rather than at the user terminal (UT), which thereby dispenses the channel parameters feedback from the UT to the BS. Moreover, the inherent channel correlation and the parsimonious parameterization properties of the KL expan- sion are respectively exploited to reduce the channel mismatch error and the computational complexity. Simulations show that the presented scheme outperforms conventional schemes in terms of both channel estimation mean square error (MSE) and ergodic capacity.
文摘Higher transmission rate is one of the technological features of promi-nently used wireless communication namely Multiple Input Multiple Output-Orthogonal Frequency Division Multiplexing(MIMO–OFDM).One among an effective solution for channel estimation in wireless communication system,spe-cifically in different environments is Deep Learning(DL)method.This research greatly utilizes channel estimator on the basis of Convolutional Neural Network Auto Encoder(CNNAE)classifier for MIMO-OFDM systems.A CNNAE classi-fier is one among Deep Learning(DL)algorithm,in which video signal is fed as input by allotting significant learnable weights and biases in various aspects/objects for video signal and capable of differentiating from one another.Improved performances are achieved by using CNNAE based channel estimation,in which extension is done for channel selection as well as achieve enhanced performances numerically,when compared with conventional estimators in quite a lot of scenar-ios.Considering reduction in number of parameters involved and re-usability of weights,CNNAE based channel estimation is quite suitable and properlyfits to the video signal.CNNAE classifier weights updation are done with minimized Sig-nal to Noise Ratio(SNR),Bit Error Rate(BER)and Mean Square Error(MSE).
基金Supported by the National Natural Science Foundation of Jiangsu province(No.08KJB510015)
文摘The channel estimation technique is investigated in OFDM communication systems with multi-antenna Amplify-and-Forward(AF) relay.The Space-Time Block Code(STBC) is applied at the transmitter of the relay to obtain diversity gain.According to the transmission characteristics of OFDM symbols on multiple antennas,a pilot-aided Linear Minimum Mean-Square-Error(LMMSE) channel estimation algorithm with low complexity is designed.Simulation results show that,the proposed LMMSE estimator outperforms least-square estimator and approaches the optimal estimator without error in the performance of Symbol Error Ratio(SER) under several modulation modes,and has a good estimation effect in the realistic relay communication scenario.
文摘In present paper, the properties of the Bayes Shrinkage estimator is studied for the measure of dispersion of an inverse Gaussian model under the Minimax estimation criteria.
文摘Our purpose is twofold: to present a prototypical example of the conditioning technique to obtain the best estimator of a parameter and to show that th</span><span style="font-family:Verdana;">is technique resides in the structure of an inner product space. Th</span><span style="font-family:Verdana;">e technique uses conditioning </span></span><span style="font-family:Verdana;">of</span><span style="font-family:Verdana;"> an unbiased estimator </span><span style="font-family:Verdana;">on</span><span style="font-family:Verdana;"> a sufficient statistic. This procedure is founded upon the conditional variance formula, which leads to an inner product space and a geometric interpretation. The example clearly illustrates the dependence on the sampling methodology. These advantages show the power and centrality of this process.
文摘In optical techniques,noise signal is a classical problem in medical image processing.Recently,there has been considerable interest in using the wavelet transform with Bayesian estimation as a powerful tool for recovering image from noisy data.In wavelet domain,if Bayesian estimator is used for denoising problem,the solution requires a prior knowledge about the distribution of wavelet coeffcients.Indeed,wavelet coeffcients might be better modeled by super Gaussian density.The super Gaussian density can be generated by Gaussian scale mixture(GSM).So,we present new minimum mean square error(MMSE)estimator for spherically-contoured GSM with Maxwell distribution in additive white Gaussian noise(AWGN).We compare our proposed method to current state-of-the-art method applied on standard test image and we quantify achieved performance improvement.
文摘The paper aims to discuss three interesting issues of statistical inferences for a common risk ratio (RR) in sparse meta-analysis data. Firstly, the conventional log-risk ratio estimator encounters a number of problems when the number of events in the experimental or control group is zero in sparse data of a 2 × 2 table. The adjusted log-risk ratio estimator with the continuity correction points based upon the minimum Bayes risk with respect to the uniform prior density over (0, 1) and the Euclidean loss function is proposed. Secondly, the interest is to find the optimal weights of the pooled estimate that minimize the mean square error (MSE) of subject to the constraint on where , , . Finally, the performance of this minimum MSE weighted estimator adjusted with various values of points is investigated to compare with other popular estimators, such as the Mantel-Haenszel (MH) estimator and the weighted least squares (WLS) estimator (also equivalently known as the inverse-variance weighted estimator) in senses of point estimation and hypothesis testing via simulation studies. The results of estimation illustrate that regardless of the true values of RR, the MH estimator achieves the best performance with the smallest MSE when the study size is rather large and the sample sizes within each study are small. The MSE of WLS estimator and the proposed-weight estimator adjusted by , or , or are close together and they are the best when the sample sizes are moderate to large (and) while the study size is rather small.
文摘Joint power control has advantages of multi-user detection and power control; and it can combat the multi-access interference and the near-far problem. A novel adaptive joint power control algorithm with channel estimation in a CDMA cellular system was designed. Simulation results show that the algorithm can control the power not only quickly but also precisely with a time change. The method is useful for increasing system capacity.
基金Supported by the National Natural Science Foundation of China (No. 61001105), the National Science and Technology Major Projects (No. 2011ZX03001- 007- 03) and Beijing Natural Science Foundation (No. 4102043).