Numerical Derivatives
|
|
- Brian Bennett
- 6 years ago
- Views:
Transcription
1 Lab 15 Numerical Derivatives Lab Objective: Understand and implement finite difference approximations of te derivative in single and multiple dimensions. Evaluate te accuracy of tese approximations. Ten use finite difference quotients to find edges in images via te Sobel filter. Derivative Approximations in One Dimension Te derivative of a function f at a point x 0 is f f(x 0 + ) f(x 0) (x 0) = lim. (15.1) 0 In tis lab, we will investigate one way a computer can calculate f (x 0). Forward Difference Quotient Suppose tat in Equation (15.1), instead of taking a limit, we just pick a small value for. Ten we would expect f (x 0) to be close to te quantity f(x 0 + ) f(x 0). (15.) Tis quotient is called te first order forward difference approximation of te derivative. Because f (x 0) is te limit of suc quotients, we expect tat wen is small, tis quotient is close to f (x 0). We can use Taylor s formula to find just ow close. By Taylor s formula, f(x 0 + ) = f(x 0) + f (x 0) + R (), ( ) 1 were R () = (1 t)f (x t)dt. (Tis is called te integral form of te remainder for Taylor s Teorem; see Volume 1 Capter 6). Wen we solve tis equation for f (x 0), we get f f(x0 + ) f(x0) (x 0) = R(). (15.3) Tus, te error in using te first order forward difference quotient to approximate f (x 0) is R () 1 1 t f (x 0 + t) dt
2 164 Lab 15. Numerical Derivatives If we assume f is continuous, ten for any δ, set M = sup x (x0 δ,x 0 +δ) f (x). Ten if < δ, we ave R() 1 Mdt = M O(). Terefore, te error in using (15.) to approximate f (x 0) grows like. Centered Difference Quotient 0 In fact, we can approximate f (x 0) to te second order wit anoter difference quotient, called te centered difference quotient. We begin by trying to find te backward difference quotient. Evaluate Taylor s formula at x 0 to derive f (x 0) = f(x0) f(x0 ) + R( ). (15.4) Te first term on te rigt and side of (15.4) is called te backward difference quotient. Tis quotient also approximates f (x 0) to first order, so it is not te quotient we are looking for. Wen we add (15.3) and (15.4) and solve for f (x 0) (by dividing by ), we get f 1 (x 0) = f(x0 + ) 1 f(x0 ) R( ) R() + (15.5) Te centered difference quotient is te first term of te rigt and side of (15.5). Let us investigate te remainder term to see ow accurate tis approximation is. Recall from te proof of Taylor s teorem tat R k = f (k) (x 0 ) R ( ) R () = 1 ( f (x 0) k + R k! k+1. Terefore, + R 3( ) f ) (x 0) R 3() = 1 (R3( ) R3()) = 1 (( 1 (1 t) f (x 0 + t)dt 0 ( 1 (1 t) = 4 0 O( ) ) ( 1 3 ) (f (x 0 + t) f (x 0 t)) 0 (1 t) ) ) f (x 0 t)dt 3 once we restrict to some δ-neigborood of 0. So te error in using te centered difference quotient to approximate f (x 0) grows like, wic is smaller tan wen < 1. Accuracy of Approximations Let us discuss wat step size we sould plug into te difference quotients to get te best approximation to f (x 0). Since f is defined as a limit as 0, you may tink tat it is best to coose as small as possible, but tis is not te case. In fact, dividing by very small numbers causes errors in floating point aritmetic. Tis means tat as we decrease, te error between f (x 0) and te difference quotient will first decrease, but ten increase wen gets too small because of floating point aritmetic. Here is an example wit te function f(x) = e x. A quick way to write f as a function in Pyton is wit te lambda keyword. >>> import numpy as np >>> from matplotlib import pyplot as plt >>> f = lambda x: np.exp(x)
3 165 1e-1 1e-3 1e-5 1e-7 1e-9 1e-11 Error 5e-3 5e-7 6e-11 6e-11 7e-9 1e-5 Table 15.1: Tis table sows tat it is best not to coose too small wen you approximate derivatives wit difference quotients. Here, Error equals te absolute value of f (1) f app (1) were f(x) = e x and f app is te centered difference approximation to f. In general, te line f = lambda <params> : <expression> is equivalent to defining a function f tat accepts te parameters params and returns expression. Next we fix a step size and define an approximation to te derivative of f using te centered difference quotient. >>> = 1e-1 >>> Df_app = lambda x:.5*(f(x+)-f(x-))/ Finally, we ceck te accuracy of tis approximation at x 0 = 1 by computing te difference between Df_app(1) and te actual derivative evaluated at 1. # Since f(x) = e^x, te derivative of f(x) is f(x) >>> np.abs( f(1)-df_app(1) ) We note tat our functions f and Df_app beave as expected wen tey are passed a NumPy array. >>> = np.array([1e-1, 1e-3, 1e-5, 1e-7, 1e-9, 1e-11]) >>> np.abs( f(1)-df_app(1) ) array([ e-03, e-07, e-11, e-11, e-09, e-05]) Tese results are summarized in Table Tus, te optimal value of is one tat is small, but not too small. A good coice is = 1e-5. Problem 1. Write a function tat accepts as input a callable function object f, an array of points pts, and a keyword argument tat defaults to 1e-5. Return an array of te centered difference quotients of f at eac point in pts wit te specified value of. You may wonder if te forward or backward difference quotients are ever used, since te centered difference quotient is a more accurate approximation of te derivative. In fact, tere are some functions tat in practice do not beave well under centered difference quotients. In tese cases, one must use te forward or backward difference quotient. Finally, we remark tat forward, backward, and centered difference quotients can be used to approximate iger-order derivatives of f. However, taking derivatives is an unstable operation. Tis means tat taking a derivative can amplify te aritmetic error in your computation. For tis reason, difference quotients are not generally used to approximate derivatives iger tan second order.
4 166 Lab 15. Numerical Derivatives Derivative Approximations in Multiple Dimensions Finite difference metods can also be used to calculate derivatives in iger dimensions. Recall tat te Jacobian of a function f : R n R m at a point x 0 R n is te m n matrix J = (J ij) defined component-wise by J ij = i x j (x 0). For example, te Jacobian for a function f : R 3 R is defined by J = ( x 1 x x 3 ) = ( 1 x 1 1 x 1 x 3 x 1 x x 3 Te Jacobian is useful in many applications. For example, te Jacobian can be used to find zeros of functions in multiple variables. Te forward difference quotient for approximating a partial derivative is x j (x 0) f(x0 + ej) f(x0), were e j is te j t standard basis vector. Similarly, te centered difference approximation is 1 (x 0) f(x0 + ej) 1 f(x0 ej). x j ). Problem. Write a function tat accepts 1. a function andle f,. an integer n tat is te dimension of te domain of f, 3. an integer m tat is te dimension of te range of f, 4. an 1 x n-dimensional NumPy array pt representing a point in R n, and 5. a keyword argument tat defaults to 1e-5. Return te approximate Jacobian matrix of f at pt using te centered difference quotient. Problem 3. Let f : R R be defined by ( ) e x sin(y) + y 3 f(x, y) = 3y cos(x) Find te error between your Jacobian function and te analytically computed derivative on te square [ 1, 1] [ 1, 1] using ten tousand grid points (100 per side). You may apply your Jacobian function to te points one at a time using a double for loop. Once you get te error matrix for a given point, calculate te Frobenius norm of tis matrix (la.norm defaults to te Frobenius norm). Tis norm will be your total error for tat point. Wat is te maximum error of your Jacobian function over all points in te square?
5 167 ( ) x Hint: Te following code defines te function f(x, y) =. x + y # f accepts a lengt- NumPy array >>> f = lambda x: np.array([x[0]**, x[0]+x[1]]) Application to Image Filters Recall tat a computer stores an image as a -D array of pixel values (i.e., a matrix of intensities). An image filter is a function tat transforms an image by operating on it locally. Tat is, to compute te ij t pixel value in te new image, an image filter uses only te pixels in a small neigborood around te ij t pixel in te original image. In tis lab, we will use a filter derived from te gradient of an image to find edges in an image. Convolutions One example of an image filter is to convolve an image wit a filter matrix. A filter matrix is a matrix wose eigt and widt are relatively small odd numbers. If te filter matrix is f 1, 1 f 1,0 f 1,1 F = f 0, 1 f 0,0 f 0,1, f 1, 1 f 1,0 f 1,1 ten te convolution of an image A wit F is A F = (C ij) were C ij = 1 k= 1 l= 1 1 f kl A i+k,j+l. (15.6) Say A is an m n matrix. Here, we take A ij = 0 wen i {1,... m} or j {1,..., n}. Te value of C ij is a linear combination of te nearby pixel values, wit coefficients given by F (see Figure 15.1). In fact, C ij equals te Frobenius inner product of F wit te 3 3 submatrix of A centered at ij. Implementation in NumPy Let us write a function tat convolves an image wit a filter. You can test tis function on te image cameraman.jpg, wic appears in Figure 15.a. Te following code loads tis image and plots it wit matplotlib. >>> image = plt.imread('cameraman.jpg') >>> plt.imsow(image, cmap = 'gray') >>> plt.sow() Here is te function definition and some setup. 1. def Filter(image, F):. m, n = image.sape 3., k = F.sape
6 168 Lab 15. Numerical Derivatives Figure 15.1: Tis diagram illustrates ow to convolve an image wit a filter. Te ligt grey rectangle represents te original image A, and te dark grey squares are te filter F. Te larger rectangle is te image padded wit zeros; i.e., all pixel values in te outer wite band are 0. To compute te entry of te convolution matrix C located at a black dot, take te inner product of F wit te submatrix of te padded image centered at te dot. To convolve image wit te filter F, we must first pad te array image wit zeros around te edges. Tis is because in (15.6), entries A ij are set to zero wen i or j is out of bounds. We do tis by creating a larger array of zeros, and ten making te interior part of te array equal to te original image (see Figure 15.1). For example, if te filter is a 3 3 matrix, ten te following code will pad te matrix wit te appropriate number of zeros. # Create a larger matrix of zeros image_pad = np.zeros((m+, n+)) # Make te interior of image_pad equal to te original image image_pad[1:1+m, 1:1+n] = image We want to do tis in general in our function. Note tat te number of zeros we need to pad our array depends on te size of te filter F. 5. image_pad = # Create an array of zeros of te appropriate size 6. # Make te interior of image_pad equal to image Finally, we iterate troug te image to compute eac entry of te convolution matrix. 7. C = np.zeros(image.sape) 8. for i in range(m): 9. for j in range(n): 10. C[i,j] = # Compute C[i, j]
7 169 Gaussian Blur A Gaussian blur is an image filter tat operates on an image by convolving wit te matrix G = Blurring an image can remove noise, or random variation tat is te visual analog of static in a radio signal (and equally undesirable). Problem 4. Finis writing te function Filter by filling in lines 5, 6, and 10. Hint: Note in 15.6, C ij was calculated by summing from -1 to 1. Tis is only te case if te filter F is 3 3. A sligt modification is needed in te general case. Test your function on te image cameraman.jpg using te Gaussian Blur. Te result is in Figure 15.b. Edge Detection Automatic detection of edges in an image can be used to segment or sarpen te image. We will find edges wit te Sobel filter, wic computes te gradient of te image at eac pixel. Te magnitude of te gradient tells us te rate of cange of te pixel values, and so large magnitudes sould correspond to edges witin te image. Te Sobel filter is not a convolution, altoug it does use convolutions. We can tink of an image as a function from a grid of points to R. Te image maps a pixel location to an intensity. It does not make sense to define te derivative of tis function as a limit because te domain is discrete a step size cannot take on arbitrarily small values. Instead, we define te derivative to be te centered difference quotient of te previous section. Tat is, we define te derivative in te x-direction at te ij t pixel to be 1 Ai+1,j 1 Ai 1,j. We can use a convolution to create a matrix A x wose ij t entry is te derivative of A at te ij t entry, in te x-direction. In fact, A x = A S, were S = Note tat tis convolution takes a weigted average of te x-derivatives at (i, j), (i, j + 1), and (i, j 1). Te derivative at (i, j) is weigted by. Using a weigted average instead of just te derivative at (i, j) makes te derivative less affected by noise. Now we can define te Sobel filter. A Sobel filter applied to an image A results in an array B = (B ij) of 0 s and 1 s, were te 1 s trace out te edges in te image. By definition, { 1 if A(ij) > M B ij = 0 oterwise. Here, A(ij) = ((A S) ij, (A S T ) ij) is te gradient of A at te ij t pixel. Te constant M sould be sufficiently large enoug to pick out tose pixels wit te largest gradient
8 170 Lab 15. Numerical Derivatives (a) Unfiltered image. (b) Image after Gaussian blur is applied. (c) Image after te Sobel filter is applied. Figure 15.: Here is an example of a Gaussian blur and te Sobel filter applied to an image. Tis poto, known as cameraman, is a standard test image in image processing. A database of suc images can be downloaded from ttp:// (i.e., tose pixels tat are part of an edge). A good coice for M is 4 times te average value of A(ij) over te wole image A. Wen te Sobel filter is applied to cameraman.jpg, we get te image in Figure 15.c. Here, te 1 s in B were mapped to wite and te 0 s were mapped to black. Problem 5. Write a function tat accepts an image as input and applies te Sobel filter to te image. Test your function on cameraman.jpg. Hint: If you want to find te average of a matrix A, use te function A.mean().
4.1 Tangent Lines. y 2 y 1 = y 2 y 1
41 Tangent Lines Introduction Recall tat te slope of a line tells us ow fast te line rises or falls Given distinct points (x 1, y 1 ) and (x 2, y 2 ), te slope of te line troug tese two points is cange
More information3.6 Directional Derivatives and the Gradient Vector
288 CHAPTER 3. FUNCTIONS OF SEVERAL VARIABLES 3.6 Directional Derivatives and te Gradient Vector 3.6.1 Functions of two Variables Directional Derivatives Let us first quickly review, one more time, te
More informationSection 2.3: Calculating Limits using the Limit Laws
Section 2.3: Calculating Limits using te Limit Laws In previous sections, we used graps and numerics to approimate te value of a it if it eists. Te problem wit tis owever is tat it does not always give
More informationMaterials: Whiteboard, TI-Nspire classroom set, quadratic tangents program, and a computer projector.
Adam Clinc Lesson: Deriving te Derivative Grade Level: 12 t grade, Calculus I class Materials: Witeboard, TI-Nspire classroom set, quadratic tangents program, and a computer projector. Goals/Objectives:
More information19.2 Surface Area of Prisms and Cylinders
Name Class Date 19 Surface Area of Prisms and Cylinders Essential Question: How can you find te surface area of a prism or cylinder? Resource Locker Explore Developing a Surface Area Formula Surface area
More informationHaar Transform CS 430 Denbigh Starkey
Haar Transform CS Denbig Starkey. Background. Computing te transform. Restoring te original image from te transform 7. Producing te transform matrix 8 5. Using Haar for lossless compression 6. Using Haar
More informationMATH 5a Spring 2018 READING ASSIGNMENTS FOR CHAPTER 2
MATH 5a Spring 2018 READING ASSIGNMENTS FOR CHAPTER 2 Note: Tere will be a very sort online reading quiz (WebWork) on eac reading assignment due one our before class on its due date. Due dates can be found
More informationMore on Functions and Their Graphs
More on Functions and Teir Graps Difference Quotient ( + ) ( ) f a f a is known as te difference quotient and is used exclusively wit functions. Te objective to keep in mind is to factor te appearing in
More information4.2 The Derivative. f(x + h) f(x) lim
4.2 Te Derivative Introduction In te previous section, it was sown tat if a function f as a nonvertical tangent line at a point (x, f(x)), ten its slope is given by te it f(x + ) f(x). (*) Tis is potentially
More informationPiecewise Polynomial Interpolation, cont d
Jim Lambers MAT 460/560 Fall Semester 2009-0 Lecture 2 Notes Tese notes correspond to Section 4 in te text Piecewise Polynomial Interpolation, cont d Constructing Cubic Splines, cont d Having determined
More informationAll truths are easy to understand once they are discovered; the point is to discover them. Galileo
Section 7. olume All truts are easy to understand once tey are discovered; te point is to discover tem. Galileo Te main topic of tis section is volume. You will specifically look at ow to find te volume
More information12.2 Techniques for Evaluating Limits
335_qd /4/5 :5 PM Page 863 Section Tecniques for Evaluating Limits 863 Tecniques for Evaluating Limits Wat ou sould learn Use te dividing out tecnique to evaluate its of functions Use te rationalizing
More informationChapter K. Geometric Optics. Blinn College - Physics Terry Honan
Capter K Geometric Optics Blinn College - Pysics 2426 - Terry Honan K. - Properties of Ligt Te Speed of Ligt Te speed of ligt in a vacuum is approximately c > 3.0µ0 8 mês. Because of its most fundamental
More information12.2 Investigate Surface Area
Investigating g Geometry ACTIVITY Use before Lesson 12.2 12.2 Investigate Surface Area MATERIALS grap paper scissors tape Q U E S T I O N How can you find te surface area of a polyedron? A net is a pattern
More information12.2 TECHNIQUES FOR EVALUATING LIMITS
Section Tecniques for Evaluating Limits 86 TECHNIQUES FOR EVALUATING LIMITS Wat ou sould learn Use te dividing out tecnique to evaluate its of functions Use te rationalizing tecnique to evaluate its of
More informationBounding Tree Cover Number and Positive Semidefinite Zero Forcing Number
Bounding Tree Cover Number and Positive Semidefinite Zero Forcing Number Sofia Burille Mentor: Micael Natanson September 15, 2014 Abstract Given a grap, G, wit a set of vertices, v, and edges, various
More informationFast Calculation of Thermodynamic Properties of Water and Steam in Process Modelling using Spline Interpolation
P R E P R N T CPWS XV Berlin, September 8, 008 Fast Calculation of Termodynamic Properties of Water and Steam in Process Modelling using Spline nterpolation Mattias Kunick a, Hans-Joacim Kretzscmar a,
More informationLinear Interpolating Splines
Jim Lambers MAT 772 Fall Semester 2010-11 Lecture 17 Notes Tese notes correspond to Sections 112, 11, and 114 in te text Linear Interpolating Splines We ave seen tat ig-degree polynomial interpolation
More information2 The Derivative. 2.0 Introduction to Derivatives. Slopes of Tangent Lines: Graphically
2 Te Derivative Te two previous capters ave laid te foundation for te study of calculus. Tey provided a review of some material you will need and started to empasize te various ways we will view and use
More informationTwo Modifications of Weight Calculation of the Non-Local Means Denoising Method
Engineering, 2013, 5, 522-526 ttp://dx.doi.org/10.4236/eng.2013.510b107 Publised Online October 2013 (ttp://www.scirp.org/journal/eng) Two Modifications of Weigt Calculation of te Non-Local Means Denoising
More information13.5 DIRECTIONAL DERIVATIVES and the GRADIENT VECTOR
13.5 Directional Derivatives and te Gradient Vector Contemporary Calculus 1 13.5 DIRECTIONAL DERIVATIVES and te GRADIENT VECTOR Directional Derivatives In Section 13.3 te partial derivatives f x and f
More informationNOTES: A quick overview of 2-D geometry
NOTES: A quick overview of 2-D geometry Wat is 2-D geometry? Also called plane geometry, it s te geometry tat deals wit two dimensional sapes flat tings tat ave lengt and widt, suc as a piece of paper.
More informationSection 1.2 The Slope of a Tangent
Section 1.2 Te Slope of a Tangent You are familiar wit te concept of a tangent to a curve. Wat geometric interpretation can be given to a tangent to te grap of a function at a point? A tangent is te straigt
More information, 1 1, A complex fraction is a quotient of rational expressions (including their sums) that result
RT. Complex Fractions Wen working wit algebraic expressions, sometimes we come across needing to simplify expressions like tese: xx 9 xx +, xx + xx + xx, yy xx + xx + +, aa Simplifying Complex Fractions
More informationMeasuring Length 11and Area
Measuring Lengt 11and Area 11.1 Areas of Triangles and Parallelograms 11.2 Areas of Trapezoids, Romuses, and Kites 11.3 Perimeter and Area of Similar Figures 11.4 Circumference and Arc Lengt 11.5 Areas
More informationMTH-112 Quiz 1 - Solutions
MTH- Quiz - Solutions Words in italics are for eplanation purposes onl (not necessar to write in te tests or. Determine weter te given relation is a function. Give te domain and range of te relation. {(,
More informationInterference and Diffraction of Light
Interference and Diffraction of Ligt References: [1] A.P. Frenc: Vibrations and Waves, Norton Publ. 1971, Capter 8, p. 280-297 [2] PASCO Interference and Diffraction EX-9918 guide (written by Ann Hanks)
More information6 Computing Derivatives the Quick and Easy Way
Jay Daigle Occiental College Mat 4: Calculus Experience 6 Computing Derivatives te Quick an Easy Way In te previous section we talke about wat te erivative is, an we compute several examples, an ten we
More information1.4 RATIONAL EXPRESSIONS
6 CHAPTER Fundamentals.4 RATIONAL EXPRESSIONS Te Domain of an Algebraic Epression Simplifying Rational Epressions Multiplying and Dividing Rational Epressions Adding and Subtracting Rational Epressions
More information2.8 The derivative as a function
CHAPTER 2. LIMITS 56 2.8 Te derivative as a function Definition. Te derivative of f(x) istefunction f (x) defined as follows f f(x + ) f(x) (x). 0 Note: tis differs from te definition in section 2.7 in
More informationCSE 332: Data Structures & Parallelism Lecture 8: AVL Trees. Ruth Anderson Winter 2019
CSE 2: Data Structures & Parallelism Lecture 8: AVL Trees Rut Anderson Winter 29 Today Dictionaries AVL Trees /25/29 2 Te AVL Balance Condition: Left and rigt subtrees of every node ave eigts differing
More informationCS 234. Module 6. October 16, CS 234 Module 6 ADT Dictionary 1 / 33
CS 234 Module 6 October 16, 2018 CS 234 Module 6 ADT Dictionary 1 / 33 Idea for an ADT Te ADT Dictionary stores pairs (key, element), were keys are distinct and elements can be any data. Notes: Tis is
More informationPYRAMID FILTERS BASED ON BILINEAR INTERPOLATION
PYRAMID FILTERS BASED ON BILINEAR INTERPOLATION Martin Kraus Computer Grapics and Visualization Group, Tecnisce Universität Müncen, Germany krausma@in.tum.de Magnus Strengert Visualization and Interactive
More informationVector Processing Contours
Vector Processing Contours Andrey Kirsanov Department of Automation and Control Processes MAMI Moscow State Tecnical University Moscow, Russia AndKirsanov@yandex.ru A.Vavilin and K-H. Jo Department of
More informationNotes: Dimensional Analysis / Conversions
Wat is a unit system? A unit system is a metod of taking a measurement. Simple as tat. We ave units for distance, time, temperature, pressure, energy, mass, and many more. Wy is it important to ave a standard?
More informationLesson 6 MA Nick Egbert
Overview From kindergarten we all know ow to find te slope of a line: rise over run, or cange in over cange in. We want to be able to determine slopes of functions wic are not lines. To do tis we use te
More informationThe Euler and trapezoidal stencils to solve d d x y x = f x, y x
restart; Te Euler and trapezoidal stencils to solve d d x y x = y x Te purpose of tis workseet is to derive te tree simplest numerical stencils to solve te first order d equation y x d x = y x, and study
More informationThe (, D) and (, N) problems in double-step digraphs with unilateral distance
Electronic Journal of Grap Teory and Applications () (), Te (, D) and (, N) problems in double-step digraps wit unilateral distance C Dalfó, MA Fiol Departament de Matemàtica Aplicada IV Universitat Politècnica
More informationYou Try: A. Dilate the following figure using a scale factor of 2 with center of dilation at the origin.
1 G.SRT.1-Some Tings To Know Dilations affect te size of te pre-image. Te pre-image will enlarge or reduce by te ratio given by te scale factor. A dilation wit a scale factor of 1> x >1enlarges it. A dilation
More information8/6/2010 Assignment Previewer
Week 4 Friday Homework (1321979) Question 1234567891011121314151617181920 1. Question DetailsSCalcET6 2.7.003. [1287988] Consider te parabola y 7x - x 2. (a) Find te slope of te tangent line to te parabola
More informationClassify solids. Find volumes of prisms and cylinders.
11.4 Volumes of Prisms and Cylinders Essential Question How can you find te volume of a prism or cylinder tat is not a rigt prism or rigt cylinder? Recall tat te volume V of a rigt prism or a rigt cylinder
More information2D transformations Homogeneous coordinates. Uses of Transformations
2D transformations omogeneous coordinates Uses of Transformations Modeling: position and resize parts of a complex model; Viewing: define and position te virtual camera Animation: define ow objects move/cange
More informationMulti-Stack Boundary Labeling Problems
Multi-Stack Boundary Labeling Problems Micael A. Bekos 1, Micael Kaufmann 2, Katerina Potika 1 Antonios Symvonis 1 1 National Tecnical University of Atens, Scool of Applied Matematical & Pysical Sciences,
More informationRECONSTRUCTING OF A GIVEN PIXEL S THREE- DIMENSIONAL COORDINATES GIVEN BY A PERSPECTIVE DIGITAL AERIAL PHOTOS BY APPLYING DIGITAL TERRAIN MODEL
IV. Évfolyam 3. szám - 2009. szeptember Horvát Zoltán orvat.zoltan@zmne.u REONSTRUTING OF GIVEN PIXEL S THREE- DIMENSIONL OORDINTES GIVEN Y PERSPETIVE DIGITL ERIL PHOTOS Y PPLYING DIGITL TERRIN MODEL bsztrakt/bstract
More informationYou should be able to visually approximate the slope of a graph. The slope m of the graph of f at the point x, f x is given by
Section. Te Tangent Line Problem 89 87. r 5 sin, e, 88. r sin sin Parabola 9 9 Hperbola e 9 9 9 89. 7,,,, 5 7 8 5 ortogonal 9. 5, 5,, 5, 5. Not multiples of eac oter; neiter parallel nor ortogonal 9.,,,
More informationCubic smoothing spline
Cubic smooting spline Menu: QCExpert Regression Cubic spline e module Cubic Spline is used to fit any functional regression curve troug data wit one independent variable x and one dependent random variable
More informationAnnouncements. Lilian s office hours rescheduled: Fri 2-4pm HW2 out tomorrow, due Thursday, 7/7. CSE373: Data Structures & Algorithms
Announcements Lilian s office ours resceduled: Fri 2-4pm HW2 out tomorrow, due Tursday, 7/7 CSE373: Data Structures & Algoritms Deletion in BST 2 5 5 2 9 20 7 0 7 30 Wy migt deletion be arder tan insertion?
More informationHash-Based Indexes. Chapter 11. Comp 521 Files and Databases Spring
Has-Based Indexes Capter 11 Comp 521 Files and Databases Spring 2010 1 Introduction As for any index, 3 alternatives for data entries k*: Data record wit key value k
More informationFault Localization Using Tarantula
Class 20 Fault localization (cont d) Test-data generation Exam review: Nov 3, after class to :30 Responsible for all material up troug Nov 3 (troug test-data generation) Send questions beforeand so all
More informationAreas of Parallelograms and Triangles. To find the area of parallelograms and triangles
10-1 reas of Parallelograms and Triangles ommon ore State Standards G-MG..1 Use geometric sapes, teir measures, and teir properties to descrie ojects. G-GPE..7 Use coordinates to compute perimeters of
More informationMAC-CPTM Situations Project
raft o not use witout permission -P ituations Project ituation 20: rea of Plane Figures Prompt teacer in a geometry class introduces formulas for te areas of parallelograms, trapezoids, and romi. e removes
More informationCHAPTER 7: TRANSCENDENTAL FUNCTIONS
7.0 Introduction and One to one Functions Contemporary Calculus 1 CHAPTER 7: TRANSCENDENTAL FUNCTIONS Introduction In te previous capters we saw ow to calculate and use te derivatives and integrals of
More informationLimits and Continuity
CHAPTER Limits and Continuit. Rates of Cange and Limits. Limits Involving Infinit.3 Continuit.4 Rates of Cange and Tangent Lines An Economic Injur Level (EIL) is a measurement of te fewest number of insect
More informationData Structures and Programming Spring 2014, Midterm Exam.
Data Structures and Programming Spring 2014, Midterm Exam. 1. (10 pts) Order te following functions 2.2 n, log(n 10 ), 2 2012, 25n log(n), 1.1 n, 2n 5.5, 4 log(n), 2 10, n 1.02, 5n 5, 76n, 8n 5 + 5n 2
More informationWhen the dimensions of a solid increase by a factor of k, how does the surface area change? How does the volume change?
8.4 Surface Areas and Volumes of Similar Solids Wen te dimensions of a solid increase by a factor of k, ow does te surface area cange? How does te volume cange? 1 ACTIVITY: Comparing Surface Areas and
More informationVOLUMES. The volume of a cylinder is determined by multiplying the cross sectional area by the height. r h V. a) 10 mm 25 mm.
OLUME OF A CYLINDER OLUMES Te volume of a cylinder is determined by multiplying te cross sectional area by te eigt. r Were: = volume r = radius = eigt Exercise 1 Complete te table ( =.14) r a) 10 mm 5
More informationDensity Estimation Over Data Stream
Density Estimation Over Data Stream Aoying Zou Dept. of Computer Science, Fudan University 22 Handan Rd. Sangai, 2433, P.R. Cina ayzou@fudan.edu.cn Ziyuan Cai Dept. of Computer Science, Fudan University
More informationMAPI Computer Vision
MAPI Computer Vision Multiple View Geometry In tis module we intend to present several tecniques in te domain of te 3D vision Manuel Joao University of Mino Dep Industrial Electronics - Applications -
More informationAreas of Triangles and Parallelograms. Bases of a parallelogram. Height of a parallelogram THEOREM 11.3: AREA OF A TRIANGLE. a and its corresponding.
11.1 Areas of Triangles and Parallelograms Goal p Find areas of triangles and parallelograms. Your Notes VOCABULARY Bases of a parallelogram Heigt of a parallelogram POSTULATE 4: AREA OF A SQUARE POSTULATE
More informationSection 3. Imaging With A Thin Lens
Section 3 Imaging Wit A Tin Lens 3- at Ininity An object at ininity produces a set o collimated set o rays entering te optical system. Consider te rays rom a inite object located on te axis. Wen te object
More informationImplementation of Integral based Digital Curvature Estimators in DGtal
Implementation of Integral based Digital Curvature Estimators in DGtal David Coeurjolly 1, Jacques-Olivier Lacaud 2, Jérémy Levallois 1,2 1 Université de Lyon, CNRS INSA-Lyon, LIRIS, UMR5205, F-69621,
More informationHash-Based Indexes. Chapter 11. Comp 521 Files and Databases Fall
Has-Based Indexes Capter 11 Comp 521 Files and Databases Fall 2012 1 Introduction Hasing maps a searc key directly to te pid of te containing page/page-overflow cain Doesn t require intermediate page fetces
More informationEXERCISES 6.1. Cross-Sectional Areas. 6.1 Volumes by Slicing and Rotation About an Axis 405
6. Volumes b Slicing and Rotation About an Ais 5 EXERCISES 6. Cross-Sectional Areas In Eercises and, find a formula for te area A() of te crosssections of te solid perpendicular to te -ais.. Te solid lies
More information5.4 Sum and Difference Formulas
380 Capter 5 Analtic Trigonometr 5. Sum and Difference Formulas Using Sum and Difference Formulas In tis section and te following section, ou will stud te uses of several trigonometric identities and formulas.
More informationImage Registration via Particle Movement
Image Registration via Particle Movement Zao Yi and Justin Wan Abstract Toug fluid model offers a good approac to nonrigid registration wit large deformations, it suffers from te blurring artifacts introduced
More information2.5 Evaluating Limits Algebraically
SECTION.5 Evaluating Limits Algebraically 3.5 Evaluating Limits Algebraically Preinary Questions. Wic of te following is indeterminate at x? x C x ; x x C ; x x C 3 ; x C x C 3 At x, x isofteform 0 xc3
More informationRead pages in the book, up to the investigation. Pay close attention to Example A and how to identify the height.
C 8 Noteseet L Key In General ON LL PROBLEMS!!. State te relationsip (or te formula).. Sustitute in known values. 3. Simplify or Solve te equation. Use te order of operations in te correct order. Order
More informationOur Calibrated Model has No Predictive Value: An Example from the Petroleum Industry
Our Calibrated Model as No Predictive Value: An Example from te Petroleum Industry J.N. Carter a, P.J. Ballester a, Z. Tavassoli a and P.R. King a a Department of Eart Sciences and Engineering, Imperial
More information15-122: Principles of Imperative Computation, Summer 2011 Assignment 6: Trees and Secret Codes
15-122: Principles of Imperative Computation, Summer 2011 Assignment 6: Trees and Secret Codes William Lovas (wlovas@cs) Karl Naden Out: Tuesday, Friday, June 10, 2011 Due: Monday, June 13, 2011 (Written
More informationCE 221 Data Structures and Algorithms
CE Data Structures and Algoritms Capter 4: Trees (AVL Trees) Text: Read Weiss, 4.4 Izmir University of Economics AVL Trees An AVL (Adelson-Velskii and Landis) tree is a binary searc tree wit a balance
More informationCSCE476/876 Spring Homework 5
CSCE476/876 Spring 2016 Assigned on: Friday, Marc 11, 2016 Due: Monday, Marc 28, 2016 Homework 5 Programming assignment sould be submitted wit andin Te report can eiter be submitted wit andin as a PDF,
More informationAVL Trees Outline and Required Reading: AVL Trees ( 11.2) CSE 2011, Winter 2017 Instructor: N. Vlajic
1 AVL Trees Outline and Required Reading: AVL Trees ( 11.2) CSE 2011, Winter 2017 Instructor: N. Vlajic AVL Trees 2 Binary Searc Trees better tan linear dictionaries; owever, te worst case performance
More informationCS211 Spring 2004 Lecture 06 Loops and their invariants. Software engineering reason for using loop invariants
CS211 Spring 2004 Lecture 06 Loops and teir invariants Reading material: Tese notes. Weiss: Noting on invariants. ProgramLive: Capter 7 and 8 O! Tou ast damnale iteration and art, indeed, ale to corrupt
More informationExperimental Studies on SMT-based Debugging
Experimental Studies on SMT-based Debugging Andre Sülflow Görscwin Fey Rolf Drecsler Institute of Computer Science University of Bremen 28359 Bremen, Germany {suelflow,fey,drecsle}@informatik.uni-bremen.de
More informationComputing geodesic paths on manifolds
Proc. Natl. Acad. Sci. USA Vol. 95, pp. 8431 8435, July 1998 Applied Matematics Computing geodesic pats on manifolds R. Kimmel* and J. A. Setian Department of Matematics and Lawrence Berkeley National
More informationJPEG Serial Camera Module. OV528 Protocol
JPEG Serial Camera Module OV528 Protocol LCF-23M1 32mmx32mm or 38mmx38mm LCF-23MA 32mm-38mm Default baudrate 9600bps~115200 bps Auto adaptive 9600bps~115200 bps Page 1 of 15 1.General Description OV528
More informationMAP MOSAICKING WITH DISSIMILAR PROJECTIONS, SPATIAL RESOLUTIONS, DATA TYPES AND NUMBER OF BANDS 1. INTRODUCTION
MP MOSICKING WITH DISSIMILR PROJECTIONS, SPTIL RESOLUTIONS, DT TYPES ND NUMBER OF BNDS Tyler J. lumbaug and Peter Bajcsy National Center for Supercomputing pplications 605 East Springfield venue, Campaign,
More informationA Novel QC-LDPC Code with Flexible Construction and Low Error Floor
A Novel QC-LDPC Code wit Flexile Construction and Low Error Floor Hanxin WANG,2, Saoping CHEN,2,CuitaoZHU,2 and Kaiyou SU Department of Electronics and Information Engineering, Sout-Central University
More informationOptimal In-Network Packet Aggregation Policy for Maximum Information Freshness
1 Optimal In-etwork Packet Aggregation Policy for Maimum Information Fresness Alper Sinan Akyurek, Tajana Simunic Rosing Electrical and Computer Engineering, University of California, San Diego aakyurek@ucsd.edu,
More informationNon-Interferometric Testing
NonInterferometric Testing.nb Optics 513 - James C. Wyant 1 Non-Interferometric Testing Introduction In tese notes four non-interferometric tests are described: (1) te Sack-Hartmann test, (2) te Foucault
More informationLecture 4: Geometry II
Lecture 4: Geometry II LPSS MATHCOUNTS 19 May 2004 Some Well-Known Pytagorean Triples A Pytagorean triple is a set of tree relatively prime 1 natural numers a,, and c satisfying a 2 + 2 = c 2 : 3 2 + 4
More informationeach node in the tree, the difference in height of its two subtrees is at the most p. AVL tree is a BST that is height-balanced-1-tree.
Data Structures CSC212 1 AVL Trees A binary tree is a eigt-balanced-p-tree if for eac node in te tree, te difference in eigt of its two subtrees is at te most p. AVL tree is a BST tat is eigt-balanced-tree.
More informationWrap up Amortized Analysis; AVL Trees. Riley Porter Winter CSE373: Data Structures & Algorithms
CSE 373: Data Structures & Wrap up Amortized Analysis; AVL Trees Riley Porter Course Logistics Symposium offered by CSE department today HW2 released, Big- O, Heaps (lecture slides ave pseudocode tat will
More information1 2 (3 + x 3) x 2 = 1 3 (3 + x 1 2x 3 ) 1. 3 ( 1 x 2) (3 + x(0) 3 ) = 1 2 (3 + 0) = 3. 2 (3 + x(0) 1 2x (0) ( ) = 1 ( 1 x(0) 2 ) = 1 3 ) = 1 3
6 Iterative Solvers Lab Objective: Many real-world problems of the form Ax = b have tens of thousands of parameters Solving such systems with Gaussian elimination or matrix factorizations could require
More informationAn Algorithm for Loopless Deflection in Photonic Packet-Switched Networks
An Algoritm for Loopless Deflection in Potonic Packet-Switced Networks Jason P. Jue Center for Advanced Telecommunications Systems and Services Te University of Texas at Dallas Ricardson, TX 75083-0688
More informationCryptanalysis of LOKI. Lars Ramkilde Knudsen. Aarhus Universitet Datalogisk Afdeling. Ny Munkegade. DK-8000 Aarhus C. Abstract
1 Cryptanalysis of LOKI Lars Ramkilde Knudsen Aarus Universitet Datalogisk Afdeling Ny Munkegade DK-8000 Aarus C. Abstract In [BPS90] Brown, Pieprzyk and Seberry proposed a new encryption primitive, wic
More informationz = x 2 xy + y 2 clf // c6.1(2)contour Change to make a contour plot of z=xy.
190 Lecture 6 3D equations formatting Open Lecture 6. See Capter 3, 10 of text for details. Draw a contour grap and a 3D grap of z = 1 x 2 y 2 = an upper emispere. For Classwork 1 and 2, you will grap
More informationProofs of Derivative Rules
Proos o Derivative Rules Ma 16010 June 2018 Altoug proos are not generally a part o tis course, it s never a ba iea to ave some sort o explanation o wy tings are te way tey are. Tis is especially relevant
More informationZernike vs. Zonal Matrix Iterative Wavefront Reconstructor. Sophia I. Panagopoulou, PhD. University of Crete Medical School Dept.
Zernie vs. Zonal Matrix terative Wavefront Reconstructor opia. Panagopoulou PD University of Crete Medical cool Dept. of Optalmology Daniel R. Neal PD Wavefront ciences nc. 480 Central.E. Albuquerque NM
More informationCESILA: Communication Circle External Square Intersection-Based WSN Localization Algorithm
Sensors & Transducers 2013 by IFSA ttp://www.sensorsportal.com CESILA: Communication Circle External Square Intersection-Based WSN Localization Algoritm Sun Hongyu, Fang Ziyi, Qu Guannan College of Computer
More informationMean Shifting Gradient Vector Flow: An Improved External Force Field for Active Surfaces in Widefield Microscopy.
Mean Sifting Gradient Vector Flow: An Improved External Force Field for Active Surfaces in Widefield Microscopy. Margret Keuper Cair of Pattern Recognition and Image Processing Computer Science Department
More informationAlternating Direction Implicit Methods for FDTD Using the Dey-Mittra Embedded Boundary Method
Te Open Plasma Pysics Journal, 2010, 3, 29-35 29 Open Access Alternating Direction Implicit Metods for FDTD Using te Dey-Mittra Embedded Boundary Metod T.M. Austin *, J.R. Cary, D.N. Smite C. Nieter Tec-X
More informationInvestigating an automated method for the sensitivity analysis of functions
Investigating an automated metod for te sensitivity analysis of functions Sibel EKER s.eker@student.tudelft.nl Jill SLINGER j..slinger@tudelft.nl Delft University of Tecnology 2628 BX, Delft, te Neterlands
More informationOn the Use of Radio Resource Tests in Wireless ad hoc Networks
Tecnical Report RT/29/2009 On te Use of Radio Resource Tests in Wireless ad oc Networks Diogo Mónica diogo.monica@gsd.inesc-id.pt João Leitão jleitao@gsd.inesc-id.pt Luis Rodrigues ler@ist.utl.pt Carlos
More informationEfficient and Not-So-Efficient Algorithms. NP-Complete Problems. Search Problems. Propositional Logic. DPV Chapter 8: NP-Complete Problems, Part 1
Efficient and Not-So-Efficient Algoritms NP-Complete Problems DPV Capter 8: NP-Complete Problems, Part 1 Jim Royer EECS November 24, 2009 Problem spaces tend to be big: A grap on n vertices can ave up
More information( )( ) ( ) MTH 95 Practice Test 1 Key = 1+ x = f x. g. ( ) ( ) The only zero of f is 7 2. The only solution to g( x ) = 4 is 2.
Mr. Simonds MTH 95 Class MTH 95 Practice Test 1 Key 1. a. g ( ) ( ) + 4( ) 4 1 c. f ( x) 7 7 7 x 14 e. + 7 + + 4 f g 1+ g. f 4 + 4 7 + 1+ i. g ( 4) ( 4) + 4( 4) k. g( x) x 16 + 16 0 x 4 + 4 4 0 x 4x+ 4
More informationwrobot k wwrobot hrobot (a) Observation area Horopter h(θ) (Virtual) horopters h(θ+ θ lim) U r U l h(θ+ θ) Base line Left camera Optical axis
Selective Acquisition of 3-D Information Enoug for Finding Passable Free Spaces Using an Active Stereo Vision System Atsusi Nisikawa, Atsusi Okubo, and Fumio Miyazaki Department of Systems and Human Science
More informationGrid Adaptation for Functional Outputs: Application to Two-Dimensional Inviscid Flows
Journal of Computational Pysics 176, 40 69 (2002) doi:10.1006/jcp.2001.6967, available online at ttp://www.idealibrary.com on Grid Adaptation for Functional Outputs: Application to Two-Dimensional Inviscid
More informationTilings of rectangles with T-tetrominoes
Tilings of rectangles wit T-tetrominoes Micael Korn and Igor Pak Department of Matematics Massacusetts Institute of Tecnology Cambridge, MA, 2139 mikekorn@mit.edu, pak@mat.mit.edu August 26, 23 Abstract
More informationAlgebra Area of Triangles
LESSON 0.3 Algera Area of Triangles FOCUS COHERENCE RIGOR LESSON AT A GLANCE F C R Focus: Common Core State Standards Learning Ojective 6.G.A. Find te area of rigt triangles, oter triangles, special quadrilaterals,
More information