estimated pulse sequene we produe estimates of the delay time struture of ripple-red events. Formulation of the Problem We model the k th seismi trae
|
|
- Brianna Jones
- 5 years ago
- Views:
Transcription
1 Bayesian Deonvolution of Seismi Array Data Using the Gibbs Sampler Eri A. Suess Robert Shumway, University of California, Davis Rong Chen, Texas A&M University Contat: Eri A. Suess, Division of Statistis, One Shields Avenue, Davis, CA 9566 ( Key Words: Bayesian analysis, Gibbs Sampling, time series, losions, ripple-ring, nulear test monitoring. Abstrat The problem of monitoring for low magnitude nulear losions, using seismi array data, under a Comprehensive Test Ban Treaty (CTBT requires a apability for distinguishing nulear losions from other seismi events. Industrial mining losions are one type of seismi event that needs to be ruled out when trying to detet nulear tests. We onsider a Bayesian approah to the problem of deteting ripple-red mining losions. Seismi array data are ressed as multidimensional time domain onvolutions of an unknown pulse funtion, representing the ripple-red delay pattern, with unknown signalpath eets sequenes on eah hannel, whih are assumed to follow independent AR proesses. Using the Gibbs Sampler a proposal of Cheng, Chen Li [] for blind deonvolution, we develop an approah to estimating the delay parameters the unknown signal-path eets sequene at eah sensor. Results for a ripple-red mining losion reorded at the Arti Experimental Seismi Station (ARCESS will be presented. Finally, the impliations of our method for monitoring a nulear test ban treaty are onsidered. Introdution The problem of monitoring seismi events for possible nulear losions is an important one that has been studied extensively is still open to further development. Muh work has been done on this problem due to the many treaties that have been signed in the past between the U.S., the former U.S.S.R., other nulear powers related to testing issues. The fous in the past has been on distinguishing possible nulear losions from earthquakes. Currently, beause the above treaties have put limitations on the permissible sizes of the nulear losions for testing other smaller seismi events suh as mining losions have beome of interest in the disrimination problem. The work we are proeeding with here is related to distinguishing low level nulear losions from ripple-red mining losions that are on the same seismi level. This tehnique of quarry mining involves the use of multiple rows of losives that are detonated with approximately equal time delays between the rows. The objetive of this paper is to develop an estimation proedure that estimates delay times in seismi events. The delay struture of a single seismi event is haraterized by the size (or amplitude of the delayed losions the delay times between the single losions. If a delay struture an be estimated, this will be useful in disriminating ripple- red mining losions from other, possibly nulear, seismi events. In addressing this problem, we have developed a model for seismi reordings olleted by an array of sensors for ripple-red events a deonvolution method based on work by Rong Chen, in the univariate ase, work by Cheng, Chen, Li [] that is used to estimate parameters in our model related to the delay time struture. Our method works by separating, or deonvolving a pulse sequene from the underlying signal-path effets sequenes. The pulse sequene is used to model the detonation design of a mining losion. The underlying signal-path eets sequenes are used to model the ommon underlying signal sent by eah detonation in a ripple-red event the eets of the path on the signal as it traes to the sensors. We approah the problem from a Bayesian perspetive. The Gibbs sampler is implemented to produe posterior estimates of the parameters in our model, using the data available prior information. From the
2 estimated pulse sequene we produe estimates of the delay time struture of ripple-red events. Formulation of the Problem We model the k th seismi trae y k (t, t ::: n,in an array of q-sensors, that is suspeted to have been produed by a ripple-red mining losion, using the following multivariate onvolution model: y k (t m j0 a j s k (t ; j+" k (t: ( We refer to the vetor of parameters a [a ::: a m ] 0 as the relative pulse sequene. This set of parameters is assumed to reet the delay pattern used in the design of the ripple-red event. We dene the k th underlying signal-path eets sequene, s k (t, t ::: n, as the ombination of the \signal" produed by eah sub-losion in a ripple-red event along with the \path eets" that are due to the rom imperfetions in the material the signal travels through on the k th hannel to the k th seismometer. We assume that the pulse sequene a is independent of the underlying signal-path eets sequene s k (t, t ::: n, on eah hannel k, k :: q. Assuming that the rst losion in the ripple-red event is used as the referene signal, for whih wex a 0 to be, that the remaining a j are alulated relative toa 0, the model an be written in a more insightful form as y k (t s k (t+ m j a j s k (t ; j+" k (t: ( This model assumes that eah delayed losion in a ripple-red event sends a signal having the same form as the rst. The only dierene is that its amplitudes a j may vary. Further modeling assumptions on the pulse sequene a, the underlying signal-path eets sequenes s k (t, k :: q, are needed to insure that the parameters in the onvolution model are identi- able. To model the pulse sequene we useavari- ant of the Bernoulli-Gaussian model ommonly used to model reetion seismology data. A disussion of the model an be found in Mendel [4]. For the pulse sequene a, wexa 0 model the relative amplitudes as follows: We dene j, j ::: m,tobe a sequene of independent Bernoulli rom variables with p( j ; p( j 0. The nonzero values of j index the possible detonation times in the ripple-red event of duration m, measured in points per seond. The parameter is the probability of observing a zero value in the pulse sequene. The parameter m is the length of the vetor a, m is assumed to be known. Next, we dene the distribution of eah a j, onditional on j, j ::: m. In ripple-red events the amplitudes will be positive beause mining losions are usually performed at shallow depths. We dene a j, j ::: m, onditional on j, j ::: m, as p(a j j j 0ifa j 0 0 otherwise. And p(a j j j ( p (3 where ( Z 0 ; (a j ; p ; (x ; I(a j > 0 dx: I(a j > 0 0 if a j > 0 0 otherwise. The trunated normal distribution is used beause the a j, the relative amplitude of the delayed detonation at time j, is assumed to be positive rom. Trunating at zero restrits the values of a j to positive values. It is reasonable to use the trunated normal distribution to model the variation in the size of a single detonation sine the distribution of the size of a single losion should be entered at a mean value with a variane. Alternatively, we write the density of the independent a j, j ::: m,as p(a j j I(a j 0+(;p(a j j j I(a j > 0 (4 whih is a mixture of a Bernoulli distribution a trunated Gaussian distribution. The mean is used to model the average size of the delayed losions in ripple-red events in the monitoring region the variane is used to model the variane of suh losions. So the a j are independent identially distributed with a density dened in (4. Note that for the Bernoulli sequene j, j ::: m, will have n 0values m ; n unit values, whih implies that the pulse sequene has m ; n nonzero values or pulses. Note that n is rom. Dene J to be the index set that inludes the values of j for nonzero values of j or thereby the nonzero values of a j. Hene, we dene the density ofa [a ::: a m ] 0 as p(aj n ( ; m;n ;n ( ( ; (a j ; J n p : (5 This modeling assumption is based on the Bernoulli- Gaussian model used by Cheng, Chen, Li [].
3 To model the underlying signal-path eets sequene s k (t, in the k th seismi trae y k (t, t ::: n, we use a low-order autoregressive, AR(p, model 0 s k (t+ s k (t ; + ::: + p s k (t ; p e k (t (6 where wex 0 let e k (t be independent identially distributed, over t k, normal with mean 0 variane. The use of this model is justied by Dargahi-Noubary [] Tjstheim [8] who gave theoretial empirial arguments for modeling seismi waveforms using low-order autoregressive models.the preision of s k (t is dened as, whih isused later in our appliation of the Gibbs sampler. We use the same set of autoregressive oeients, [ ::: p ] 0, to model the k th underlying signal-path eets sequene s k (t, sine the same event is assumed to generate the trae y k (t reorded at the k th seismometer. Hene, we believe the same autoregressive model will be appliable for eah underlying signal-path eet sequene s k (t, k ::: q, t ::: n, although the funtions generated will be dierent. Finally, we assume that the additive observation noise term " k (t in ( is distributed independent identially, over t k, normal with mean 0 variane, where >0 is xed. The parameter an be thought ofastheinverse of the signalto-noise ratio, dened as SNR e ", where e is the variane of the underlying signal-path eets sequene on eah hannel " is the variane of the noise term. Reall that e ", whih implies that SNR. The use of the xed parameter alleviates a saling problem that exists if two separate variane terms are inluded in our model. 3 The Bayesian Approah In this paper the deonvolution of seismi array data for ommon delay patterns is ahieved using Bayesian methods. Speially in the model under onsideration in (, the parameter set is f a Sg: (7 The set ontains +m+p++nq elements, where we dene a [a a ::: a m ] 0, [ ::: p ] 0, S [s s ::: s q ] 0, suh that eah row of S is s k [s k ( s k ( ::: s k (n] 0. The given data onsists of the olletion of q seismi traes reorded for a duration of n points. Here we dene Y [y y ::: y q ] 0, suh that eah row of Y is y k [y k ( y k ( ::: y k (n] 0. Hene in this problem only qn data points are available to estimate all of the parameters in. To develop a Bayesian deonvolution method we follow the Bayesian approahto statistial data analysis. We develop a model likelihood p(yj. Prior distributions for the parameters are hosen. The overall prior is p(. Lastly, the posterior distributions of the parameters in, given the data Y, is alulated. By using Bayes' Rule the posterior distributions an be alulated as p(jy / p(yjp(. Then, to draw inferenes about the unknown parameters in our model, we alulate the marginal posteriors. For point estimates we alulate the posterior means. In our deonvolution problem, is a vetor of many parameters, due to our modeling assumptions the marginal posterior distributions are diult to alulate. However, the onditional marginal distributions are available. From the latter we an implement the Gibbs sampler, a Markov Chain Monte Carlo tehnique that simulates \rom samples" from the onditional marginal distributions, whih an be used to alulate posterior estimates of the parameters in our model. Point estimates are alulated by taking the means of the Markov hains samples after a suient \burn in" period. See Gelf Smith [3] Tanner [7] for desriptions of the Gibbs sampler. 3. Prior Distributions To perform our Bayesian analysis we hoose prior distributions for eah unknown parameter in the parameter set. We speify priors on,, sine they are unknown. We have modeling assumptions that determine the distributions of the rom parameters a S. The probability of seeing a zero at any point in the pulse sequene a j, j ::: m, is assumed to have a beta prior, with predetermined hyperparameters > 0 > 0. For the AR(p oeients [ ::: p ] 0, we assume a p-variate normal prior distribution where 0 0 are predetermined hyperparameters. To speify the prior for the ommon preision of the underlying signal-path eets, for eah hannel k, we assume that has a gamma prior where > 0 > 0 are predetermined hyperparameters.
4 3. Overall Prior, Likelihood, Conditional Posterior Distributions The overall prior distribution on the parameter set is dened using independene assumptions is p( p(p(p( my j p(a jj qy k p(s k j :(8 We speify the onditional likelihood use it as an approximation of the full likelihood. Due to the dependene of y k (t on the past m values of s k (t, whih an be seen in (, the dependene of s k (t on its past p values, whih an be seen in (6, we dene the likelihood onditional on y k ( ::: y k (l, k ::: q, where l m + p. We now redene the data matrix as Y [y y ::: y q ], where eah y k, k ::: q, starts at l +(m + p +, so y k [y k (l + ::: y k (n] 0. Using the model assumption that the noise, " k (t, is independent identially distributed normal with mean 0 variane, the onditional likelihood is p(yj ( ;q(n;l ( q(n;l ; " k(t P m where " k (t y k (t ; j0 a js k (t ; j. The joint density of the unknown parameters the data Y is p(y p(yjp(, where the prior is given in (8 the likelihood is given in (9. The joint posterior of the unknown parameters, given the data, an be ressed using Bayes' Rule as proportional to the likelihood times the prior, here p(jy / p(yjp(. For our deonvolution, alulations of p(jy, p(ajy, p(jy, p( jy, p(sjy are not easily alulated analytially. So we now proeed to develop the onditional marginal posterior distributions of the parameters in our model, whih are used in the Gibbs sampler. It is through the Gibbs sampler that we produe estimates of the parameters in our model. We proeed to develop the onditional marginal posterior distributions of the sets of parameters in our model. The following onditional marginal posteriors distributions are alulated: p(jy rest, p(ajy rest, p(jy rest, p( jy rest, p(sjy rest, where rest refers to the remaining parameters in for a spei parameter or set of parameters of interest.. The onditional posterior of is (9 jy rest beta( (0 with parameters n + m ; n +.. The onditional posterior of a j, j ::: m is p(a jjy rest / ji(aj 0 +( ; j (aj a j p aj ; where j [ +(; d j ], d j (aj a j ( " aj a j + (a j ; aj a j I(a j > 0 ( aj aj a j ; " k(t[;j]s k (t ; j P where " k (t[;j] y m k(t ; i0 i6j a is k (t ; i, a j " + s k(t ; j : 3. The onditional posterior of is jy rest N p ( ( where ~s k (t [s k (t ; ::: s k (t ; p] 0, "; s k (t~s k (t ~s k (t~s 0 k(t + 0 : 4. The onditional posterior of is jy rest gamma( (3 where q(n ; l ; p+ " k(t+ k tp+ e k(t+ : 5. For k :: q i (p + ::: n, s k (i has a Normal posterior distribution with s k (ijy rest N( sk (i s k (i (4
5 with sk (i s k 4 min((i+m n (i " 0 k(t[;i]a t;i tmax(i (l+ min((i+p n ; e 0 k(t[;i] t;i where " 0 k(t[;i] y k (t ; e 0 k(t[;i] s k (i ti m j0 j6t;i p j0 j6t;i 4 min((i+m n tmax(i (l+ + a j s k (t ; j j s k (t ; j a t;i min((i+p n ti t;i 6. For k :: q the onditional posterior of ~s k (p +[s k (p ::: s k (] 0 : We dene t [ t ::: p 0 ::: 0] 0 for t ::: p also dene k (t s k (p + t + s k (p + t ; + ::: + ts k (p + for t f ::: pg. Therefore, we have ~s k (p + N p ( ~sk (p+ ~s 0 k (p+ (5 where " p ~sk (p+ ~sk (p+ ; k (t t " ~s 0 k (p+ p t t t 0 t + 4 Results for Real Data We present in this setion the results produed by applying our Bayesian deonvolution tehnique to real data. The array data we analyze is Event 054, a mine blast reorded at the Arti Experimental Seismi Station (ARCESS in northern Norway, previously analyzed by Shumway, Baumgart, Der [5]. The seismi soure of this data is known to be a ripple-red mining losion that was reorded at a regional distane with a sampling rate of 40 points per seond. : : 4. Hyperparameters Fixed Model Parameters In this setion we rst present our hoie of hyperparameters for the prior distributions on the parameters in our hoie of the xed model parameters. For, thehyperparameters we use in our beta prior are For the autoregressive oeients themeanvalue we use for our multivariate normal prior on is 0 [;0:9 0:7 ;0:] 0 the ovariane matrix is 0 4 (0:45 ;(0:5 (0:5 ;(0:5 (0:30 ;(0:5 (0:5 ;(0:5 (0:5 3 5 : (6 The hyperparameters we hoose for the gamma prior distribution of the preision of the underlying signal-path eets sequene s k on eah hannel k are 5 4. The hyperparameters that determine the normal distribution used to model the independent nonzero elements of the pulse sequene a j are 0:7 0:5. The value of the parameter is hosen to reet the signal-to-noise ratio dened as SNR e ". We x 0:0. And we xm Results for Event 054 Here we present the results produed by applying our Bayesian deonvolution program to the P n -phase of Event 054, see Figure for plots of the q 5 reordings of the P n -phase. Estimates of the model parameters are alulated from a typial run of our Gibbs sampling program. For the results presented here the program was run 40,000 iterations a \burn in" of 35,000 iterations was used. The estimates were omputed by alulating the means of the Markov hains after \burn in" for eah parameter. The peak values of a estimated from the P n -phase our at the times 4, 8,, 4, 8, 5, 8, 35, whih an be seen in Figure. The estimated delay struture ontains multiple delays, indiating that the P n -phase of Event 054was generated by ripple- ring. 5 Conlusions To summarize our work, we have developed a Bayesisan deonvolution method for seismi array data that detets a ommon delay pattern. This work was based on the univariate work of Rong Chen Cheng, Chen, Li []. The priors we hose were derived from the previous work
6 Pn054 Channel One - Observed Trae Pn054 Channel Two - Observed Trae Pn054 Channel There - Observed Trae the statistial analysis ert opinion related to the parameters in our model or researh related to the seismi ativity of the region being monitored. Seondly, we see the possibility of automating the analysis so a quik rst analysis an be run to hek seismi events for ommon ripple-red mining losions. 6 Aknowledgments This work was supported by the Center for Statistis in Siene Tehnology, U.C. Davis, under grant DMS Pn054 Channel Four - Observed Trae Pn054 Channel Five - Observed Trae Figure : P n -phase Event 054: Saled data, y k, k ::: q, whereq 5 n 0. aest Pn054 - Delay Amplitude Estimates Figure : P n -phase Event 054: Posterior estimates of the pulse sequene, ^a j, j ::: m, where m 40. of Tjstheim [8], Smith [6], Shumway, Baumgart, Der [5]. We tested our program extensively on simulated data we used the lessons learned from the simulated data analysis to implement our program on the data from Event 054 reorded at ARCESS. We foresee two main impliations of our method for monitoring seismi ativity for low level nulear losions. First, our method an inorporate into j Referenes [] Cheng, Q., Chen, R., IL, T.-H. (996 Simultaneous Wavelet Estimation Deonvolution of Reetion Seismi Signals IEEE Transations On Geosiene Remote Sensing 34, [] Dargah-Noubary, G.R. (995 Stohasti modeling identiation of seismi reords based on established deterministi formulations. J. Series Analysis 6, 0-9. [3] Gelf, A.E. Smith, A.F.M. (990 Sampling-Based Approahes to Calulating Marginal Densities Journal of the Amerian Statistial Assoiation 85, [4] Mendel, J.M. (990 Maximum-Likelihood Deonvolution: A Journey into Model-Based Signal Proessing. New York: Springer-Verlag. [5] Shumway, R.H., Baumgart, D.R. Der, Z.A. (998 A Cepstral F Statisti for Deteting Delay-Fired Seismi Signals. Tehnometris 40, [6] Smith, A.T. (993 Disrimination of Explosions from Simultaneous Mining Blasts Bulletin of the Seismologial Soiety of Ameria 83, [7] Tanner, M.A. (996 Tools for Statistial Inferene, Methods for the Exploration of Posterior Distributions Likelihood Funtions, 3rd edition. New York: Springer-Verlag. [8] Tjstheim, D. (975 Autoregressive Representation of Seismi P-wave Signals with an Appliation to the Problem of Short-Period Disriminants Geophys. J. R. Astr. So. 43, 69-9.
Complexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationExploring the feasibility of on-site earthquake early warning using close-in records of the 2007 Noto Hanto earthquake
Exploring the feasibility of on-site earthquake early warning using lose-in reords of the 2007 Noto Hanto earthquake Yih-Min Wu 1 and Hiroo Kanamori 2 1. Department of Geosienes, National Taiwan University,
More informationOptimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach
Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)
More informationDanielle Maddix AA238 Final Project December 9, 2016
Struture and Parameter Learning in Bayesian Networks with Appliations to Prediting Breast Caner Tumor Malignany in a Lower Dimension Feature Spae Danielle Maddix AA238 Final Projet Deember 9, 2016 Abstrat
More informationA Spatiotemporal Approach to Passive Sound Source Localization
A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,
More informationModel-based mixture discriminant analysis an experimental study
Model-based mixture disriminant analysis an experimental study Zohar Halbe and Mayer Aladjem Department of Eletrial and Computer Engineering, Ben-Gurion University of the Negev P.O.Box 653, Beer-Sheva,
More informationMeasuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach
Measuring & Induing Neural Ativity Using Extraellular Fields I: Inverse systems approah Keith Dillon Department of Eletrial and Computer Engineering University of California San Diego 9500 Gilman Dr. La
More informationProbabilistic Graphical Models
Probabilisti Graphial Models David Sontag New York University Leture 12, April 19, 2012 Aknowledgement: Partially based on slides by Eri Xing at CMU and Andrew MCallum at UMass Amherst David Sontag (NYU)
More informationNonreversibility of Multiple Unicast Networks
Nonreversibility of Multiple Uniast Networks Randall Dougherty and Kenneth Zeger September 27, 2005 Abstrat We prove that for any finite direted ayli network, there exists a orresponding multiple uniast
More informationSection 3. Interstellar absorption lines. 3.1 Equivalent width
Setion 3 Interstellar absorption lines 3.1 Equivalent width We an study diuse interstellar louds through the absorption lines they produe in the spetra of bakground stars. Beause of the low temperatures
More informationCase I: 2 users In case of 2 users, the probability of error for user 1 was earlier derived to be 2 A1
MUTLIUSER DETECTION (Letures 9 and 0) 6:33:546 Wireless Communiations Tehnologies Instrutor: Dr. Narayan Mandayam Summary By Shweta Shrivastava (shwetash@winlab.rutgers.edu) bstrat This artile ontinues
More informationComputer Science 786S - Statistical Methods in Natural Language Processing and Data Analysis Page 1
Computer Siene 786S - Statistial Methods in Natural Language Proessing and Data Analysis Page 1 Hypothesis Testing A statistial hypothesis is a statement about the nature of the distribution of a random
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationChapter 8 Hypothesis Testing
Leture 5 for BST 63: Statistial Theory II Kui Zhang, Spring Chapter 8 Hypothesis Testing Setion 8 Introdution Definition 8 A hypothesis is a statement about a population parameter Definition 8 The two
More informationEstimation for Unknown Parameters of the Extended Burr Type-XII Distribution Based on Type-I Hybrid Progressive Censoring Scheme
Estimation for Unknown Parameters of the Extended Burr Type-XII Distribution Based on Type-I Hybrid Progressive Censoring Sheme M.M. Amein Department of Mathematis Statistis, Faulty of Siene, Taif University,
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationProbabilistic and nondeterministic aspects of Anonymity 1
MFPS XX1 Preliminary Version Probabilisti and nondeterministi aspets of Anonymity 1 Catusia Palamidessi 2 INRIA and LIX Éole Polytehnique, Rue de Salay, 91128 Palaiseau Cedex, FRANCE Abstrat Anonymity
More informationSensor management for PRF selection in the track-before-detect context
Sensor management for PRF seletion in the tra-before-detet ontext Fotios Katsilieris, Yvo Boers, and Hans Driessen Thales Nederland B.V. Haasbergerstraat 49, 7554 PA Hengelo, the Netherlands Email: {Fotios.Katsilieris,
More informationSingular Event Detection
Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate
More informationGeweke, J. (1989). Bayesian inference in econometric models using Monte Carlo integration.
Geweke, J (1989) Bayesian inferene in eonometri models using Monte Carlo integration Eonometria, 57, 1317{1339 Geyer, C J (1992) On the onvergene of Monte Carlo maximum likelihood alulations Tehnial Report
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationSensitivity Analysis in Markov Networks
Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationAnalysis of discretization in the direct simulation Monte Carlo
PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,
More information10.5 Unsupervised Bayesian Learning
The Bayes Classifier Maximum-likelihood methods: Li Yu Hongda Mao Joan Wang parameter vetor is a fixed but unknown value Bayes methods: parameter vetor is a random variable with known prior distribution
More informationAn I-Vector Backend for Speaker Verification
An I-Vetor Bakend for Speaker Verifiation Patrik Kenny, 1 Themos Stafylakis, 1 Jahangir Alam, 1 and Marel Kokmann 2 1 CRIM, Canada, {patrik.kenny, themos.stafylakis, jahangir.alam}@rim.a 2 VoieTrust, Canada,
More informationarxiv: v2 [math.pr] 9 Dec 2016
Omnithermal Perfet Simulation for Multi-server Queues Stephen B. Connor 3th Deember 206 arxiv:60.0602v2 [math.pr] 9 De 206 Abstrat A number of perfet simulation algorithms for multi-server First Come First
More informationAn Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances
An aptive Optimization Approah to Ative Canellation of Repeated Transient Vibration Disturbanes David L. Bowen RH Lyon Corp / Aenteh, 33 Moulton St., Cambridge, MA 138, U.S.A., owen@lyonorp.om J. Gregory
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More information2 The Bayesian Perspective of Distributions Viewed as Information
A PRIMER ON BAYESIAN INFERENCE For the next few assignments, we are going to fous on the Bayesian way of thinking and learn how a Bayesian approahes the problem of statistial modeling and inferene. The
More informationArray Design for Superresolution Direction-Finding Algorithms
Array Design for Superresolution Diretion-Finding Algorithms Naushad Hussein Dowlut BEng, ACGI, AMIEE Athanassios Manikas PhD, DIC, AMIEE, MIEEE Department of Eletrial Eletroni Engineering Imperial College
More informationA Queueing Model for Call Blending in Call Centers
A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl
More informationOn the Bit Error Probability of Noisy Channel Networks With Intermediate Node Encoding I. INTRODUCTION
5188 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 54, NO. 11, NOVEMBER 2008 [8] A. P. Dempster, N. M. Laird, and D. B. Rubin, Maximum likelihood estimation from inomplete data via the EM algorithm, J.
More informationMultiPhysics Analysis of Trapped Field in Multi-Layer YBCO Plates
Exerpt from the Proeedings of the COMSOL Conferene 9 Boston MultiPhysis Analysis of Trapped Field in Multi-Layer YBCO Plates Philippe. Masson Advaned Magnet Lab *7 Main Street, Bldg. #4, Palm Bay, Fl-95,
More informationQCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines
DOI.56/sensoren6/P3. QLAS Sensor for Purity Monitoring in Medial Gas Supply Lines Henrik Zimmermann, Mathias Wiese, Alessandro Ragnoni neoplas ontrol GmbH, Walther-Rathenau-Str. 49a, 7489 Greifswald, Germany
More informationDIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS
CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install
More information15.12 Applications of Suffix Trees
248 Algorithms in Bioinformatis II, SoSe 07, ZBIT, D. Huson, May 14, 2007 15.12 Appliations of Suffix Trees 1. Searhing for exat patterns 2. Minimal unique substrings 3. Maximum unique mathes 4. Maximum
More informationOn the Designs and Challenges of Practical Binary Dirty Paper Coding
On the Designs and Challenges of Pratial Binary Dirty Paper Coding 04 / 08 / 2009 Gyu Bum Kyung and Chih-Chun Wang Center for Wireless Systems and Appliations Shool of Eletrial and Computer Eng. Outline
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationMethods of evaluating tests
Methods of evaluating tests Let X,, 1 Xn be i.i.d. Bernoulli( p ). Then 5 j= 1 j ( 5, ) T = X Binomial p. We test 1 H : p vs. 1 1 H : p>. We saw that a LRT is 1 if t k* φ ( x ) =. otherwise (t is the observed
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationMOMENT TENSOR AND SOURCE PROCESS OF EARTHQUAKES IN FIJI REGION OBTAINED BY WAVEFORM INVERSION
MOMENT TENSOR AND SOURCE PROCESS OF EARTHQUAKES IN FIJI REGION OBTAINED BY WAVEFORM INVERSION Seru Sefanaia* Supervisor: Yui YAGI** MEE0767 ABSTRACT We evaluated the uality of the loal seismi waveform
More informationAssessing the Performance of a BCI: A Task-Oriented Approach
Assessing the Performane of a BCI: A Task-Oriented Approah B. Dal Seno, L. Mainardi 2, M. Matteui Department of Eletronis and Information, IIT-Unit, Politenio di Milano, Italy 2 Department of Bioengineering,
More informationSupplementary Material
Supplementary Material Joint Disriminative Bayesian Ditionary and Classifier Learning Joint probability distribution Aording to the proposed model, the joint probability distribution over the data of the
More informationSome recent developments in probability distributions
Proeedings 59th ISI World Statistis Congress, 25-30 August 2013, Hong Kong (Session STS084) p.2873 Some reent developments in probability distributions Felix Famoye *1, Carl Lee 1, and Ayman Alzaatreh
More informationA new method of measuring similarity between two neutrosophic soft sets and its application in pattern recognition problems
Neutrosophi Sets and Systems, Vol. 8, 05 63 A new method of measuring similarity between two neutrosophi soft sets and its appliation in pattern reognition problems Anjan Mukherjee, Sadhan Sarkar, Department
More informationApplication of the Dyson-type boson mapping for low-lying electron excited states in molecules
Prog. Theor. Exp. Phys. 05, 063I0 ( pages DOI: 0.093/ptep/ptv068 Appliation of the Dyson-type boson mapping for low-lying eletron exited states in moleules adao Ohkido, and Makoto Takahashi Teaher-training
More informationSupplementary Materials
Supplementary Materials Neural population partitioning and a onurrent brain-mahine interfae for sequential motor funtion Maryam M. Shanehi, Rollin C. Hu, Marissa Powers, Gregory W. Wornell, Emery N. Brown
More informationThe universal model of error of active power measuring channel
7 th Symposium EKO TC 4 3 rd Symposium EKO TC 9 and 5 th WADC Workshop nstrumentation for the CT Era Sept. 8-2 Kosie Slovakia The universal model of error of ative power measuring hannel Boris Stogny Evgeny
More informationTemperature Control of Batch Suspension Polyvinyl Chloride Reactors
1285 A publiation of CHEMICAL ENGINEERING TRANSACTIONS VOL. 39, 2014 Guest Editors: Petar Sabev Varbanov, Jiří Jaromír Klemeš, Peng Yen Liew, Jun Yow Yong Copyright 2014, AIDIC Servizi S.r.l., ISBN 978-88-95608-30-3;
More informationContact Block Reduction Method for Ballistic Quantum Transport with Semi-empirical sp3d5s* Tight Binding band models
Purdue University Purdue e-pubs Other Nanotehnology Publiations Birk Nanotehnology Center -2-28 Contat Redution Method for Ballisti Quantum Transport with Semi-empirial sp3d5s* Tight Binding band models
More informationSHIELDING MATERIALS FOR HIGH-ENERGY NEUTRONS
SHELDNG MATERALS FOR HGH-ENERGY NEUTRONS Hsiao-Hua Hsu Health Physis Measurements Group Los Alamos National Laboratory Los Alamos, New Mexio, 87545 USA Abstrat We used the Monte Carlo transport ode Los
More information7 Max-Flow Problems. Business Computing and Operations Research 608
7 Max-Flow Problems Business Computing and Operations Researh 68 7. Max-Flow Problems In what follows, we onsider a somewhat modified problem onstellation Instead of osts of transmission, vetor now indiates
More informationRelativistic Addition of Velocities *
OpenStax-CNX module: m42540 1 Relativisti Addition of Veloities * OpenStax This work is produed by OpenStax-CNX and liensed under the Creative Commons Attribution Liense 3.0 Abstrat Calulate relativisti
More informationTHE TWIN PARADOX A RELATIVISTIC DOMAIN RESOLUTION
THE TWIN PARADOX A RELATIVISTIC DOMAIN RESOLUTION Peter G.Bass P.G.Bass www.relativitydomains.om January 0 ABSTRACT This short paper shows that the so alled "Twin Paradox" of Speial Relativity, is in fat
More informationCommon Trends in European School Populations
Common rends in European Shool Populations P. Sebastiani 1 (1) and M. Ramoni (2) (1) Department of Mathematis and Statistis, University of Massahusetts. (2) Children's Hospital Informatis Program, Harvard
More informationSINCE Zadeh s compositional rule of fuzzy inference
IEEE TRANSACTIONS ON FUZZY SYSTEMS, VOL. 14, NO. 6, DECEMBER 2006 709 Error Estimation of Perturbations Under CRI Guosheng Cheng Yuxi Fu Abstrat The analysis of stability robustness of fuzzy reasoning
More informationA Heuristic Approach for Design and Calculation of Pressure Distribution over Naca 4 Digit Airfoil
IOSR Journal of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 PP 11-15 www.iosrjen.org A Heuristi Approah for Design and Calulation of Pressure Distribution over Naa 4 Digit Airfoil G.
More informationAyan Kumar Bandyopadhyay
Charaterization of radiating apertures using Multiple Multipole Method And Modeling and Optimization of a Spiral Antenna for Ground Penetrating Radar Appliations Ayan Kumar Bandyopadhyay FET-IESK, Otto-von-Guerike-University,
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 4, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume, No 4, 01 Copyright 010 All rights reserved Integrated Publishing servies Researh artile ISSN 0976 4399 Strutural Modelling of Stability
More informationDevelopment of Fuzzy Extreme Value Theory. Populations
Applied Mathematial Sienes, Vol. 6, 0, no. 7, 58 5834 Development of Fuzzy Extreme Value Theory Control Charts Using α -uts for Sewed Populations Rungsarit Intaramo Department of Mathematis, Faulty of
More informationData-Recursive Algorithms for Blind Channel Identication in. Direct-Sequence Systems. November 30, Abstract
Data-Reursive Algorithms for Blind Channel Identiation in Diret-Sequene Systems Dennis L. Goekel 1, Alfred O. Hero III 2, and Wayne E. Stark 3 November 30, 1995 Abstrat Algorithms for performing blind
More information7.1 Roots of a Polynomial
7.1 Roots of a Polynomial A. Purpose Given the oeffiients a i of a polynomial of degree n = NDEG > 0, a 1 z n + a 2 z n 1 +... + a n z + a n+1 with a 1 0, this subroutine omputes the NDEG roots of the
More informationA GENERATION METHOD OF SIMULATED EARTHQUAKE GROUND MOTION CONSIDERING PHASE DIFFERENCE CHARACTERISTICS
Otober 1-17, 8, Beijing, China A GENERATION METHOD OF SIMULATED EARTHQUAKE GROUND MOTION CONSIDERING PHASE DIFFERENCE CHARACTERISTICS T. Yamane 1 and S. Nagahashi 1 Senior Strutural Engineer, Strutural
More informationDept. of Computer Science. Raleigh, NC 27695, USA. May 14, Abstract. 1, u 2 q i+1 :
Anti-Leture Hall Compositions Sylvie Corteel CNRS PRiSM, UVSQ 45 Avenue des Etats-Unis 78035 Versailles, Frane syl@prism.uvsq.fr Carla D. Savage Dept. of Computer Siene N. C. State University, Box 8206
More informationTests of fit for symmetric variance gamma distributions
Tests of fit for symmetri variane gamma distributions Fragiadakis Kostas UADPhilEon, National and Kapodistrian University of Athens, 4 Euripidou Street, 05 59 Athens, Greee. Keywords: Variane Gamma Distribution,
More informationRESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED
The Seventh Asia-Paifi Conferene on Wind Engineering, November 8-1, 9, Taipei, Taiwan RESEARCH ON RANDOM FORIER WAVE-NMBER SPECTRM OF FLCTATING WIND SPEED Qi Yan 1, Jie Li 1 Ph D. andidate, Department
More informationVariation Based Online Travel Time Prediction Using Clustered Neural Networks
Variation Based Online Travel Time Predition Using lustered Neural Networks Jie Yu, Gang-Len hang, H.W. Ho and Yue Liu Abstrat-This paper proposes a variation-based online travel time predition approah
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 16 Aug 2004
Computational omplexity and fundamental limitations to fermioni quantum Monte Carlo simulations arxiv:ond-mat/0408370v1 [ond-mat.stat-meh] 16 Aug 2004 Matthias Troyer, 1 Uwe-Jens Wiese 2 1 Theoretishe
More informationWeighted K-Nearest Neighbor Revisited
Weighted -Nearest Neighbor Revisited M. Biego University of Verona Verona, Italy Email: manuele.biego@univr.it M. Loog Delft University of Tehnology Delft, The Netherlands Email: m.loog@tudelft.nl Abstrat
More informationTextual Document Indexing and Retrieval via Knowledge Sources and Data Mining
Textual Doument Indexing and Retrieval via Knowledge Soures and Data Mining Wesley. W. Chu, Zhenyu Liu and Wenlei Mao Computer Siene Department, University of California, Los Angeles 90095 {ww, viliu,
More informationA Characterization of Wavelet Convergence in Sobolev Spaces
A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat
More informationTaste for variety and optimum product diversity in an open economy
Taste for variety and optimum produt diversity in an open eonomy Javier Coto-Martínez City University Paul Levine University of Surrey Otober 0, 005 María D.C. Garía-Alonso University of Kent Abstrat We
More informationAharonov-Bohm effect. Dan Solomon.
Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that
More informationA NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS
Vietnam Journal of Mehanis, VAST, Vol. 4, No. (), pp. A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Le Thanh Tung Hanoi University of Siene and Tehnology, Vietnam Abstrat. Conventional ship autopilots are
More informationSufficient Conditions for a Flexible Manufacturing System to be Deadlocked
Paper 0, INT 0 Suffiient Conditions for a Flexile Manufaturing System to e Deadloked Paul E Deering, PhD Department of Engineering Tehnology and Management Ohio University deering@ohioedu Astrat In reent
More informationLightpath routing for maximum reliability in optical mesh networks
Vol. 7, No. 5 / May 2008 / JOURNAL OF OPTICAL NETWORKING 449 Lightpath routing for maximum reliability in optial mesh networks Shengli Yuan, 1, * Saket Varma, 2 and Jason P. Jue 2 1 Department of Computer
More informationFailure Assessment Diagram Analysis of Creep Crack Initiation in 316H Stainless Steel
Failure Assessment Diagram Analysis of Creep Crak Initiation in 316H Stainless Steel C. M. Davies *, N. P. O Dowd, D. W. Dean, K. M. Nikbin, R. A. Ainsworth Department of Mehanial Engineering, Imperial
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More informationConformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices
AAS 03-190 Conformal Mapping among Orthogonal, Symmetri, and Skew-Symmetri Matries Daniele Mortari Department of Aerospae Engineering, Texas A&M University, College Station, TX 77843-3141 Abstrat This
More information3.2 Gaussian (Normal) Random Numbers and Vectors
3.2 Gaussian (Normal) Random Numbers and Vetors A. Purpose Generate pseudorandom numbers or vetors from the Gaussian (normal) distribution. B. Usage B.1 Generating Gaussian (normal) pseudorandom numbers
More informationA NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM
NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationEffect of magnetization process on levitation force between a superconducting. disk and a permanent magnet
Effet of magnetization proess on levitation fore between a superonduting disk and a permanent magnet L. Liu, Y. Hou, C.Y. He, Z.X. Gao Department of Physis, State Key Laboratory for Artifiial Mirostruture
More informationFrequency hopping does not increase anti-jamming resilience of wireless channels
Frequeny hopping does not inrease anti-jamming resiliene of wireless hannels Moritz Wiese and Panos Papadimitratos Networed Systems Seurity Group KTH Royal Institute of Tehnology, Stoholm, Sweden {moritzw,
More informationarxiv:cond-mat/ v1 [cond-mat.str-el] 3 Aug 2006
arxiv:ond-mat/0608083v1 [ond-mat.str-el] 3 Aug 006 Raman sattering for triangular latties spin- 1 Heisenberg antiferromagnets 1. Introdution F. Vernay 1,, T. P. Devereaux 1, and M. J. P. Gingras 1,3 1
More informationParallel disrete-event simulation is an attempt to speed-up the simulation proess through the use of multiple proessors. In some sense parallel disret
Exploiting intra-objet dependenies in parallel simulation Franeso Quaglia a;1 Roberto Baldoni a;2 a Dipartimento di Informatia e Sistemistia Universita \La Sapienza" Via Salaria 113, 198 Roma, Italy Abstrat
More informationCapacity-achieving Input Covariance for Correlated Multi-Antenna Channels
Capaity-ahieving Input Covariane for Correlated Multi-Antenna Channels Antonia M. Tulino Universita Degli Studi di Napoli Federio II 85 Napoli, Italy atulino@ee.prineton.edu Angel Lozano Bell Labs (Luent
More informationThe Laws of Acceleration
The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the
More informationProduct Policy in Markets with Word-of-Mouth Communication. Technical Appendix
rodut oliy in Markets with Word-of-Mouth Communiation Tehnial Appendix August 05 Miro-Model for Inreasing Awareness In the paper, we make the assumption that awareness is inreasing in ustomer type. I.e.,
More informationOptimal Control of Air Pollution
Punjab University Journal of Mathematis (ISSN 1016-2526) Vol. 49(1)(2017) pp. 139-148 Optimal Control of Air Pollution Y. O. Aderinto and O. M. Bamigbola Mathematis Department, University of Ilorin, Ilorin,
More informationEstimating the probability law of the codelength as a function of the approximation error in image compression
Estimating the probability law of the odelength as a funtion of the approximation error in image ompression François Malgouyres Marh 7, 2007 Abstrat After some reolletions on ompression of images using
More informationThe experimental plan of displacement- and frequency-noise free laser interferometer
7th Edoardo Amaldi Conferene on Gravitational Waves (Amaldi7) Journal of Physis: Conferene Series 122 (2008) 012022 The experimental plan of displaement- and frequeny-noise free laser interferometer K
More informationCentroid Detection by Gaussian Pattern Matching In Adaptive Optics
010 International Journal of Computer Appliations (0975-8887) Centroid Detetion by Gaussian Pattern Mathing In Adaptive Optis Akondi Vyas Indian Institute of Astrophysis Indian Institute of Siene M B Roopashree
More informationGeneralized Neutrosophic Soft Set
International Journal of Computer Siene, Engineering and Information Tehnology (IJCSEIT), Vol.3, No.2,April2013 Generalized Neutrosophi Soft Set Said Broumi Faulty of Arts and Humanities, Hay El Baraka
More informationContact State Estimation using Multiple Model Estimation and Hidden Markov Models
Thomas J. Debus Pierre E. Dupont Department of Aerospae and Mehanial Engineering Boston University Boston, MA 02445, USA Robert D. Howe Division of Engineering and Applied Sienes Harvard University Cambridge,
More informationModeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers
A. E. Romanov et al.: Threading Disloation Density Redution in Layers (II) 33 phys. stat. sol. (b) 99, 33 (997) Subjet lassifiation: 6.72.C; 68.55.Ln; S5.; S5.2; S7.; S7.2 Modeling of Threading Disloation
More information