This paper addresses the issues of channel estimation in a Multiple-Input/Multiple-Output (MIMO) system. Markov Chain Monte Carlo (MCMC) method is employed to jointly estimate the Channel State Information (CSI) and t...This paper addresses the issues of channel estimation in a Multiple-Input/Multiple-Output (MIMO) system. Markov Chain Monte Carlo (MCMC) method is employed to jointly estimate the Channel State Information (CSI) and the transmitted signals. The deduced algorithms can work well under circumstances of low Signal-to-Noise Ratio (SNR). Simulation results are presented to demonstrate their effectiveness.展开更多
The application of Monte Carlo method in estimating rate constants for polymerization was described, A general program for Monte Carlo simulation was determined first according to the elementary reactions, after which...The application of Monte Carlo method in estimating rate constants for polymerization was described, A general program for Monte Carlo simulation was determined first according to the elementary reactions, after which the rate constants could be automatically adjusted and optimized through comparing of experimental and simulated data with an error expression that meeted a given minimum criterion. Such a process made the rate constants to be estimated without kinetic model in advance. The technique was applied to estimate the rate constants of the bulk polymerization of styrene catalyzed by the rare earth catalyst. The esthnated results showed the Monte Carlo method was feasible and effective for estimating rate constants in polymerization engineering.展开更多
Optimal detection of liquid ionization calorimeter signal in experimental particle physics is considered. A few linear and nonlinear approaches for amplitude and arrival time estimation based on the χ2 function are c...Optimal detection of liquid ionization calorimeter signal in experimental particle physics is considered. A few linear and nonlinear approaches for amplitude and arrival time estimation based on the χ2 function are compared in simulation considering the noise sample correlation introduced by the analog pulse shaper. The estimation bias of the first-order approximation, a.k.a linear optimal filtering, is studied and contrasted to those of the second-order as well as the exhaustive search. A gradient-descent technique is presented as an alternative to the exhaustive search with significantly reduced search time and computation complexity. Results from various pulse shapers including the CR-RC2, CR-RC3, and CR2-RC2 are also compared.展开更多
In this paper, an adaptive estimation algorithm is proposed for non-linear dynamic systems with unknown static parameters based on combination of particle filtering and Simultaneous Perturbation Stochastic Approxi- ma...In this paper, an adaptive estimation algorithm is proposed for non-linear dynamic systems with unknown static parameters based on combination of particle filtering and Simultaneous Perturbation Stochastic Approxi- mation (SPSA) technique. The estimations of parameters are obtained by maximum-likelihood estimation and sampling within particle filtering framework, and the SPSA is used for stochastic optimization and to approximate the gradient of the cost function. The proposed algorithm achieves combined estimation of dynamic state and static parameters of nonlinear systems. Simulation result demonstrates the feasibilitv and efficiency of the proposed algorithm展开更多
In this paper, we use the lower record values from the inverse Weibull distribution (IWD) to develop and discuss different methods of estimation in two different cases, 1) when the shape parameter is known and 2) when...In this paper, we use the lower record values from the inverse Weibull distribution (IWD) to develop and discuss different methods of estimation in two different cases, 1) when the shape parameter is known and 2) when both of the shape and scale parameters are unknown. First, we derive the best linear unbiased estimate (BLUE) of the scale parameter of the IWD. To compare the different methods of estimation, we present the results of Sultan (2007) for calculating the best linear unbiased estimates (BLUEs) of the location and scale parameters of IWD. Second, we derive the maximum likelihood estimates (MLEs) of the location and scale parameters. Further, we discuss some properties of the MLEs of the location and scale parameters. To compare the different estimates we calculate the relative efficiency between the obtained estimates. Finally, we propose some numerical illustrations by using Monte Carlo simulations and apply the findings of the paper to some simulated data.展开更多
Hydrocracking is a catalytic reaction process in the petroleum refineries for converting the higher boiling temperature residue of crude oil into a lighter fraction of hydrocarbons such as gasoline and diesel. In this...Hydrocracking is a catalytic reaction process in the petroleum refineries for converting the higher boiling temperature residue of crude oil into a lighter fraction of hydrocarbons such as gasoline and diesel. In this study, a modified continuous lumping kinetic approach is applied to model the hydro-cracking of vacuum gas oil. The model is modified to take into consideration the reactor temperature on the reaction yield distribution. The model is calibrated by maximizing the likelihood function between the modeled and measured data at four different reactor temperatures. Bayesian approach parameter estimation is also applied to obtain the confidence interval of model parameters by considering the uncertainty associated with the measured errors and the model structural errors. Then Monte Carlo simulation is applied to the posterior range of the model parameters to obtain the 95% confidence interval of the model outputs for each individual fraction of the hydrocracking products. A good agreement is observed between the output of the calibrated model and the measured data points. The Bayesian approach based on the Markov Chain Monte Carlo simulation is shown to be efficient to quantify the uncertainty associated with the parameter values of the continuous lumping model.展开更多
In this paper, an importance sampling maximum likelihood(ISML) estimator for direction-of-arrival(DOA) of incoherently distributed(ID) sources is proposed. Starting from the maximum likelihood estimation description o...In this paper, an importance sampling maximum likelihood(ISML) estimator for direction-of-arrival(DOA) of incoherently distributed(ID) sources is proposed. Starting from the maximum likelihood estimation description of the uniform linear array(ULA), a decoupled concentrated likelihood function(CLF) is presented. A new objective function based on CLF which can obtain a closed-form solution of global maximum is constructed according to Pincus theorem. To obtain the optimal value of the objective function which is a complex high-dimensional integral,we propose an importance sampling approach based on Monte Carlo random calculation. Next, an importance function is derived, which can simplify the problem of generating random vector from a high-dimensional probability density function(PDF) to generate random variable from a one-dimensional PDF. Compared with the existing maximum likelihood(ML) algorithms for DOA estimation of ID sources, the proposed algorithm does not require initial estimates, and its performance is closer to CramerRao lower bound(CRLB). The proposed algorithm performs better than the existing methods when the interval between sources to be estimated is small and in low signal to noise ratio(SNR)scenarios.展开更多
A diagnostic procedure based on maximum likelihood estimation, to study the convergence of the Markov chain produced by Gibbs sampler, is presented. The unbiasedness, consistent and asymptotic normality are considered...A diagnostic procedure based on maximum likelihood estimation, to study the convergence of the Markov chain produced by Gibbs sampler, is presented. The unbiasedness, consistent and asymptotic normality are considered for the estimation of the parameters produced by the procedure. An example is provided to illustrate the procedure, and the numerical result is consistent with the theoretical one.展开更多
This paper presents augmented input estimation(AIE)for multiple maneuvering target tracking.Multi-target tracking(MTT)is based on two main parts,data association and estimation.In data association(DA),the best observa...This paper presents augmented input estimation(AIE)for multiple maneuvering target tracking.Multi-target tracking(MTT)is based on two main parts,data association and estimation.In data association(DA),the best observations are assigned to the considered tracks.In real conditions,the number of observations is more than targets and also locations of observations are often so scattered that the association between targets and observations cannot be done simply.In this case,for general MTT problems with unknown numbers of targets,we present a Markov chain Monte-Carlo DA(MCMCDA)algorithm that approximates the optimal Bayesian filter with low complexity in computations.After DA,estimation and tracking should be done.Since in general cases,many targets can have maneuvering motions,then AIE is proposed to cover both the non-maneuvering and maneuvering parts of motion and the maneuver detection procedure is eliminated.This model with an input estimation(IE)approach is a special augmentation in the state space model which considers both the state vector and the unknown input vector as a new augmented state vector.Some comparisons based on the Monte-Carlo simulations are also made to evaluate the performances of the proposed method and other older methods in MTT.展开更多
A particle filter is proposed to perform joint estimation of the carrier frequency offset (CFO) and the channel in multiple-input multiple-output orthogonal frequency division multiplexing (MIMO-OFDM) wireless com...A particle filter is proposed to perform joint estimation of the carrier frequency offset (CFO) and the channel in multiple-input multiple-output orthogonal frequency division multiplexing (MIMO-OFDM) wireless communication systems. It marginalizes out the channel parameters from the sampling space in sequential importance sampling (SIS), and propagates them with the Kalman filter. Then the importance weights of the CFO particles are evaluated according to the imaginary part of the error between measurement and estimation. The varieties of particles are maintained by sequential importance resampling (SIR). Simulation results demonstrate this algorithm can estimate the CFO and the channel parameters with high accuracy. At the same time, some robustness is kept when the channel model has small variations.展开更多
This paper deals with the Bayesian inferences of unknown parameters of the progressively Type II censored Weibull-geometric (WG) distribution. The Bayes estimators cannot be obtained in explicit forms of the unknown p...This paper deals with the Bayesian inferences of unknown parameters of the progressively Type II censored Weibull-geometric (WG) distribution. The Bayes estimators cannot be obtained in explicit forms of the unknown parameters under a squared error loss function. The approximate Bayes estimators will be computed using the idea of Markov Chain Monte Carlo (MCMC) method to generate from the posterior distributions. Also the point estimation and confidence intervals based on maximum likelihood and bootstrap technique are also proposed. The approximate Bayes estimators will be obtained under the assumptions of informative and non-informative priors are compared with the maximum likelihood estimators. A numerical example is provided to illustrate the proposed estimation methods here. Maximum likelihood, bootstrap and the different Bayes estimates are compared via a Monte Carlo Simulation展开更多
A multi-constituent water quality model is presented,Which relates carbonaceous biochemical oxygen demand (CBOD),amonia (NH3-N), nitrite(NO2-N), nitrate(NO3-N) and dissolvedoxygen(DO). The parameters are solved by Mar...A multi-constituent water quality model is presented,Which relates carbonaceous biochemical oxygen demand (CBOD),amonia (NH3-N), nitrite(NO2-N), nitrate(NO3-N) and dissolvedoxygen(DO). The parameters are solved by Marquardt Method (i. e.,Dampled Least Square Method) while initial values inoptimization are produced by Monte-Carlo Method. The Potential ofthe method as a parameter estimation aid is demonstrated for theapplication to the Liangyi Rver, JiangSu Province of China and by aspecial comparison with Gauss Method.展开更多
The Generalized Markov Fluid Model(GMFM)is assumed for modeling sources in the network because it is versatile to describe the traffic fluctuations.In order to estimate resources allocations or in other words the chan...The Generalized Markov Fluid Model(GMFM)is assumed for modeling sources in the network because it is versatile to describe the traffic fluctuations.In order to estimate resources allocations or in other words the channel occupation of each source,the concept of effective bandwidth(EB)proposed by Kelly is used.In this paper we use an expression to determine the EB for this model which is of particular interest because it allows expressing said magnitude depending on the parameters of the model.This paper provides EB estimates for this model applying Kernel Estimation techniques in data networking.In particular we will study two differentiated cases:dispatches following a Gaussian and Exponential distribution.The performance of the proposed method is analyzed using simulated traffic traces generated by Monte Carlo Markov Chain algorithms.The estimation process worked much better in the Gaussian distribution case than in the Exponential one.展开更多
Cone penetration testing (CPT) is an extensively utilized and cost effective tool for geotechnical site characterization. CPT consists of pushing at a constant rate an electronic cone into penetrable soils and recordi...Cone penetration testing (CPT) is an extensively utilized and cost effective tool for geotechnical site characterization. CPT consists of pushing at a constant rate an electronic cone into penetrable soils and recording the resistance to the cone tip (q<sub>c</sub> value). The measured q<sub>c</sub> values (after correction for the pore water pressure) are utilized to estimate soil type and associated soil properties based predominantly on empirical correlations. The most common cone tips have associated areas of 10 cm<sup>2</sup> and 15 cm<sup>2</sup>. Investigators also utilized significantly larger cone tips (33 cm<sup>2</sup> and 40 cm<sup>2</sup>) so that gravelly soils can be penetrated. Small cone tips (2 cm<sup>2</sup> and 5 cm<sup>2</sup>) are utilized for shallow soil investigations. The cone tip resistance measured at a particular depth is affected by the values above and below the depth of interest which results in a smoothing or blurring of the true bearing values. Extensive work has been carried out in mathematically modelling the smoothing function which results in the blurred cone bearing measurements. This paper outlines a technique which facilitates estimating the dominant parameters of the cone smoothing function from processing real cone bearing data sets. This cone calibration technique is referred to as the so-called CPSPE algorithm. The mathematical details of the CPSPE algorithm are outlined in this paper along with the results from a challenging test bed simulation.展开更多
文摘This paper addresses the issues of channel estimation in a Multiple-Input/Multiple-Output (MIMO) system. Markov Chain Monte Carlo (MCMC) method is employed to jointly estimate the Channel State Information (CSI) and the transmitted signals. The deduced algorithms can work well under circumstances of low Signal-to-Noise Ratio (SNR). Simulation results are presented to demonstrate their effectiveness.
基金Project supported by National Natural Science Foundation of China (Grant No, 20406016)
文摘The application of Monte Carlo method in estimating rate constants for polymerization was described, A general program for Monte Carlo simulation was determined first according to the elementary reactions, after which the rate constants could be automatically adjusted and optimized through comparing of experimental and simulated data with an error expression that meeted a given minimum criterion. Such a process made the rate constants to be estimated without kinetic model in advance. The technique was applied to estimate the rate constants of the bulk polymerization of styrene catalyzed by the rare earth catalyst. The esthnated results showed the Monte Carlo method was feasible and effective for estimating rate constants in polymerization engineering.
文摘Optimal detection of liquid ionization calorimeter signal in experimental particle physics is considered. A few linear and nonlinear approaches for amplitude and arrival time estimation based on the χ2 function are compared in simulation considering the noise sample correlation introduced by the analog pulse shaper. The estimation bias of the first-order approximation, a.k.a linear optimal filtering, is studied and contrasted to those of the second-order as well as the exhaustive search. A gradient-descent technique is presented as an alternative to the exhaustive search with significantly reduced search time and computation complexity. Results from various pulse shapers including the CR-RC2, CR-RC3, and CR2-RC2 are also compared.
基金the National Natural Science Foundation of China (No. 60404011)
文摘In this paper, an adaptive estimation algorithm is proposed for non-linear dynamic systems with unknown static parameters based on combination of particle filtering and Simultaneous Perturbation Stochastic Approxi- mation (SPSA) technique. The estimations of parameters are obtained by maximum-likelihood estimation and sampling within particle filtering framework, and the SPSA is used for stochastic optimization and to approximate the gradient of the cost function. The proposed algorithm achieves combined estimation of dynamic state and static parameters of nonlinear systems. Simulation result demonstrates the feasibilitv and efficiency of the proposed algorithm
文摘In this paper, we use the lower record values from the inverse Weibull distribution (IWD) to develop and discuss different methods of estimation in two different cases, 1) when the shape parameter is known and 2) when both of the shape and scale parameters are unknown. First, we derive the best linear unbiased estimate (BLUE) of the scale parameter of the IWD. To compare the different methods of estimation, we present the results of Sultan (2007) for calculating the best linear unbiased estimates (BLUEs) of the location and scale parameters of IWD. Second, we derive the maximum likelihood estimates (MLEs) of the location and scale parameters. Further, we discuss some properties of the MLEs of the location and scale parameters. To compare the different estimates we calculate the relative efficiency between the obtained estimates. Finally, we propose some numerical illustrations by using Monte Carlo simulations and apply the findings of the paper to some simulated data.
文摘Hydrocracking is a catalytic reaction process in the petroleum refineries for converting the higher boiling temperature residue of crude oil into a lighter fraction of hydrocarbons such as gasoline and diesel. In this study, a modified continuous lumping kinetic approach is applied to model the hydro-cracking of vacuum gas oil. The model is modified to take into consideration the reactor temperature on the reaction yield distribution. The model is calibrated by maximizing the likelihood function between the modeled and measured data at four different reactor temperatures. Bayesian approach parameter estimation is also applied to obtain the confidence interval of model parameters by considering the uncertainty associated with the measured errors and the model structural errors. Then Monte Carlo simulation is applied to the posterior range of the model parameters to obtain the 95% confidence interval of the model outputs for each individual fraction of the hydrocracking products. A good agreement is observed between the output of the calibrated model and the measured data points. The Bayesian approach based on the Markov Chain Monte Carlo simulation is shown to be efficient to quantify the uncertainty associated with the parameter values of the continuous lumping model.
基金supported by the basic research program of Natural Science in Shannxi province of China (2021JQ-369)。
文摘In this paper, an importance sampling maximum likelihood(ISML) estimator for direction-of-arrival(DOA) of incoherently distributed(ID) sources is proposed. Starting from the maximum likelihood estimation description of the uniform linear array(ULA), a decoupled concentrated likelihood function(CLF) is presented. A new objective function based on CLF which can obtain a closed-form solution of global maximum is constructed according to Pincus theorem. To obtain the optimal value of the objective function which is a complex high-dimensional integral,we propose an importance sampling approach based on Monte Carlo random calculation. Next, an importance function is derived, which can simplify the problem of generating random vector from a high-dimensional probability density function(PDF) to generate random variable from a one-dimensional PDF. Compared with the existing maximum likelihood(ML) algorithms for DOA estimation of ID sources, the proposed algorithm does not require initial estimates, and its performance is closer to CramerRao lower bound(CRLB). The proposed algorithm performs better than the existing methods when the interval between sources to be estimated is small and in low signal to noise ratio(SNR)scenarios.
文摘A diagnostic procedure based on maximum likelihood estimation, to study the convergence of the Markov chain produced by Gibbs sampler, is presented. The unbiasedness, consistent and asymptotic normality are considered for the estimation of the parameters produced by the procedure. An example is provided to illustrate the procedure, and the numerical result is consistent with the theoretical one.
文摘This paper presents augmented input estimation(AIE)for multiple maneuvering target tracking.Multi-target tracking(MTT)is based on two main parts,data association and estimation.In data association(DA),the best observations are assigned to the considered tracks.In real conditions,the number of observations is more than targets and also locations of observations are often so scattered that the association between targets and observations cannot be done simply.In this case,for general MTT problems with unknown numbers of targets,we present a Markov chain Monte-Carlo DA(MCMCDA)algorithm that approximates the optimal Bayesian filter with low complexity in computations.After DA,estimation and tracking should be done.Since in general cases,many targets can have maneuvering motions,then AIE is proposed to cover both the non-maneuvering and maneuvering parts of motion and the maneuver detection procedure is eliminated.This model with an input estimation(IE)approach is a special augmentation in the state space model which considers both the state vector and the unknown input vector as a new augmented state vector.Some comparisons based on the Monte-Carlo simulations are also made to evaluate the performances of the proposed method and other older methods in MTT.
基金Project supported by the National Natural Science Foundation of China (Grant No.60572157)the International Cooper-ation Foundation (Grant No.2008DFA11950)
文摘A particle filter is proposed to perform joint estimation of the carrier frequency offset (CFO) and the channel in multiple-input multiple-output orthogonal frequency division multiplexing (MIMO-OFDM) wireless communication systems. It marginalizes out the channel parameters from the sampling space in sequential importance sampling (SIS), and propagates them with the Kalman filter. Then the importance weights of the CFO particles are evaluated according to the imaginary part of the error between measurement and estimation. The varieties of particles are maintained by sequential importance resampling (SIR). Simulation results demonstrate this algorithm can estimate the CFO and the channel parameters with high accuracy. At the same time, some robustness is kept when the channel model has small variations.
文摘This paper deals with the Bayesian inferences of unknown parameters of the progressively Type II censored Weibull-geometric (WG) distribution. The Bayes estimators cannot be obtained in explicit forms of the unknown parameters under a squared error loss function. The approximate Bayes estimators will be computed using the idea of Markov Chain Monte Carlo (MCMC) method to generate from the posterior distributions. Also the point estimation and confidence intervals based on maximum likelihood and bootstrap technique are also proposed. The approximate Bayes estimators will be obtained under the assumptions of informative and non-informative priors are compared with the maximum likelihood estimators. A numerical example is provided to illustrate the proposed estimation methods here. Maximum likelihood, bootstrap and the different Bayes estimates are compared via a Monte Carlo Simulation
文摘A multi-constituent water quality model is presented,Which relates carbonaceous biochemical oxygen demand (CBOD),amonia (NH3-N), nitrite(NO2-N), nitrate(NO3-N) and dissolvedoxygen(DO). The parameters are solved by Marquardt Method (i. e.,Dampled Least Square Method) while initial values inoptimization are produced by Monte-Carlo Method. The Potential ofthe method as a parameter estimation aid is demonstrated for theapplication to the Liangyi Rver, JiangSu Province of China and by aspecial comparison with Gauss Method.
文摘The Generalized Markov Fluid Model(GMFM)is assumed for modeling sources in the network because it is versatile to describe the traffic fluctuations.In order to estimate resources allocations or in other words the channel occupation of each source,the concept of effective bandwidth(EB)proposed by Kelly is used.In this paper we use an expression to determine the EB for this model which is of particular interest because it allows expressing said magnitude depending on the parameters of the model.This paper provides EB estimates for this model applying Kernel Estimation techniques in data networking.In particular we will study two differentiated cases:dispatches following a Gaussian and Exponential distribution.The performance of the proposed method is analyzed using simulated traffic traces generated by Monte Carlo Markov Chain algorithms.The estimation process worked much better in the Gaussian distribution case than in the Exponential one.
文摘Cone penetration testing (CPT) is an extensively utilized and cost effective tool for geotechnical site characterization. CPT consists of pushing at a constant rate an electronic cone into penetrable soils and recording the resistance to the cone tip (q<sub>c</sub> value). The measured q<sub>c</sub> values (after correction for the pore water pressure) are utilized to estimate soil type and associated soil properties based predominantly on empirical correlations. The most common cone tips have associated areas of 10 cm<sup>2</sup> and 15 cm<sup>2</sup>. Investigators also utilized significantly larger cone tips (33 cm<sup>2</sup> and 40 cm<sup>2</sup>) so that gravelly soils can be penetrated. Small cone tips (2 cm<sup>2</sup> and 5 cm<sup>2</sup>) are utilized for shallow soil investigations. The cone tip resistance measured at a particular depth is affected by the values above and below the depth of interest which results in a smoothing or blurring of the true bearing values. Extensive work has been carried out in mathematically modelling the smoothing function which results in the blurred cone bearing measurements. This paper outlines a technique which facilitates estimating the dominant parameters of the cone smoothing function from processing real cone bearing data sets. This cone calibration technique is referred to as the so-called CPSPE algorithm. The mathematical details of the CPSPE algorithm are outlined in this paper along with the results from a challenging test bed simulation.