3.2 Gaussian (Normal) Random Numbers and Vectors
|
|
- Angelica Casey
- 6 years ago
- Views:
Transcription
1 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 B.1.a REAL B.1.b Program Prototype, Single Preision SRANG, X X = SRANG() Argument Definitions SRANG [out] The funtion returns a pseudorandom number from the Gaussian (normal) distribution with mean zero and unit standard deviation. B.2 Generating Gaussian (normal) pseudorandom vetors Given an N-vetor, µ, and a symmetri positive-definite N N matrix, A, the objetive is to ompute pseudorandom N-vetors, x, from the N-dimensional Gaussian (normal) distribution having mean vetor, µ, and ovariane matrix, A. On the first all to SRANGV with a new ovariane matrix, A, the user must set HAVEC =.false. to indiate that the Cholesky fator of A has not yet been omputed. SRANGV will replae A in storage by its Cholesky fator, C, and set HAVEC =.true. to indiate that the array A() now ontains C rather than A. On eah all to SRANGV an N-vetor, x, is returned as a pseudorandom sample from the speified distribution. B.2.a INTEGER REAL LOGICAL Program Prototype, Single Preision NDIM, N, IERR A(NDIM N, N), U( N), X( N) HAVEC Assign values to A(,), NDIM, N, U(), and HAVEC. CALL SRANGV(A, NDIM, N, U, X, HAVEC, IERR) Computed values will be returned in X() and IERR, and the ontents of A(,) and HAVEC may be hanged Calif. Inst. of Tehnology, 2015 Math à la Carte, In. B.2.b Argument Definitions A(,) [inout] When HAVEC is.false., A(,) ontains a user-supplied N N symmetri positive-definite ovariane matrix, A. Only the diagonal and subdiagonal elements of A need be given. This subroutine will not aess the super-diagonal loations in the array, A(,). This subroutine will replae A by its lower-triangular Cholesky fator, C, and set HAVEC =.true. When HAVEC is.true., A(,) is assumed to ontain the Cholesky fator, C. NDIM [in] First dimension of the array A(,). Require NDIM N. N [in] Order of the ovariane matrix A and dimension of the vetors U and X. Require N 1. U() [in] Contains the N-dimensional mean vetor, µ. X() [out] On return will ontain the N-dimensional generated random vetor. HAVEC [inout] See desription above for A(,). IERR [out] IERR will only be referened when the subroutine is entered with HAVEC =.false. If the Cholesky fatorization is suessful, the subroutine will set HAVEC =.true. and IERR = 0. Otherwise, it will leave HAVEC =.false. and set IERR to the index of the row of A in whih failure was noted. In this latter ase the results returned in A(,) and X() will not be useful. B.3 Modifiations for Double Preision For double-preision usage, hange the REAL statements to DOUBLE PRECISION and hange the initial S of the funtion and subroutine names to D. Note partiularly that if the funtion name, DRANG, is used it must be typed DOUBLE PRECISION either expliitly or via an IMPLICIT statement. C. Examples and Remarks The program DRSRANG demonstrates the use of SRANG to ompute Gaussian random numbers and uses SSTAT1 and SSTAT2 to ompute and print statistis and a histogram based on a sample of numbers delivered by SRANG. To ompute Gaussian random numbers with mean XMEAN and standard deviation STDDEV, one an use the statement X = XMEAN + STDDEV SRANG() July 11, 2015 Gaussian (Normal) Random Numbers and Vetors 3.2 1
2 The program DRSRANGV demonstrates the use of SRANGV to ompute pseudorandom vetors from a 3- dimensional Gaussian distribution with mean vetor and ovariane matrix speified as µ = and A = To feth or set the seed used in the underlying pseudorandom integer sequene use the subroutines desribed in Chapter 3.1. D. Funtional Desription Method The algorithm for generation of Gaussian random numbers is based on [1] [3]. This method draws pairs of uniform random numbers x from [0, 1] and y from [ 1, 1] until a pair is obtained that satisfies x 2 + y 2 1. The probability of satisfying this onstraint is π/ It then draws another uniform random number u from [0, 1] and omputes s = x 2 + y 2 t = 2 log u /s g 1 = t(x 2 y 2 ) g 2 = 2xyt The numbers g 1 and g 2 are independent samples from the Gaussian distribution with mean zero and unit standard deviation. The number g 1 is returned when it is omputed and g 2 is saved and returned the next time a Gaussian random number is requested. The saved value will be disarded if the underlying uniform sequene is reinitialized by a all to RAN1 or RANPUT of Chapter 3.1. This method is a mathematially exat transformation from the uniform distribution to the Gaussian distribution so the statistial quality of the delivered numbers depends entirely on the quality of the uniform pseudorandom numbers used. The uniform numbers are obtained by alling SRANUA or DRANUA, using the array in ommon blok /RANCMS/ or /RANCMD/ as a buffer as desribed in Chapter 3.1. For the generation of Gaussian random vetors we are given a mean vetor µ and a symmetri positive-definite ovariane matrix A. The Cholesky method is used to fator the given ovariane matrix A as A = CC t where C is a lower triangular matrix. Then for eah vetor to be delivered the method first onstruts a vetor g whose omponents are independent samples from the Gaussian distribution with mean zero and unit standard deviation. It then omputes x = µ + Cg. whih is a sample from the N-dimensional Gaussian distribution with mean µ and ovariane matrix A. Values returned as double-preision random numbers will have random bits throughout the word, however the quality of randomness should not be expeted to be as good in a low-order segment of the word as in a highorder part. Referenes 1. James R. Bell, Algorithm 334 Normal random deviates, Comm. ACM 11, 7 (July 1968) G. Box and M. Muller, Note on the generation of normal deviates, Ann. Math. Stat. 28 (1958) John von Neumann, Various tehniques used in onnetion with random digits in the National Bureau of Standards, Applied Mathematis Series 12, National Bureau of Standards (1959) 36 pages. E. Error Proedures and Restritions While omputing the Cholesky fatorization of A, subroutine SRANGV or DRANGV may determine that the matrix is not positive-definite. In that ase it will return with IERR set to the index of the row of the matrix at whih the problem was deteted. In this non-positivedefinite ase the ontents of the array A(,) will have been altered, HAVEC will still have the value.false., and X() will not ontain useful results. When the Cholesky fatorization is suessful, IERR will be set to zero. The onditions that N 1, and NDIM N are required but not heked. F. Supporting Information The soure language is ANSI Fortran 77. Entry Required Files DRANG DRANG, ERFIN, ERMSG, RANPK1, RANPK2 DRANGV DRANG, DRANGV, ERFIN, ERMSG, RANPK1, RANPK2 SRANG ERFIN, ERMSG, RANPK1, RANPK2, SRANG SRANGV ERFIN, ERMSG, RANPK1, RANPK2, SRANG, SRANGV Based on subprograms written for JPL by Carl Pitts, Heliodyne Corp., April, Adapted to Fortran 77 for the JPL MATH77 library by C. L. Lawson and S. Y. Chiu, JPL, April November 1991: Lawson reorganized and renamed ommon bloks. Improved oordination between Gaussian and uniform sequenes on reinitializations Gaussian (Normal) Random Numbers and Vetors July 11, 2015
3 DRSRANG program DRSRANG >> DRSRANG Krogh Minor hange f o r making. f90 v e r s i o n. >> DRSRANG Krogh Added e x t e r n a l statement. >> DRSRANG Krogh Changes t o use M77CON >> DRSRANG Lawson I n i t i a l Code. S r e p l a e s? : DR?RANG,?RANG,?STAT1,?STAT2 Demonstration d r i v e r f o r SRANG. Generates random numbers from t h e Gaussian d i s t r i b u t i o n. integer NCELLS parameter (NCELLS = 12+2) external SRANG real SRANG, STATS( 5 ), Y1, Y2, YTAB( 1 ), ZERO integer I, IHIST (NCELLS), N parameter (ZERO = 0. 0 E0) data N / 10000/ data Y1, Y2 / 3.0E0, +3.0E0/ STATS( 1 ) = ZERO do 20 I =1,N GET RANDOM NUMBER YTAB( 1 ) = SRANG( ) Aumulate s t a t i s t i s and histogram. a l l SSTAT1(YTAB( 1 ), 1, STATS, IHIST, NCELLS, Y1, Y2) 20 ontinue Print t h e s t a t i s t i s and histogram. write (, (13 x, a //) ) Gaussian random numbers from SRANG a l l SSTAT2(STATS, IHIST, NCELLS, Y1, Y2) stop end July 11, 2015 Gaussian (Normal) Random Numbers and Vetors 3.2 3
4 ODSRANG Gaussian random numbers from SRANG BREAK PT COUNT PLOT OF COUNT Count Minimum Maximum Mean Std. Deviation E Gaussian (Normal) Random Numbers and Vetors July 11, 2015
5 DRSRANGV program DRSRANGV >> DRSRANGV Krogh Minor hange f o r making. f90 v e r s i o n. >> DRSRANGV Krogh Changes t o use M77CON >> DRSRANGV Lawson I n i t i a l Code. S r e p l a e s? : DR?RANGV,?RANGV Demonstration d r i v e r f o r t h e Gaussian random v e t o r generator, SRANGV. C. L. Lawson & S. Y. Chiu, JPL, Apr integer I, J, N, IERR, NMAX, NSAMPL parameter (NMAX = 3) real A(NMAX,NMAX), U(NMAX), X(NMAX) l o g i a l HAVEC data (U( J ), J=1,NMAX)/ 1. 0 E0, 2. 0 E0, 3. 0 E0 / data (A( 1, J ), J=1,NMAX)/ E0, E0, E0 / data (A( 2, J ), J=1,NMAX)/ E0, E0, 0.03E0 / data (A( 3, J ), J=1,NMAX)/ E0, 0.03E0, E0 / N = NMAX NSAMPL = 20 HAVEC =. f a l s e. print (1 x, a/1x ), Pseudorandom v e t o r s omputed by SRANGV do 10 I = 1,NSAMPL a l l SRANGV(A, NMAX, N, U, X, HAVEC, IERR) i f (IERR. ne. 0) then print (1 x, i3, i 6 ),N, IERR else print (1 x, i3, 5 x, 3 f ), I,X endif 10 ontinue stop end July 11, 2015 Gaussian (Normal) Random Numbers and Vetors 3.2 5
6 ODSRANGV Pseudorandom v e t o r s omputed by SRANGV Gaussian (Normal) Random Numbers and Vetors July 11, 2015
11.1 Polynomial Least-Squares Curve Fit
11.1 Polynomial Least-Squares Curve Fit A. Purpose This subroutine determines a univariate polynomial that fits a given disrete set of data in the sense of minimizing the weighted sum of squares of residuals.
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 information4.3 Singular Value Decomposition and Analysis
4.3 Singular Value Deomposition and Analysis A. Purpose Any M N matrix, A, has a Singular Value Deomposition (SVD) of the form A = USV t where U is an M M orthogonal matrix, V is an N N orthogonal matrix,
More information2.9 Incomplete Elliptic Integrals
.9 Inomplete Ellipti Integrals A. Purpose An integral of the form R t, P t) 1/) dt 1) in whih Pt) is a polynomial of the third or fourth degree that has no multiple roots, and R is a rational funtion of
More information6.3 Basic Linear Algebra Subprograms (BLAS1)
6.3 Basi Linear Algebra Subprograms (BLAS1) A. Purpose This is a set of subroutine and funtion subprograms for basi mathematial operations on a single vetor or a pair of vetors. The operations provided
More informationComplexity 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 information8.2 Solve System of Nonlinear Equations
8.2 Solve System of Nonlinear Equations A. Purpose This subroutine solves a system of n nonlinear equations in n unknowns. The problem can be expressed as f(x) = 0 where x denotes an n-vector of unknowns
More informationmax min z i i=1 x j k s.t. j=1 x j j:i T j
AM 221: Advaned Optimization Spring 2016 Prof. Yaron Singer Leture 22 April 18th 1 Overview In this leture, we will study the pipage rounding tehnique whih is a deterministi rounding proedure that an be
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 informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationSOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS
SOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS Prepared by S. Broverman e-mail 2brove@rogers.om website http://members.rogers.om/2brove 1. We identify the following events:. - wathed gymnastis, ) - wathed baseball,
More informationQ2. [40 points] Bishop-Hill Model: Calculation of Taylor Factors for Multiple Slip
27-750, A.D. Rollett Due: 20 th Ot., 2011. Homework 5, Volume Frations, Single and Multiple Slip Crystal Plastiity Note the 2 extra redit questions (at the end). Q1. [40 points] Single Slip: Calulating
More informationA Unified View on Multi-class Support Vector Classification Supplement
Journal of Mahine Learning Researh??) Submitted 7/15; Published?/?? A Unified View on Multi-lass Support Vetor Classifiation Supplement Ürün Doğan Mirosoft Researh Tobias Glasmahers Institut für Neuroinformatik
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 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 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 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 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 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 informationAn iterative least-square method suitable for solving large sparse matrices
An iteratie least-square method suitable for soling large sparse matries By I. M. Khabaza The purpose of this paper is to report on the results of numerial experiments with an iteratie least-square method
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematis A NEW ARRANGEMENT INEQUALITY MOHAMMAD JAVAHERI University of Oregon Department of Mathematis Fenton Hall, Eugene, OR 97403. EMail: javaheri@uoregon.edu
More informationLecture 7: Sampling/Projections for Least-squares Approximation, Cont. 7 Sampling/Projections for Least-squares Approximation, Cont.
Stat60/CS94: Randomized Algorithms for Matries and Data Leture 7-09/5/013 Leture 7: Sampling/Projetions for Least-squares Approximation, Cont. Leturer: Mihael Mahoney Sribe: Mihael Mahoney Warning: these
More informationNAG Fortran Library Routine Document F04CFF.1
F04 Simultaneous Linear Equations NAG Fortran Library Routine Document Note: before using this routine, please read the Users Note for your implementation to check the interpretation of bold italicised
More informationProbability Distributions
STATGRAPHICS Re. 6805 Probability Distributions Summary... Data Input... Analysis Summary... 3 Analysis Options... 3 Cumulatie Distribution... 4 Inerse CDF... 5 Random Numbers... 6 DensityMass Funtion...
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 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 informationChapter 3 Church-Turing Thesis. CS 341: Foundations of CS II
CS 341: Foundations of CS II Marvin K. Nakayama Computer Siene Department New Jersey Institute of Tehnology Newark, NJ 07102 CS 341: Chapter 3 3-2 Contents Turing Mahines Turing-reognizable Turing-deidable
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 informationIntegration of the Finite Toda Lattice with Complex-Valued Initial Data
Integration of the Finite Toda Lattie with Complex-Valued Initial Data Aydin Huseynov* and Gusein Sh Guseinov** *Institute of Mathematis and Mehanis, Azerbaijan National Aademy of Sienes, AZ4 Baku, Azerbaijan
More informationQUANTUM MECHANICS II PHYS 517. Solutions to Problem Set # 1
QUANTUM MECHANICS II PHYS 57 Solutions to Problem Set #. The hamiltonian for a lassial harmoni osillator an be written in many different forms, suh as use ω = k/m H = p m + kx H = P + Q hω a. Find a anonial
More informationEE 321 Project Spring 2018
EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,
More informationOn maximal inequalities via comparison principle
Makasu Journal of Inequalities and Appliations (2015 2015:348 DOI 10.1186/s13660-015-0873-3 R E S E A R C H Open Aess On maximal inequalities via omparison priniple Cloud Makasu * * Correspondene: makasu@uw.a.za
More informationMODELING MATTER AT NANOSCALES. 4. Introduction to quantum treatments Eigenvectors and eigenvalues of a matrix
MODELING MATTER AT NANOSCALES 4 Introdution to quantum treatments 403 Eigenvetors and eigenvalues of a matrix Simultaneous equations in the variational method The problem of simultaneous equations in the
More informationThe transition between quasi-static and fully dynamic for interfaces
Physia D 198 (24) 136 147 The transition between quasi-stati and fully dynami for interfaes G. Caginalp, H. Merdan Department of Mathematis, University of Pittsburgh, Pittsburgh, PA 1526, USA Reeived 6
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 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 informationParticle-wave symmetry in Quantum Mechanics And Special Relativity Theory
Partile-wave symmetry in Quantum Mehanis And Speial Relativity Theory Author one: XiaoLin Li,Chongqing,China,hidebrain@hotmail.om Corresponding author: XiaoLin Li, Chongqing,China,hidebrain@hotmail.om
More informationA Small Footprint i-vector Extractor
A Small Footprint i-vetor Extrator Patrik Kenny Centre de reherhe informatique de Montréal (CRIM) Patrik.Kenny@rim.a Abstrat Both the memory and omputational requirements of algorithms traditionally used
More informationThe Tetrahedron Quality Factors of CSDS
MAX PLANCK INSTITUT FÜR AERONOMIE D 37191 Katlenburg-Lindau, Federal Republi of Germany MPAE W 100 94 27 The Tetrahedron Quality Fators of CSDS PATRICK W. DALY 1994 June 7 This report available from http://www.mpae.gwdg.de/~daly/tetra.html
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationA New Version of Flusser Moment Set for Pattern Feature Extraction
A New Version of Flusser Moment Set for Pattern Feature Extration Constantin-Iulian VIZITIU, Doru MUNTEANU, Cristian MOLDER Communiations and Eletroni Systems Department Military Tehnial Aademy George
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 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 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 informationestimated 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
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,
More informationETSI EN V8.0.1 ( )
EN 300 728 V8.0.1 (2000-11) European Standard (Teleommuniations series) Digital ellular teleommuniations system (Phase 2+); Comfort noise aspets for Enhaned Full Rate (EFR) speeh traffi hannels (GSM 06.62
More informationELG 5372 Error Control Coding. Claude D Amours Lecture 2: Introduction to Coding 2
ELG 5372 Error Control Coding Claude D Amours Leture 2: Introdution to Coding 2 Deoding Tehniques Hard Deision Reeiver detets data before deoding Soft Deision Reeiver quantizes reeived data and deoder
More informationCSC2515 Winter 2015 Introduc3on to Machine Learning. Lecture 5: Clustering, mixture models, and EM
CSC2515 Winter 2015 Introdu3on to Mahine Learning Leture 5: Clustering, mixture models, and EM All leture slides will be available as.pdf on the ourse website: http://www.s.toronto.edu/~urtasun/ourses/csc2515/
More informationTransformation to approximate independence for locally stationary Gaussian processes
ransformation to approximate independene for loally stationary Gaussian proesses Joseph Guinness, Mihael L. Stein We provide new approximations for the likelihood of a time series under the loally stationary
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 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 informationFormal Specification for Transportation Cyber Physical Systems
Formal Speifiation for Transportation Cyber Physial Systems ihen Zhang, Jifeng He and Wensheng Yu Shanghai Key aboratory of Trustworthy Computing East China Normal University Shanghai 200062, China Zhanglihen1962@163.om
More informationMATHEMATICAL AND NUMERICAL BASIS OF BINARY ALLOY SOLIDIFICATION MODELS WITH SUBSTITUTE THERMAL CAPACITY. PART II
Journal of Applied Mathematis and Computational Mehanis 2014, 13(2), 141-147 MATHEMATICA AND NUMERICA BAI OF BINARY AOY OIDIFICATION MODE WITH UBTITUTE THERMA CAPACITY. PART II Ewa Węgrzyn-krzypzak 1,
More informationFitting a Covariance Selection Model to a Matrix
Algorithm AS 105 Fitting a ovariance Selection Model to a Matrix By NANNY WERMUTH and EBERHARD SHEIDT Institutfur Medizinische Statistik und Dokumentation, Johannes-Gutenberg-Universitdt, Mainz, Federal
More informationCMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017
CMSC 451: Leture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017 Reading: Chapt 11 of KT and Set 54 of DPV Set Cover: An important lass of optimization problems involves overing a ertain domain,
More informationAn Integer Solution of Fractional Programming Problem
Gen. Math. Notes, Vol. 4, No., June 0, pp. -9 ISSN 9-784; Copyright ICSRS Publiation, 0 www.i-srs.org Available free online at http://www.geman.in An Integer Solution of Frational Programming Problem S.C.
More informationMath 151 Introduction to Eigenvectors
Math 151 Introdution to Eigenvetors The motivating example we used to desrie matrixes was landsape hange and vegetation suession. We hose the simple example of Bare Soil (B), eing replaed y Grasses (G)
More informationThe Lorenz Transform
The Lorenz Transform Flameno Chuk Keyser Part I The Einstein/Bergmann deriation of the Lorentz Transform I follow the deriation of the Lorentz Transform, following Peter S Bergmann in Introdution to the
More informationF04JGF NAG Fortran Library Routine Document
F4 Simultaneous Linear Equations F4JGF NAG Fortran Library Routine Document Note. Before using this routine, please read the Users Note for your implementation to check the interpretation of bold italicised
More informationChapter IX. Roots of Equations
14 Chapter IX Roots of Equations Root solving typially onsists of finding values of x whih satisfy a relationship suh as: 3 3 x + x = 1 + e x The problem an be restated as: Find the root of 3 f ( x) =
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 informationRelativity in Classical Physics
Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of
More informationLOGISTIC REGRESSION IN DEPRESSION CLASSIFICATION
LOGISIC REGRESSIO I DEPRESSIO CLASSIFICAIO J. Kual,. V. ran, M. Bareš KSE, FJFI, CVU v Praze PCP, CS, 3LF UK v Praze Abstrat Well nown logisti regression and the other binary response models an be used
More informationLCPFCT Flux-Corrected Transport Algorithm for Solving Generalized Continuity Equations
NAVAL RESEARCH LABORATORY Washington, DC 0375-530 NRL/MR/6410-93-719 LCPFCT Flux-Correted Transport Algorithm for Solving Generalized Continuity Equations JAY P. BORIS ALEXANDRA M. LANDSBERG ELAINE S.
More information(c) Calculate the tensile yield stress based on a critical resolved shear stress that we will (arbitrarily) set at 100 MPa. (c) MPa.
27-750, A.D. Rollett Due: 20 th Ot., 2011. Homework 5, Volume Frations, Single and Multiple Slip Crystal Plastiity Note the 2 extra redit questions (at the end). Q1. [40 points] Single Slip: Calulating
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 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 informationHamiltonian with z as the Independent Variable
Hamiltonian with z as the Independent Variable 1 Problem Kirk T. MDonald Joseph Henry Laboratories, Prineton University, Prineton, NJ 08544 Marh 19, 2011; updated June 19, 2015) Dedue the form of the Hamiltonian
More informationKAMILLA OLIVER AND HELMUT PRODINGER
THE CONTINUED RACTION EXPANSION O GAUSS HYPERGEOMETRIC UNCTION AND A NEW APPLICATION TO THE TANGENT UNCTION KAMILLA OLIVER AND HELMUT PRODINGER Abstrat Starting from a formula for tan(nx in the elebrated
More informationSome Math of the PAR Problem
Some Math of the PAR Problem Peter Oswald, Jaobs University Bremen http://www.faulty.jaobs-university.de/poswald Summer Shool Math & Comm July 2007 1. Introdution A pilot tone design problem and PAR optimization
More informationGrasp Planning: How to Choose a Suitable Task Wrench Space
Grasp Planning: How to Choose a Suitable Task Wrenh Spae Ch. Borst, M. Fisher and G. Hirzinger German Aerospae Center - DLR Institute for Robotis and Mehatronis 8223 Wessling, Germany Email: [Christoph.Borst,
More informationStructural Reconfiguration of Systems under Behavioral Adaptation
Strutural Reonfiguration of Systems under Behavioral Adaptation Carlos Canal a, Javier Cámara b, Gwen Salaün a Department of Computer Siene, University of Málaga, Spain b Department of Informatis Engineering,
More informationThe simulation analysis of the bridge rectifier continuous operation in AC circuit
Computer Appliations in Eletrial Engineering Vol. 4 6 DOI 8/j.8-448.6. The simulation analysis of the bridge retifier ontinuous operation in AC iruit Mirosław Wiślik, Paweł Strząbała Kiele University of
More informationON THE MOVING BOUNDARY HITTING PROBABILITY FOR THE BROWNIAN MOTION. Dobromir P. Kralchev
Pliska Stud. Math. Bulgar. 8 2007, 83 94 STUDIA MATHEMATICA BULGARICA ON THE MOVING BOUNDARY HITTING PROBABILITY FOR THE BROWNIAN MOTION Dobromir P. Kralhev Consider the probability that the Brownian motion
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 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 informationCommunicating Special Relativity Theory s Mathematical Inconsistencies
Communiating Speial Relatiity Theory s Mathematial Inonsistenies Steen B Bryant Primitie Logi, In, 704 Sansome Street, San Franiso, California 94111 Stee.Bryant@RelatiityChallenge.Com Einstein s Speial
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 informationExperiment 3: Basic Electronic Circuits II (tbc 1/7/2007)
Experiment 3: Basi Eletroni iruits II (tb /7/007) Objetive: a) To study the first-order dynamis of a apaitive iruits with the appliation of Kirhoff s law, Ohm s law and apaitane formula. b) To learn how
More informationPHYSICS 104B, WINTER 2010 ASSIGNMENT FIVE SOLUTIONS N { N 2 + 1, N 2 + 2, N
PHYSICS 4B, WINTER 2 ASSIGNMENT FIVE SOLUTIONS [.] For uniform m =. and springs g = 2. the eigenvalues have the analyti form, The ode jaobi test yields λ k = 2g ( os k) m k = k n = 2πn N = 2π N { N 2 +,
More informationReconstruction of the Dirichlet to Neumann Map for a Resistor Network
Reconstruction of the Dirichlet to Neumann Map for a Resistor Network Adrian V. Mariano University of Washington Miriam A. Myjak Seattle University June 21, 23 Abstract We describe methods of reconstructing
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 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 informationSOME CONSIDERATIONS IN THE ADJUSTMENT OF GPS-DERIVED BASELINES IN THE NETWORK MODE
SOME CONSIDERATIONS IN THE ADJUSTMENT OF GPS-DERIVED BASELINES IN THE NETWORK MODE A Thesis Presented in Partial Fulfillment of the Requirements for The degree Master of Siene in the Graduate Shool of
More informationMath 220A - Fall 2002 Homework 8 Solutions
Math A - Fall Homework 8 Solutions 1. Consider u tt u = x R 3, t > u(x, ) = φ(x) u t (x, ) = ψ(x). Suppose φ, ψ are supported in the annular region a < x < b. (a) Find the time T 1 > suh that u(x, t) is
More informationarxiv: v1 [physics.gen-ph] 5 Jan 2018
The Real Quaternion Relativity Viktor Ariel arxiv:1801.03393v1 [physis.gen-ph] 5 Jan 2018 In this work, we use real quaternions and the basi onept of the final speed of light in an attempt to enhane the
More informationeappendix for: SAS macro for causal mediation analysis with survival data
eappendix for: SAS maro for ausal mediation analysis with survival data Linda Valeri and Tyler J. VanderWeele 1 Causal effets under the ounterfatual framework and their estimators We let T a and M a denote
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 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 informationAnomaly cancellation and modularity, II: The E 8 E 8 case
SCIENCE CHINA Mathematis. ARTICLES. June 07 Vol. 60 No. 6: 985 994 doi: 0.007/s45-06-9034- Anomaly anellation and modularity, II: The E 8 E 8 ase In memory of Professor LU QiKeng 97 05 HAN Fei, LIU KeFeng,3
More informationA EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM.
A EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM. S. Kanagaraj Eulidean Relativity s.kana.raj@gmail.om (1 August 009) Abstrat By re-interpreting the speial relativity (SR) postulates based on Eulidean
More informationUse of prior information in the form of interval constraints for the improved estimation of linear regression models with some missing responses
Journal of Statistial Planning and Inferene 136 (2006) 2430 2445 www.elsevier.om/loate/jspi Use of prior information in the form of interval onstraints for the improved estimation of linear regression
More informationModel for a Distributed Radio Telescope. Patrick Fleckenstein
Model for a Distributed Radio Telesope Patrik Flekenstein Department of Mathematis, Rohester Institute of Tehnology, Rohester, NY 14623, U.S.A. E-mail address: pat@sh.rit.edu Abstrat. In this paper, a
More informationLecture 44. Better and successive approximations x2, x3,, xn to the root are obtained from
Lecture 44 Solution of Non-Linear Equations Regula-Falsi Method Method of iteration Newton - Raphson Method Muller s Method Graeffe s Root Squaring Method Newton -Raphson Method An approximation to the
More information9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry
9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived
More informationChapter 9. The excitation process
Chapter 9 The exitation proess qualitative explanation of the formation of negative ion states Ne and He in He-Ne ollisions an be given by using a state orrelation diagram. state orrelation diagram is
More informationA Functional Representation of Fuzzy Preferences
Theoretial Eonomis Letters, 017, 7, 13- http://wwwsirporg/journal/tel ISSN Online: 16-086 ISSN Print: 16-078 A Funtional Representation of Fuzzy Preferenes Susheng Wang Department of Eonomis, Hong Kong
More informationNAG Library Routine Document F08UBF (DSBGVX)
NAG Library Routine Document (DSBGVX) Note: before using this routine, please read the Users Note for your implementation to check the interpretation of bold italicised terms and other implementation-dependent
More informationCovariances for Bivariate Selection Model Second-Step Estimates
Eon 240B, Fall 2001 Daniel MFadden Covarianes for Bivariate Seletion Model Seond-Step Estimates Consider the bivariate latent variable model with normal disturbanes, (1) y * = + C, w * = zα + σν, where
More informationPerturbation Analyses for the Cholesky Factorization with Backward Rounding Errors
Perturbation Analyses for the holesky Fatorization with Bakward Rounding Errors Xiao-Wen hang Shool of omputer Siene, MGill University, Montreal, Quebe, anada, H3A A7 Abstrat. This paper gives perturbation
More information