Uncertainty and Parameter Space Analysis in Visualization -
|
|
- Muriel McCoy
- 5 years ago
- Views:
Transcription
1 Uncertaint and Parameter Space Analsis in Visualiation - Session 4: Structural Uncertaint Analing the effect of uncertaint on the appearance of structures in scalar fields Rüdiger Westermann and Tobias Pfaffelmoser Computer Graphics & Visualiation Specific literature [1] C. Johnson and A. Sanderson, A net step: Visualiing errors and uncertaint, Computer Graphics and Applications, IEEE, vol. 23, no. 5, pp. 6 10, [2] A. Pang, C. Wittenbrink, and S. Lodha, Approaches to uncertaint visualiation, The Visual Computer, vol. 13, no. 8, pp ,1997. [3] T. Pfaffelmoser, M. Reitinger, and R. Westermann, Visualiing the positional and geometrical variabilit of isosurfaces in uncertain scalar fields, in Computer Graphics Forum, vol. 30, no. 3. Wile Online Librar, 2011, pp [4] T. Pfaffelmoser and R. Westermann, Visualiation of global correlation structures in uncertain 2d scalar fields, in Computer Graphics Forum, vol. 31, no. 3. Wile Online Librar, 2012, pp [5] T. Pfaffelmoser and R. Westermann, Correlation Visualiation for Structural Uncertaint Analsis, in Journal of Uncertaint Quantification, 2012 Facult of Informatics Session 4: Structural uncertaint Part 1: Structural uncertaint Part 1: Structural uncertaint [3] [5] Stochastic modelling of uncertaint Variance as an uncertaint indicator Limitations of variance as an uncertaint indicator Definition of structural uncertaint and correlation Part 2: Local and global correlation visualiation [4] [5] Requirements and challenges Local correlation tensor Glph based correlation visualiation Global correlation clustering Concluding remarks Uncertain scalar fields Uncertain scalar field Y: M R on compact domain M R n Discretel sampled on nodes i I with values Y i I We assume a stochastic uncertaint model: Y i I are (correlated) Gaussian distributed random variables 1, 2, 3... Described b a multivariate Gaussian distribution 1 p = ep ( 0.5 μ T Σ 1 ( μ)) 2π detσ with mean μ = (μ1,, μn) and covariance matri Σ Y i = μ i is the most likel configuration Multivariate Gaussian distribution Uncertaint parameters in the multivariate case The degree of uncertaint is modeled b the Σ 1 p = 2π detσ ep ( 0.5 μ T Σ 1 ( μ)) For 3 grid points: Σ = 2 σ 1 σ 12 σ 13 σ 12 2 σ 2 σ 23 σ 13 σ 23 2 σ 3 σ i 2 = var Y i = E Y i μ i 2 σ ij = Cov Y i, Y j = E( Y i μ i (Y j μ j )) 1
2 Multivariate Gaussian distribution Compute μ i, σ ij from a given ensemble of data sets Ensemble members are treated as realiations of a corresponding multivariate random variable ensemble μ i σ i Specif μ i, σ ij and map independent normall distributed random number vectors to a corresponding realiation Σ Gaussian covariance matri Variance information is often used as primar indicator of the uncertaint (square of the standard deviation σ) Σ = 2 σ 1 σ 12 σ 13 σ 12 2 σ 2 σ 23 σ 13 σ 23 2 σ 3 Measure of the amount of variation of the values of a random variable Often visualied directl via confidence regions, uncertaint glphs, or specific color or opacit mappings Variance as uncertaint indicator Eample: an ensemble of 2D seismic tomograph wave velocities Relative velocities are color coded from blue (negative) to red (positive)) ensemble μ i σ i Variance as uncertaint indicator Uncertaint mapping in 3D via Stochastic Distance Functions μ i θ Ψ θ i ma σ i, σ min Mahalanobis distance: in numbers of standard deviation to the level- surface in the mean field ( Two prominent circular features are observed in the mean values Mapping the standard deviations to colors shows low, and almost constant uncertaint in both regions SDF surfaces are level-sets in the Mahalanobis distance field, enclosing the volume which contains the uncertain surface with a certain probabilit For Gaussian distributions, the less standard deviations an observation is from the mean, the higher the probabilit of this observation Variance as uncertaint indicator Uncertaint visualiation for an iso-surface in the mean values of a 3D temperature ensemble Variance as uncertaint indicator Uncertaint visualiation for an iso-surface in the mean values of a 3D temperature ensemble Color coding of stochastic distance Color coding of Euclidian distance along normal curves [3] Iso-surface in the mean values of a 3D temperature field Positional variabilit of the surface due to uncertaint 2
3 The standard deviation describes the local uncertaint but does not allow inferring on possible variations at different positions relative to each other e.g., if mean values and standard deviations at two adjacent points are identical, but the realiations of the random variables at both points behave independentl, it cannot be predicted whether there is a positive, ero, or negative derivative of the data between the two points In general, the effect of uncertaint on structures which depend on random values at multiple points cannot be predicted from the mean values and standard deviations alone The effect is to a large etent arbitrar if the realiations of the random values are stochasticall independent, while a structure s shape can be assumed stable if the realiations are stochasticall dependent Variations of an uncertain curve due to uncertaint Curve points Y i = Y i are modelled via a multivariate Gaussian random variable with smoothl varing mean values (green curve) and a constant standard deviation (blue curves show corresponding confidence interval) realiations realiations In a) and b), high and low stochastic between the values at adjacent i was modelled between scalar values was modelled eplicitl : between scalar values was modelled eplicitl : (5) along (5) along between scalar values was modelled eplicitl : between scalar values was modelled eplicitl : (5) along (5) along 3
4 between scalar values was modelled eplicitl : between scalar values was modelled eplicitl : (5) along (5) along between scalar values was modelled eplicitl : between scalar values was modelled eplicitl : (5) along (5) along between scalar values was modelled eplicitl : between scalar values was modelled eplicitl : (5) along (5) along 4
5 Random values in the 3D scalar field shown before were generated using a multivariate Gaussian distribution with constant standard deviation and linearl increasing mean (along ) Conclusion: just b looking at the standard deviation it is impossible to infer on the stabilit (or confidence) of structures in the mean values A structure is ver likel to change its shape due to uncertaint, if the stochastic between the data values making up the region is low A structure is rather unlikel to change (it is stochasticall stable), if the stochastic between the data values making up the structure is high, even if the standard deviations are high (a) Iso-surface in the mean values (b) containing all points that belong to the surface with a certain probabilit (c) Occurrence of the surface in (a) for one possible realiation of the random values Thus, one important goal in uncertaint visualiation is to conve the stochastic and conclude on the stabilit of structures in the data Structural uncertaint We call the effect of uncertaint in the data values on structures depending on the values at multiple points a structural uncertaint It is associated with the occurrence of particular structures in the data which are affected b the degree of between the values at two or more data points Indicators for structural uncertaint are given b the subdiagonals of Σ, which contain covariance information on relative uncertainties between random variables Structural uncertaint For Gaussian distributed random variables, the stochastic of random values at two points is given b the correlation: ρ ij = Σ ij Σ ii Σ jj = σ ij σ i σ j, 1 ρ ij 1 Positive correlation: ρ ij > 0 Σ = 2 σ 1 σ 12 σ 13 σ 12 2 σ 2 σ 23 σ 13 σ 23 2 σ 3 Stochastic in: ρ ij = 0 Negative/ correlation: ρ ij < 0 Variabilit of a scalar value around its mean position: Variabilit of a scalar value around its mean position: 5
6 Variabilit of a scalar value around its mean position: Variabilit of a scalar value around its mean position: Variabilit of a scalar value around its mean position: Variabilit of a scalar value around its mean position: Variabilit of a scalar value around its mean position: Correlation describes relative variabilit of one or more random variables Standard Deviation Confidence Interval Mean Position Positive Correlation Negative Correlation 6
7 Correlation describes relative variabilit of one or more random variables Correlation describes relative variabilit of one or more random variables Positive Correlation Negative Correlation Positive Correlation Negative Correlation Correlation describes relative variabilit of one or more random variables Part 2: Local and global correlation visualiation 2.1 Local Correlation Visualiation Positive Correlation Negative Correlation Structural uncertaint Structural uncertaint Remember the eample in Part 1: locall varing structural variations of an iso-surface in different realiations of an uncertain 3D scalar field Caused b different correlation structures Resulting in stochasticall stable and unstable surface parts Cannot be revealed b mean and standard deviation values High structural uncertaint Low structural uncertaint Problems addressed in this part of the tutorial: 1. How can one visualie correlations to anale the structural uncertaint of particular features in 2D and 3D scalar fields, given a set of realiations or a (Gaussian based) stochastic uncertaint model Two different approaches will be discussed, both aim at conveing correlation structures in the data: 1. Local correlation analsis using characteristic correlation tensors 2. Global correlation analsis showing longe-range dependencies such as inverse correlations via correlation clusters 7
8 Correlation Correlation Local correlation analsis Assumption of a distance dependent correlation model: higher correlations between the data values at points with shorter Euclidean distance Local correlation analsis Assumption of a distance dependent correlation model implies that correlation strength depends on grid resolution Weak correlation Strong correlation Domain of a 2D scalar data set Low correlation to 1-ring neighborhood High correlation to 1-ring neighborhood Distance dependent correlation model The Gaussian correlation function relates correlation to spatial distance Distance dependent correlation model Correlation strength parameter τ relates correlation to distance and is independent of grid resolution for local correlation analsis 1 τ 1 τ 2 τ 3 correlation strength parameter 1 τ 1 τ 2 τ 3 0 Euclidean distance 0 Euclidean distance Correlation anisotrop Problem: Correlation values and thus correlation strength parameters depend on direction, ielding a directional correlation distribution at each grid point τ 2 τ 3 τ 1 Correlation strength tensor Idea: The (anisotropic) distribution of the correlation strength parameter is modelled b a rank-2 tensor T The correlation strength tensor T models the correlation strength to 8 (2D) or 26 (3D) neighbors Given T, correlation strength τ at position into direction r can be computed via tensor-vector multiplication A single correlation value can not represent the stochastic s! 8
9 Correlation strength tensor Given the tensor T, the distance and direction dependent correlation becomes Correlation strength tensor computation T is smmetric and has 3 (2D) or 6 (3D) free values to be determined T models the correlation strength to 8 (2D) or 26 (3D) neighbors T is obtained for ever grid point b solving an over-determined linear sstem using a least squares approach At ever grid point a correlation tensor is computed i.e., correlation data is transformed into a distance dependent tensor model Reduces memor requirments because it avoids storing (global) correlation values eplicitl T is independent of the grid resolution, because correlation is set in relation to Euclidean distances Correlation visualiation on iso-surfaces Goal: Visualiing the structural uncertaint of specific features in the mean data, i.e. an iso-surface Requires correlation visualiation with respect to the mean surface s geometric shape Correlation visualiation on iso-surfaces Approach: Correlation values are etracted for the two principal correlation directions in the surface s tangent plane (highest and lowest strength) and normal direction [4] Basis transformation is necessar to obtain correlation strength parameters in principal tangential directions (τtan _ma) and (τtan _min), and normal direction (τnormal) at ever surface point B interactivel specifing an Euclidean radius d, correlation values along the three principal directions are computed as T Idea: Distinguish between correlations in the surface s local tangent planes and normal direction, and visualie these correlations via glphs and used to model the appearance of circular correlation glphs Iso-surface Correlation glphs Correlation glphs encode the correlation values in the first and second pricipal tangent and surface normal direction Correlation glphs Correlation ratio between first and second tangent direction is encoded in the shape of the elliptic one (2) Glph coloring wrt correlation in surface Glph coloring wrt correlation in surface (1) normal direction (2) first principal tangent direction (3) second principal tangent direction (1) normal direction (2) first principal tangent direction (3) second principal tangent direction 0 Correlation 1 0 Correlation 1 9
10 Correlation glphs Absolute and relative correlation anisotrop in the surface tangent plane are encoded via color differences Structural uncertaint analsis Surface realiation and glph based correlation visualiation on mean surface 0 Correlation 1 Surface realiation Correlation glphs Structural uncertaint analsis Structural uncertaint analsis High and low correlation in vertical and horiontal direction High and low correlation in horiontal and vertical direction Structural uncertaint analsis Structural uncertaint analsis Low correlation in both vertical and horiontal direction Low correlation especiall in normal direction 10
11 Structural uncertaint analsis Structural uncertaint analsis Mean surface in a simulated temperature ensemble Glph based correlation visualiation Structural uncertaint analsis Local correlation analsis summar Local distance dependent correlation tensor model allows memor reduction to O(n) for n grid points Visualiation of local correlation anisotrop possible Simultaneous visualiation of correlation ratios as well as absolute correlation values Region with high structural uncertaint Interactive analsis of anisotrop in surface tangent plane and between tangent plane and surface normal direction Local correlation analsis summar No integration of absolute uncertaint information (e.g. standard deviations) Part 2: Local and global correlation visualiation 2.2 Global correlation visualiation Global s are ignored correlation cannot be visualied as correlation model onl accounts for correlation strengths (magnitudes) 11
12 Cardinal number Global correlation analsis To anale long-range spatial dependencies between the data values in uncertain scalar fields, in principle the full covariance matri needs to be visualied This is problematic due to the following reasons For n spatial grid points it requires O n 2 memor Entries in Σ are not in a spatial conte, making it difficult to infer on spatial dependencies Uncertaint information has to be visuall separated from correlation information Σ Global correlation analsis Requirements due to the aforementioned problems The correlation information first has to be condensed Requires to define and seek for the most prominent short- and long-range dependencies Correlation structures have to be embedded into standard visualiations providing spatial relationships Approach Definition of a measure for the degree of dependenc of a random variable to its local and global spatial surrounding Spatial clustering of random variables based on this measure Color coding of data points wrt. cluster IDs Eample data set 2D temperature ensemble simulated b the European Centre for Medium-Range Weather Forecasts Correlation strength indicator Assumption: the more correlation partners a particular point in the domain has to which the correlation is larger than a given level, the more important this point is The sets of partners spread anisotropicall across the domain The epansion in different directions is directl related to the correlation distribution in the respective region Mean values Mean values as heightfield Standard deviations For a given level and point i, the number of partners indicates the degree of dependenc between the random variable Y( i ) to its local and global spatial surrounding It counts the most prominent partners of i, independent of their position in the domain The number of partners is used as the correlation strength indicator Correlation strength indicator Correlation strength indicator Definitions Correlation neighborhood: η ρ + i j D ρ Y i, Y j ρ + } D: data domain ρ 1 : correlation level Eample showing relation between correlation neighborhoods and cardinal numbers high Cardinal number: η ρ + i Correlation neighborhoods of 2 grid points Cardinal numbers low 12
13 Cardinal number Cardinal number Spatial correlation clustering Correlation clustering algorithm Basic idea: use correlation neighborhoods as clusters Problem: correlation neighborhoods of different points can overlap Results in ambiguities in the assignment of points to clusters Solution: process correlation neighborhoods sequentiall and eclude those neighborhoods which overlap an previousl processed neighborhood The strateg avoids overlapping clusters but, depending on selection order, will eclude large clusters which overlap small cluster Thus, neighborhoods are selected in descending order of cardinalit such that the largest clusters are alwas selected first Step 1: Select correlation level ρ + Step 2: Compute cardinal numbers for all points in the 2D domain Step 3: Select domain point with largest cardinal number as first cluster center Step 4: Create cluster (correlation neighborhood) high Mean values Cardinal numbers low Correlation clustering algorithm continued Correlation clustering results Step 5: In the set of points not et assigned to a cluster, search for point p with highest cardinal number whose correlation neighborhood c does not overlap an eisting cluster Step 6: Create cluster using p and c and repeat with Step 5 generate custers for different levels ρ + (multilevel clustering) high Visualiation of correlation clusters and cluster centers on the mean height surface for a particular correlation level ρ + Mean values Cardinal numbers low Correlation clustering results for varing levels dependencies stochastic dependencies tpicall eist between spatiall separated regions; ie. pairs which are inversel correlated to each other correlation structures are global features and cannot be modeled b distant dependent correlation models
14 dependencies stochastic dependencies tpicall eist between spatiall separated regions; ie. pairs which are inversel correlated to each other correlation structures are global features and cannot be modeled b distant dependent correlation models dependencies stochastic dependencies tpicall eist between spatiall separated regions; ie. pairs which are inversel correlated to each other correlation structures are global features and cannot be modeled b distant dependent correlation models dependencies stochastic dependencies tpicall eist between spatiall separated regions; ie. pairs which are inversel correlated to each other correlation structures are global features and cannot be modeled b distant dependent correlation models correlation clustering algorithm Step 1: Select negative correlation level ρ Step 2: Compute correlation neighborhoods (ρ Y i, Y j ρ ) and corresponding cardinal numbers for all points Step 3: Select point with largest cardinal number as first cluster center correlation clustering algorithm Step 1: Select negative correlation level ρ Step 2: Compute correlation neighborhoods (ρ Y i, Y j ρ ) and corresponding cardinal numbers for all points Step 3: Select point with largest cardinal number as first cluster center Step 4: Assign negativel correlated points to first cluster correlation clustering algorithm Step 1: Select negative correlation level ρ Step 2: Compute correlation neighborhoods (ρ Y i, Y j ρ ) and corresponding cardinal numbers for all points Step 3: Select point with largest cardinal number as first cluster center Step 4: Assign negativel correlated points to first cluster Step 5: Select point in cluster with largest number of inversel correlated points 14
15 correlation clustering algorithm Step 1: Select negative correlation level ρ Step 2: Compute correlation neighborhoods (ρ Y i, Y j ρ ) and corresponding cardinal numbers for all points Step 3: Select point with largest cardinal number as first cluster center Step 4: Assign negativel correlated points to first cluster Step 5: Select point in cluster with largest number of inversel correlated points correlation clustering algorithm Step 1: Select negative correlation level ρ Step 2: Compute correlation neighborhoods (ρ Y i, Y j ρ ) and corresponding cardinal numbers for all points Step 3: Select point with largest cardinal number as first cluster center Step 4: Assign negativel correlated points to first cluster Step 5: Select point in cluster with largest number of inversel correlated points correlation clustering algorithm Step 1: Select negative correlation level ρ Step 2: Compute correlation neighborhoods (ρ Y i, Y j ρ ) and corresponding cardinal numbers for all points Step 3: Select point with largest cardinal number as first cluster center Step 4: Assign negativel correlated points to first cluster Step 5: Select point in cluster with largest number of inversel correlated points Step 6: Assign both clusters to one pair correlation clustering algorithm Step 1: Select negative correlation level ρ Step 2: Compute correlation neighborhoods (ρ Y i, Y j ρ ) and corresponding cardinal numbers for all points Step 3: Select point with largest cardinal number as first cluster center Step 4: Assign negativel correlated points to first cluster Step 5: Select point in cluster with largest number of inversel correlated points Step 6: Assign both clusters to one pair Step 7: Proceed with Step 3, taking into account that clusters do not overlap dependencies Each pair of inversel correlated spatial regions gets assigned one distinct color Clusters in each pair are distinguished b differentl oriented stripe patterns correlation dependencies for varing negative levels ρ
16 Eample: An uncertain 2D scalar field with constant means and standard deviations. Specific orrelation structures have been enforced: Strong positive correlation within the sets of points colored blue. Zero correlation between these sets correlation between the two sets in the upper left region Mean surface Random realiation Random realiation Random realiation 16
17 Random realiation Random realiation Random realiation Random realiation Random realiation Random realiation 17
18 Correlation clustering algorithm identifies positivel correlated regions correctl Correlation clustering algorithm identifies inversel correlated regions correctl Validation Validation Global correlation analsis summar Correlation clustering allows analing global correlation structures such as inverse correlations, requiring an amount of memor that is linear in the number of data points The clusters distributions reveal anisotropic correlation effects Uncertaint information can be integrated easil, for instance, b etruding clusters depending on standard deviation [5] Etension to 3D is possible, but special projection or restriction schemes are required when used to anale correlation structures on iso-surfaces Future work and challenges in correlation visualiation Approaches for visualiing correlation structures in 3D Correlation analsis/visualiation for other data tpes (e.g. vector fields) Integration of correlation information in eisting uncertaint visualiation approaches (e.g. positional uncertaint of features) Quantification of the effect of correlation on the occurence of differential quantities and higher order features, like critical points Modeling and interpretation of stochastic dependencies in non- Gaussian distributed random fields 18
CONTINUOUS SPATIAL DATA ANALYSIS
CONTINUOUS SPATIAL DATA ANALSIS 1. Overview of Spatial Stochastic Processes The ke difference between continuous spatial data and point patterns is that there is now assumed to be a meaningful value, s
More informationINF Introduction to classifiction Anne Solberg Based on Chapter 2 ( ) in Duda and Hart: Pattern Classification
INF 4300 151014 Introduction to classifiction Anne Solberg anne@ifiuiono Based on Chapter 1-6 in Duda and Hart: Pattern Classification 151014 INF 4300 1 Introduction to classification One of the most challenging
More informationINF Introduction to classifiction Anne Solberg
INF 4300 8.09.17 Introduction to classifiction Anne Solberg anne@ifi.uio.no Introduction to classification Based on handout from Pattern Recognition b Theodoridis, available after the lecture INF 4300
More informationStrain Transformation and Rosette Gage Theory
Strain Transformation and Rosette Gage Theor It is often desired to measure the full state of strain on the surface of a part, that is to measure not onl the two etensional strains, and, but also the shear
More informationINF Anne Solberg One of the most challenging topics in image analysis is recognizing a specific object in an image.
INF 4300 700 Introduction to classifiction Anne Solberg anne@ifiuiono Based on Chapter -6 6inDuda and Hart: attern Classification 303 INF 4300 Introduction to classification One of the most challenging
More informationRandom Vectors. 1 Joint distribution of a random vector. 1 Joint distribution of a random vector
Random Vectors Joint distribution of a random vector Joint distributionof of a random vector Marginal and conditional distributions Previousl, we studied probabilit distributions of a random variable.
More information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
More informationMAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation
The Stress Equilibrium Equation As we mentioned in Chapter 2, using the Galerkin formulation and a choice of shape functions, we can derive a discretized form of most differential equations. In Structural
More information2: Distributions of Several Variables, Error Propagation
: Distributions of Several Variables, Error Propagation Distribution of several variables. variables The joint probabilit distribution function of two variables and can be genericall written f(, with the
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationScatter Plot Quadrants. Setting. Data pairs of two attributes X & Y, measured at N sampling units:
Geog 20C: Phaedon C Kriakidis Setting Data pairs of two attributes X & Y, measured at sampling units: ṇ and ṇ there are pairs of attribute values {( n, n ),,,} Scatter plot: graph of - versus -values in
More informationParameterized Joint Densities with Gaussian and Gaussian Mixture Marginals
Parameterized Joint Densities with Gaussian and Gaussian Miture Marginals Feli Sawo, Dietrich Brunn, and Uwe D. Hanebeck Intelligent Sensor-Actuator-Sstems Laborator Institute of Computer Science and Engineering
More informationHigher order method for non linear equations resolution: application to mobile robot control
Higher order method for non linear equations resolution: application to mobile robot control Aldo Balestrino and Lucia Pallottino Abstract In this paper a novel higher order method for the resolution of
More information6. Vector Random Variables
6. Vector Random Variables In the previous chapter we presented methods for dealing with two random variables. In this chapter we etend these methods to the case of n random variables in the following
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationStability Analysis of a Geometrically Imperfect Structure using a Random Field Model
Stabilit Analsis of a Geometricall Imperfect Structure using a Random Field Model JAN VALEŠ, ZDENĚK KALA Department of Structural Mechanics Brno Universit of Technolog, Facult of Civil Engineering Veveří
More informationMATRIX TRANSFORMATIONS
CHAPTER 5. MATRIX TRANSFORMATIONS INSTITIÚID TEICNEOLAÍOCHTA CHEATHARLACH INSTITUTE OF TECHNOLOGY CARLOW MATRIX TRANSFORMATIONS Matri Transformations Definition Let A and B be sets. A function f : A B
More informationFORCES AND MOMENTS IN CFD ANALYSIS
ORCES AND OENS IN CD ANALYSIS Authors Ing. Zdeněk Říha, PhD., Institute of Geonics ASCR v.v.i., Email: denek.riha@ugn.cas.c Ing. Josef oldna, CSc., Institute of Geonics ASCR v.v.i. Anotace V příspěvku
More informationComputer Graphics: 2D Transformations. Course Website:
Computer Graphics: D Transformations Course Website: http://www.comp.dit.ie/bmacnamee 5 Contents Wh transformations Transformations Translation Scaling Rotation Homogeneous coordinates Matri multiplications
More informationReview of Probability
Review of robabilit robabilit Theor: Man techniques in speech processing require the manipulation of probabilities and statistics. The two principal application areas we will encounter are: Statistical
More informationCH.7. PLANE LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics
CH.7. PLANE LINEAR ELASTICITY Multimedia Course on Continuum Mechanics Overview Plane Linear Elasticit Theor Plane Stress Simplifing Hpothesis Strain Field Constitutive Equation Displacement Field The
More informationResearch Design - - Topic 15a Introduction to Multivariate Analyses 2009 R.C. Gardner, Ph.D.
Research Design - - Topic 15a Introduction to Multivariate Analses 009 R.C. Gardner, Ph.D. Major Characteristics of Multivariate Procedures Overview of Multivariate Techniques Bivariate Regression and
More informationPolynomial and Rational Functions
Polnomial and Rational Functions Figure -mm film, once the standard for capturing photographic images, has been made largel obsolete b digital photograph. (credit film : modification of ork b Horia Varlan;
More informationBiostatistics in Research Practice - Regression I
Biostatistics in Research Practice - Regression I Simon Crouch 30th Januar 2007 In scientific studies, we often wish to model the relationships between observed variables over a sample of different subjects.
More informationPerturbation Theory for Variational Inference
Perturbation heor for Variational Inference Manfred Opper U Berlin Marco Fraccaro echnical Universit of Denmark Ulrich Paquet Apple Ale Susemihl U Berlin Ole Winther echnical Universit of Denmark Abstract
More informationME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites
ME 75 Lecture Effective Properties of Particulate and Unidirectional Composites Concepts from Elasticit Theor Statistical Homogeneit, Representative Volume Element, Composite Material Effective Stress-
More informationIntroduction to Differential Equations. National Chiao Tung University Chun-Jen Tsai 9/14/2011
Introduction to Differential Equations National Chiao Tung Universit Chun-Jen Tsai 9/14/011 Differential Equations Definition: An equation containing the derivatives of one or more dependent variables,
More informationChapter 3. Theory of measurement
Chapter. Introduction An energetic He + -ion beam is incident on thermal sodium atoms. Figure. shows the configuration in which the interaction one is determined b the crossing of the laser-, sodium- and
More informationHamiltonicity and Fault Tolerance
Hamiltonicit and Fault Tolerance in the k-ar n-cube B Clifford R. Haithcock Portland State Universit Department of Mathematics and Statistics 006 In partial fulfillment of the requirements of the degree
More informationCopyright, 2008, R.E. Kass, E.N. Brown, and U. Eden REPRODUCTION OR CIRCULATION REQUIRES PERMISSION OF THE AUTHORS
Copright, 8, RE Kass, EN Brown, and U Eden REPRODUCTION OR CIRCULATION REQUIRES PERMISSION OF THE AUTHORS Chapter 6 Random Vectors and Multivariate Distributions 6 Random Vectors In Section?? we etended
More informationChapter 9 BIAXIAL SHEARING
9. DEFNTON Chapter 9 BAXAL SHEARNG As we have seen in the previous chapter, biaial (oblique) shearing produced b the shear forces and, appears in a bar onl accompanied b biaial bending (we ma discuss about
More informationBifurcations of the Controlled Escape Equation
Bifurcations of the Controlled Escape Equation Tobias Gaer Institut für Mathematik, Universität Augsburg 86135 Augsburg, German gaer@math.uni-augsburg.de Abstract In this paper we present numerical methods
More informationMath 214 Spring problem set (a) Consider these two first order equations. (I) dy dx = x + 1 dy
Math 4 Spring 08 problem set. (a) Consider these two first order equations. (I) d d = + d (II) d = Below are four direction fields. Match the differential equations above to their direction fields. Provide
More information* τσ σκ. Supporting Text. A. Stability Analysis of System 2
Supporting Tet A. Stabilit Analsis of Sstem In this Appendi, we stud the stabilit of the equilibria of sstem. If we redefine the sstem as, T when -, then there are at most three equilibria: E,, E κ -,,
More informationwe must pay attention to the role of the coordinate system w.r.t. which we perform a tform
linear SO... we will want to represent the geometr of points in space we will often want to perform (rigid) transformations to these objects to position them translate rotate or move them in an animation
More informationImaging Metrics. Frequency response Coherent systems Incoherent systems MTF OTF Strehl ratio Other Zemax Metrics. ECE 5616 Curtis
Imaging Metrics Frequenc response Coherent sstems Incoherent sstems MTF OTF Strehl ratio Other Zema Metrics Where we are going with this Use linear sstems concept of transfer function to characterize sstem
More informationE(x i ) = µ i. 2 d. + sin 1 d θ 2. for d < θ 2 0 for d θ 2
1 Gaussian Processes Definition 1.1 A Gaussian process { i } over sites i is defined by its mean function and its covariance function E( i ) = µ i c ij = Cov( i, j ) plus joint normality of the finite
More informationCONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM
Chapter 4 CONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM Figure 4.1: 4.1 Conservation of Angular Momentum Angular momentum is defined as the moment of the linear momentum about some spatial reference
More informationThe first change comes in how we associate operators with classical observables. In one dimension, we had. p p ˆ
VI. Angular momentum Up to this point, we have been dealing primaril with one dimensional sstems. In practice, of course, most of the sstems we deal with live in three dimensions and 1D quantum mechanics
More information1.6 ELECTRONIC STRUCTURE OF THE HYDROGEN ATOM
1.6 ELECTRONIC STRUCTURE OF THE HYDROGEN ATOM 23 How does this wave-particle dualit require us to alter our thinking about the electron? In our everda lives, we re accustomed to a deterministic world.
More information1 GSW Gaussian Elimination
Gaussian elimination is probabl the simplest technique for solving a set of simultaneous linear equations, such as: = A x + A x + A x +... + A x,,,, n n = A x + A x + A x +... + A x,,,, n n... m = Am,x
More informationCovariance Tracking Algorithm on Bilateral Filtering under Lie Group Structure Yinghong Xie 1,2,a Chengdong Wu 1,b
Applied Mechanics and Materials Online: 014-0-06 ISSN: 166-748, Vols. 519-50, pp 684-688 doi:10.408/www.scientific.net/amm.519-50.684 014 Trans Tech Publications, Switzerland Covariance Tracking Algorithm
More informationInference about the Slope and Intercept
Inference about the Slope and Intercept Recall, we have established that the least square estimates and 0 are linear combinations of the Y i s. Further, we have showed that the are unbiased and have the
More informationAnd similarly in the other directions, so the overall result is expressed compactly as,
SQEP Tutorial Session 5: T7S0 Relates to Knowledge & Skills.5,.8 Last Update: //3 Force on an element of area; Definition of principal stresses and strains; Definition of Tresca and Mises equivalent stresses;
More informationCH. 1 FUNDAMENTAL PRINCIPLES OF MECHANICS
446.201 (Solid echanics) Professor Youn, eng Dong CH. 1 FUNDENTL PRINCIPLES OF ECHNICS Ch. 1 Fundamental Principles of echanics 1 / 14 446.201 (Solid echanics) Professor Youn, eng Dong 1.2 Generalied Procedure
More informationFrom the help desk: It s all about the sampling
The Stata Journal (2002) 2, Number 2, pp. 90 20 From the help desk: It s all about the sampling Allen McDowell Stata Corporation amcdowell@stata.com Jeff Pitblado Stata Corporation jsp@stata.com Abstract.
More informationPolynomial and Rational Functions
Polnomial and Rational Functions Figure -mm film, once the standard for capturing photographic images, has been made largel obsolete b digital photograph. (credit film : modification of ork b Horia Varlan;
More informationExercise solutions: concepts from chapter 5
1) Stud the oöids depicted in Figure 1a and 1b. a) Assume that the thin sections of Figure 1 lie in a principal plane of the deformation. Measure and record the lengths and orientations of the principal
More informationVector Calculus Review
Course Instructor Dr. Ramond C. Rumpf Office: A-337 Phone: (915) 747-6958 E-Mail: rcrumpf@utep.edu Vector Calculus Review EE3321 Electromagnetic Field Theor Outline Mathematical Preliminaries Phasors,
More informationLECTURE 14 Strength of a Bar in Transverse Bending. 1 Introduction. As we have seen, only normal stresses occur at cross sections of a rod in pure
V. DEMENKO MECHNCS OF MTERLS 015 1 LECTURE 14 Strength of a Bar in Transverse Bending 1 ntroduction s we have seen, onl normal stresses occur at cross sections of a rod in pure bending. The corresponding
More informationwe must pay attention to the role of the coordinate system w.r.t. which we perform a tform
linear SO... we will want to represent the geometr of points in space we will often want to perform (rigid) transformations to these objects to position them translate rotate or move them in an animation
More information520 Chapter 9. Nonlinear Differential Equations and Stability. dt =
5 Chapter 9. Nonlinear Differential Equations and Stabilit dt L dθ. g cos θ cos α Wh was the negative square root chosen in the last equation? (b) If T is the natural period of oscillation, derive the
More informationDemonstrate solution methods for systems of linear equations. Show that a system of equations can be represented in matrix-vector form.
Chapter Linear lgebra Objective Demonstrate solution methods for sstems of linear equations. Show that a sstem of equations can be represented in matri-vector form. 4 Flowrates in kmol/hr Figure.: Two
More informationOn the Stability of a Differential-Difference Analogue of a Two-Dimensional Problem of Integral Geometry
Filomat 3:3 18, 933 938 https://doi.org/1.98/fil183933b Published b Facult of Sciences and Mathematics, Universit of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat On the Stabilit of a Differential-Difference
More informationUnit 12 Study Notes 1 Systems of Equations
You should learn to: Unit Stud Notes Sstems of Equations. Solve sstems of equations b substitution.. Solve sstems of equations b graphing (calculator). 3. Solve sstems of equations b elimination. 4. Solve
More informationAn Information Theory For Preferences
An Information Theor For Preferences Ali E. Abbas Department of Management Science and Engineering, Stanford Universit, Stanford, Ca, 94305 Abstract. Recent literature in the last Maimum Entrop workshop
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
More informationChapter 4 Analytic Trigonometry
Analtic Trigonometr Chapter Analtic Trigonometr Inverse Trigonometric Functions The trigonometric functions act as an operator on the variable (angle, resulting in an output value Suppose this process
More informationMACHINE LEARNING 2012 MACHINE LEARNING
1 MACHINE LEARNING kernel CCA, kernel Kmeans Spectral Clustering 2 Change in timetable: We have practical session net week! 3 Structure of toda s and net week s class 1) Briefl go through some etension
More informationComputer Vision & Digital Image Processing
Computer Vision & Digital Image Processing Image Segmentation Dr. D. J. Jackson Lecture 6- Image segmentation Segmentation divides an image into its constituent parts or objects Level of subdivision depends
More informationKey words: Analysis of dependence, chi-plot, confidence intervals.
Journal of Data Science 10(2012), 711-722 The Chi-plot and Its Asmptotic Confidence Interval for Analzing Bivariate Dependence: An Application to the Average Intelligence and Atheism Rates across Nations
More informationMetric-based classifiers. Nuno Vasconcelos UCSD
Metric-based classifiers Nuno Vasconcelos UCSD Statistical learning goal: given a function f. y f and a collection of eample data-points, learn what the function f. is. this is called training. two major
More informationFATIGUE ANALYSIS: THE SUPER-NEUBER TECHNIQUE FOR CORRECTION OF LINEAR ELASTIC FE RESULTS
6 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES FATIGUE ANALYSIS: THE SUPER-NEUBER TECHNIQUE FOR CORRECTION OF LINEAR ELASTIC FE RESULTS L. Samuelsson * * Volvo Aero Corporation, Sweden Kewords:
More informationOn Maximizing the Second Smallest Eigenvalue of a State-dependent Graph Laplacian
25 American Control Conference June 8-, 25. Portland, OR, USA WeA3.6 On Maimizing the Second Smallest Eigenvalue of a State-dependent Graph Laplacian Yoonsoo Kim and Mehran Mesbahi Abstract We consider
More informationLINEAR REGRESSION ANALYSIS
LINEAR REGRESSION ANALYSIS MODULE V Lecture - 2 Correcting Model Inadequacies Through Transformation and Weighting Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technolog Kanpur
More informationLaurie s Notes. Overview of Section 3.5
Overview of Section.5 Introduction Sstems of linear equations were solved in Algebra using substitution, elimination, and graphing. These same techniques are applied to nonlinear sstems in this lesson.
More informationReading. 4. Affine transformations. Required: Watt, Section 1.1. Further reading:
Reading Required: Watt, Section.. Further reading: 4. Affine transformations Fole, et al, Chapter 5.-5.5. David F. Rogers and J. Alan Adams, Mathematical Elements for Computer Graphics, 2 nd Ed., McGraw-Hill,
More information5.3.3 The general solution for plane waves incident on a layered halfspace. The general solution to the Helmholz equation in rectangular coordinates
5.3.3 The general solution for plane waves incident on a laered halfspace The general solution to the elmhol equation in rectangular coordinates The vector propagation constant Vector relationships between
More informationNotes 7 Analytic Continuation
ECE 6382 Fall 27 David R. Jackson Notes 7 Analtic Continuation Notes are from D. R. Wilton, Dept. of ECE Analtic Continuation of Functions We define analtic continuation as the process of continuing a
More informationSection 3.1. ; X = (0, 1]. (i) f : R R R, f (x, y) = x y
Paul J. Bruillard MATH 0.970 Problem Set 6 An Introduction to Abstract Mathematics R. Bond and W. Keane Section 3.1: 3b,c,e,i, 4bd, 6, 9, 15, 16, 18c,e, 19a, 0, 1b Section 3.: 1f,i, e, 6, 1e,f,h, 13e,
More information2.4 Orthogonal Coordinate Systems (pp.16-33)
8/26/2004 sec 2_4 blank.doc 1/6 2.4 Orthogonal Coordinate Sstems (pp.16-33) 1) 2) Q: A: 1. 2. 3. Definition: ). 8/26/2004 sec 2_4 blank.doc 2/6 A. Coordinates * * * Point P(0,0,0) is alwas the origin.
More informationMobile Robot Localization
Mobile Robot Localization 1 The Problem of Robot Localization Given a map of the environment, how can a robot determine its pose (planar coordinates + orientation)? Two sources of uncertainty: - observations
More informationManeuvering Target Tracking Method based on Unknown but Bounded Uncertainties
Preprints of the 8th IFAC World Congress Milano (Ital) ust 8 - Septemer, Maneuvering arget racing Method ased on Unnown ut Bounded Uncertainties Hodjat Rahmati*, Hamid Khaloozadeh**, Moosa Aati*** * Facult
More informationCamera calibration. Outline. Pinhole camera. Camera projection models. Nonlinear least square methods A camera calibration tool
Outline Camera calibration Camera projection models Camera calibration i Nonlinear least square methods A camera calibration tool Applications Digital Visual Effects Yung-Yu Chuang with slides b Richard
More informationFigure 1: Visualising the input features. Figure 2: Visualising the input-output data pairs.
Regression. Data visualisation (a) To plot the data in x trn and x tst, we could use MATLAB function scatter. For example, the code for plotting x trn is scatter(x_trn(:,), x_trn(:,));........... x (a)
More informationExploiting Correlation in Stochastic Circuit Design
Eploiting Correlation in Stochastic Circuit Design Armin Alaghi and John P. Haes Advanced Computer Architecture Laborator Department of Electrical Engineering and Computer Science Universit of Michigan,
More informationf x, y x 2 y 2 2x 6y 14. Then
SECTION 11.7 MAXIMUM AND MINIMUM VALUES 645 absolute minimum FIGURE 1 local maimum local minimum absolute maimum Look at the hills and valles in the graph of f shown in Figure 1. There are two points a,
More informationEVALUATION OF STRESS IN BMI-CARBON FIBER LAMINATE TO DETERMINE THE ONSET OF MICROCRACKING
EVALUATION OF STRESS IN BMI-CARBON FIBER LAMINATE TO DETERMINE THE ONSET OF MICROCRACKING A Thesis b BRENT DURRELL PICKLE Submitted to the Office of Graduate Studies of Teas A&M Universit in partial fulfillment
More information3.7 InveRSe FUnCTIOnS
CHAPTER functions learning ObjeCTIveS In this section, ou will: Verif inverse functions. Determine the domain and range of an inverse function, and restrict the domain of a function to make it one-to-one.
More informationProbability Densities in Data Mining
Probabilit Densities in Data Mining Note to other teachers and users of these slides. Andrew would be delighted if ou found this source material useful in giving our own lectures. Feel free to use these
More informationMultiscale Principal Components Analysis for Image Local Orientation Estimation
Multiscale Principal Components Analsis for Image Local Orientation Estimation XiaoGuang Feng Peman Milanfar Computer Engineering Electrical Engineering Universit of California, Santa Cruz Universit of
More informationParticular Solutions
Particular Solutions Our eamples so far in this section have involved some constant of integration, K. We now move on to see particular solutions, where we know some boundar conditions and we substitute
More informationLévy stable distribution and [0,2] power law dependence of. acoustic absorption on frequency
Lév stable distribution and [,] power law dependence of acoustic absorption on frequenc W. Chen Institute of Applied Phsics and Computational Mathematics, P.O. Box 89, Division Box 6, Beijing 88, China
More informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationDecision tables and decision spaces
Abstract Decision tables and decision spaces Z. Pawlak 1 Abstract. In this paper an Euclidean space, called a decision space is associated with ever decision table. This can be viewed as a generalization
More informationTENSOR TRANSFORMATION OF STRESSES
GG303 Lecture 18 9/4/01 1 TENSOR TRANSFORMATION OF STRESSES Transformation of stresses between planes of arbitrar orientation In the 2-D eample of lecture 16, the normal and shear stresses (tractions)
More informationSample-based Optimal Transport and Barycenter Problems
Sample-based Optimal Transport and Barcenter Problems MAX UANG New York Universit, Courant Institute of Mathematical Sciences AND ESTEBAN G. TABA New York Universit, Courant Institute of Mathematical Sciences
More informationRoberto s Notes on Integral Calculus Chapter 3: Basics of differential equations Section 3. Separable ODE s
Roberto s Notes on Integral Calculus Chapter 3: Basics of differential equations Section 3 Separable ODE s What ou need to know alread: What an ODE is and how to solve an eponential ODE. What ou can learn
More informationChapter 13 TORSION OF THIN-WALLED BARS WHICH HAVE THE CROSS SECTIONS PREVENTED FROM WARPING (Prevented or non-uniform torsion)
Chapter 13 TORSION OF THIN-WALLED BARS WHICH HAVE THE CROSS SECTIONS PREVENTED FROM WARPING (Prevented or non-uniform torsion) 13.1 GENERALS In our previous chapter named Pure (uniform) Torsion, it was
More informationx = 1 n (x 1 + x 2 + +x n )
3 PARTIAL DERIVATIVES Science Photo Librar Three-dimensional surfaces have high points and low points that are analogous to the peaks and valles of a mountain range. In this chapter we will use derivatives
More information1 Differential Equations for Solid Mechanics
1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to
More informationStructure-integrated Active Damping System: Integral Strain-based Design Strategy for the Optimal Placement of Functional Elements
International Journal of Composite Materials 213, 3(6B): 53-58 DOI: 1.5923/s.cmaterials.2131.6 Structure-egrated Active Damping Sstem: Integral Strain-based Design Strateg for the Optimal Placement of
More informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
More informationSchool of Computer and Communication Sciences. Information Theory and Coding Notes on Random Coding December 12, 2003.
ÉCOLE POLYTECHNIQUE FÉDÉRALE DE LAUSANNE School of Computer and Communication Sciences Handout 8 Information Theor and Coding Notes on Random Coding December 2, 2003 Random Coding In this note we prove
More informationNon-parametric estimation of geometric anisotropy from environmental sensor network measurements
Non-parametric estimation of geometric anisotrop from environmental sensor network measurements Dionissios T. Hristopulos 1, M. P. Petrakis 1, G. Spiliopoulos 1, Arsenia Chorti 2 1 Technical Universit
More informationKernel PCA, clustering and canonical correlation analysis
ernel PCA, clustering and canonical correlation analsis Le Song Machine Learning II: Advanced opics CSE 8803ML, Spring 2012 Support Vector Machines (SVM) 1 min w 2 w w + C j ξ j s. t. w j + b j 1 ξ j,
More informationInterspecific Segregation and Phase Transition in a Lattice Ecosystem with Intraspecific Competition
Interspecific Segregation and Phase Transition in a Lattice Ecosstem with Intraspecific Competition K. Tainaka a, M. Kushida a, Y. Ito a and J. Yoshimura a,b,c a Department of Sstems Engineering, Shizuoka
More informationy R T However, the calculations are easier, when carried out using the polar set of co-ordinates ϕ,r. The relations between the co-ordinates are:
Curved beams. Introduction Curved beams also called arches were invented about ears ago. he purpose was to form such a structure that would transfer loads, mainl the dead weight, to the ground b the elements
More informationGreen s Theorem Jeremy Orloff
Green s Theorem Jerem Orloff Line integrals and Green s theorem. Vector Fields Vector notation. In 8.4 we will mostl use the notation (v) = (a, b) for vectors. The other common notation (v) = ai + bj runs
More information