Using a Projected Subgradient Method to Solve a Constrained Optimization Problem for Separating an Arbitrary Set of Points into Uniform Segments
|
|
- Bertram Barnett
- 5 years ago
- Views:
Transcription
1 Using a Projected Subgradient Method to Solve a Constrained Optimization Problem or Separating an Arbitrary Set o Points into Uniorm Segments Michael Johnson May 31, Background Inormation The Airborne Topographic Mapper [Kra09] is a circular scanning LIDAR that measures surace height along a track as the aircrat lies overhead. We desire to run analysis on how the topography is changing along the track, using bins o constant along-track length to do so. As the data is stored in long vectors o latitude, longitude, and height sorted only in time, the circular pattern in which the data is recorded does not allow simple sorting methods to be used to bin the data. We propose to use optimization methods to locate the points on the boundary o the track that yield the bins o constant length. Figure 1 illustrates a sample track o data being considered, which was taken over sea ice above Northern Greenland in Please reer to the project proposal and progress report or more background inormation on the problem and data. 2 Model Figure 2 illustrates a simpliied instance o this problem which assumes the track boundaries, B 1 and B 2, are straight, parallel lines. A track bin is deined as the area enclosed by a starting line segment connecting B 1 and B 2 and an ending line segment connecting B 1 and B 2. The segment length is deined as the distance between the center o the starting line segment, x s, and the center o the ending line segment, x. The points x 1 and x 2 are known, and we wish to ind the points x 3, x 4, and x (which are all related) that enclose the track bin. The circle o radius α deines the possible points or the center o the next line segment that is a distance o α away rom our starting point. In this problem instance, the obvious choice o x 3, x 4, and x are shown, as they are the points that create the most uniorm segment o length α. The length d o our optimal ending line segment has the minimum 1
2 83.7 Sample Track o ATM Measurements Latitude Longitude Figure 1: Sample track o ATM measurements length o any line segment ormed by the intersection o the boundaries with a line tangent to our circle. Thereore, solving this problem involves minimizing the length o this line segment. Figure 2: Simpliied Instance This simpliied instance will be used to solve this problem or the actual data. The nonconvex nature o the boundaries will be addressed by using an aine approximation which is updated as the solution progresses. A mathematical model o this problem is: minimize x 3 x 4 2 subject to x 3 = A 1 1 b 1 x 4 = A 1 2 b 2 x x s 2 = α The equations or x 3 and x 4 are ormed as the solution to the problem o inding the 2 (1)
3 intersection between the boundaries and the line tangent to the circle at point x. For simplicities sake, we will consider the point x s to be the origin o our rame so can be ignored in uture equations. Likewise, we will use x = [x y] T. In solving or the intersection o the two lines, we get: [ ] [ ] x y x T A 1 =, b m B1 1 1 = x (2) b B1 In equation (2), m B1 reers to the slope o our aine approximation to the track boundaries, and b B1 reers to the intercept. Because A and b are dependent on the point x, there is not a simple closed orm solution to this problem. The problem can be progressed by inputting the deinitions or x 3 and x 4 into the objective unction and splitting the vector into its constituents to give us F (x, y). The gradient can then be ormed as F (x, y) = [δf/δx δf/δy] T. Although in general F (x, y) is a non-convex unction, it can be shown to be convex when x is constrained to lie on the circle. The derivation or F (x, y) using (1) and (2) is as ollows. Using the deinitions or x 3 and x 4, the objective unction becomes: F (x, y) = x 3 x 4 2 = A 1 1 b 1 A 1 2 b 2 2 (3) Using the deinitions o A and b and x = [x y] T gives us an equation o the orm: Where: F (x, y) = 1 (x, y) + g(x, y) (4) h(x, y) h(x, y) = (x + m 1 y)(x + m 2 y) (x, y) = (x 2 + y 2 b 1 y)(x + m 2 y) (x 2 + y 2 b 2 y)(x + m 1 y) g(x, y) = (m 1 x 2 + b 1 x + m 1 y 2 )(x + m 2 y) (m 2 x 2 + b 2 x + m 2 y 2 )(x + m 1 y) Using these equations, we similarly get expressions or δf/δx and δf/δy. Dropping the (x, y) notation or simplicity, we have: δf/δx = (1/h 2 )(h[2(δ/δx) + 2g(δg/δx)] (δh/δx)[ 2 + g 2 ] δf/δy = (1/h 2 )(h[2(δ/δy) + 2g(δg/δy)] (δh/δy)[ 2 + g 2 ] Using the equations in (5), we get: δh/δx = 2x + (m 1 + m 2 )y δh/δy = 2m 1 m 2 y + (m 1 + m 2 )x δ/δx = 2x(x + m 2 y) + (x 2 + y 2 b 1 y) 2x(x + m 1 y) + (x 2 + y 2 b 2 y) δ/δy = (2y b 1 )(x + m 2 y) + m 2 (x 2 + y 2 b 1 y) (2y b 2 )(x + m 1 y) + m 1 (x 2 + y 2 b 2 y) δg/δx = (2m 1 x + b 1 )(x + m 2 y) + (m 1 x 2 + b 1 x + m 1 y 2 )... (2m 2 x + b 2 )(x + m 1 y) + (m 2 x 2 + b 2 x + m 2 y 2 ) δg/δy = 2m 1 y(x + m 2 y) + m 2 (m 1 x 2 + b 1 x + m 1 y 2 )... 2m 2 y(x + m 1 y) + m 1 (m 2 x 2 + b 2 x + m 2 y 2 ) (7) 3 (5) (6)
4 The equations presented in (5), (6), and (7) allow us to orm the gradient o the objective unction. Figure 3: Unconstrained objective unction Figure 3 shows a plot o the unconstrained objective unction. The large increase in value is attributed to when the line tangent to the circle is parallel to the boundaries. The unction is convex on each halplane and using the gradient will lower the objective value and ensures that x stays on the proper side o x s ; we do not want the segment to go to a section o the track that is already segmented. 3 Algorithm It is proposed to solve this problem using a projected subgradient method [Boy11] because o the constraint that the solution lie on the circle centered around x s. Do while not converged: Generate linear approximation o boundaries Do while not converged: Calculate F (x (k) ), the gradient at F (x(k) Take x (k+1) = P(x (k) γ k F (x (k) )) ) where P is the projection onto the circle centered around x s, γ k is the step length (could be 1/k) Update linear approximation o boundaries based on solution 4
5 As The objective is dierentiable, the subgradient is the gradient o the unction which simpliies the procedure. The constraint that x lie on the circle centered at x s with radius α requires the projection, as the gradient o the unconstrained objective unction will not keep the value o x (k) constrained to the circle. Projection onto the circle can be done in the ollowing manner. Using x s = [0 0] T, we calculate θ = tan 1 (y (k) /x (k) ) to determine the angle that the current value o x (k) makes with x s. Projecting the point onto the circle is then: [ ] x (k) αcosθ = (8) αsinθ Equation (8) will give the point on the circle that is closest to the x calculated using the gradient o the objective. 4 Results and Conclusion The method and algorithm described in this paper was applied to the set o data shown in igure 1. This dataset was chosen as it has multiple traits that are to be tested. The track contains near-straight sections, curved sections, kinks, sections where the track width narrows, and sections o missing measurements. The latter traits introduce non-uniormity to the boundaries which will test the robustness o the model and algorithm Figure 4: The entire track segmented. 5
6 (a) Detail (b) Detail (c) Detail (d) Detail 4 Figure 5: Details o segmented track Figure 4 shows the entire track ater segmentation to show the method worked on the entire track. The details in igure 5 show close-ups o several segments o the track to illustrate the perormance o the algorithm. In the igures, the dashed lines connect x 3 and x 4 that are calculated or each segment, which are shown as stars along the track boundary. Figure 5a shows a relatively straight section o track ater segmentation. Figure 5b shows the portion o the track where the boundary narrows due to bad measurements. The algorithm is able to successully segment and pass through this region. Likewise, igure 5c shows the portion o the track that has a large gap in both boundaries, again due to bad measurements. Even though the boundaries are missing, the algorithm makes a visually good approximation to what the segments may be like. Figure 5d shows the area o the track with modest kinks. The algorithm segments this region without problem, and the pattern o the segments at the kinks visually makes: adjacent boundary points are closer together on the inside o the kink and arther apart on the outside o the kink. The results presented in this section suggest the model and algorithm described is adequate at solving this problem. Ater careul selection o problem parameters, the method is shown to be robust, eicient, and accurate. 6
7 Reerences [Boy11] S. Boyd. Subgradient Methods or Constrained Problems, EE364B Lecture Notes. Stanord Univeristy, [Kra09] William B Krabill. IceBridge ATM L1B Qit Elevation and Return Strength. Boulder, Colorado USA: National Snow and Ice Data Center. Digital Media,
Reflection and Refraction
Relection and Reraction Object To determine ocal lengths o lenses and mirrors and to determine the index o reraction o glass. Apparatus Lenses, optical bench, mirrors, light source, screen, plastic or
More information5.2 Properties of Rational functions
5. Properties o Rational unctions A rational unction is a unction o the orm n n1 polynomial p an an 1 a1 a0 k k1 polynomial q bk bk 1 b1 b0 Eample 3 5 1 The domain o a rational unction is the set o all
More informationLarger K-maps. So far we have only discussed 2 and 3-variable K-maps. We can now create a 4-variable map in the
EET 3 Chapter 3 7/3/2 PAGE - 23 Larger K-maps The -variable K-map So ar we have only discussed 2 and 3-variable K-maps. We can now create a -variable map in the same way that we created the 3-variable
More information9.8 Graphing Rational Functions
9. Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm P where P and Q are polynomials. Q An eample o a simple rational unction
More informationL ENSES. Lenses Spherical refracting surfaces. n 1 n 2
Lenses 2 L ENSES 2. Sherical reracting suraces In order to start discussing lenses uantitatively, it is useul to consider a simle sherical surace, as shown in Fig. 2.. Our lens is a semi-ininte rod with
More informationTrigonometric Functions. Copyright Cengage Learning. All rights reserved.
4 Trigonometric Functions Copyright Cengage Learning. All rights reserved. 4.7 Inverse Trigonometric Functions Copyright Cengage Learning. All rights reserved. What You Should Learn Evaluate and graph
More informationSystems of Equations and Inequalities. Copyright Cengage Learning. All rights reserved.
5 Systems of Equations and Inequalities Copyright Cengage Learning. All rights reserved. 5.5 Systems of Inequalities Copyright Cengage Learning. All rights reserved. Objectives Graphing an Inequality Systems
More informationPATH PLANNING OF UNMANNED AERIAL VEHICLE USING DUBINS GEOMETRY WITH AN OBSTACLE
Proceeding o International Conerence On Research, Implementation And Education O Mathematics And Sciences 2015, Yogyakarta State University, 17-19 May 2015 PATH PLANNING OF UNMANNED AERIAL VEHICLE USING
More informationPiecewise polynomial interpolation
Chapter 2 Piecewise polynomial interpolation In ection.6., and in Lab, we learned that it is not a good idea to interpolate unctions by a highorder polynomials at equally spaced points. However, it transpires
More informationConcavity. Notice the location of the tangents to each type of curve.
Concavity We ve seen how knowing where a unction is increasing and decreasing gives a us a good sense o the shape o its graph We can reine that sense o shape by determining which way the unction bends
More informationSYSTEMS OF NONLINEAR EQUATIONS
SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two
More informationRational Functions. Definition A rational function can be written in the form. where N(x) and D(x) are
Rational Functions Deinition A rational unction can be written in the orm () N() where N() and D() are D() polynomials and D() is not the zero polynomial. *To ind the domain o a rational unction we must
More informationQuiz 6 Practice Problems
Quiz 6 Practice Problems Practice problems are similar, both in difficulty and in scope, to the type of problems you will see on the quiz. Problems marked with a are for your entertainment and are not
More informationDirectional Derivatives. Directional Derivatives. Directional Derivatives. Directional Derivatives. Directional Derivatives. Directional Derivatives
Recall that if z = f(x, y), then the partial derivatives f x and f y are defined as and represent the rates of change of z in the x- and y-directions, that is, in the directions of the unit vectors i and
More informationGeology and Calculus
GEOL 452 - Mathematical Tools in Geology Lab Assignment # 6 - Feb 25, 2010 (Due March 9, 2010) Name: Geology and Calculus A. Volume of San Nicolas Island San Nicolas Island is one of the remote and smaller
More informationGrad operator, triple and line integrals. Notice: this material must not be used as a substitute for attending the lectures
Grad operator, triple and line integrals Notice: this material must not be used as a substitute for attending the lectures 1 .1 The grad operator Let f(x 1, x,..., x n ) be a function of the n variables
More informationScalable Test Problems for Evolutionary Multi-Objective Optimization
Scalable Test Problems or Evolutionary Multi-Objective Optimization Kalyanmoy Deb Kanpur Genetic Algorithms Laboratory Indian Institute o Technology Kanpur PIN 8 6, India deb@iitk.ac.in Lothar Thiele,
More information6, 1 0, f x x 1 2 x h x x x 3, f x sin x cos x, f x x 2 6x 5 f x 4x 3 5x 30. g x x3 8x 31. f x x f x x2 3x 4 33.
Chapter Applications o Dierentiation Review Eercises See CalcChat.com or tutorial help and worked-out solutions to odd-numbered eercises. Finding Etrema on a Closed Interval In Eercises, ind the absolute
More informationComputer Graphics : Bresenham Line Drawing Algorithm, Circle Drawing & Polygon Filling
Computer Graphics : Bresenham Line Drawing Algorithm, Circle Drawing & Polygon Filling Downloaded from :www.comp.dit.ie/bmacnamee/materials/graphics/006- Contents In today s lecture we ll have a loo at:
More informationSkill Sets Chapter 5 Functions
Skill Sets Chapter 5 Functions No. Skills Examples o questions involving the skills. Sketch the graph o the (Lecture Notes Example (b)) unction according to the g : x x x, domain. x, x - Students tend
More informationIn today s lecture we ll have a look at: A simple technique The mid-point circle algorithm
Drawing Circles In today s lecture we ll have a look at: Circle drawing algorithms A simple technique The mid-point circle algorithm Polygon fill algorithms Summary raster drawing algorithms A Simple Circle
More informationJim Lambers MAT 169 Fall Semester Lecture 33 Notes
Jim Lambers MAT 169 Fall Semester 2009-10 Lecture 33 Notes These notes correspond to Section 9.3 in the text. Polar Coordinates Throughout this course, we have denoted a point in the plane by an ordered
More informationChapter 3 Image Enhancement in the Spatial Domain
Chapter 3 Image Enhancement in the Spatial Domain Yinghua He School o Computer Science and Technology Tianjin University Image enhancement approaches Spatial domain image plane itsel Spatial domain methods
More information13.1. Functions of Several Variables. Introduction to Functions of Several Variables. Functions of Several Variables. Objectives. Example 1 Solution
13 Functions of Several Variables 13.1 Introduction to Functions of Several Variables Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. Objectives Understand
More informationKANGAL REPORT
Individual Penalty Based Constraint handling Using a Hybrid Bi-Objective and Penalty Function Approach Rituparna Datta Kalyanmoy Deb Mechanical Engineering IIT Kanpur, India KANGAL REPORT 2013005 Abstract
More informationMAPI Computer Vision. Multiple View Geometry
MAPI Computer Vision Multiple View Geometry Geometry o Multiple Views 2- and 3- view geometry p p Kpˆ [ K R t]p Geometry o Multiple Views 2- and 3- view geometry Epipolar Geometry The epipolar geometry
More information521466S Machine Vision Exercise #1 Camera models
52466S Machine Vision Exercise # Camera models. Pinhole camera. The perspective projection equations or a pinhole camera are x n = x c, = y c, where x n = [x n, ] are the normalized image coordinates,
More informationWHAT YOU SHOULD LEARN
GRAPHS OF EQUATIONS WHAT YOU SHOULD LEARN Sketch graphs of equations. Find x- and y-intercepts of graphs of equations. Use symmetry to sketch graphs of equations. Find equations of and sketch graphs of
More informationPolar Coordinates. Calculus 2 Lia Vas. If P = (x, y) is a point in the xy-plane and O denotes the origin, let
Calculus Lia Vas Polar Coordinates If P = (x, y) is a point in the xy-plane and O denotes the origin, let r denote the distance from the origin O to the point P = (x, y). Thus, x + y = r ; θ be the angle
More informationLagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers. Lagrange Multipliers
In this section we present Lagrange s method for maximizing or minimizing a general function f(x, y, z) subject to a constraint (or side condition) of the form g(x, y, z) = k. Figure 1 shows this curve
More informationA function: A mathematical relationship between two variables (x and y), where every input value (usually x) has one output value (usually y)
SESSION 9: FUNCTIONS KEY CONCEPTS: Definitions & Terminology Graphs of Functions - Straight line - Parabola - Hyperbola - Exponential Sketching graphs Finding Equations Combinations of graphs TERMINOLOGY
More informationThe Graph of an Equation Graph the following by using a table of values and plotting points.
Calculus Preparation - Section 1 Graphs and Models Success in math as well as Calculus is to use a multiple perspective -- graphical, analytical, and numerical. Thanks to Rene Descartes we can represent
More information3.3 Optimizing Functions of Several Variables 3.4 Lagrange Multipliers
3.3 Optimizing Functions of Several Variables 3.4 Lagrange Multipliers Prof. Tesler Math 20C Fall 2018 Prof. Tesler 3.3 3.4 Optimization Math 20C / Fall 2018 1 / 56 Optimizing y = f (x) In Math 20A, we
More informationUNIT #2 TRANSFORMATIONS OF FUNCTIONS
Name: Date: UNIT # TRANSFORMATIONS OF FUNCTIONS Part I Questions. The quadratic unction ollowing does,, () has a turning point at have a turning point? 7, 3, 5 5, 8. I g 7 3, then at which o the The structure
More informationLevel Curves and Contour Maps (sec 14.1) & Partial Derivatives (sec 14.3) Today, We Will: Level Curves. Notes. Notes. Notes.
Level Curves and Contour Maps (sec 14.1) & Partial Derivatives (sec 14.3) 11 February 2015 Today, We Will: Wrap-up our discussion of level curves and contour maps. Define partial derivatives of a function
More informationCurves, Tangent Planes, and Differentials ( ) Feb. 26, 2012 (Sun) Lecture 9. Partial Derivatives: Signs on Level Curves, Tangent
Lecture 9. Partial Derivatives: Signs on Level Curves, Tangent Planes, and Differentials ( 11.3-11.4) Feb. 26, 2012 (Sun) Signs of Partial Derivatives on Level Curves Level curves are shown for a function
More information14.5 Directional Derivatives and the Gradient Vector
14.5 Directional Derivatives and the Gradient Vector 1. Directional Derivatives. Recall z = f (x, y) and the partial derivatives f x and f y are defined as f (x 0 + h, y 0 ) f (x 0, y 0 ) f x (x 0, y 0
More informationAppendix E Calculating Normal Vectors
OpenGL Programming Guide (Addison-Wesley Publishing Company) Appendix E Calculating Normal Vectors This appendix describes how to calculate normal vectors for surfaces. You need to define normals to use
More informationFinding the slope to base angle of the virtual pyramid Case 1
inding the slope to base angle of the virtual pyramid ase 1 igure 1 What we are seeking to find is the measure of angle, or conversely, as triangle is isosceles the two angles at the base will be equal.
More informationwhere ~n = ( ) t is the normal o the plane (ie, Q ~ t ~n =,8 Q ~ ), =( ~ X Y Z ) t and ~t = (t X t Y t Z ) t are the camera rotation and translation,
Multi-Frame Alignment o Planes Lihi Zelnik-Manor Michal Irani Dept o Computer Science and Applied Math The Weizmann Institute o Science Rehovot, Israel Abstract Traditional plane alignment techniques are
More informationLecture 15: Segmentation (Edge Based, Hough Transform)
Lecture 15: Segmentation (Edge Based, Hough Transform) c Bryan S. Morse, Brigham Young University, 1998 000 Last modified on February 3, 000 at :00 PM Contents 15.1 Introduction..............................................
More informationPartial Derivatives. Partial Derivatives. Partial Derivatives. Partial Derivatives. Partial Derivatives. Partial Derivatives
In general, if f is a function of two variables x and y, suppose we let only x vary while keeping y fixed, say y = b, where b is a constant. By the definition of a derivative, we have Then we are really
More information14.6 Directional Derivatives and the Gradient Vector
14 Partial Derivatives 14.6 and the Gradient Vector Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. and the Gradient Vector In this section we introduce
More informationOn the Geometry of Visual Correspondence
International Journal o Computer Vision 21(3), 223 247 (1997) c 1997 Kluwer Academic Publishers. Manuactured in The Netherlands. On the Geometry o Visual Correspondence CORNELIA FERMÜLLER AND YIANNIS ALOIMONOS
More informationHOUGH TRANSFORM CS 6350 C V
HOUGH TRANSFORM CS 6350 C V HOUGH TRANSFORM The problem: Given a set of points in 2-D, find if a sub-set of these points, fall on a LINE. Hough Transform One powerful global method for detecting edges
More informationGeometric Registration for Deformable Shapes 2.2 Deformable Registration
Geometric Registration or Deormable Shapes 2.2 Deormable Registration Variational Model Deormable ICP Variational Model What is deormable shape matching? Example? What are the Correspondences? Eurographics
More informationRepresentative Fraction (RF) is a numerical statement of the scale relationship. e.g., 1:24,000 or 1/24,000
Scale and Measurement Scale can be deined as a statement o the relationship between the distance o a map or image in relation to the distance on the Earth surace. This association may be displayed using
More informationCalypso Construction Features. Construction Features 1
Calypso 1 The Construction dropdown menu contains several useful construction features that can be used to compare two other features or perform special calculations. Construction features will show up
More informationCS485/685 Computer Vision Spring 2012 Dr. George Bebis Programming Assignment 2 Due Date: 3/27/2012
CS8/68 Computer Vision Spring 0 Dr. George Bebis Programming Assignment Due Date: /7/0 In this assignment, you will implement an algorithm or normalizing ace image using SVD. Face normalization is a required
More informationNeighbourhood Operations
Neighbourhood Operations Neighbourhood operations simply operate on a larger neighbourhood o piels than point operations Origin Neighbourhoods are mostly a rectangle around a central piel Any size rectangle
More informationTest Name: Chapter 3 Review
Test Name: Chapter 3 Review 1. For the following equation, determine the values of the missing entries. If needed, write your answer as a fraction reduced to lowest terms. 10x - 8y = 18 Note: Each column
More informationUnit 10 Reflection. Grading: Show all work, keeping it neat and organized. Show equations used and include all units.
Name: Hr: Unit 0 Relection Grading: Show all work, keeping it neat and organized. Show equations used and include all units. REFLECTION Vocabulary Relection: The bouncing o light. The angle a beam o light
More informationAdvanced Algebra. Equation of a Circle
Advanced Algebra Equation of a Circle Task on Entry Plotting Equations Using the table and axis below, plot the graph for - x 2 + y 2 = 25 x -5-4 -3 0 3 4 5 y 1 4 y 2-4 3 2 + y 2 = 25 9 + y 2 = 25 y 2
More informationSPM Add Math Form 5 Chapter 3 Integration
SPM Add Math Form Chapter Integration INDEFINITE INTEGRAL CHAPTER : INTEGRATION Integration as the reverse process of differentiation ) y if dy = x. Given that d Integral of ax n x + c = x, where c is
More informationA Cylindrical Surface Model to Rectify the Bound Document Image
A Cylindrical Surace Model to Rectiy the Bound Document Image Huaigu Cao, Xiaoqing Ding, Changsong Liu Department o Electronic Engineering, Tsinghua University State Key Laboratory o Intelligent Technology
More informationComputational Geometry. Geometry Cross Product Convex Hull Problem Sweep Line Algorithm
GEOMETRY COMP 321 McGill University These slides are mainly compiled from the following resources. - Professor Jaehyun Park slides CS 97SI - Top-coder tutorials. - Programming Challenges books. Computational
More informationComputer Data Analysis and Use of Excel
Computer Data Analysis and Use o Excel I. Theory In this lab we will use Microsot EXCEL to do our calculations and error analysis. This program was written primarily or use by the business community, so
More informationFall 2016 Semester METR 3113 Atmospheric Dynamics I: Introduction to Atmospheric Kinematics and Dynamics
Fall 2016 Semester METR 3113 Atmospheric Dynamics I: Introduction to Atmospheric Kinematics and Dynamics Lecture 5 August 31 2016 Topics: Polar coordinate system Conversion of polar coordinates to 2-D
More informationES 240: Scientific and Engineering Computation. a function f(x) that can be written as a finite series of power functions like
Polynomial Deinition a unction () that can be written as a inite series o power unctions like n is a polynomial o order n n ( ) = A polynomial is represented by coeicient vector rom highest power. p=[3-5
More informationThe Differential df, the Gradient f, & the Directional Derivative Dû f sec 14.4 (cont), Goals. Warm-up: Differentiability. Notes. Notes.
The Differential df, the Gradient f, & the Directional Derivative Dû f sec 14.4 (cont), 14.5 10 March 2014 Goals. We will: Define the differential df and use it to approximate changes in a function s value.
More information16.6. Parametric Surfaces. Parametric Surfaces. Parametric Surfaces. Vector Calculus. Parametric Surfaces and Their Areas
16 Vector Calculus 16.6 and Their Areas Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. and Their Areas Here we use vector functions to describe more general
More informationSection 18-1: Graphical Representation of Linear Equations and Functions
Section 18-1: Graphical Representation of Linear Equations and Functions Prepare a table of solutions and locate the solutions on a coordinate system: f(x) = 2x 5 Learning Outcome 2 Write x + 3 = 5 as
More information2. Recommended Design Flow
2. Recommended Design Flow This chapter describes the Altera-recommended design low or successully implementing external memory interaces in Altera devices. Altera recommends that you create an example
More informationLab 21.1 The Tangent Galvanometer
Name School Date Lab 21.1 The Tangent Galvanometer Purpose To investigate the magnetic field at the center of a current-carrying loop of wire. To verify the right-hand rule for the field inside a current
More informationComputer Data Analysis and Plotting
Phys 122 February 6, 2006 quark%//~bland/docs/manuals/ph122/pcintro/pcintro.doc Computer Data Analysis and Plotting In this lab we will use Microsot EXCEL to do our calculations. This program has been
More informationQuestion 2: How do you solve a linear programming problem with a graph?
Question : How do you solve a linear programming problem with a graph? Now that we have several linear programming problems, let s look at how we can solve them using the graph of the system of inequalities.
More informationChapter 34. Images. Two Types of Images. A Common Mirage. Plane Mirrors, Extended Object. Plane Mirrors, Point Object
Capter Images One o te most important uses o te basic laws governing ligt is te production o images. Images are critical to a variety o ields and industries ranging rom entertainment, security, and medicine
More informationd f(g(t), h(t)) = x dt + f ( y dt = 0. Notice that we can rewrite the relationship on the left hand side of the equality using the dot product: ( f
Gradients and the Directional Derivative In 14.3, we discussed the partial derivatives f f and, which tell us the rate of change of the x y height of the surface defined by f in the x direction and the
More informationReflection and Refraction
Reflection and Refraction Theory: Whenever a wave traveling in some medium encounters an interface or boundary with another medium either (or both) of the processes of (1) reflection and (2) refraction
More informationReducing the Bandwidth of a Sparse Matrix with Tabu Search
Reducing the Bandwidth o a Sparse Matrix with Tabu Search Raael Martí a, Manuel Laguna b, Fred Glover b and Vicente Campos a a b Dpto. de Estadística e Investigación Operativa, Facultad de Matemáticas,
More informationBinary Morphological Model in Refining Local Fitting Active Contour in Segmenting Weak/Missing Edges
0 International Conerence on Advanced Computer Science Applications and Technologies Binary Morphological Model in Reining Local Fitting Active Contour in Segmenting Weak/Missing Edges Norshaliza Kamaruddin,
More informationOutput Primitives Lecture: 4. Lecture 4
Lecture 4 Circle Generating Algorithms Since the circle is a frequently used component in pictures and graphs, a procedure for generating either full circles or circular arcs is included in most graphics
More informationA Proposed Approach for Solving Rough Bi-Level. Programming Problems by Genetic Algorithm
Int J Contemp Math Sciences, Vol 6, 0, no 0, 45 465 A Proposed Approach or Solving Rough Bi-Level Programming Problems by Genetic Algorithm M S Osman Department o Basic Science, Higher Technological Institute
More informationArea and Volume. where x right and x left are written in terms of y.
Area and Volume Area between two curves Sketch the region and determine the points of intersection. Draw a small strip either as dx or dy slicing. Use the following templates to set up a definite integral:
More information2.3 Additional Relations
3 2.3 Additional Relations Figure 2.3 identiies additional relations, indicating te locations o te object and image, and te ratio o teir eigts (magniication) and orientations. Ray enters te lens parallel
More informationChapter 3. Sukhwinder Singh
Chapter 3 Sukhwinder Singh PIXEL ADDRESSING AND OBJECT GEOMETRY Object descriptions are given in a world reference frame, chosen to suit a particular application, and input world coordinates are ultimately
More informationEXTREME POINTS AND AFFINE EQUIVALENCE
EXTREME POINTS AND AFFINE EQUIVALENCE The purpose of this note is to use the notions of extreme points and affine transformations which are studied in the file affine-convex.pdf to prove that certain standard
More informationThe Rational Zero Theorem
The Rational Zero Theorem Our goal in this section is to learn how we can ind the rational zeros o the polynomials. For example: x = x 4 + x x x + ( ) We could randomly try some actors and use synthetic
More informationLifting Transform, Voronoi, Delaunay, Convex Hulls
Lifting Transform, Voronoi, Delaunay, Convex Hulls Subhash Suri Department of Computer Science University of California Santa Barbara, CA 93106 1 Lifting Transform (A combination of Pless notes and my
More informationRelaxing the 3L algorithm for an accurate implicit polynomial fitting
Relaxing the 3L algorithm or an accurate implicit polynomial itting Mohammad Rouhani Computer Vision Center Ediici O, Campus UAB 08193 Bellaterra, Barcelona, Spain rouhani@cvc.uab.es Angel D. Sappa Computer
More information(c) 0 (d) (a) 27 (b) (e) x 2 3x2
1. Sarah the architect is designing a modern building. The base of the building is the region in the xy-plane bounded by x =, y =, and y = 3 x. The building itself has a height bounded between z = and
More informationHimanshu Jain and Kalyanmoy Deb, Fellow, IEEE
An Evolutionary Many-Objective Optimization Algorithm Using Reerence-point Based Non-dominated Sorting Approach, Part II: Handling Constraints and Extending to an Adaptive Approach Himanshu Jain and Kalyanmoy
More informationThe diagram above shows a sketch of the curve C with parametric equations
1. The diagram above shows a sketch of the curve C with parametric equations x = 5t 4, y = t(9 t ) The curve C cuts the x-axis at the points A and B. (a) Find the x-coordinate at the point A and the x-coordinate
More informationMath 2A Vector Calculus Chapter 11 Test Fall 07 Name Show your work. Don t use a calculator. Write responses on separate paper.
Math A Vector Calculus Chapter Test Fall 7 Name Show our work. Don t use a calculator. Write responses on separate paper.. Consider the nice, smooth unction z, whose contour map is shown at right. a. Estimate
More informationChapter 5 Partial Differentiation
Chapter 5 Partial Differentiation For functions of one variable, y = f (x), the rate of change of the dependent variable can dy be found unambiguously by differentiation: f x. In this chapter we explore
More informationLog1 Contest Round 1 Theta Circles & Polygons. 4 points each. 5 points each
014 015 Log1 Contest Round 1 Theta Circles & Polygons 1 Find the area, in square inches, enclosed by a circle whose diameter is 8 inches. A rectangle has sides of length 4 and 6. Find the area enclosed
More informationResearch on Image Splicing Based on Weighted POISSON Fusion
Research on Image Splicing Based on Weighted POISSO Fusion Dan Li, Ling Yuan*, Song Hu, Zeqi Wang School o Computer Science & Technology HuaZhong University o Science & Technology Wuhan, 430074, China
More informationThe Cut Locus and the Jordan Curve Theorem
The Cut Locus and the Jordan Curve Theorem Rich Schwartz November 19, 2015 1 Introduction A Jordan curve is a subset of R 2 which is homeomorphic to the circle, S 1. The famous Jordan Curve Theorem says
More information6. Find the equation of the plane that passes through the point (-1,2,1) and contains the line x = y = z.
Week 1 Worksheet Sections from Thomas 13 th edition: 12.4, 12.5, 12.6, 13.1 1. A plane is a set of points that satisfies an equation of the form c 1 x + c 2 y + c 3 z = c 4. (a) Find any three distinct
More informationMultivariate Calculus Review Problems for Examination Two
Multivariate Calculus Review Problems for Examination Two Note: Exam Two is on Thursday, February 28, class time. The coverage is multivariate differential calculus and double integration: sections 13.3,
More informationPolar Coordinates. 2, π and ( )
Polar Coordinates Up to this point we ve dealt exclusively with the Cartesian (or Rectangular, or x-y) coordinate system. However, as we will see, this is not always the easiest coordinate system to work
More informationA Singular Example for the Averaged Mean Curvature Flow
To appear in Experimental Mathematics Preprint Vol. No. () pp. 3 7 February 9, A Singular Example for the Averaged Mean Curvature Flow Uwe F. Mayer Abstract An embedded curve is presented which under numerical
More informationComputer Vision. Image Segmentation. 10. Segmentation. Computer Engineering, Sejong University. Dongil Han
Computer Vision 10. Segmentation Computer Engineering, Sejong University Dongil Han Image Segmentation Image segmentation Subdivides an image into its constituent regions or objects - After an image has
More informationRoad Sign Analysis Using Multisensory Data
Road Sign Analysis Using Multisensory Data R.J. López-Sastre, S. Lauente-Arroyo, P. Gil-Jiménez, P. Siegmann, and S. Maldonado-Bascón University o Alcalá, Department o Signal Theory and Communications
More informationBalanced Multiresolution for Symmetric/Antisymmetric Filters
Balanced Multiresolution or Symmetric/Antisymmetric Filters Mahmudul Hasan Faramarz F. Samavati Mario C. Sousa Department o Computer Science University o Calgary Alberta Canada {mhasan samavati smcosta}@ucalgary.ca
More informationMath 20C. Lecture Examples.
Math 20C. Lecture Eamples. (8/7/0) Section 4.4. Linear approimations and tangent planes A function z = f(,) of two variables is linear if its graph in z-space is a plane. Equations of planes were found
More informationUNIT 2 GRAPHIC PRIMITIVES
UNIT 2 GRAPHIC PRIMITIVES Structure Page Nos. 2.1 Introduction 46 2.2 Objectives 46 2.3 Points and Lines 46 2.4 Line Generation Algorithms 48 2.4.1 DDA Algorithm 49 2.4.2 Bresenhams Line Generation Algorithm
More informationRenderer Implementation: Basics and Clipping. Overview. Preliminaries. David Carr Virtual Environments, Fundamentals Spring 2005
INSTITUTIONEN FÖR SYSTEMTEKNIK LULEÅ TEKNISKA UNIVERSITET Renderer Implementation: Basics and Clipping David Carr Virtual Environments, Fundamentals Spring 2005 Feb-28-05 SMM009, Basics and Clipping 1
More informationClassifier Evasion: Models and Open Problems
Classiier Evasion: Models and Open Problems Blaine Nelson 1, Benjamin I. P. Rubinstein 2, Ling Huang 3, Anthony D. Joseph 1,3, and J. D. Tygar 1 1 UC Berkeley 2 Microsot Research 3 Intel Labs Berkeley
More informationRasterization. COMP 575/770 Spring 2013
Rasterization COMP 575/770 Spring 2013 The Rasterization Pipeline you are here APPLICATION COMMAND STREAM 3D transformations; shading VERTEX PROCESSING TRANSFORMED GEOMETRY conversion of primitives to
More information