DOCSLIB.ORG

Delta Functions and All That1 We Can Define the Step Function

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Delta Functions and All That1 We Can Define the Step Function Delta functions and all that1 We can de¯ne the step function (sometimes called the Heaviside step function), θ(x), by θ(x) = 0 for x < 0 (1) = 1 for x > 0 : (2) It is strictly not de¯ned for x = 0 but physicists like to de¯ne it to be symmetric, θ(0) = 1=2. One can think of this discontinuous function as a limit: 1 + tanh(x=a) θ(x) = Lim + ; a!0 2 this is easy to check if one remembers that the hyperbolic tangent assumes the values §1 as x ! §1 respectively. a goes to zero from the positive direction. This limiting procedure is consistent with θ(0) = 1=2. The step function is useful when something is switched on abruptly, say a voltage which is 0 for t < 0 and a constant above. In physical situations the step function is a convenient mathematical approximation. The (Dirac) delta function is strictly not a function but is a so-called distribution or a generalized function. The theory behind this was worked out by Laurent Schwartz who received the Fields Medal, the mathematics equivalent of the Nobel Prize. We will use the delta function carefully, con¯dent that the rigorous aspects have been checked. The delta function, ±(x ¡ a), is de¯ned by the property Z c2 dx f(x) ±(x ¡ a) = f(a) if a 2 (c1; c2) (3) c1 = 0 otherwise (4) for any suitably continuous function.2 Note that if a does not lie in the range of integration then the argument of the delta function does not vanish and therefore 1For a rigorous yet clear exposition see M. J. Lighthill, `Introduction to Fourier Analysis and Generalized Functions. 2From the point of view of physics, this de¯nition implies that the dimensions(units) of the delta function are the inverse of the dimensions of its argument. 1 it is always zero. In words, the delta function, when multiplied by a function and integrated, picks out the value of the function at the point where the argument of the delta function vanishes inside the range of the integral. If viewed naively as a function, the delta function x is zero when its argument is nonzero, i.e, when x 6= a and equal to 1 when x = a. We can think of the delta function as the derivative of the step function: d ±(x) = θ(x) ; (5) dx clearly the derivative vanishes for x 6= 0 and because of the ¯nite jump within an in¯nitesimal interval the derivative is in¯nite at x = 0. The delta function can also be viewed as a limit of a sequence of functions or the limit of a function; for example, : 1 ² ±(x) = lim (6) ²!0 ¼ x2 + ²2 1 x2 ¡ 2 = limσ!0 p e 2σ (7) 2¼σ and other such representations. Check in each case that when x 6= 0 taking the spec- i¯ed limit yields 0 and when x = 0 the limit yields in¯nity. How are the constants in the above limits chosen ? From the fundamental property of the delta function we have Z 1 dx ±(x) = 1 ¡1 and the constants in the limits above are chosen to obey this integral constraint. For example, (changing variables from x to y = x/σ in the integral below) we have Z 1 2 Z 1 1 ¡ x 1 ¡y2=2 dx f(x) p e 2σ2 = p dy f(σy) e ¡1 2¼σ 2¼ ¡1 and upon expanding f(σy) for small σ as f(0) + O(σ) we ¯nd that all the terms proportional to powers of σ vanish and we are left behind with f(0) which is precisely the de¯nition of a delta function. We also note that 1 sin ºx lim = ±(x) : (8) º!1 ¼ x From the study of Fourier integrals we ¯nd another of the more useful representa- tions of the delta function: Z 1 dk 0 ±(x ¡ x0) = eik(x¡x ) : (9) ¡1 2¼ 2 Equivalently, since the variables in question are dummy variables, Z 1 dx ei(k¡k0)x = 2¼ ±(k ¡ k0) : (10) ¡1 The material presented may be used in this course. I list below some properties of the delta function which are useful. All the results in which a delta function appears without being integrated over should be viewed as being multiplied (on both sides of the equation) by a a good function f(x) being integrated over a range of x, say (¡1; +1). The delta function is even : ±(x) = ±(¡x) ; (11) For a real constant c, we have 1 ±( c(x ¡ a) ) = ±(x ¡ a) : (12) j c j We also have, as can be shown with a bit of work, X 0 ±( f(x) ) = [ ±(x ¡ xi)= j f (xi) j ] (13) i 0 0 where f (x) = df=dx, f (xi) = df=dxjx=xi , and the fxig are the roots of the equation f(x) = 0 and the sum is over all the roots. For example, 1 ±(x2 ¡ a2) = [ ±(x ¡ a) + ±(x + a)] 2a for a > 0. One can (blithely) de¯ne the derivative of the delta function, ±0(x), even though it is not a function! Z c2 0 df dx ± (x) f(x) = ¡ jx=0 if 0 2 (c1; c2) (14) c1 dx i.e., the derivative of the delta function picks out the value of the derivative of the function at the point where the argument of the delta function vanishes apart from a negative sign. This is easily veri¯ed by using integration by parts. 3 Exercise: Generalize the result to the nth derivative of the delta function by re- peated integration by parts. Z 1 Answer: dx ±(n)(x) f(x) = (¡1)n f (n)(0) ¡1 where the superscript n denotes the nth derivative. Another obvious result is f(x) ±(x ¡ a) = f(a) ±(x ¡ a) : One can generalize delta functions de¯ned in one dimension above to two and three dimensions in Cartesian coordinates easily and in non-Cartesian coordinates with some care and work. The Jacobian appears naturally as is shown below. In two dimensions let ±(~r ¡ ~r0) = ±(x ¡ x0) ±(y ¡ y0) = D(r; θ) ±(r ¡ r0) ±(θ ¡ θ0) in polar coordinates. Note that D has dimensions of 1=L. We have Z d2r g(~r) ±(~r ¡ ~r0) = g(~r0) = g(r0; θ0) ­ ZZ = rdr dθ D(r; θ) ±(r ¡ r0) ±(θ ¡ θ0) g(r; θ) (15) and for this to be true for all g we have D = 1=r, i.e., it is the inverse of the Jacobian. In three dimensions we have 1 1 ±(~r ¡ ~r0) = ±(r ¡r0) ±(θ ¡θ0) ±(Á¡Á0) = ±(r ¡r0) ±(cos θ ¡ cos θ0) ±(Á¡Á0) r2 sin θ r2 where we have used ±(θ ¡ θ0) ±(cos θ ¡ cos θ0) = sin θ observing that in spherical polar coordinates sin θ is non-negative. For cylindrical coordinates 1 ±(~r ¡ ~r0) = ±(½ ¡ ½0) ±(Á ¡ Á0) ±(z ¡ z0) : ½ 4.
Load more

Recommended publications
  • Discontinuous Forcing Functions
    Math 135A, Winter 2012 Discontinuous forcing functions 1 Preliminaries If f(t) is defined on the interval [0; 1), then its Laplace transform is defined to be Z 1 F (s) = L (f(t)) = e−stf(t) dt; 0 as long as this integral is defined and converges. In particular, if f is of exponential order and is piecewise continuous, the Laplace transform of f(t) will be defined. • f is of exponential order if there are constants M and c so that jf(t)j ≤ Mect: Z 1 Since the integral e−stMect dt converges if s > c, then by a comparison test (like (11.7.2) 0 in Salas-Hille-Etgen), the integral defining the Laplace transform of f(t) will converge. • f is piecewise continuous if over each interval [0; b], f(t) has only finitely many discontinuities, and at each point a in [0; b], both of the limits lim f(t) and lim f(t) t!a− t!a+ exist { they need not be equal, but they must exist. (At the endpoints 0 and b, the appropriate one-sided limits must exist.) 2 Step functions Define u(t) to be the function ( 0 if t < 0; u(t) = 1 if t ≥ 0: Then u(t) is called the step function, or sometimes the Heaviside step function: it jumps from 0 to 1 at t = 0. Note that for any number a > 0, the graph of the function u(t − a) is the same as the graph of u(t), but translated right by a: u(t − a) jumps from 0 to 1 at t = a.
    [Show full text]
  • Laplace Transforms: Theory, Problems, and Solutions
    Laplace Transforms: Theory, Problems, and Solutions Marcel B. Finan Arkansas Tech University c All Rights Reserved 1 Contents 43 The Laplace Transform: Basic Definitions and Results 3 44 Further Studies of Laplace Transform 15 45 The Laplace Transform and the Method of Partial Fractions 28 46 Laplace Transforms of Periodic Functions 35 47 Convolution Integrals 45 48 The Dirac Delta Function and Impulse Response 53 49 Solving Systems of Differential Equations Using Laplace Trans- form 61 50 Solutions to Problems 68 2 43 The Laplace Transform: Basic Definitions and Results Laplace transform is yet another operational tool for solving constant coeffi- cients linear differential equations. The process of solution consists of three main steps: • The given \hard" problem is transformed into a \simple" equation. • This simple equation is solved by purely algebraic manipulations. • The solution of the simple equation is transformed back to obtain the so- lution of the given problem. In this way the Laplace transformation reduces the problem of solving a dif- ferential equation to an algebraic problem. The third step is made easier by tables, whose role is similar to that of integral tables in integration. The above procedure can be summarized by Figure 43.1 Figure 43.1 In this section we introduce the concept of Laplace transform and discuss some of its properties. The Laplace transform is defined in the following way. Let f(t) be defined for t ≥ 0: Then the Laplace transform of f; which is denoted by L[f(t)] or by F (s), is defined by the following equation Z T Z 1 L[f(t)] = F (s) = lim f(t)e−stdt = f(t)e−stdt T !1 0 0 The integral which defined a Laplace transform is an improper integral.
    [Show full text]
  • 10 Heat Equation: Interpretation of the Solution
    Math 124A { October 26, 2011 «Viktor Grigoryan 10 Heat equation: interpretation of the solution Last time we considered the IVP for the heat equation on the whole line u − ku = 0 (−∞ < x < 1; 0 < t < 1); t xx (1) u(x; 0) = φ(x); and derived the solution formula Z 1 u(x; t) = S(x − y; t)φ(y) dy; for t > 0; (2) −∞ where S(x; t) is the heat kernel, 1 2 S(x; t) = p e−x =4kt: (3) 4πkt Substituting this expression into (2), we can rewrite the solution as 1 1 Z 2 u(x; t) = p e−(x−y) =4ktφ(y) dy; for t > 0: (4) 4πkt −∞ Recall that to derive the solution formula we first considered the heat IVP with the following particular initial data 1; x > 0; Q(x; 0) = H(x) = (5) 0; x < 0: Then using dilation invariance of the Heaviside step function H(x), and the uniquenessp of solutions to the heat IVP on the whole line, we deduced that Q depends only on the ratio x= t, which lead to a reduction of the heat equation to an ODE. Solving the ODE and checking the initial condition (5), we arrived at the following explicit solution p x= 4kt 1 1 Z 2 Q(x; t) = + p e−p dp; for t > 0: (6) 2 π 0 The heat kernel S(x; t) was then defined as the spatial derivative of this particular solution Q(x; t), i.e. @Q S(x; t) = (x; t); (7) @x and hence it also solves the heat equation by the differentiation property.
    [Show full text]
  • Delta Functions and Distributions
    When functions have no value(s): Delta functions and distributions Steven G. Johnson, MIT course 18.303 notes Created October 2010, updated March 8, 2017. Abstract x = 0. That is, one would like the function δ(x) = 0 for all x 6= 0, but with R δ(x)dx = 1 for any in- These notes give a brief introduction to the mo- tegration region that includes x = 0; this concept tivations, concepts, and properties of distributions, is called a “Dirac delta function” or simply a “delta which generalize the notion of functions f(x) to al- function.” δ(x) is usually the simplest right-hand- low derivatives of discontinuities, “delta” functions, side for which to solve differential equations, yielding and other nice things. This generalization is in- a Green’s function. It is also the simplest way to creasingly important the more you work with linear consider physical effects that are concentrated within PDEs, as we do in 18.303. For example, Green’s func- very small volumes or times, for which you don’t ac- tions are extremely cumbersome if one does not al- tually want to worry about the microscopic details low delta functions. Moreover, solving PDEs with in this volume—for example, think of the concepts of functions that are not classically differentiable is of a “point charge,” a “point mass,” a force plucking a great practical importance (e.g. a plucked string with string at “one point,” a “kick” that “suddenly” imparts a triangle shape is not twice differentiable, making some momentum to an object, and so on.
    [Show full text]
  • An Analytic Exact Form of the Unit Step Function
    Mathematics and Statistics 2(7): 235-237, 2014 http://www.hrpub.org DOI: 10.13189/ms.2014.020702 An Analytic Exact Form of the Unit Step Function J. Venetis Section of Mechanics, Faculty of Applied Mathematics and Physical Sciences, National Technical University of Athens *Corresponding Author: [email protected] Copyright © 2014 Horizon Research Publishing All rights reserved. Abstract In this paper, the author obtains an analytic Meanwhile, there are many smooth analytic exact form of the unit step function, which is also known as approximations to the unit step function as it can be seen in Heaviside function and constitutes a fundamental concept of the literature [4,5,6]. Besides, Sullivan et al [7] obtained a the Operational Calculus. Particularly, this function is linear algebraic approximation to this function by means of a equivalently expressed in a closed form as the summation of linear combination of exponential functions. two inverse trigonometric functions. The novelty of this However, the majority of all these approaches lead to work is that the exact representation which is proposed here closed – form representations consisting of non - elementary is not performed in terms of non – elementary special special functions, e.g. Logistic function, Hyperfunction, or functions, e.g. Dirac delta function or Error function and Error function and also most of its algebraic exact forms are also is neither the limit of a function, nor the limit of a expressed in terms generalized integrals or infinitesimal sequence of functions with point wise or uniform terms, something that complicates the related computational convergence. Therefore it may be much more appropriate in procedures.
    [Show full text]
  • Appendix I Integrating Step Functions
    Appendix I Integrating Step Functions In Chapter I, section 1, we were considering step functions of the form (1) where II>"" In were brick functions and AI> ... ,An were real coeffi­ cients. The integral of I was defined by the formula (2) In case of a single real variable, it is almost obvious that the value JI does not depend on the representation (2). In Chapter I, we admitted formula (2) for several variables, relying on analogy only. There are various methods to prove that formula precisely. In this Appendix we are going to give a proof which uses the concept of Heaviside functions. Since the geometrical intuition, above all in multidimensional case, might be decep­ tive, all the given arguments will be purely algebraic. 1. Heaviside functions The Heaviside function H(x) on the real line Rl admits the value 1 for x;;;o °and vanishes elsewhere; it thus can be defined on writing H(x) = {I for x ;;;0O, ° for x<O. By a Heaviside function H(x) in the plane R2, where x actually denotes points x = (~I> ~2) of R2, we understand the function which admits the value 1 for x ;;;0 0, i.e., for ~l ~ 0, ~2;;;o 0, and vanishes elsewhere. In symbols, H(x) = {I for x ;;;0O, (1.1) ° for x to. Note that we cannot use here the notation x < °(which would mean ~l <0, ~2<0) instead of x~O, because then H(x) would not be defined in the whole plane R 2 • Similarly, in Euclidean q-dimensional space R'l, the Heaviside function H(x), where x = (~I>' .
    [Show full text]
  • On the Approximation of the Step Function by Raised-Cosine and Laplace Cumulative Distribution Functions
    European International Journal of Science and Technology Vol. 4 No. 9 December, 2015 On the Approximation of the Step function by Raised-Cosine and Laplace Cumulative Distribution Functions Vesselin Kyurkchiev 1 Faculty of Mathematics and Informatics Institute of Mathematics and Informatics, Paisii Hilendarski University of Plovdiv, Bulgarian Academy of Sciences, 24 Tsar Assen Str., 4000 Plovdiv, Bulgaria Email: [email protected] Nikolay Kyurkchiev Faculty of Mathematics and Informatics Institute of Mathematics and Informatics, Paisii Hilendarski University of Plovdiv, Bulgarian Academy of Sciences, Acad.G.Bonchev Str.,Bl.8 1113Sofia, Bulgaria Email: [email protected] 1 Corresponding Author Abstract In this note the Hausdorff approximation of the Heaviside step function by Raised-Cosine and Laplace Cumulative Distribution Functions arising from lifetime analysis, financial mathematics, signal theory and communications systems is considered and precise upper and lower bounds for the Hausdorff distance are obtained. Numerical examples, that illustrate our results are given, too. Key words: Heaviside step function, Raised-Cosine Cumulative Distribution Function, Laplace Cumulative Distribution Function, Alpha–Skew–Laplace cumulative distribution function, Hausdorff distance, upper and lower bounds 2010 Mathematics Subject Classification: 41A46 75 European International Journal of Science and Technology ISSN: 2304-9693 www.eijst.org.uk 1. Introduction The Cosine Distribution is sometimes used as a simple, and more computationally tractable, approximation to the Normal Distribution. The Raised-Cosine Distribution Function (RC.pdf) and Raised-Cosine Cumulative Distribution Function (RC.cdf) are functions commonly used to avoid inter symbol interference in communications systems [1], [2], [3]. The Laplace distribution function (L.pdf) and Laplace Cumulative Distribution function (L.cdf) is used for modeling in signal processing, various biological processes, finance, and economics.
    [Show full text]
  • On the Approximation of the Cut and Step Functions by Logistic and Gompertz Functions
    www.biomathforum.org/biomath/index.php/biomath ORIGINAL ARTICLE On the Approximation of the Cut and Step Functions by Logistic and Gompertz Functions Anton Iliev∗y, Nikolay Kyurkchievy, Svetoslav Markovy ∗Faculty of Mathematics and Informatics Paisii Hilendarski University of Plovdiv, Plovdiv, Bulgaria Email: [email protected] yInstitute of Mathematics and Informatics Bulgarian Academy of Sciences, Sofia, Bulgaria Emails: [email protected], [email protected] Received: , accepted: , published: will be added later Abstract—We study the uniform approximation practically important classes of sigmoid functions. of the sigmoid cut function by smooth sigmoid Two of them are the cut (or ramp) functions and functions such as the logistic and the Gompertz the step functions. Cut functions are continuous functions. The limiting case of the interval-valued but they are not smooth (differentiable) at the two step function is discussed using Hausdorff metric. endpoints of the interval where they increase. Step Various expressions for the error estimates of the corresponding uniform and Hausdorff approxima- functions can be viewed as limiting case of cut tions are obtained. Numerical examples are pre- functions; they are not continuous but they are sented using CAS MATHEMATICA. Hausdorff continuous (H-continuous) [4], [43]. In some applications smooth sigmoid functions are Keywords-cut function; step function; sigmoid preferred, some authors even require smoothness function; logistic function; Gompertz function; squashing function; Hausdorff approximation. in the definition of sigmoid functions. Two famil- iar classes of smooth sigmoid functions are the logistic and the Gompertz functions. There are I. INTRODUCTION situations when one needs to pass from nonsmooth In this paper we discuss some computational, sigmoid functions (e.
    [Show full text]
  • Lecture 2 ELE 301: Signals and Systems
    Lecture 2 ELE 301: Signals and Systems Prof. Paul Cuff Princeton University Fall 2011-12 Cuff (Lecture 2) ELE 301: Signals and Systems Fall 2011-12 1 / 70 Models of Continuous Time Signals Today's topics: Signals I Sinuoidal signals I Exponential signals I Complex exponential signals I Unit step and unit ramp I Impulse functions Systems I Memory I Invertibility I Causality I Stability I Time invariance I Linearity Cuff (Lecture 2) ELE 301: Signals and Systems Fall 2011-12 2 / 70 Sinusoidal Signals A sinusoidal signal is of the form x(t) = cos(!t + θ): where the radian frequency is !, which has the units of radians/s. Also very commonly written as x(t) = A cos(2πft + θ): where f is the frequency in Hertz. We will often refer to ! as the frequency, but it must be kept in mind that it is really the radian frequency, and the frequency is actually f . Cuff (Lecture 2) ELE 301: Signals and Systems Fall 2011-12 3 / 70 The period of the sinuoid is 1 2π T = = f ! with the units of seconds. The phase or phase angle of the signal is θ, given in radians. cos(ωt) cos(ωt θ) − 0 0 -2T -T T 2T t -2T -T T 2T t Cuff (Lecture 2) ELE 301: Signals and Systems Fall 2011-12 4 / 70 Complex Sinusoids The Euler relation defines ejφ = cos φ + j sin φ. A complex sinusoid is Aej(!t+θ) = A cos(!t + θ) + jA sin(!t + θ): e jωt -2T -T 0 T 2T Real sinusoid can be represented as the real part of a complex sinusoid Aej(!t+θ) = A cos(!t + θ) <f g Cuff (Lecture 2) ELE 301: Signals and Systems Fall 2011-12 5 / 70 Exponential Signals An exponential signal is given by x(t) = eσt If σ < 0 this is exponential decay.
    [Show full text]
  • Maximizing Acquisition Functions for Bayesian Optimization
    Maximizing acquisition functions for Bayesian optimization James T. Wilson∗ Frank Hutter Marc Peter Deisenroth Imperial College London University of Freiburg Imperial College London PROWLER.io Abstract Bayesian optimization is a sample-efficient approach to global optimization that relies on theoretically motivated value heuristics (acquisition functions) to guide its search process. Fully maximizing acquisition functions produces the Bayes’ decision rule, but this ideal is difficult to achieve since these functions are fre- quently non-trivial to optimize. This statement is especially true when evaluating queries in parallel, where acquisition functions are routinely non-convex, high- dimensional, and intractable. We first show that acquisition functions estimated via Monte Carlo integration are consistently amenable to gradient-based optimiza- tion. Subsequently, we identify a common family of acquisition functions, includ- ing EI and UCB, whose properties not only facilitate but justify use of greedy approaches for their maximization. 1 Introduction Bayesian optimization (BO) is a powerful framework for tackling complicated global optimization problems [32, 40, 44]. Given a black-box function f : , BO seeks to identify a maximizer X!Y x∗ arg maxx f(x) while simultaneously minimizing incurred costs. Recently, these strategies have2 demonstrated2X state-of-the-art results on many important, real-world problems ranging from material sciences [17, 57], to robotics [3, 7], to algorithm tuning and configuration [16, 29, 53, 56]. From a high-level perspective, BO can be understood as the application of Bayesian decision theory to optimization problems [11, 14, 45]. One first specifies a belief over possible explanations for f using a probabilistic surrogate model and then combines this belief with an acquisition function to convey the expected utility for evaluating a set of queries X.
    [Show full text]
  • Step and Delta Functions 18.031 Haynes Miller and Jeremy Orloff 1
    Step and Delta Functions 18.031 Haynes Miller and Jeremy Orloff 1 The unit step function 1.1 Definition Let's start with the definition of the unit step function, u(t): ( 0 for t < 0 u(t) = 1 for t > 0 We do not define u(t) at t = 0. Rather, at t = 0 we think of it as in transition between 0 and 1. It is called the unit step function because it takes a unit step at t = 0. It is sometimes called the Heaviside function. The graph of u(t) is simple. u(t) 1 t We will use u(t) as an idealized model of a natural system that goes from 0 to 1 very quickly. In reality it will make a smooth transition, such as the following. 1 t Figure 1. u(t) is an idealized version of this curve But, if the transition happens on a time scale much smaller than the time scale of the phenomenon we care about then the function u(t) is a good approximation. It is also much easier to deal with mathematically. One of our main uses for u(t) will be as a switch. It is clear that multiplying a function f(t) by u(t) gives ( 0 for t < 0 u(t)f(t) = f(t) for t > 0: We say the effect of multiplying by u(t) is that for t < 0 the function f(t) is switched off and for t > 0 it is switched on. 1 18.031 Step and Delta Functions 2 1.2 Integrals of u0(t) From calculus we know that Z Z b u0(t) dt = u(t) + c and u0(t) dt = u(b) − u(a): a For example: Z 5 u0(t) dt = u(5) − u(−2) = 1; −2 Z 3 u0(t) dt = u(3) − u(1) = 0; 1 Z −3 u0(t) dt = u(−3) − u(−5) = 0: −5 In fact, the following rule for the integral of u0(t) over any interval is obvious ( Z b 1 if 0 is inside the interval (a; b) u0(t) = (1) a 0 if 0 is outside the interval [a; b].
    [Show full text]
  • Expanded Use of Discontinuity and Singularity Functions in Mechanics
    AC 2010-1666: EXPANDED USE OF DISCONTINUITY AND SINGULARITY FUNCTIONS IN MECHANICS Robert Soutas-Little, Michigan State University Professor Soutas-Little received his B.S. in Mechanical Engineering 1955 from Duke University and his MS 1959 and PhD 1962 in Mechanics from University of Wisconsin. His name until 1982 when he married Patricia Soutas was Robert William Little. Under this name, he wrote an Elasticity book and published over 20 articles. Since 1982, he has written over 100 papers for journals and conference proceedings and written books in Statics and Dynamics. He has taught courses in Statics, Dynamics, Mechanics of Materials, Elasticity, Continuum Mechanics, Viscoelasticity, Advanced Dynamics, Advanced Elasticity, Tissue Biomechanics and Biodynamics. He has won teaching excellence awards and the Distinguished Faculty Award. During his tenure at Michigan State University, he chaired the Department of Mechanical Engineering for 5 years and the Department of Biomechanics for 13 years. He directed the Biomechanics Evaluation Laboratory from 1990 until he retired in 2002. He served as Major Professor for 22 PhD students and over 100 MS students. He has received numerous research grants and consulted with engineering companies. He now is Professor Emeritus of Mechanical Engineering at Michigan State University. Page 15.549.1 Page © American Society for Engineering Education, 2010 Expanded Use of Discontinuity and Singularity Functions in Mechanics Abstract W. H. Macauley published Notes on the Deflection of Beams 1 in 1919 introducing the use of discontinuity functions into the calculation of the deflection of beams. In particular, he introduced the singularity functions, the unit doublet to model a concentrated moment, the Dirac delta function to model a concentrated load and the Heaviside step function to start a uniform load at any point on the beam.
    [Show full text]