ANGLE OF OF ARRIVAL ESTIMATION WITH A POLARIZATION DIVERSE ARRAY
|
|
- Timothy Evans
- 6 years ago
- Views:
Transcription
1 ANGLE OF OF ARRIVAL ESTIMATION WIT A POLARIZATION DIVERSE ARRAY DR. DR. WILLIAM ILLIAM P. P. BALLANCE ** RALP ALP A. A. COAN COAN MARC * No it Te Boeing Company WPB-
2 Problem formulation OUTLINE Obseration model is linear in polarization Maximum Likeliood Estimator MLE) Te MLE is identical to te Beamform, Witen, and Energy Detect BWE) tecnique of Allan Steinardt Problem is a minor ariation on an already soled MLE problem, te single snapsot -source problem Example to proide intuitie feel for solution Cramer-Rao Bound CRB) on estimation performance Agreement beteen MLE and CRB asymptotic efficiency) Performance impact due to required estimation of polarization Performance implications of assumed polarization Angle of arrial AOA) estimate generally as reduced standard deiation but is potentially biased WPB-
3 PROBLEM FORMULATION: ARRAY OBSERVATION MODEL Standard array obseration model for an N-element 3-dimensional array is x d )a + n ere [φ;]. Te direction ector is gien by jωoτ n ) [ d )] f )e ; n,, N n n ere f n ) is te complex-alued element pattern and te pase term is due to element separation. Te unknon model parameters are and a, for a total of 4 real-alued unknons. Te generalization for a polarization dierse array is simply x d, )a + n ere te polarization state is gien by [ ; ]. Te direction ector is gien by jωoτ n d, ) f, )e z c [ ] ) n n x φ c y Azimut φ 9 o - φ c, Eleation 9 o - c Time sift τ n ) n n + n c x sinφ cos + y cosφ cos z sin ) ere x n,y n, z n ) is te location of te n t element. WPB-3
4 PROBLEM FORMULATION: POLARIZATION SPECIFICS Electric field of a polarized ae traces an ellipse as time progresses E φˆ ˆ Transerse Plane ˆ β φˆ ere + Te real-alued quantities γ and η determine te ellipse orientation β and axial ratio r minor axis)/major axis). E cosγ E sinγ e jη Examples orizontal Vertical Linear RCP LCP Elliptical WPB-4
5 PROBLEM FORMULATION: POLARIZATION SPECIFICS CONT.) Just as electric field can be represented in terms of a basis set of polarizations e.g., orizontal and ertical, RCP and LCP, an elliptical polarization and its ortogonal polarization pol and cross-pol )), so can te element pattern. f n, ) f ) +,n f, n ) So tat d, ) d ) + d ) ere jωoτ n ) jωoτ n ) [ d )] f )e, [ d )] f ) e n,n Te unknon model parameters are, γ, η, and a, for a total of 6 real-alued unknons. Te model is linear in polarization en ieed in terms of te parameters and. n,n WPB-5
6 Rater tan begin it te likeliood function, begin it te resultant standard ML performance index and ten proide te requisite incremental deriation Maximizing te likeliood function is equialent to maximizing te folloing performance index Te key is tat te direction ector is linear in te unknon polarization parameters d, ) d ) + d ) D ) So tat ere J, ), MLE DERIVATION: AN INCREMENTALIST APPROAC, ) x, ) R, ) D ) ere R d, ) [ d ) d )], W ) [ ) )] ), ) ) + ) W ) R d ), ) R d ) R D ), WPB-6
7 MLE DERIVATION: AN INCREMENTALIST APPROAC CONT.) Define te by ector output of a orizontal and ertical polarization beamformer as y, and te corresponding by coariance matrix by R yy. Making use of, ) W ), y W ) x, R yy ) W ) R W ) Te performance index may be ritten J, ), ) x, ) R, ) R y ) yy ) ic is a Rayleig quotient in te unknon polarization state. Te closed-form ML estimate of polarization is ten R yy y Substituting tis into J results in te simple expression J ) y ) R- yy y ) ence, te performance index need only be searced oer. Note tat J is identical to te BWE performance index and so te BWE is te Maximum Likeliood Estimator. Tis assertion as first made by Roland Stougton of SAIC.) WPB-7
8 MLE DERIVATION: SUCCINCT PROOF FROM FIRST PRINCIPLES Te model may be ritten as x d, )a + n a d ) + a d ) + n D )a + n ere a a a. a Tis model is analogous to te single snapsot -source model x a d ) + a d ) + n D, + )a n except tat D [d ) d )] and a [a ; a ] are replaced by D [d ) d )] and a [a ; a ], respectiely. Te deriation readily follos. Te likeliood function is gien by N L π det R ) exp{ x D )a) R x D )a) } Maximizing te likeliood function is equialent to minimizing te nonlinear least squares cost function C x D )a) R x D )a) Witen for notational conenience C x D )a WPB-8
9 MLE DERIVATION: SUCCINCT PROOF FROM FIRST PRINCIPLES CONT.) Te problem is a separable nonlinear least squares problem. It is nonlinear in and linear in a, alloing a closed form solution for a. Substitution yields C I P) x ere P D D D ) D is a projection matrix. Minimizing C is equialent to maximizing or C a D J x D a) a D D ) D x x J P x R D D R D ) Q.E.D. Te BWE is te MLE only for te deterministic model it a single snapsot. For multiple snapsots, te model is [ d ) + d )]a k, ic is akin to d )a,k + d )a,k ere a,k β a,k, te fully correlated multipat problem toug for many snapsots, a stocastic-model ML solution may be preferred). x D P x R x y R yy y WPB-9
10 MLE: COERENT VS. NONCOERENT COMBINATION Te MLE for knon polarization coerently combines in polarization., ) ) + ) W ) α ) x N α ) Te MLE for unknon polarization noncoerently combines. * * J x N y Witener R yy z J Te polarization estimate is gien by R But because te pase must be estimated, it cannot be exploited to aciee te same leel of gain as in coerent combination for knon polarization. WPB- γ tan / ), yy y η / )
11 MLE PROCESSING FLOW If adaptie beamforming is performed before oter coerent processing stages, ten cannels per beam must be carried. N Pre- Beamforming Coerent Processing N Beamforming - it R yy itening zl Ver Post-Beamforming Coerent Processing Post-Beamforming Coerent Processing z If an explicit polarization estimate is desired, it sould be computed after te last stage in te coerent processing cain. Te improed SNR ill improe te quality of te estimate. In te example aboe, te by data z after te post-beamforming coerent processing stage is used. Te estimate is R yy z WPB-
12 SIMPLE EXAMPLE: BASIC ARRAY CONFIGURATION Underlying teory olds for a general array, it arbitrary element polarization responses and element locations. Tis includes elements tat are separated by many aelengts, as often found in space telescope applications. To proide an intuitie feel for te concepts inoled, oeer, consider a ery simple example. linear segments of dipoles, 3 elements eac, alf-ae spaced. 45 degree slope. Dipoles oriented along segment axes. Various segment separations considered. Direction ectors for broadside o azimut, o eleation) are T T d )[ ], d )[ ] 8 Element Locations 6 z aelengts) x aelengts) WPB-
13 SIMPLE EXAMPLE: RESPONSE OF ORIZONTAL & VERTICAL BEAMFORMERS Separate te to segments orizontally by 4 aelengts default config.). T orizontally polarized source at broadside, x )[ ] At broadside, d x peaks and d x. As scan off in azimut, ertical beamformer output is nonzero een toug source is orizontally polarized. Witener accounts for cross-coupling of beamformer. For tis example, te itener is simply te identity matrix for o eleation cut. orizontal -) and Vertical - -) Witened Beamformer Output 5 5 db Azimut degrees) WPB-3
14 SIMPLE EXAMPLE: MLE PERFORMANCE INDEX Solid plot is MLE performance index. Dased plots are responses of itened) orizontal and ertical beamformers. Te MLE performance index is te enelope of te to responses. Te mainlobe is idened, resulting in larger RMSE angle estimation tan for knon orizontal) polarization. 5 5 db Azimut degrees) WPB-4
15 SIMPLE EXAMPLE: POLARIZATION ESTIMATE VS. AZIMUT MLE polarization estimate aries it angle to maximize performance index. Witin te mainlobe te polarization estimate takes on many alues. For example, en ertical output is zero, estimate is orizontal polarization. 3 5 ZL RCP VER LCP ZL Pol. 5 db Azimut degrees) WPB-5
16 WPB-6 CRB ON AOA ESTIMATION CRB ON AOA ESTIMATION For te single snapsot) obseration model x ξ) + n, ere ξ is te M by unknon real-alued parameter ector and n is complex-gaussian it zero mean and coariance R, te M by M Fiser information matrix is gien by ere te m t column of te N by M matrix is ξ) /ξ m. For te unknon-polarization problem ae Te 6 columns of are gien by ] R [ Re al Γ )a ) d e sin ) d cos )a, d j γ γ η + )a d e sin j )a ) d e cos ) d sin - )a, d j a )e, d a )a ) d e sin ) d cos )a ) d e sin ) d cos j j a j j j γ η γ γ γ γ γ φ γ φ γ φ η η η η For te knon-polarization problem, only te first four columns are used.)
17 CRB ON AOA ESTIMATION CONT.) Te CRB is gien by te elements of B Γ - -. Denoting [ B ] m,n by b mn, one as b and b are te CRB on RMSE for φ and estimation, respectiely. b is te correlation coefficient beteen te φ and estimation. / b b ) Example equi-probablilty contour Scatter plot example Eleation Eleation Azimut Azimut To obtain CRB for te rotated coordinate system, compute te eigenalues and eigenectors ordered so tat λ λ ) of b b M b b λ and λ are te CRB on RMSE along te major and minor axes, respectiely. β tan [ u ] / [ u ] ) is te ellipse orientation, ere M u λ u. Direct numeric computation is preferred. Wile analytic expressions can be obtained using a partitioned matrix inerse approac, te resultant expressions do not lend great insigt due to teir complexity. To understand AOA result, study sape of corresponding MLE perf index and its indiidual components i.e., itened zl and er beamformer responses). WPB-7
18 AGREEMENT BETWEEN MLE AND CRB orizontally polarized source at broadside. No noise. Contours at -.5 db for eac performance index bot indices ae te same peak). Ellipse orientation and axial ratio agree it te respectie CRBs. Monte Carlo MLE RMSE performance agrees it te respectie CRBs. Knon Polarization: MLE Perf Index 5 4 db Unknon Polarization: MLE Perf Index 5 4 db Eleation degrees) Eleation degrees) Azimut degrees) Azimut degrees) -3 WPB-8
19 CRB RMSE ELLIPSES: UNKNOWN VS. KNOWN POLARIZATION orizontally polarized source. Source at broadside for first plot, at 3 o azimut and 3 o eleation for second plot. Solid plot is for unknon polarization, dased plot is for knon polarization. Te unknon polarization RMSE ellipse circumscribes te knon polarization RMSE ellipse..6 Source at Broadside Source at 3 o Azimut, 3 o Eleation Eleation degrees).. Eleation degrees) Azimut degrees) Azimut degrees) WPB-9
20 PERFORMANCE IMPACT OF ESTIMATING POLARIZATION: RMSE VS. SEPARATION FOR LINEAR SEGMENTS Circularly polarized source at broadside. Vary te orizontal separation beteen te segments. For knon polarization, te RMSE decreases as separation increases. For unknon polarization, te RMSE is independent of separation. For zero separation, te RMSE is te same. Minor Axis RMSE Major Axis RMSE RMSE degrees) Unknon Pol Knon Pol RMSE degrees) Unknon Pol Knon Pol Segment Separation aelengts) Segment Separation aelengts) WPB-
21 PERFORMANCE IMPACT OF ESTIMATING POLARIZATION: RMSE VS. SEP FOR MORE TAN LINEAR SEGMENTS Circularly polarized source at broadside. Add anoter linear segment, identical to one of te oters, but it fixed location miday beteen te oter to. Vary te orizontal separation beteen te outer segments. Te RMSE for unknon polarization no decreases it increasing array separation, but is still larger tan te RMSE for knon polarization. For linear segments of crossed dipoles 4 segments total), te RMSE is te same. Minor Axis RMSE Major Axis RMSE RMSE degrees) Unknon Pol Knon Pol RMSE degrees) Unknon Pol Knon Pol Segment Separation aelengts) Segment Separation aelengts) WPB-
22 Eleation degrees) Azimut degrees) WPB- IMPACT OF ASSUMED POLARIZATION: CIRCULARLY POLARIZED BEAMFORMER, ZL SOURCE Assumed-polarization performance using default config..5 db mismatc loss.8 o bias.6 of peak-to-null idt of unknon polarization performance index) Standard deiation commensurate it tat of correctly assumed polarization, but increased by loss due to mismatc 6% increase) Unknon-Pol Perf Index db Eleation degrees) Assumed-Pol Perf Index Azimut degrees) Contours at -.5 db from on peak x true AOA, o mean of estimated AOA Polarization Plot Beamformer polarization True source polarization db
23 Vertical pol component is 6dB don it respect to zl component, it 9º pase offset γ 6.6 o, η 9 o ) Assumed-polarization performance. db mismatc loss.49 o bias. of peak-to-null idt of unknon pol perf index) Standard deiation commensurate it tat of correctly assumed polarization, but increased by loss due to mismatc % increase) Eleation degrees) Azimut degrees) WPB-3 IMPACT OF ASSUMED POLARIZATION: ZL BEAMFORMER, ELLIPTICALLY POLARIZED SOURCE Unknon-Pol Perf Index db Eleation degrees) Contours at -.5 db from on peak x true AOA, o mean of estimated AOA Polarization Plot Beamformer polarization Assumed-Pol Perf Index Azimut degrees) True source polarization db
24 SUMMARY Model General 3-dimensional array consisting of elements it arbitrary polarization responses Found to be linear in polarization Maximum Likeliood Estimator MLE) Performance index need only be searced oer azimut and eleation) Problem analogous to te single-snapsot to-source problem Te BWE is te MLE for a single snapsot, toug not for multiple snapsots To deriations proided incremental and from first principles) Te MLE uses noncoerent combination in polarization, as compared to coerent for knon polarization Processing flo exemplified and explicit polarization estimation proided Example to proide intuitie feel for MLE estimator WPB-4
25 SUMMARY CONT.) Cramer-Rao Bound CRB) Agreement beteen MLE and CRB asymptotic efficiency), for unknon and knon polarization - Monte Carlo - Contour cut of performance indices Performance impact due to estimation of polarization - RMSE for unknon pol. case > RMSE for knon pol. case - Obsere RMSE ellipses circumscribed - For linear segments of dipoles, RMSE is independent of segment separation for unknon polarization Assumed Polarization Angle of arrial AOA) estimate standard deiation commensurate it tat for correctly assumed polarization, but increased by SNR loss due to polarization mismatc Potentially biased AOA estimate WPB-5
qwertyuiopasdfghjklzxcvbnmqwerty uiopasdfghjklzxcvbnmqwertyuiopasd fghjklzxcvbnmqwertyuiopasdfghjklzx cvbnmqwertyuiopasdfghjklzxcvbnmq
qwertyuiopasdfgjklzxcbnmqwerty uiopasdfgjklzxcbnmqwertyuiopasd fgjklzxcbnmqwertyuiopasdfgjklzx cbnmqwertyuiopasdfgjklzxcbnmq Projectile Motion Quick concepts regarding Projectile Motion wertyuiopasdfgjklzxcbnmqwertyui
More informationNotes on Linear Minimum Mean Square Error Estimators
Notes on Linear Minimum Mean Square Error Estimators Ça gatay Candan January, 0 Abstract Some connections between linear minimum mean square error estimators, maximum output SNR filters and the least square
More informationAstrometric Errors Correlated Strongly Across Multiple SIRTF Images
Astrometric Errors Correlated Strongly Across Multiple SIRTF Images John Fowler 28 March 23 The possibility exists that after pointing transfer has been performed for each BCD (i.e. a calibrated image
More informationPredictive Coding. Lossy or lossless. Feedforward or feedback. Intraframe or interframe. Fixed or Adaptive
Predictie Coding Predictie coding is a compression tecnique based on te difference between te original and predicted alues. It is also called DPCM Differential Pulse Code Modulation Lossy or lossless Feedforward
More informationLINEAR MULTISCALE AUTOREGRESSIVE MODEL FOR FORECASTING SEASONAL TIME SERIES DATA
INEAR UTISCAE AUTOREGRESSIVE ODE FOR FORECASTING SEASONA TIE SERIES DATA Brodjol S.S.U. Suartono and A.J. Endarta Department of Statistics Faculty of atematic and Natural Sciences Sepulu Nopember Institute
More informationDrive power of a discharge belt conveyor comparison of calculation, simulation and measurement
Martin Egger Professor Uniersity of pplied Sciences Upper ustria Scool of Engineering and Enironmental Sciences Klaus Decker Senior Scientist Vienna Uniersity of Tecnology aculty of Mecanical and Industrial
More informationEstimation of Efficiency with the Stochastic Frontier Cost. Function and Heteroscedasticity: A Monte Carlo Study
Estimation of Efficiency ith the Stochastic Frontier Cost Function and Heteroscedasticity: A Monte Carlo Study By Taeyoon Kim Graduate Student Oklahoma State Uniersity Department of Agricultural Economics
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationERAD THE SEVENTH EUROPEAN CONFERENCE ON RADAR IN METEOROLOGY AND HYDROLOGY
Multi-beam raindrop size distribution retrieals on the oppler spectra Christine Unal Geoscience and Remote Sensing, TU-elft Climate Institute, Steinweg 1, 68 CN elft, Netherlands, c.m.h.unal@tudelft.nl
More informationTangent Lines-1. Tangent Lines
Tangent Lines- Tangent Lines In geometry, te tangent line to a circle wit centre O at a point A on te circle is defined to be te perpendicular line at A to te line OA. Te tangent lines ave te special property
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More informationA.P. CALCULUS (AB) Outline Chapter 3 (Derivatives)
A.P. CALCULUS (AB) Outline Capter 3 (Derivatives) NAME Date Previously in Capter 2 we determined te slope of a tangent line to a curve at a point as te limit of te slopes of secant lines using tat point
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationEdge Detection Based on the Newton Interpolation s Fractional Differentiation
Te International Arab Journal of Information Tecnology, Vol. 11, No. 3, May 014 3 Edge Detection Based on te Newton Interpolation s Fractional Differentiation Caobang Gao 1,, Jiliu Zou, 3, and Weiua Zang
More informationv represents the uncertainties of the model. The observations are related
HARMO13-1-4 June 010, Paris, France - 13t Conference on Harmonisation witin Atmosperic Dispersion Modelling for Regulatory Purposes H13-137 DATA ASSIMILATION IN AIR QUALITY MODELLING OVER PO VALLEY REGION
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationMath 262 Exam 1 - Practice Problems. 1. Find the area between the given curves:
Mat 6 Exam - Practice Problems. Find te area between te given curves: (a) = x + and = x First notice tat tese curves intersect wen x + = x, or wen x x+ =. Tat is, wen (x )(x ) =, or wen x = and x =. Next,
More informationTHE LANDAUER LIMIT AND THERMODYNAMICS OF BIOLOGICAL SYSTEMS
THE LANDAUER LIMIT AND THERMODYNAMICS OF BIOLOGICAL SYSTEMS Daid H. Wolpert Santa Fe Institute, MIT ttp://daidwolpert.weebly.com Jos Grocow Santa Fe Institute HEAT OF ERASING A BIT a) b) c) R 1 2 2R 1
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
More informationKinematics on oblique axes
Bolina 1 Kinematics on oblique axes arxi:physics/01111951 [physics.ed-ph] 27 No 2001 Oscar Bolina Departamento de Física-Matemática Uniersidade de São Paulo Caixa Postal 66318 São Paulo 05315-970 Brasil
More informationMEASUREMENTS OF POLARIMETRIC PARAMETERS AT LOW SIGNAL-TO-NOISE RATIOS
P9. MEAUEMENT OF POLAIMETIC PAAMETE AT LOW IGNAL-TO-NOIE ATIO V. M. Melniko and D.. rnic Cooperatie Institute for Mesoscale Meteorological tudies, Uniersity of Oklaoa and NOAA/National eere tors Laboratory,
More informationINTRODUCTION TO CALCULUS LIMITS
Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and
More informationLast lecture (#4): J vortex. J tr
Last lecture (#4): We completed te discussion of te B-T pase diagram of type- and type- superconductors. n contrast to type-, te type- state as finite resistance unless vortices are pinned by defects.
More informationASYMPTOTIC MEAN SQUARED ERROR PERFORMANCE OF DIAGONALLY LOADED CAPON-MVDR PROCESSOR
ASYMPTOTIC MEAN SQUARED ERROR PERFORMANCE OF DIAGONALLY LOADED CAPON-MVDR PROCESSOR Crist D. Ricmond Lincoln Laboratory Massacusetts Institute of Tecnology Lexington, MA 02420 email: crist@ll.mit.edu Raj
More informationSIMG Solution Set #5
SIMG-303-0033 Solution Set #5. Describe completely te state of polarization of eac of te following waves: (a) E [z,t] =ˆxE 0 cos [k 0 z ω 0 t] ŷe 0 cos [k 0 z ω 0 t] Bot components are traveling down te
More informationExcluded Volume Effects in Gene Stretching. Pui-Man Lam Physics Department, Southern University Baton Rouge, Louisiana
Excluded Volume Effects in Gene Stretcing Pui-Man Lam Pysics Department, Soutern University Baton Rouge, Louisiana 7083 Abstract We investigate te effects excluded volume on te stretcing of a single DNA
More informationON THE CONSTRUCTION OF REGULAR A-OPTIMAL SPRING BALANCE WEIGHING DESIGNS
Colloquium Biometricum 43 3, 3 9 ON THE CONSTRUCTION OF REGUAR A-OPTIMA SPRING BAANCE WEIGHING DESIGNS Bronisław Ceranka, Małgorzata Graczyk Department of Matematical and Statistical Metods Poznań Uniersity
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
More informationOn seismic landslide hazard assessment
Yang, J. (27). Géotecnique 57, No. 8, 77 713 doi: 1.168/geot.27.57.8.77 TECHNICAL NOTE On seismic landslide azard assessment J. YANG* KEYWORDS: eartquakes; landslides; slopes INTRODUCTION Seismic landslides
More informationThe Importance of Anisotropy for Prestack Imaging
The Importance of Anisotropy for Prestack Imaging Xiaogui Miao*, Daid Wilkinson and Scott Cheadle Veritas DGC Inc., 715-5th Ae. S.W., Calgary, AB, T2P 5A2 Calgary, AB Xiaogui_Miao@eritasdgc.com ABSTRACT
More informationNumerical analysis of the influence of porosity and void size on soil stiffness using random fields
Numerical analysis of te influence of porosity and oid size on soil stiffness using random fields D.V. Griffits, Jumpol Paiboon, Jinsong Huang Colorado Scool of Mines, Golden, CO, USA Gordon A. Fenton
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationGiant Squid - Hidden Canard: the 3D Geometry of the Hodgkin-Huxley Model
Noname manuscript No. (will be inserted by te editor) Jonatan Rubin Martin Wecselberger Giant Squid - Hidden Canard: te 3D Geometry of te Hodgkin-Huxley Model te date of receipt and acceptance sould be
More informationBayesian ML Sequence Detection for ISI Channels
Bayesian ML Sequence Detection for ISI Cannels Jill K. Nelson Department of Electrical and Computer Engineering George Mason University Fairfax, VA 030 Email: jnelson@gmu.edu Andrew C. Singer Department
More informationContinuous Stochastic Processes
Continuous Stocastic Processes Te term stocastic is often applied to penomena tat vary in time, wile te word random is reserved for penomena tat vary in space. Apart from tis distinction, te modelling
More information[db]
Blind Source Separation based on Second-Order Statistics wit Asymptotically Optimal Weigting Arie Yeredor Department of EE-Systems, el-aviv University P.O.Box 3900, el-aviv 69978, Israel Abstract Blind
More informationTeaching Differentiation: A Rare Case for the Problem of the Slope of the Tangent Line
Teacing Differentiation: A Rare Case for te Problem of te Slope of te Tangent Line arxiv:1805.00343v1 [mat.ho] 29 Apr 2018 Roman Kvasov Department of Matematics University of Puerto Rico at Aguadilla Aguadilla,
More information1 :: Mathematical notation
1 :: Mathematical notation x A means x is a member of the set A. A B means the set A is contained in the set B. {a 1,..., a n } means the set hose elements are a 1,..., a n. {x A : P } means the set of
More informationENGI Gradient, Divergence, Curl Page 5.01
ENGI 940 5.0 - Gradient, Divergence, Curl Page 5.0 5. e Gradient Operator A brief review is provided ere for te gradient operator in bot Cartesian and ortogonal non-cartesian coordinate systems. Sections
More information1 2 x Solution. The function f x is only defined when x 0, so we will assume that x 0 for the remainder of the solution. f x. f x h f x.
Problem. Let f x x. Using te definition of te derivative prove tat f x x Solution. Te function f x is only defined wen x 0, so we will assume tat x 0 for te remainder of te solution. By te definition of
More informationFundamentals of Concept Learning
Aims 09s: COMP947 Macine Learning and Data Mining Fundamentals of Concept Learning Marc, 009 Acknowledgement: Material derived from slides for te book Macine Learning, Tom Mitcell, McGraw-Hill, 997 ttp://www-.cs.cmu.edu/~tom/mlbook.tml
More information(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?
Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is
More information& Tech Note: Jin Huang (M.I.T.), Yi Qiang (JLab) Hall A Analysis Workshop Dec 8, JLab
& Tec Note: ttp://www.jlab.org/~jinuang/transversity/mle.pdf Jin Huang (M.I.T., Yi Qiang (JLab Hall A Analysis Worksop Dec 8, 010 @ JLab Wat are Azimutal Asymmetries TMD & SIDIS Definition Related Experiments
More informationNew Distribution Theory for the Estimation of Structural Break Point in Mean
New Distribution Teory for te Estimation of Structural Break Point in Mean Liang Jiang Singapore Management University Xiaou Wang Te Cinese University of Hong Kong Jun Yu Singapore Management University
More informationWeek #15 - Word Problems & Differential Equations Section 8.2
Week #1 - Word Problems & Differential Equations Section 8. From Calculus, Single Variable by Huges-Hallett, Gleason, McCallum et. al. Copyrigt 00 by Jon Wiley & Sons, Inc. Tis material is used by permission
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More informationTe comparison of dierent models M i is based on teir relative probabilities, wic can be expressed, again using Bayes' teorem, in terms of prior probab
To appear in: Advances in Neural Information Processing Systems 9, eds. M. C. Mozer, M. I. Jordan and T. Petsce. MIT Press, 997 Bayesian Model Comparison by Monte Carlo Caining David Barber D.Barber@aston.ac.uk
More information(a) Taking the derivative of the position vector with respect to time, we have, in SI units (m/s),
Chapter 4 Student Solutions Manual. We apply Eq. 4- and Eq. 4-6. (a) Taking the deriatie of the position ector with respect to time, we hae, in SI units (m/s), d ˆ = (i + 4t ˆj + tk) ˆ = 8tˆj + k ˆ. dt
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationMA 266 Review Topics - Exam # 2 (updated)
MA 66 Reiew Topics - Exam # updated Spring First Order Differential Equations Separable, st Order Linear, Homogeneous, Exact Second Order Linear Homogeneous with Equations Constant Coefficients The differential
More informationCOLLEGE WB3413 DREDGING PROCESSES
COLLEGE WB343 DREDGING PROCESSES Te Breacing Process (September 003) Prof. ir. W.J. Vlasblom . Intoduction Breacing is te occurrence of continues and/or local instabilities on a sandy slope causing a density
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationPredicting the Intelligibility of Noisy and Non-Linearly Processed Binaural Speech
1 Predicting te Intelligibility of Noisy and Non-Linearly Processed Binaural Speec Asger Heidemann Andersen, Jan Mark de Haan, Zeng-Hua Tan, and Jesper Jensen Abstract Objective speec intelligibility measures
More informationNotes on Planetary Motion
(1) Te motion is planar Notes on Planetary Motion Use 3-dimensional coordinates wit te sun at te origin. Since F = ma and te gravitational pull is in towards te sun, te acceleration A is parallel to te
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationChapter 1 Functions and Graphs. Section 1.5 = = = 4. Check Point Exercises The slope of the line y = 3x+ 1 is 3.
Capter Functions and Graps Section. Ceck Point Exercises. Te slope of te line y x+ is. y y m( x x y ( x ( y ( x+ point-slope y x+ 6 y x+ slope-intercept. a. Write te equation in slope-intercept form: x+
More informationChapter 1D - Rational Expressions
- Capter 1D Capter 1D - Rational Expressions Definition of a Rational Expression A rational expression is te quotient of two polynomials. (Recall: A function px is a polynomial in x of degree n, if tere
More informationDifferential Geometry of Surfaces
Differential Geometry of urfaces Jordan mith and Carlo équin C Diision, UC Berkeley Introduction These are notes on differential geometry of surfaces ased on reading Greiner et al. n. d.. Differential
More informationRepresentational Power of Restricted Boltzmann Machines
1 Representational Power of Restricted Boltzmann Macines and Deep Belief Networks Nicolas Le Roux and Yosua Bengio Dept. IRO, Uniersité de Montréal C.P. 6128, Montreal, Qc, H3C 3J7, Canada {lerouxni,bengioy}@iro.umontreal.ca
More informationEIGENVALUES AND EIGENVECTORS
CHAPTER 6 EIGENVALUES AND EIGENVECTORS SECTION 6. INTRODUCTION TO EIGENVALUES In each of Problems we first list the characteristic polynomial p( λ) A λi of the gien matrix A and then the roots of p( λ
More informationTrajectory Estimation for Tactical Ballistic Missiles in Terminal Phase Using On-line Input Estimator
Proc. Natl. Sci. Counc. ROC(A) Vol. 23, No. 5, 1999. pp. 644-653 Trajectory Estimation for Tactical Ballistic Missiles in Terminal Phase Using On-line Input Estimator SOU-CHEN LEE, YU-CHAO HUANG, AND CHENG-YU
More informationLines, Conics, Tangents, Limits and the Derivative
Lines, Conics, Tangents, Limits and te Derivative Te Straigt Line An two points on te (,) plane wen joined form a line segment. If te line segment is etended beond te two points ten it is called a straigt
More informationarxiv: v1 [hep-ex] 23 Mar 2018
EUROPEAN ORGANIZAION FOR NUCLEAR RESEARCH (CERN CERN-EP-08-039 08/03/6 -OP7-00 arxiv:803.08856v [ep-ex 3 Mar 08 Measurement of differential cross sections for te production of top quark pairs and of additional
More information1 Introduction. Control 2:
Introduction Kinematics of Wheeled Mobile Robots (No legs here. They ork like arms. We e seen arms already) For control purposes, the kinematics of heeled mobile robots (WMRs) that e care about are the
More informationFinancial Econometrics Prof. Massimo Guidolin
CLEFIN A.A. 2010/2011 Financial Econometrics Prof. Massimo Guidolin A Quick Review of Basic Estimation Metods 1. Were te OLS World Ends... Consider two time series 1: = { 1 2 } and 1: = { 1 2 }. At tis
More information130 Chapter 3 Differentiation
0 Capter Differentiation 20. (a) (b) 2. C position, A velocity, an B acceleration. Neiter A nor C can be te erivative of B because B's erivative is constant. Grap C cannot be te erivative of A eiter, because
More informationPre-lab Quiz/PHYS 224 Earth s Magnetic Field. Your name Lab section
Pre-lab Quiz/PHYS 4 Eart s Magnetic Field Your name Lab section 1. Wat do you investigate in tis lab?. For a pair of Helmoltz coils described in tis manual and sown in Figure, r=.15 m, N=13, I =.4 A, wat
More informationMath 312 Lecture Notes Modeling
Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a
More informationCalculus I Practice Exam 1A
Calculus I Practice Exam A Calculus I Practice Exam A Tis practice exam empasizes conceptual connections and understanding to a greater degree tan te exams tat are usually administered in introductory
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More information7 Semiparametric Methods and Partially Linear Regression
7 Semiparametric Metods and Partially Linear Regression 7. Overview A model is called semiparametric if it is described by and were is nite-dimensional (e.g. parametric) and is in nite-dimensional (nonparametric).
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More information6.1.1 Angle between Two Lines Intersection of Two lines Shortest Distance from a Point to a Line
CHAPTER 6 : VECTORS 6. Lines in Space 6.. Angle between Two Lines 6.. Intersection of Two lines 6..3 Shortest Distance from a Point to a Line 6. Planes in Space 6.. Intersection of Two Planes 6.. Angle
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationFORCED CONVECTION FILM CONDENSTION OF DOWNWARD-FLOWING VAPOR ON HORIZONTAL TUBE WITH WALL SUCTION EFFECT
Journal of Marine Science and Tecnology DOI: 1.6119/JMST-11-15-1 Tis article as been peer reiewed and accepted for publication in JMST but as not yet been copyediting, typesetting, pagination and proofreading
More informationRules of Differentiation
LECTURE 2 Rules of Differentiation At te en of Capter 2, we finally arrive at te following efinition of te erivative of a function f f x + f x x := x 0 oing so only after an extene iscussion as wat te
More informationFast Computation of Capacitance Matrix and Potential Distribution for Multiconductor in Non-Homogenous Multilayered Dielectric Media
Excerpt from te Proceedings of te OMSOL onference 2009 Boston Fast omputation of apacitance Matrix and Potential Distribution for Multiconductor in Non-Homogenous Multilayered Dielectric Media S. M. Musa
More informationResidual migration in VTI media using anisotropy continuation
Stanford Exploration Project, Report SERGEY, Noember 9, 2000, pages 671?? Residual migration in VTI media using anisotropy continuation Tariq Alkhalifah Sergey Fomel 1 ABSTRACT We introduce anisotropy
More informationRegularized Regression
Regularized Regression David M. Blei Columbia University December 5, 205 Modern regression problems are ig dimensional, wic means tat te number of covariates p is large. In practice statisticians regularize
More informationA Geometric Review of Linear Algebra
A Geometric Reiew of Linear Algebra The following is a compact reiew of the primary concepts of linear algebra. The order of presentation is unconentional, with emphasis on geometric intuition rather than
More informationName: Answer Key No calculators. Show your work! 1. (21 points) All answers should either be,, a (finite) real number, or DNE ( does not exist ).
Mat - Final Exam August 3 rd, Name: Answer Key No calculators. Sow your work!. points) All answers sould eiter be,, a finite) real number, or DNE does not exist ). a) Use te grap of te function to evaluate
More informationKernel Density Based Linear Regression Estimate
Kernel Density Based Linear Regression Estimate Weixin Yao and Zibiao Zao Abstract For linear regression models wit non-normally distributed errors, te least squares estimate (LSE will lose some efficiency
More informationMultilag Correlation Estimators for Polarimetric Radar Measurements in the Presence of Noise
77 J O U R N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y VOLUME 9 Multilag Correlation Estimators for Polarimetric Radar Measurements in te Presence of Noise LEI LEI,*,# GUIFU
More information64 IX. The Exceptional Lie Algebras
64 IX. Te Exceptional Lie Algebras IX. Te Exceptional Lie Algebras We ave displayed te four series of classical Lie algebras and teir Dynkin diagrams. How many more simple Lie algebras are tere? Surprisingly,
More informationEconometrics II - EXAM Outline Solutions All questions have 25pts Answer each question in separate sheets
Econometrics II - EXAM Outline Solutions All questions hae 5pts Answer each question in separate sheets. Consider the two linear simultaneous equations G with two exogeneous ariables K, y γ + y γ + x δ
More informationDiscrete Monogenic Functions in Image processing
Discrete Function Teory Discrete Monogenic Functions in Image processing U. Käler CIDMA and Departamento de Matemática Universidade de Aveiro ukaeler@ua.pt MOIMA Scloss Herrenausen Hannover, June 20-24,
More informationSpace Probe and Relative Motion of Orbiting Bodies
Space robe and Relatie Motion of Orbiting Bodies Eugene I. Butiko Saint etersburg State Uniersity, Saint etersburg, Russia E-mail: e.butiko@phys.spbu.ru bstract. Seeral possibilities to launch a space
More informationAssignment 4 (Solutions) NPTEL MOOC (Bayesian/ MMSE Estimation for MIMO/OFDM Wireless Communications)
Assignment 4 Solutions NPTEL MOOC Bayesian/ MMSE Estimation for MIMO/OFDM Wireless Communications The system model can be written as, y hx + The MSE of the MMSE estimate ĥ of the aboe mentioned system
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationSEISMIC PASSIVE EARTH PRESSURE WITH VARYING SHEAR MODULUS: PSEUDO-DYNAMIC APPROACH
IGC 29, Guntur, INDIA SEISMIC PASSIE EART PRESSURE WIT ARYING SEAR MODULUS: PSEUDO-DYNAMIC APPROAC P. Gos Assistant Professor, Deartment of Ciil Engineering, Indian Institute of Tecnology Kanur, Kanur
More information