Top Banner
BHARATH EDUCATIONAL SOCIETY’S GROUP OF INSTITUTIONS: MADANAPALLI GOLDEN VALLEY COLLEGE OF ENGINEERING [ GVIC ] Branch: ECE-A ROLL NO-109M1A0416 ROLL NO-109M1A0435 Laplace Transform & Fourier transforms & It’s Applications in Engineering ABSTRACT:
157
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
  • BHARATH EDUCATIONAL SOCIETYS

    GROUP OF INSTITUTIONS:

    MADANAPALLI

    GOLDEN

    VALLEY COLLEGE OF

    ENGINEERING [GVIC]

    Branch: ECE-A

    ROLL NO-109M1A0416 ROLL NO-109M1A0435

    Laplace Transform & Fourier transforms & Its Applications in Engineering ABSTRACT:

  • Laplace transform is useful in all branches of Engineering, in probability theory ,in

    electrical and mechanical systems. These transform we are using often.

    Inverse Laplace transform is also useful in many domains

    In periodic function problems and also to deal with differential equations, they

    are all very useful.

    Fourier transforms are useful for electronics and communications students

    and to solve many problems in their subjects. Particularly Z-transforms have

    properties similar to that of Laplace .The main difference is Laplace transforms is

    defined for the functions of continuous variables and Z-transform defined for

    sequences.

    The partial differential equation in a boundary value problem gets transform to

    an ordinary differential equation by the application of Fourier transform.

    Introduction To Laplace Transform:-

    Laplace transform is useful since

    1. Particular solution is obtained without first determining the general solution.

    2. Non homogenous equations are solved without obtaining the complementary integral.

    3. Solution of mechanical or electrical problems involving discontinuous force function (RHS function

    F(x)) (or) periodic functions other than cosine and sine are obtained easily.

    4. System of DE, PDE and integral functions

    Before the advent of calculators and computers logarithms were extensively used to replace

    multiplication (or division) of two large numbers. The crucial idea which made the Laplace transform

    (L.T.). A very powerful technique is that it replaces operations of calculus by operations of algebra. For

    example, with application of Laplace transform to an initial value problem, consisting of an ordinary (or

    partial) differential equation (O.D.E) together with initial conditions (I.C.) is reduced to problem of

    solving an algebraic equation (with any given initial conditions automatically taken care) as shown in

    fig. below.

    The z-transform plays an important role in the field of communication engineering and control

    engineering at the stage of

    analy

    sis

    and

    repre

    senta

    tion of discrete time linear

    shift invariance system.

    Original space

    t-shape

    O.D.E

    + (particular) solution of D.E

    IC

  • :

    :

    The z-transform plays an important role in the field of communication engineering and control

    engineering at the stage of analysis and representation of discrete time linear shift invariance system.

    In mathematics ,the Laplace transform is widely used integral transform denoted L{f(t)}, it

    is a linear operator of a function f(t) with a real argument t(t0) that transforms to a function F(s) with

    a complex argument S. this transformation is essentially bijective for the majority of practical uses ; the

    respective pairs of f(t) and f(s) are matched in the tables. The Laplace transform has the useful property

    that many relationships and operations over the originals of f(t) correspond to similar relationships and

    operations over the images f(s)[1] The Laplace transform has many important applications in his work on

    probability theory.

    The Laplace transform is related to the Fourier transform, but whereas the Fourier

    transform resolves a function or signal into its modes of vibrations, the Laplace resolves a function into

    its moments like the Fourier transform, the Laplace transform is used for solving differential and integral

    equations. In Physics and engineering, it is used for analysis of linear time-in variant systems such as

    electrical circuits, harmonic oscillations, optical devices, and mechanical systems.

    Given a simple mathematical or functional description of an input or output of a system, the Laplace

    transform provides an alternative functional description that often simplifies the process of analyzing

    the behavior of the system, or in synthesizing a new system based on a set of qualification

    INTRODUCTION OF FOURIER TRANSFORMS :-

    It may seem unusual that we begin a course on Geodynamics by reviewing Fourier transforms and

    Fourier series. However you will that Fourier analysis is used in almost every aspect of the subject

    ranging from solving linear differential equation to depending computer modules, to the processing and

    analysis of data .we wont be using Fourier analysis for first few lectures, but Ill introduce the concepts

    today so that people who are use familiar with the topic can have time for review. In the first few

    lectures, Ill also discuss plate tectonics. I imagine that students with physics and math backgrounds may

    Laplace

    Transform

    Image space

    s-space

    Algebraic solution of AE

    (or) subsidiary (by pure algebraic manipulations)

    Equation (AE)

  • have to spend some time reviewing plate tectonics. Hopefully everyone will busy the first two weeks

    next class Ill give homework assignment involving Fourier transforms. Later well have a short quiz on

    plate tectonics

    FOURIER TRANSFORMS:-

    Definitions of Fourier transforms:-

    The 1-dimensional Fourier transform is defined as:

    f(k) = ()2

    Forward Transform

    f(k) = ()2

    Inverse Transform

    Where x is distance &

    K is wave number

    Where k=1

    and is wavelength.

    These equations are commonly written in terms of time t and frequency v where = 1/ and T is the

    period.

    The 2-dimensional Fourier transform is defined as:-

    F(k) = ()2( .)2

    F(x) = ()2( .)2

    Where X=(x,y) is the position vector

    K = (Kx,Ky) is the wave number vector

    & (k.x) = kxx+kyy

    Why use Fourier transforms on a nearly special earth?

    It you have taken Geomagnetism or global seismology, we were taught to expand a function of

    latitude and longitude in spherical harmonics. Later in the course we will also use spherical harmonics to

    represent large scale variations in the gravity field and to represent viscous mantle flow. However,

    throughout the course we will dealing with problems related to the crust and lithosphere. In these cases

    a flat earth approximation is both adequate and practical for the following reasons:

  • Cartesian geometry is good approximation consider a small patch of crust or litho sphere on the surface

    of a sphere. If the area of patch is A is much less than the area of the earth and thickness L of the patch

    is much less than the radius of the earth Re, then the Cartesian geometry will be adequate

    A

  • Fourier sine and cosine transforms:- Any function f(x) can be decomposed into add O(x) and even E(x) components

    f(x)=E(x)+O(x)

    E(x =1

    2[f(x)+f(-x)]

    O(x) =1

    2[f(x)-f(-x)]

    f(k)= ()2

    =

    f(k) = cos 2

    s in 2

    odd part cancels even part cancels

    Fourier Transform :- Fourier series can be generalized to complex numbers, and further transform

    Forward Fourier transform

    f(k) = ()2

    *inverse Fourier transform

    f(x) = ()2

    Note = cos(x)+isin(x)

    Fourier Transform:- Fourier transform maps a time series.( e.g. : Audio samples) into the series of frequencies(their

    amplitudes and phases) that composed the time series.

    Inverse Fourier transform maps the series of frequencies (their amplitudes and phases) back

    into the corresponding time series.

    The two functions and inverses of each other.

    Discrete Fourier transform If we wish to find the frequency spectrum of a function that we have sampled, the

    continuous Fourier transform is not so useful.

    We need a discreet version.

    Discreet Fourier transforms.

  • Forward DFT:

    =

    The complex numbers f0.fn are transformed into complex numbers f0.fn

    Inverse DFT :

    =

    The complex numbers f0..fn are transformed into complex numbers f0fn.

    DFT Example Interpreting a DFT can be slightly difficult, because the DFT for really data includes complex

    numbers.

    Basically :The magnitude of the complex number for a DFT component is the power at that

    frequency.

    The phase of the wave form can be determined from the relatives values of the real

    and imaginary coefficients

    Also both POSITIVE and NEGATIVE frequencies show up.

  • Fourier Transform :- The continuous Fourier transform is evaluating the bilateral Laplace transform with

    complex argument.

    S = iw (or) S = 2fi:

    f(w) =F{f(t)}

    = =

    = f(t)dt

    This expression excludes the scaling factor , which is often included in definitions of the Fourier

    transform. This relationship between the Laplace and Fourier transform is often used to determine the

    frequency spectrum of a single or dynamical system.

    Applications:-

    The solution of a IBVP consisting of a partial differential equation together with if the

    boundary & initial conditions can be solved by the Fourier transform method. If the boundary

    conditions are of the Dirichlet type where the function value is prescribed on the boundary,

    then the Fourier sine transform is used .If the boundary conditions are of the Neumann type

    where the derivative of function is prescribed on the boundary, then Fourier cosine transform is

    applied . In either case, the P.D.E reduces to an O.D.E in Fourier transform which is solved then

    the inverse Fourier sine (or cosine) transforms will give the solution to the problem.

    Fourier transforms are widely used to solve various boundary values problems

    of engineering such as

    1. Vibrations of a string,

    2. Conductions of heat ,

    3. Oscillations of electric beam ,

    4. Transmission lines etc.,

    conclusion:- The partial differential equation in a boundary value problem gets transformed to an ordinary differential equation by the application of Fourier sine and cosine transform are useful when

    f(x) is defined ina finite interval for say x in o

  • transformation or Fourier transform is an integral transform which is used to solve the partial

    differential equations. The advantages of applying Fourier transform or any integral transform is to

    reduce the number of independent variables by one. Fourier transforms are widely used to solve various

    boundary value problems of engineeringsuch as vibrations of a string, conductions of heat, oscillations

    of an electric beam, transmission lines.

    Laplace Transform:-

    Definition:- Let f(t) be the given function defined for all t 0. Laplace transforms of f(t) denoted by L{f(t)}

    or simply L{f} is defined as

    L{f(t)} =

    0 = ()

    L is known as Laplace transform operator. The original given function f(t) known as determining function depends on t, while the new function to be determined f(s), called as generating function to be depends only on S (because the improper integral on the R.H.S of (1) is integrated with respect to T ). f(s) in equation (1) is known as the Laplace transform of f(t). Eq. (1) is known as direct transform or simply transform, in which f(t) is given and f(s) is to be determined. Thus, Laplace transform transforms one class of complicated function f(t) to produce another class of simpler functions F(s).

    APPLICATIONS:- Laplace transforms is very useful in obtaining solutions of linear differential equations , both ordinary and partial , solutions of system of simultaneous differential equations , solution of integral equations , solutions of linear difference equations and in the evaluation of definite integrals. In the PHYSICS and ENGINEERINGLaplace transforms are widely used for analysis of linear time variantsystems such as

    1. Electric circuits, 2. Harmonic oscillators, 3. Optical devices,& 4. Mechanical systems. 5. In the analysis, the Laplace transform is often interpreted as a transformation from the

    TIME-DOMAIN , in which inputs and outputs are functions of time , to the frequency domain , where the same inputs and outputs are functions of complex angular frequency , in radians per unit time

    6. Laplace transforms are used in the Area of digital signal processing and digital filters

    APPLICATION OF THE LAPLACE TRANSFORM IN CIRCUIT ANALYSIS:- The Laplace transform is an attractive tool in circuit analysis . it transforms a set of linear constant - coefficient differential equations into a set of linear polynomial equations . It automatically

  • introduces into the polynomial equations the initial values of the current and voltage variables. In the circuit analysis, we can develop the s-domain circuit models for various elements and s-domain equations can be written directly. CIRCUIT ELEMENTS IN THE S-DOMAIN :-

    We take the Laplace transform of the time domain equation . This gives an algebraic relation between s-domain current and voltage . The dimensions of a transformed voltage is volt-seconds and the dimensions of a transformed current is ampere-seconds .

    ADVANTAGES:- 1. With the application of L.T. particular solution of differential equation(D.E.) is obtained directly

    without the necessity of first determining general solution and then obtaining the particular solution (by substitution of initial condition ).

    2. L.T. solves non-homogenous D.E. without the necessity of first solving the corresponding homogenous D.E.

    3. L.T. is applicable not only to continuous functions but also to piecewise continuous functions, complicated periodic functions, step functions and impulse functions.

    SUFFICIENT CONDITIONS FOR THE EXISTENCE OF LAPLACE

    TRANSFORMS OF f(t):- The L.T. of f(t) exists i.e., the improper integral in the R.H.S. of (1) converges (has a finite valve ) when the following sufficient conditions are satisfied :- (a). f(t) is piecewise (or section ally ) continuous

    i.e., f(t) is continuous in every sub interval and has finite limits at end points of each of these sub intervals and

    (b). f(t) is of exponential order of i.e., there exists M, such that |f(t)| 2

    3 = finite,

    f(t) = t2 is of exponential order say 3.

    ex:- since lim>

    2

    = not finite,

    f(t) = 2 is not of exponential order.

    NOTE:- Above conditions (a) and (b) are not necessary conditions.

    Conclusion:- The partial differential equation in a boundary value problem gets transformed to an ordinary differential equation by the application. Laplace Transform can be used to solve linear differential equations with constant coefficients. The advantages by using Laplace transform is that the particular solution can be obtained for given initial condition with obtaining the general solution. The formulae to be used in applying Laplace transform to solve differential equations are L [ (t)] =s (s) f(0)

  • L [ (t)] = 2 (s) s f(0) - (0)

    L [ (t)] = 3 (s) 2f (0) s (0) - (0)

    While apply the above formulae, replace , replace f(t) by Y(x), (s) by (s).

  • 2. Laplace Transforms and their Applications

    To

    Differential Equations

    Authors:-

    K.Ravi kiran, E.E.E, MITS

    K.B.Venu gopal, E.E.E, MITS

    Email id [email protected]

    [email protected]

    Contents:-

    1. Definition of Laplace Transforms

    2.Conditions for Existance

    3. Applications of Differential Equations

    4. Problem solving

    Abstract:-

    The main idea behind Laplace Transformation is ordinary differential equations with constant coefficients can be easily solved by the laplace transform method, with out

    the necessity of first finding the general solution and then evaluating the arbitrary

    constants. This method is, in general, shorter than our earlier method and is especially

    suitable to obtain the solution of linear non-homogeneous ordinary differential equations

    with constant coefficients.

  • INTRODUCTION TO THE LAPLACE TRANSFORM METHOD

    The Laplace Transform method is a technique for solving linear differential

    equations with initial conditions. It is commonly used to solve electrical circuit

    and systems problems.

    What is a Transform Method?

    The simplest way to describe a transform method is to consider an example.

    Suppose we wish to compute the product of VI and XIV, both Roman numerals,

    and express the answer as a Roman numeral. Unless you are a Roman(!), the

    first thing to do is transform the Roman numerals to Arabic numerals. VI is 6

    and XIV is 14. The transformed problem is: compute the product of 6 and 14. We

    can all do this! The solution to the transformed problem is 84. We then convert

    the solution of the transformed problem to the solution to original problem. 84 in

    Roman numerals is LXXXIV. This last step is called the inverse transformation.

    The following diagram summarizes what we have done.

    Why use a transform method? Some problems are difficult to solve directly.

    With a transform method, the hope is that the transformed problem is easy

    to solve. That is certainly the case for the simple example above. One must

    also take into account the difficulty of transforming the original problem

    and inverse transforming the solution to the transformed problem.

    DEFINITION OF LAPLACE TRANSFORM

  • CONDITIONS FOR THE EXISTENCE OF LT

    While finding the laplace transforms of elementary functions, it can be

    noticed that the integral exists under certain conditions, such as s > 0 or s > a etc.

    in general, the functions f(t) must satisfy the following conditions for the

    existence of the laplace transform.

    (i) The function f(t) must be piece- wise continuous or sectionally

    continuous in any limited interval 0 < a t b.

    (ii) The function f(t) is of exponential order.

    APPLICATION TO DIFFERENTIAL EQUATIONS

    Consider the linear differential equation with constant coefficients

    under the initial conditions

    The Laplace transform directly gives the solution without going through the general solution.

    The steps to follow are:

    (1)

    Evaluate the Laplace transform of the two sides of the equation (C);

  • (2)

    ;

    (3)

    After algebraic manipulation, write down

    ;

    (4)

    Make use of the properties of the inverse Laplace transform , to find the solution y(t).

    Example: Find the solution of the IVP

    ,

    where

    .

    Solution: Let us follow these steps:

    (1)

    We have

    ;

    (2)

    Using properties of Laplace transform, we get

    ,

    where . Since , we get

    ;

  • (3)

    Inverse Laplace:

    Using partial decomposition technique we get

    ,

    which implies

    Since

    ,

    which gives

    ,

    and

    Hence,

    PROBLEM SOLVING

    Solve the following initial value problem by using laplace transform 4y"+2 y = 0, y (0) =2, y'(0) =0. Solution: Given equation is 4y"+2 y = 0 Taking laplace transform on both sides, we have 4L,y"-+2L{y}=0

    i.e., 4[s2L{y}-sy(0)- y'(0)++2L[y]=0

    => 4[s2 L{y}-2s] + 2L{y}=0

  • i.e., L{y}[4s2+2] = 2s or L{y}= 2s

    4s2+2

    Y=(1/2) L-1 s =(1/2) cos t/2 S2+ 2/4

  • 3. FOURIER TRANSFORMS AND ITS APPLICATIONS IN ENGINEERING

    FIELDS

    G.Sneha Geetha,

    I B.Tech, EEE,

    SDIT, Nandyal,

    E-mail: [email protected]

    S.Talat Misba,

    I B.Tech, EEE,

    SDIT, Nandyal,

    E-mail: [email protected]

    Abstract:

    The Fourier transform is a mathematical operation that decomposes a signal into its

    constituent frequencies. Thus the Fourier transform of a musical chord is a mathematical

    representation of the amplitudes of the individual notes that make it up. The original signal

    depends on time, and therefore is called the time domain representation of the signal, whereas

    the Fourier transform depends on frequency and is called the frequency domain representation

    of the signal. The term Fourier transform refers both to the frequency domain representation of

    the signal and the process that transforms the signal to its frequency domain representation .In

    this paper Fourier transform and its applications in a step by step procedure.

    Introduction:

    In mathematical terms, the Fourier transform transforms one complex-valued function of

    a real variable into another. In effect, the Fourier transform decomposes a function into

    oscillatory functions. The Fourier transform and its generalizations are the subject of Fourier

    analysis. In this specific case, both the time and frequency domains are unbounded linear

    continua. It is possible to define the Fourier transform of a function of several variables, which is

    important for instance in the physical study of wave motion and optics. It is also possible to

    generalize the Fourier transform on discrete structures such as finite groups. The efficient

    computation of such structures, by fast Fourier transform, is essential for high-speed computing.

    The motivation for the Fourier transform comes from the study of Fourier series. In the

    study of Fourier series, complicated functions are written as the sum of simple waves

    mathematically represented by sines and cosines. Due to the properties of sine and cosine it is

  • possible to recover the amount of each wave in the sum by an integral. In many cases it is

    desirable to use Euler's formula, which states that e2i

    = cos 2 + i sin 2, to write Fourier

    series in terms of the basic waves e2i

    . This has the advantage of simplifying many of the

    formulas involved and providing a formulation for Fourier series that more closely resembles the

    definition followed in this article. This passage from sines and cosines to complex exponentials

    makes it necessary for the Fourier coefficients to be complex valued. The usual interpretation of

    this complex number is that it gives both the amplitude (or size) of the wave present in the

    function and the phase (or the initial angle) of the wave. This passage also introduces the need

    for negative "frequencies". If were measured in seconds then the waves e2i and e2i would

    both complete one cycle per second, but they represent different frequencies in the Fourier

    transform. Hence, frequency no longer measures the number of cycles per unit time, but is

    closely related.

    There is a close connection between the definition of Fourier series and the Fourier transform for

    functions which are zero outside of an interval. For such a function we can calculate its Fourier

    series on any interval that includes the interval where is not identically zero. The Fourier

    transform is also defined for such a function. As we increase the length of the interval on which

    we calculate the Fourier series, then the Fourier series coefficients begin to look like the Fourier

    transform and the sum of the Fourier series of begins to look like the inverse Fourier transform.

    To explain this more precisely, suppose that T is large enough so that the interval [T/2,T/2]

    contains the interval on which is not identically zero. Then the n-th series coefficient cn is

    given by:

    Comparing this to the definition of the Fourier transform it follows that since

    (x) is zero outside [T/2,T/2]. Thus the Fourier coefficients are just the values of the Fourier

    transform sampled on a grid of width 1/T. As T increases the Fourier coefficients more closely

    represent the Fourier transform of the function.

    Under appropriate conditions the sum of the Fourier series of will equal the function . In other

    words can be written:

    where the last sum is simply the first sum rewritten using the definitions n = n/T, and

    = (n + 1)/T n/T = 1/T.

    This second sum is a Riemann sum, and so by letting T it will converge to the integral for

    the inverse Fourier transform given in the definition section. Under suitable conditions this

    argument may be made precise (Stein & Shakarchi 2003).

  • In the study of Fourier series the numbers cn could be thought of as the "amount" of the wave in

    the Fourier series of . Similarly, as seen above, the Fourier transform can be thought of as a

    function that measures how much of each individual frequency is present in our function , and

    we can recombine these waves by using an integral (or "continuous sum") to reproduce the

    original function.

    The following images provide a visual illustration of how the Fourier transform measures

    whether a frequency is present in a particular function. The function depicted

    oscillates at 3 hertz (if t measures seconds) and tends quickly to 0.

    This function was specially chosen to have a real Fourier transform which can easily be plotted.

    The first image contains its graph. In order to calculate we must integrate e2i(3t)(t). The

    second image shows the plot of the real and imaginary parts of this function. The real part of the

    integrand is almost always positive, this is because when (t) is negative, then the real part of

    e2i(3t)

    is negative as well. Because they oscillate at the same rate, when (t) is positive, so is the

    real part of e2i(3t)

    . The result is that when you integrate the real part of the integrand you get a

    relatively large number (in this case 0.5). On the other hand, when you try to measure a

    frequency that is not present, as in the case when we look at , the integrand oscillates

    enough so that the integral is very small. The general situation may be a bit more complicated

    than this, but this in spirit is how the Fourier transform measures how much of an individual

    frequency is present in a function (t).

    Applications in Engineering Fields:

    Analysis of differential equations

    Fourier transforms and the closely related Laplace transforms are widely used in solving

    differential equations. The Fourier transform is compatible with differentiation in the following

    sense: if f(x) is a differentiable function with Fourier transform , then the Fourier transform

    of its derivative is given by . This can be used to transform differential equations into

    algebraic equations. Note that this technique only applies to problems whose domain is the

    whole set of real numbers. By extending the Fourier transform to functions of several variables

    partial differential equations with domain Rn can also be translated into algebraic equations.

    Fourier transform spectroscopy

    The Fourier transform is also used in nuclear magnetic resonance (NMR) and in other

    kinds of spectroscopy, e.g. infrared (FTIR). In NMR an exponentially-shaped free induction

    decay (FID) signal is acquired in the time domain and Fourier-transformed to a Lorentzian line-

    shape in the frequency domain. The Fourier transform is also used in magnetic resonance

    imaging (MRI) and mass spectrometry.

  • One of the most basic tasks in spectroscopy is to characterize the spectrum of a light

    source: How much light is emitted at each different wavelength. The most straightforward way

    to measure a spectrum is to pass the light through a monochromator, an instrument that blocks all

    of the light except the light at a certain wavelength (the un-blocked wavelength is set by a knob

    on the monochromator). Then the intensity of this remaining (single-wavelength) light is

    measured. The measured intensity directly indicates how much light is emitted at that

    wavelength. By varying the monochromator's wavelength setting, the full spectrum can be

    measured. This simple scheme in fact describes how some spectrometers work.

    Fourier transform spectroscopy is a less intuitive way to get the same information. Rather

    than allowing only one wavelength at a time to pass through to the detector, this technique lets

    through a beam containing many different wavelengths of light at once, and measures the total

    beam intensity. Next, the beam is modified to contain a different combination of wavelengths,

    giving a second data point. This process is repeated many times. Afterwards, a computer takes all

    this data and works backwards to infer how much light there is at each wavelength.

    To be more specific, between the light source and the detector, there is a certain

    configuration of mirrors that allows some wavelengths to pass through but blocks others (due to

    wave interference). The beam is modified for each new data point by moving one of the mirrors;

    this changes the set of wavelengths that can pass through.

    Conclusion:

    This paper mainly deals with the introduction of Fourier Transforms and its applications

    to Engineering Fields.

  • APPLICATIONS OF LAPLACE TRANSFORMS

    1. V.SRAVANTHI 1stECE, AITS, [email protected]

    2.D.VIJITHA 1stECE, AITS,Rajampet [email protected]

    ABSTRACT

    In mathematics the Laplace transformation is a widely used integral transform. Laplace

    Transform has many important applications throughout the sciences. It is named for pierre-simon

    Laplace who introduced the transform in his work on theory. The Laplace transform has the

    useful property that many relationships and operations over the originals correspond to simpler

    relationships and operations over the images. In this paper the initial value problem has been

    arrived and solved by showing pictorial representations in a procedural way.

    INTRODUCTION

    The knowledge of Laplace transform is an essential part of mathematics required by engineers

    and scientists. The Laplace Transforms is an excellent tool for solving linear differential

    equations with given initial values of an unknown function and its derivatives without the

    necessity of first finding the general solution (complementary function +particular integral) and

    then evaluating from it the particular solution satisfying the given conditions. The technique is

    useful to solve some partial differential equations as well. This is a powerful tool in diverse

    fields of engineering.

    What are Laplace Transforms?

    A Laplace Transform is a type of integral transform.

    Let f(t) be a function defined for all positive values of t. Then Laplace Transform of f(t),

    denoted by L{f(t)} is defined by

  • 0f(t)dt = F(s)

    Plug one function in and get another function out.

    The new function is in a different domain when

    0f(t)dt = F(s)

    F(s) is the Laplace Transform of f(t)

    Write L {f (t)} =F(s),

    L{y(t)}=Y(s),

    L{x(t)}=X(s)etc

    L{f(t)}=

    0f(t)dt it is a Laplace Transform and it can be written as

    f(t)=L-1

    0f(t)dt So that the function f(t) is said to be Inverse Laplace Transform.

    STANDARD FORMULAE:

    1. L{1}=1/s

    2. L{}=1/s-a

    3. L{ }=1/s+a

    4. L{coshat}=s/s2-a2

    5. L{sinhat}=a/s2-a2

    6. L{cosat}=s/ s2+a2

    7. L{sinat}=a/ s2+a2

    8. L{t}=1/s2

  • To what end does one use Laplace Transforms?

    We can use Laplace Transforms to turn an initial value problem.

    solve for Y(t) into an algebraic equation

    solve for Y(s)

    Laplace Transforms are particularly effective on differential equation with forcing functions that

    are piece-wise like the head wise function, and other functions that turn on and of

    " 3 ' 4 ( 1)

    (0) 1, '(0) 2

    y y y t u t

    y y

    2

    2 1( )*( 3 4) ( 1) ss

    s eY s s s s

  • Then

    If you solve the algebraic equation

    and find the inverse Laplace Transform of the solution Y(s), you have the solution to the I.V.P

    The Inverse Laplace Transform of

    is

    2

    2 2

    ( 1) ( 1)( )

    ( 3 4)

    s ss s e eY s

    s s s

    2

    2 2

    ( 1) ( 1)( )

    ( 3 4)

    s ss s e eY s

    s s s

    4 43 32 15 80 4 16

    4325 5

    ( ) ( 1)( + ( ) )

    ( )( ( ) )

    t tee

    t t

    y t u t e e t

    u t e e

  • Thus

    is the solution of Initial value problem

    REFERENCES

    1. Engineering Mathematics volume 1by Dr.T.K.Iyengar, Dr.B.KrishnaGandhi, S.Ranganatham,

    M.V.S.S.N.Prasad.

    2. Text Bookof Engineering Mathematics B.V.Ramana.

    3. TextBookof Engineering Mathematics Thomson Book collection.

    4 43 32 15 80 4 16

    4325 5

    ( ) ( 1)( + ( ) )

    ( )( ( ) )

    t tee

    t t

    y t u t e e t

    u t e e

  • 5. APPLICATIONS OF LAPLACE TRANSFORMS

    Y. SILPA, I B.TECH (EEE)

    Sri Sai Institute of Science & Technology, Rayachoti.

    ABSTRACT Laplace transform or Laplace transformation is a method for solving linear differential

    equation arising in physics and engineering. It reduces the problem of solving a differential

    equation to an algebraic problem.

    Let K(s,t) be a function of two variables s, and t, where s is a parameter (may be real or

    complex) independent of t. The function F(s) defined by the integral (assumed to be

    convergent).

    sFdttftsK

    ),

    is called the Integral transform of the function f(t) and denoted by T{f(t)}. The function

    K(s, t) is called the kernal of the transformation.

    If the kernal K(s, t) is defined as

    0

    00,

    tfore

    tfortsK

    st

    then sFdttfe st

    0

    (1)

    The function F(s) defined by the integral (1) is called the Laplace transform of the

    function f(t) and is also denoted by

    L{f(t)} or F(s)

    Thus Laplace transform is a function of a new variable s given by (1)

    A semi-infinite insulated bar (x>0) which is initial at constant temperature (To > 0) and in

    which the end is held at a temperature of 0oC is considered for presentation of the solution. The

    solution for the problem deals with the determination of the temperature at any point of the semi-

    infinite insulated bar.

    INTRODUCTION

    They are many partial differential equations problems in engineering cases, which their

    quantities vary with the time. To solve these problems Laplace transforms, as a powerful

    technique can be used to transform the original differential equation into integral algebraic

  • expression. For mechanical engineering problem, Laplace transform is addressed to analyze the

    linear time invariant such as harmonic oscillators in vibration analysis and various problems in

    mechanical systems. Basically Laplace transforms changes integral and differential equations

    into polynomial equations.

    FORMATION OF THE PROBLEM

    A semi-infinite insulated bar (x > 0) which is initially at constant temperature (To> 0) and

    in which the end is hold at a temperature at 0oC is considered.

    We are to solve the diffusion equation

    2

    22

    x

    wc

    t

    w

    (2)

    Subject to the initial and boundary conditions

    w(x, 0) = To

    w(0, t) = 0 (3)

    w(x, t) 0 as x

    Applying the Laplace transform on both sides of equation 92)

    2

    22

    x

    wcL

    t

    wL

    t

    wL

    x

    wcL

    2

    22

    t

    wL

    x

    wLc

    2

    22

    oTsxwssxwdx

    dc ,,

    2

    22

    Equation (4) in ordinary differential equation

    O x

  • We have

    x

    c

    sx

    c

    s

    esBesAsxwFC

    ,:.

    s

    TsxwPI o,:

    Solution is s

    TesBesAsxw o

    xc

    sx

    c

    s

    , (5)

    Evidently B(s) = 0 from the third condition

    s

    TesAsxw o

    xc

    s

    ,

    Since w(0, t) = 0, we have w(0, s) = 0

    0s

    TsA o

    and so s

    c

    x

    oo es

    T

    s

    Tsxw

    , (6)

    Inverse Laplace Transforms:

    111

    sL

    t

    aerf

    t

    aerfc

    s

    eL

    sa

    21

    2

    1

    where erf is the error function.

    Applying the Inverse Laplace Transforms on both sides of equation (6)

    s

    c

    x

    oo es

    T

    s

    TLsxwL 11 ,

    s

    eTL

    s

    TLtxw

    sc

    x

    oo 11,

    s

    eLT

    sLTtxw

    sc

    x

    oo

    11 1,

    t

    cxerfTTtxw oo

    2

    /11,

  • tc

    xerfTTTtxw ooo

    2,

    tc

    xerfTtxw o

    2,

    CONCLUSION:

    The temperature of any the semi-infinite insulated at different points by calculated with the help

    of Laplace transformations.

  • 6. LAPLACE TRANSFORMS

    NAME :D.SUSHMITHA

    BRANCH:EEE

    DEPARTMENT OF HUMANITIES

    Email: [email protected]

    Phone no: 9490059399

    2) NAME :P. RAJYALAKSHMI

    BRANCH: EEE

    DEPARTMENT OF HUMANITIES

    Email: [email protected]

    ABSTRACT

    In this paper we here analyzed different electrical circuits here we have used laplace transform

    technique to save these electrical circuits. We have analyzed different parameters namely

    capacitor,resistor,battery,transistors, in this paper the solved problems are represented graphically. In

    this the applications of laplace transforms also discussed. The terms explained in this are inverse laplace

    transformations, applications of laplace transforms, table of laplace transformations advantage of

    laplace transformations,sufficient conditions for yjr exostamce pf ;a[;ace tramsfpr, pf f(t). general

    properties of laplace transform. Linear properties, applications of simulataneous differential equations,

    methods of solution to system of differential equations. In this the graphs involved the increasing curve

    and alternativecurves. In this mainly discussed about of electrical circuits, with examples and

    explanations. In this each expression is explained with an example circuit problems and with diagramic

    representation, and explained graphically. These are the terms explained in this paper.

    KEYWORDS: Capacitors, switch, registors, transistors, inductors etc.

    Pierre Simon Laplace, after whom the Laplace Transform is named, lived from 1749 to 1827.

    The Laplace Transform is a powerful tool that is very useful in Electrical Engineering. The

    transform allows equations in the "time domain" to be transformed into an equivalent equation in

    the Complex S Domain. The laplace transform is an integral transform, although the reader does

    not need to have a knowledge of integral calculus because all results will be provided. This page

    will discuss the Laplace transform as being simply a tool for solving and manipulating ordinary

    differential equations.

    Laplace transformations of circuit elements are similar to phasor representations, but they are not

    the same. Laplace transformations are more general than phasors, and can be easier to use in

  • some instances. Also, do not confuse the term "Complex S Domain" with the complex power

    ideas that we have been talking about earlier. Complex power uses the variable , while the

    Laplace transform uses the variable s. The Laplace variable s has nothing to do with power.

    The transform is named after the mathematician Pierre Simon Laplace, who lived in the 18th

    century. The transform itself did not become popular until Oliver Heaviside, a famous electrical

    engineer, began using a variation of it to solve electrical circuits.

    The Transform

    The mathematical definition of the Laplace transform is as follows:

    [The Laplace Transform]

    Note:

    The letter s has no special significance, and is used with the Laplace Transform as a matter of common

    convention.

    The transform, by virtue of the definite integral, removes all t from the resulting equation,

    leaving instead the new variable s, a complex number that is normally written as s = + j. In essence, this transform takes the function f(t), and "transforms it" into a function in terms of s,

    F(s). As a general rule the transform of a function f(t) is written as F(s). Time-domain functions

    are written in lower-case, and the resultant s-domain functions are written in upper-case.

    There is a table of Laplace Transform pairs in

    we will use the following notation to show the transform of a function:

    We use this notation, because we can convert F(s) back into f(t) using the inverse Laplace

    transform.

    The Inverse Transform

    The inverse laplace transform converts a function in the complex S-domain to its counterpart in

    the time-domain. Its mathematical definition is as follows:

    [Inverse Laplace Transform]

  • where c is a real constant such that all of the poles s1,s2,...,sn of F(s) fall in the region

    . In other words, c is chosen so that all of the poles of F(s) are to the left of the

    vertical line intersecting the real axis at s = c

    Before the advent of calculators and computers, logarithms were extensively used to

    replace multiplication (or) division of two large numbers by addition (or substraction) of

    two numbers. The cucial idea which made the laplace transformation a very powerful

    technique is that it replaces operations of calculus by operations of algebra for example

    with the applications of laplace transform to an initial value problem, consisting of an

    ordinary (or partial) differential equation. Together with initial conditions is reduced to a

    problem of solving an algebraic equation (with any given initial conditions automatically

    taken care).

    APPLICATIONS OF A LAPLACETRANSFORMATIONS:

    Laplace transform is very seful in obtaining solution of linear differential equations

    both ordinary and partial, solution of system of simultaneous differential equations. And in

    the evaluation of definite integrals.

    ADVANTAGES:

    1 . with the application of laplae transform, particular solution of differential equation is

    obtained directly without the necessity of first determining general solution and then

    obtaining the particular solution (by substitution of initial conditions).

    2 . laplace transformations solves non-homogeneous differential equation without the

    necessity of first solving the corresponding homogeneous differential equations.

    3 . laplace transforms is applicable not only to continuous functions but also to piecewise

    continuous functions, complicated periodic functions, step functions and impulse function.

    4 . laplace transforms of various functions are readily available (in tabulated form). In

    section used functions are tabulated.

    SUFFICIENT CONDITIONS FOR THE EXISTENCE OF LAPLACE TRANSFORM OF

    f(t):-

    The laplacetransform of f(t) exists that is the improper integral in the R.H.S. of (1)

    convages (has a finite value) when the following sufficient conditions are satisfied.

    1 . f(t) is piecewise (or sectionally) continuous that is f(t) is continuous in every subinterval

    and has finite limits at end points of each of these sub interval.

    GENERAL PROPERTIES OF LAPLACE TRANSFORM:

    Although theoretically F(s), the laplace transforms of f(t) is obtained from the

    definition in practice most of the time laplace transforms are obtained by the judical

    application of some of the following important properties. In a nutshell, they are:

  • 1 . Linearity property states that laplace transform of a linear combination(sum) of laplace

    transforms.

    2 . in change of scale, where the argument t of f is multiplied by a constant a,s is replaced

    by s/a in F(s) and then multiplied by 1/a.

    3. First shift theorem proves that multiplication of f(t) by eat

    amounts to replacement of s by

    s-a in F(s).

    4 . Laplace transfomation of a first derivative amounts to multiplication of F(s) by s

    (approximately but for the constant f(o).

    5 . Laplace transform of an integral of f amount to division of F(s) bt s.

    6 . Multiplication of f(t) by t power n amounts to differentiations of F(s) n times with respect to s

    (with(-1)power n as sign).

    APPLICATIONS OF LAPLAE TRANSFORMS TO SYSTEM OF SIMULTANEOUS

    DIFFERENTIAL EQUATIONS:

    Lplace transform can also be used to solve a system (or) family of m simultaneous

    ordinary differential equations in m dependent variables which ar functions of independent

    variable t.

    EXAMPLE 1

    In the circuit shown below, the capacitor is uncharged at time t = 0. If the switch is then closed, find the currents i1 and i2, and the charge on

    C at time t greater than zero.

    Answer

    It is easier in this example to do the second method. In many examples, it is easier to do the first method.

    For the first loop, we have:

  • For the second loop, we have:

    Substituting (2) into (1) gives:

    Next we take the Laplace Transform of both sides.

    Note:

  • In this example, . So

    Now taking Inverse Laplace:

    And using result (2) from above, we have:

    For charge on the capacitor, we first need voltage across the capacitor:

  • So, since , we have:

    Graph of q(t):

    EXAMPLE 2

    In the circuit shown, the capacitor has an initial charge of 1 mC and the switch is in position 1 long enough to establish the steady state.

    The switch is moved from position 1 to 2 at t = 0. Obtain the transient

    current i(t) for t > 0.

  • Answer

    Quiescent implies i1, i2 and their derivatives are zero for t = 0, ie

    i1(0) = i2(0) = i1'(0) = i2'(0) = 0.

    For loop 1:

    For loop 2:

  • Substituting our result from (1) gives:

    Taking Laplace transform:

    Let

    So

    So

    Taking Inverse Laplace:

  • So

    Alternative answer using Scientific Notebook. (.tex file)

    EXAMPLE 6

    Consider a series RLC circuit where R = 20 W, L = 0.05 H and C = 10-

    4 F and is driven by an alternating emf given by E = 100 cos 200t.

    Given that both the circuit current i and the capacitor charge q are

    zero at time t = 0, find an expression for i(t) in the region t > 0.

    Answer

    EXAMPLE 6

    Consider a series RLC circuit where R = 20 W, L = 0.05 H and C = 10-

    4 F and is driven by an alternating emf given by E = 100 cos 200t.

    Given that both the circuit current i and the capacitor charge q are

    zero at time t = 0, find an expression for i(t) in the region t > 0.

    Answer

    We use the following:

  • and obtain:

    After multiplying throughout by 20, we have:

    Taking Laplace transform and using the fact that i(0) = 0:

    Using Scientific Notebook to find the partial fractions:

    So

  • So

    + cos200t 2 sin 200t

    NOTE: Scientific Notebook can do all this for us very easily. In one step, we have:

    + cos200t 2 sin 200t

  • Transient part:

    Steady state part:

  • RESULTS AND DISCUSSION:

    The solutions of problems are finally represented graphically. From example 1 problem the graph shows the graph of the q(t) . q(t) is a function of t. the graph is taken between q and t. there are directly independt on each other. There fpre ot the va;ie pf q os omcreased tjem the t os a;sp omcreased. This is about example 1. From exapmple2 we observed that the gra[h os taklenm between i and t if I increases t also increases.

    CONCLUSION:

    Finally we got the solution for the ordinary differential equations . and in the example 2 the value of current is also got by graph the value of current is also finding in the example by using laplace transforms solutions are easily got this is the conclusion of laplace transformations.

    REFERENCES:

    E.J. Berg, Hevisides operational calculus, end ed.

    R.V. Chur chill,modern operational mathematics in engineering.

    D.V. Widder, advances calculus, end ed.

    Engineering mathematics B.V. Ramana.

  • 7. LAPLACE TRANSFORMS

    VAAGDEVI INSTITUTE OF TECHNOLOGY

    AND SCIENCE

    PRODDATUR

    KADAPA(DIST)

    PRESENTED BY:

    P.KAVYA, N.VARA LAKSHMI,

    I B-Tech, I B-Tech,

    E-mail ID:[email protected] , E-Mail ID:[email protected]

  • Abstract:

    In mathematics, the Laplace transform is a widely used integral transform. Denoted

    , it is a linear operator of a function f(t) with a real argument t (t 0) that transforms it to a

    function F(s) with a complex argument s. This transformation is essentially bijective for the

    majority of practical uses; the respective pairs of f(t) and F(s) are matched in tables. The Laplace

    transform has the useful property that many relationships and operations over the originals f(t)

    correspond to simpler relationships and operations over the images F(s).[1]

    The Laplace

    transform has many important applications throughout the sciences. It is named for Pierre-Simon

    Laplace who introduced the transform in his work on probability theory.

    The Laplace transform is related to the Fourier transform, but whereas the Fourier transform

    resolves a function or signal into its modes of vibration, the Laplace transform resolves a

    function into its moments. Like the Fourier transform, the Laplace transform is used for solving

    differential and integral equations. In physics and engineering, it is used for analysis oflinear

    time-invariant systems such as electrical circuits, harmonic oscillators, optical devices, and

    mechanical systems. In this analysis, the Laplace transform is often interpreted as a

    transformation from the time-domain, in which inputs and outputs are functions of time, to

    the frequency-domain, where the same inputs and outputs are functions of complex angular

    frequency, in radians

  • Introduction:

    The knowledge of Laplace transforms is essential part of mathematics required by the

    engineers and scientists. The Laplace transform is an excellent tool for solving linear differential

    equations with given initial values of an unknown function and its Derivatives without the

    necessity of first finding the general solution (complementary function + particular integral) and

    then evaluating from its the particular solution satisfying the given conditions. This technique is

    useful to solve from fractional differential equations as well. This is a powerful tool in diverse

    fields of engineering.

    The Laplace Transformation

    Pierre-Simon Laplace (1749-1827)

    Laplace was a French mathematician, astronomer, and physicist who applied the

    Newtonian theory of gravitation to the solar system (an important problem of his day).

    He played a leading role in the development of the metric system.

    The Laplace Transform is widely used in engineering applications (mechanical and

    electronic), especially where the driving force is discontinuous. It is also used in process

    control.

  • Definition of Laplace Transform of f(t)

    The Laplace transform of a function f(t) for t> 0 is defined by the following integral defined over

    0 to :

    { f(t)} =

    The resulting expression is a function of s, which we write as F(s). In words we say

    "The Laplace Transform of f(t) equals function F of s"

    and write:

    {f(t)} = F(s)

    Similarly, the Laplace transform of a function g(t) would be written:

    {g(t)} = G(s)

    .

    PROPERTIES:

    PROPERTY 1:

    If a is a constant and f(t) is a function of t, then

    {a f(t)} = a {f(t)}

    {7 sin t} = 7 {sin t}

    PROPERTY 2: Linearity property

    If a and b are constants while f(t) and g(t) are functions of t, then

    {a f(t) + b g(t)} = a {f(t)} + b {g(t)}

    PROPERTY 3: Change of scale property

  • If {f(t)} = F(s) then

    PROPERTY:4:Shifting property (Shift theorem)

    {eat

    f(t)} = F(s a)

    {e3tf(t)} = F(s 3)

    PROPERTY 5:

    Property 6:

    The Laplace transforms of the real (or imaginary) part of a complex function is equal to the real

    (or imaginary) part of the transform of the complex function.

    Let Re denote the real part of a complex function C(t) and Im denote the imaginary part of C(t),

    then

    {Re[C(t)]} = Re {C(t)} and {Im[C(t)]} = Im {C(t)}

    Properties of the unilateral Laplace transform

    Time domain 's' domain Comment

    Linearity

    Can be proved using

    basic rules of

    integration.

    Frequency

    differentia

    tion

    is the first derivative

    of .

    Frequency

    differentia

    tion

    More general form,

    nthderivative of F(s).

  • Differenti

    ation

    is assumed to be a

    differentiable function,

    and its derivative is

    assumed to be of

    exponential type. This

    can then be obtained

    by integration by parts

    Second

    Differenti

    ation

    is assumed twice

    differentiable and the

    second derivative to be

    of exponential type.

    Follows by applying the

    Differentiation property

    to .

    General

    Differenti

    ation

    is assumed to be n-

    times differentiable,

    with nth derivative of

    exponential type.

    Follow by mathematical

    induction.

    Frequency

    integratio

    n

    Integratio

    n

    u(t) is the Heaviside

    step function. Note (u *

    f)(t) is the convolution

    of u(t) and f(t).

    Scaling

    wherea is positive.

    Frequency

    shifting

    Time

    shifting u(t) is the Heaviside

    step function

    Multiplica

    tion

    the integration is done

    along the vertical line

    Re() = c that lies

    entirely within the

    region of convergence

  • of F.[12]

    Convoluti

    on

    (t) and g(t) are

    extended by zero for

    t < 0 in the definition of

    the convolution.

    z

    f(t) is a periodic

    function of periodT so

    that

    . This is the result of the

    time shifting property

    and the geometric

    series

    APPLICATIONS :

    There are two (related) approaches:

    1. Derive the circuit (differential) equations in the time domain, then transform these ODEs

    to the s-domain;

    2. Transform the circuit to the s-domain, then derive the circuit equations in the s-domain

    (using the concept of "impedance").

    We will use the first approach. We will derive the system equations(s) in the t-plane, then

    transform the equations to the s-plane. We will usually then transform back to the t-plane.

    EXAMPLE 1:

    Consider the circuit when the switch is closed at t = 0 with VC(0) = 1.0 V. Solve for the current

    i(t) in the circuit.

  • EXAMPLE 2:Solve for i(t) for the circuit, given that V(t) = 10 sin5t V, R =

    4W and L=2H

    Conclusion

    In this article, the basic development of the NumericalLaplace Transform has been

    presented. This techniquehasproven to be efficient for the analysis of electromagnetic transients

    in power systems. The main advantages of the NLT aresummarized below:

    The modeling of components with distributed and frequencydependentparameters can be

    done in a straightforwardmanner.

    Since its basic principles are different from those oftime domain methods, the NLT is

    very useful to verifyingtime domain methods, as well as in the developmentof new

    time domain models and techniques.

    The application of the NLT can be very important whena high accuracy of results is

    mandatory. The examplesgiven show that time domain methods may require amuch smaller

    discretization step.

  • 8.

    MATHEMATICAL MODELLING TO ESTIMATE THE TIME OF

    DEATH OF A MURDERED PERSON

    Harshitha T and Madhavi P

    I B.Tech, Electronics and Communication in Engineering

    Madanapalle Institute of Technology and Science, Madanapalle, India.

    Email:[email protected]

    [email protected]

    ABSTRACT

    This paper is concerned with an exposition of the methods of solving some classes of ordinary

    differential equations. An analysis is presented for solving ordinary differential equations of first order

    and first degree. A mathematical model is presented to investigate the time of death of a murdered

    person using the Newtons law of cooling. The time of death of a murdered person can be determined

    with the help of modeling through differential equation. To formulate this process mathematically,

    solved analytically and results are shown in graphical representation. It is noticed that the object

    cools, the temperature difference gets smaller, and the cooling rate decreases; thus, the object

  • cools more and more slowly as time passes. Differential equations arise whenever we want to

    represent mathematically a problem involving rate measure. Many real world phenomena can be

    described through either differential equation involving ordinary derivatives/partial derivatives.

    Differential equations play an important role in many applications in the fields of science and

    engineering, such as (i) problems relating to motion of particles (ii) problems involving bending of beams

    (iii) problems related to stability of electric system, chemical process, Economics, Anthropology and

    diverse branches.

    1. INTRODUCTION

    A differential equation is a mathematical equation for an unknown function of one or several

    variables that relates the values of the function itself and its derivatives of various orders. An example of

    modelling a real world problem using differential equations is determination of the velocity of a ball

    falling through the air, considering only gravity and air resistance. The ball's acceleration towards the

    ground is the acceleration due to gravity minus the deceleration due to air resistance. Gravity is

    constant but air resistance may be modelled as proportional to the ball's velocity. This means the ball's

    acceleration, which is the derivative of its velocity, depends on the velocity. Finding the velocity as a

    function of time requires solving a differential equation.

    The theory of differential equations is quite developed and the methods used to study them

    vary significantly with the type of the equation.

    An ordinary differential equation (ODE) is a differential equation in which the unknown function

    (also known as the dependent variable) is a function of a single independent variable. In the simplest

    form, the unknown function is a real or complex valued function, but more generally, it may be vector-

    valued or

    matrix-valued: this corresponds to considering a system of ordinary differential

    equations for a single function. Ordinary differential equations are further

    classified according to the order of the highest derivative with respect to the

  • dependent variable appearing in the equation. The most important cases for

    applications are first order and second order differential equations. In the classical

    literature also distinction is made between differential equations explicitly solved

    with respect to the highest derivative and differential equations in an implicit

    form.

    A partial differential equation (PDE) is a differential equation in which the

    unknown function is a function of multiple independent variables and the equation

    involves its partial derivatives. The order is defined similarly to the case of

    ordinary differential equations, but further classification into elliptic, hyperbolic,

    and parabolic equations, especially for second order linear equations, is of utmost

    importance. Some partial differential equations do not fall into any of these

    categories over the whole domain of the independent variables and they are said

    to be of mixed type.

    Both ordinary and partial differential equations are broadly classified as linear and nonlinear. A

    differential equation is linear if the unknown function and its derivatives appear to the power 1

    (products are not allowed) and non linear otherwise. The

    characteristic property of linear equations is that their solutions form an affine subspace

    of an appropriate function space, which results in much more developed theory of linear

    differential equations. Homogeneous linear differential equations are a further subclass

    for which the space of solutions is a linear subspace i.e. the sum of any set of solutions or

    multiples of solutions is also a solution. The coefficients of the unknown function and its

    derivatives in a linear differential equation are allowed to be (known) functions of the

    independent variable or variables; if these coefficients are constants then one speaks of a constant

    coefficient linear differential equation.

  • Differential equations are mathematically studied from several different perspectives, mostly

    concerned with their solutions, the set of functions that satisfy the equation. Only the simplest

    differential equations admit solutions given by explicit formulas; however, some properties of solutions

    of a given differential equation may be determined without finding their exact form. If a self-contained

    formula for the solution is not available, the solution may be numerically approximated using

    computers. The theory of dynamical systems puts emphasis on qualitative analysis of systems described

    by differential equations, while many numerical methods have been developed to determine solutions

    with a given degree of accuracy.

    Many of the general laws in all branches of engineering, physics, chemistry, biology, astronomy,

    population studies, geometry, economics, etc. can be expressed through Mathematical Modelling in

    the form of an equation connecting the variables and their rate of change, resulting in differential

    equations. The time of death of a murdered person can be determined with the help of modeling

    through differential equation.

    The objective of the present study is to investigate the time of death of a murdered person

    using the Newtons law of cooling.

    2. MATHEMATICAL MODEL

    When a hot object is placed in a cool room, the object dissipates heat to the

    surroundings, and its temperature decreases. Newton's Law of Cooling states that the rate

    at which the object's temperature decreases is proport ional to the difference between the

    temperature of the object and the ambient temperature. At the beginning of the cooling

    process, the difference between these temperatures is greatest, so this is when the rate of

    temperature decrease is greatest. However, as the object cools, the temperature difference

    gets smaller, and the cooling rate decreases; thus, the object cools more and more slowly as

    time passes. To formulate this process mathematically, let T(t) denote the temperature of the

  • object at time t and let Ts denote the (essentially constant) temperature of the surroundings.

    Newton's Law of Cooling then says,

    dT

    (T T )sdt (1)

    i.e. K dT

    (T T ) (K > 0)sdt (2)

    Where K is the unknown proportionality constant, the equation (2) governing the Newtons law of

    cooling, is a first order and first degree linear separable differential equation.

    Since Ts

  • 3. METHOD OF SOLUTION

    We therefore present a more detailed exposition here to solve the above problem

    analytically. Essentially 3 phases are central to the solution. These are:

    I. Identify Ts, the temperature of the surrounding medium, so that the general

    solution is given by (3).

    II. Use two conditions given to determine the constant of integration C and the

    unknown proportionality constant K.

    III. Substituting C and K, obtained from step II, in equation (3) (a) the value of T

    for a given time t or (b) the value of time t for a given temperature T can

    be determined from equation (3)

    Let T(t) denote the temperature of the body at time t.

    Then the ( ) KtsT t T Ce since Ts =

    022 C

    Constants K and C can be determined provided the following information is available:

    Time of arrival of the person, the temperature of the body just after his arrival, temperature of the

    body after certain interval of time.

    Let the officer arrived at 8.30 a.m. and the body temperature was 30 degrees. This means

    that if the officer considers 8:30 a.m. as t = 0, i.e. T(0)=30 then

    k*0

    30 22 C e

    8C

    Hence k*tT t 22 8e (4)

  • Let the officer makes another measurement of the temperature say after 1 hour (60

    minutes), that is, at 9.30 a.m. and temperature was 28 degrees. This means that T(1) = 028

    ,

    k*128 22 8e

    k*1 6

    8e

    ln 4 ln 3k

    k = 0.2877

    Now it only remains to find out when the murder occurred. At the time of death the body

    was 037 C (i.e. normal human body temperature).

    So, equation (4) k*t

    37 22 8e

    k*t 15

    8e

    k*t ln 1.875

    1

    t *ln 1.875k

    2.19t Hours

    Here sign minus represents the time before found the victim.

  • From this, we can conclude that the murder occurred about 2 hours and 19 minutes before

    the body was found.

    The death occurred approximately 139 minutes before the first measurement at 8.30a.m that is

    at 6:19a.m approximately

    4. RESULTS AND DISCUSSIONS

    The present analysis integrates the equations by the analytical method. The details of the

    solution method is shown and the variation of the temperature T of the body before identifying

    and after for different time t are also shown graphically.

    From experimental observations it is known that (up to a satisfactory approximation) the

    surface temperature of an object changes at a rate proportional to its relative temperature. That

    is, the difference between its temperature and the temperature of the surrounding environment.

    This is what is known as Newton's law of cooling. The time of death of a murdered person is

    determined with the help of modeling through differential equation.

    According to Newtons law of cooling, the body will radiate heat energy into the room at

    a rate proportional to the difference in temperature between the body and the room.

    For some time after the death, the body will radiate heat into the cooler room, causing the

    bodys temperature to decrease assuming that the victims temperature was normal 37oC (98.6oF) at the

    time of death. Forensic expert will try to estimate this time from bodys current temperature and

    calculating how long it would have had to lose heat to reach this point.

    The following figure display the temperature of the body of victim for different time. It is

    found that temperature decrease monotonically from normal human body temperature (37oC) for

    increasing time because the body will radiate heat into the cooler room, causing the bodys

  • temperature to decrease. This is owing to the fact that the body temperature suddenly falls in the

    first 10 minutes after murdered and then slowly decrease with increasing time.

    The death occurred approximately 139 minutes before the first measurement at 8.30a.m

    that is at 6:19a.m approximately

    TABLE AND GRAPH:

    t (time in minutes) T(t) (temperature in

    degrees) t (time in minutes)

    T(t) (temperature in

    degrees)

    -130 37 0 30

    -120 32.8 10 29.8353

    -110 32.1 20 29.674

    -100 31.8497 30 29.516

    -90 31.6469 40 29.361

    -80 31.4483 50 29.209

    -70 31.2539 60 29.0613

    -60 31.0634 70 28.916

    -50 30.9124 80 28.7736

    -40 30.6941 90 28.634

    -30 30.5151 100 28.497

    -20 30.3398 110 28.364

    -10 30.1681 120 28.2329

    130 28.104

  • 5. CONCLUSIONS

    The present paper analytically studied to investigate the time of death using the Newtons law of cooling. The time of death of a murdered person is determined with the help of modeling through differential equation. To formulate this process mathematically, solved analytically and results are shown in graphical representation. It is noticed that the murdered body cools, the temperature difference gets smaller, and the cooling rate decreases; thus, the object cools more and more slowly as time passes. The death occurred approximately 139 minutes before the first measurement at 8.30a.m that is at 6:19a.m approximately. The problem of solving differential equations is a natural goal of differential and integral calculus. Further many of the general laws of nature in Physics, Chemistry, Biology and Astronomy can be expressed in the language of differential equations and hence the theory of differential equations is the most important part of mathematics for understanding Physical sciences. 6. REFERENCES 1. Engineering mathematics volume-1, Dr.T.K.V.Iyengar, Dr.B.krishnagandhi, S.Ranganatham,M.V.S.S.N.Prasad 2. AText book of EngineringMathematics, B.V.Ramana. 3. A Text Book of Engineering Mathematics, Thomson Book Collection.

  • 9.

    CONCEPT ON

    MATRICES MATHEMATICS

    MOULA ALI COLLEGE OF ENGINEERING AND TECHNOLOGY

    PREPARED BY

    T.GOWSIA B.PRIYANKA

    REG.NO:10F51A0424 10F51A0460

    E.C.E E.C.E

    I-BTECH I-BTECH

    Email:[email protected]

    ABSTRACT

    DEFINATION

    NOTATION

    BASIC OPERATORS

    MATRICES MULTIPLICATION,LINEAR EQUATIONS

    RANK OF MATRICES

    SQUARE MATRICES

    INVERSE MATRICES

    TRACE

    TRIANGULAR MATRICES

    DIAGONAL MATRICES

    DETERMINENT

  • EIGENVECTORS AND EIGENVALUES

    SYMMETRY

    MATRICES DECOMPOSITION METHODS(LU-decomposition method)

    MATRICES:

    Definition

    A matrix is a rectangular arrangement of numbers.

    For example A= 43

    21

    The horizontal and vertical lines in a matrix are called rows and columns, respectively.

    The numbers in the matrix are called its entries or its elements.

    To specify the size of a matrix, a matrix with m rows and n columns is called an m-by-n matrix or m n matrix, while m and n are called its dimensions.

    ROW MATRICES:

    A matrix with one row is called a row matrix.

    For example: A= 21

    COLUMN MATRICES:

    A matrix with one column is called acolumn matrix.

    For example: A=3

    2

    TRANSPOSE MATRICES:

    The rows of a matrix are equal to the corresponding columns of itsmatrix, the matrix is transpose

    matrix.

    Most of this article focuses on real and complex matrices, i.e., matrices whose entries are real or complex numbers.

  • Notation

    The entry that lies in the i-th row and the j-th column of a matrix is typically referred to as the

    (i,j), or (i,j)th entry of the matrix.

    For example, the (2, 2) entry of the above matrix A is 4.

    The (i, j) th

    entry of a matrix A is most commonly written as ai,j.

    Alternative notations for that entry are A [i,j] or Ai,j.

    Letai, refers to the ith row of A, and a,j refers to the j

    th column of A. The set of all m-by-n

    matrices is denoted (m, n).

    A = [a, j]i=1,...,m; j=1,...,n or more briefly A = [ai,j]mn

    To define an m n matrix A. Usually the entries ai,j are defined separately for all integers 1 i m and 1 j n.

    Basic operations

    Matrix addition

    Scalar multiplication

    Transpose

    Row operations

    There are a number of operations that can be applied to modify matrices called matrix addition,

    scalar multiplication and transposition. These form the basic techniques to deal with matrices.

    Operation Definition

    Addition

    The sumA+B of two m-by-n matrices A and B is calculated entrywise:

    (A + B)i,j = Ai,j + Bi,j, where 1 i m and 1 j n.

    Scalar

    multiplication

    The scalar multiplication cA of a matrix A and a number c is given by

    multiplying every entry of A by c:

    (cA)i,j = c Ai,j.

    Transpose

    The transpose of an m-by-n matrix A is the n-by-m matrix AT (also denoted

    Atr or

    tA) formed by turning rows into columns and vice versa:

    (AT)i,j = Aj,i.

    Addition of matrices iscommutative,

  • A + B = B + A.

    The transpose is compatible with addition and scalar multiplication

    (cA)T = c (A

    T) and (A + B)

    T = A

    T + B

    T. Finally, (A

    T)

    T = A.

    Row operations

    There are three types of row operations: row switching, that is interchanging two rows of a

    matrix; row multiplication, multiplying all entries of a row by a non-zero constant; and finally

    row addition, which means adding a multiple of a row to another row. These row operations are

    used in a number of ways including solving linear equations and finding inverses

    Matrix multiplication, linear equations

    Matrix multiplication

    Schematic depiction of the matrix product AB of two matrices A and B.

    Multiplication of two matrices is defined only if the number of columns of the left matrix is the

    same as the number of rows of the right matrix.

    If A is an m-by-n matrix and B is an n-by-p matrix, then their matrix productAB is the m-by-p

    matrix whose entries are given by dot-product of the corresponding row of A and the

    corresponding column of B:

    Matrix multiplication satisfies the rules (AB) C = A (BC) (associativity),

    And (A+B) C = AC+BC as well as C(A+B) = CA+CB (left and right distributivity),

    Whenever the size of the matrices is such that the various products are defined.

    IfA and B are m-by-n and n-by-k matrices, respectively, and m k. Even if both products are defined, they need not be equal, i.e. generally one has

    AB BA,

    i.e., matrix multiplication is not commutative,

    The identity matrixIn of size n is the n-by-n matrix in which all the elements on the main

    diagonal are equal to 1 and all other elements are equal to 0,

    It is called identity matrix because multiplication with it leaves a matrix unchanged: MIn = ImM

    = M for any m-by-n matrix M.

    Linear equations

    Linear equation and System of linear equations

  • A particular case of matrix multiplication is tightly linked to linear equations: if x designates a

    column vector (i.e. n1-matrix) of n variables x1, x2, ..., xn, and A is an m-by-n matrix, then the

    matrix equation

    Ax = b,

    Whereb is some m1-column vector, is equivalent to the system of linear equations

    A1,1x1 + A1, 2x2 + ... + A1,nxn = b1

    ...

    Am, 1x1 + Am, 2x2 + ... + Am, nxn = bm.

    This way, matrices can be used to compactly write and deal with multiple linear equations, i.e.

    systems of linear equations.

    RANK OF MATICES:

    The rank of a matrixA is the maximum number of linearly independent row vectors of the

    matrix, which is the same as the maximum number of linearly independent column vectors.

    Square matrices

    A square matrix is a matrix which has the same number of rows and columns.

    For example: A= 32

    21

    An n-by-n matrix is known as a square matrix of order n. Any two square matrices of the same

    order can be added and multiplied.

    A square matrix A is called invertible or non-singular if there exists a matrix B such that

    AB = In.

    INVERSE MATRICES:

    This is equivalent to BA = In. Moreover, if B exists, it is unique and is called the inverse matrix

    of A, denoted A1

    .

    TRACE:

    The entries Ai,i form the main diagonal of a matrix. The trace, tr(A) of a square matrix A is the

    sum of its diagonal entries,

    the trace of the product of two matrices is independent of the order of the factors:

    tr(AB) = tr(BA).

    Also, the trace of a matrix is equal to that of its transpose, i.e. tr(A) = tr(AT).

  • DIAGONAL MATRICES:

    If all entries outside the main diagonal are zero, A is called a diagonal matrix.

    TRIANGULAR MATRICES:

    If only all entries above (below) the main diagonal are zero, A is called a lower triangular matrix

    (upper triangular matrix, respectively). For example, if n = 3, they look like (Diagonal), (lower)

    and (upper triangular matrix).

    Lower triangular matrix: A= 30

    21

    Upper triangular matrix: A= 43

    01

    Determinant

    The determinantdet (A) or |A| of a square matrix A is a number encoding certain properties of the

    matrix. A matrix is invertible if and only if its determinant is nonzero.

    The determinant of a product of square matrices equals the product of their determinants: det

    (AB) = det (A) det (B).

    Adding a multiple of any row to another row, or a multiple of any column to another column,

    does not change the determinant.

    To solve linear systems using Cramer's rule, where the division of the determinants of two

    related square matrices equates to the value of each of the system's variables.

    Eigenvalues and eigenvectors

    A number and a non-zero vector v satisfying

    AX = X

    are called an eigenvalue and an eigenvector of A, respectively. The number is an eigenvalue of an nn-matrix A if and only if AIn is not invertible.

    The function pA(t) = det(AtI) is called the characteristic polynomial of A, its degree is n. Therefore pA(t) has at most n different roots, i.e., eigenvalues of the matrix

    They may be complex even if the entries of A are real. According to the Cayley-Hamilton

    theorem, pA(A) = 0, that is to say, the characteristic polynomial applied to the matrix itself yields

    the zero matrix.

  • SYMMETRY:

    SYMMETRIC MATRICES:

    A square matrix A that is equal to its transpose, i.e. A = AT, is a symmetric matrix.

    SKEW-SYMMETRICE MATRICES:

    If instead, A was equal to the negative of its transpose, i.e. A = AT, then A is a skew-symmetric matrix.

    HERMITIAN MATRICES:

    In complex matrices, symmetry is often replaced by the concept of Hermitian matrices, which

    satisfy A = A, where the star or asterisk denotes the conjugate transpose of the matrix, i.e. the

    transpose of the complex conjugate of A.

    Matrix decomposition methods

    Matrix decomposition

    Matrix diagonalization

    Gaussian elimination

    Montante's method

    There are several methods to render matrices into a more easily accessible form. They are

    generally referred to as matrix transformation or matrix decomposition techniques.

    LU decomposition:

    The LU decomposition factors matrices as a product of lower (L) and an upper triangular

    matrices (U).

    Problem on LU-decomposition method:

    Using LU-decomposition method solve x+y+z=1, 3x+y-3z=5, x-y-5z=10.

    Solution:

    521

    313

    111

    z

    y

    x

    =

    10

    5

    1

    i.e., AX=B

    Let A=LU

  • Where L=

    13231

    012

    001

    ll

    l

    U=

    3300

    23220

    131211

    u

    uu

    uuu

    A=LU

    =

    13231

    0121

    001

    ll

    l

    =

    3323321331223212311131

    2313212212211121

    131211

    uululululul

    uuluulul

    uuu

    Compare the corresponding elements of LU and the elements of A

    u11=1 ,

    u12=1 ,

    u13=1,

    l21u11=3 ,

    l21=3

    l21u12+u22=1

    u22=-2

    l21u13+u23=-3

    u23=-6

    l31u11=1

  • l31=1

    l31u12+l32u22=-2

    l32=3/2

    l31u13+l32u23+u33=-5

    u33=3

    12/31

    013

    001

    300

    620

    111

    AX=B

    LUX=B

    Let UX=Y

    LY=B

    12/31

    013

    001

    3

    2

    1

    Y

    Y

    Y

    =

    10

    5

    1

    322/31

    213

    1

    YYY

    YY

    Y

    =

    10

    5

    1

    Y1=1

    3Y1+Y2=5

    Y2=2

    Y1+3/2Y2+Y3=10

    Y3=6

  • 32

    1

    Y

    Y

    Y

    =

    6

    2

    1

    Since UX=Y

    300

    620

    111

    Z

    Y

    X

    6

    2

    1

    Z

    ZY

    ZYX

    3

    62

    =

    6

    2

    1

    X+Y+Z=1

    -2Y-6Z=2

    3Z=6

    Z=2

    -2Y-6Z=2

    Y=-7

    X=6

    X =

    z

    y

    x

    =

    2

    7

    6

    Therefore the problem is sloved.

  • 10. Moula ali college of engineering and

    technology,

    Anantapur.

    Prepared by

    B.Peddy reddy B.Aravind kumar

    ECE(I B.TECH) ECE(I B.TECH)

    REG.NO 10F51A0458 REG.NO 10F51A0408

    EMAIL.ID [email protected]

    ABSTRACT

  • Applications of Differential Equations

    Application 1 : Exponential Growth - Population

    Application 2 : Exponential Decay - Radioactive Material

    Application 3 : Falling Object

    MATRICES

    Definition of a Matrix

    Introduction to Matrix:

    Other Types of Matrix

    Transpose Matrix:

    Definition

    Adjoint Matrix

    Minors

    Co-factors

    Properties of Determinants

  • Applications of Differential Equations

    We present examples where differential equations are widely applied to model natural phenomena, engineering systems and many other situations.

    Application 1 : Exponential Growth - Population

    Let P(t) be a quantity that increases with time t and the rate of increase is proportional to the same quantity P as follows

    d P / d t = k P where d p / d t is the first derivative of P, k > 0 and t is the time. The solution to the above first order differential equation is given by

    P(t) = A ek t where A is a constant not equal to 0. If P = P0 at t = 0, then P0 = A e

    0 which gives A = P0 The final form of the solution is given by

    P(t) = P0 ek t

    Assuming P0 is positive and since k is positive, P(t) is an increasing exponential. d P / d t = k P is also called an exponential growth model.

    Application 2 : Exponential Decay - Radioactive Material

    Let M(t) be the amount of a product that decreases with time t and the rate of decrease is proportional to the amount M as follows

    d M / d t = - k M where d M / d t is the first derivative of M, k > 0 and t is the time.

  • Solve the above first order differential equation to obtain

    M(t) = A e- k t where A is non zero constant. It we assume that M = M0 at t = 0, then M0 = A e

    0 which gives A = M0 The solution may be written as follows

    M(t) = M0 e- k t

    Assuming M0 is positive and since k is positive, M(t) is an decreasing exponential. d M / d t = - k M is also called an exponential decay model.

    Application 3 : Falling Object

    An object is dropped from a height at time t = 0. If h(t) is the height of the object at time t, a(t) the acceleration and v(t) the velocity. The relationships between a, v and h are as follows: a(t) = dv / dt , v(t) = dh / dt. For a falling object, a(t) is constant and is equal to g = -9.8 m/s. Combining the above differential equations, we can easily deduce the follwoing equation d 2h / dt 2 = g Integrate both sides of the above equation to obtain dh / dt = g t + v0 Integrate one more time to obtain h(t) = (1/2) g t + v0 t + h0 The above equation describes the height of a falling object, from an initial height h0 at an initial velocity v0, as a function of time. Substitute c in the solution - (L / R) ln(E - R i) = t + (-L/R) ln (E)

  • which may be written (L/R) ln (E)- (L / R) ln(E - R i) = t ln[E/(E - Ri)] = t(R/L) Change into exponential form [E/(E - Ri)] = et(R/L) Solve for i to obtain i = (E/R) (1-e-Rt/L) The starting model for the circuit is a differential equation which when solved, gives an expression of the current in the circuit as a function of time.

  • Definition of a Matrix

    A rectangular array of entries is called a Matrix. The entries may be real, complex or functions.

    The entries are also called as the elements of the matrix.

    The rectangular array of entries are enclosed in an ordinary bracket or in square bracket.

    Matrices are denoted by capital letters.

    Example:

    (i)

    A matrix having m rows and n columns is called as matrix of order mxn. Such a matrix has mn

    elements.

    In general, an mxn matrix is in the form

    Where aij represents the element in ith column.

    The above matrix may be denoted as A = [aij]mxn.

    Introduction to Matrix:

    A Matrix is a rectangular arrangement of number in rows and columns and enclosed by

    Parenthesis (or) Brackets. Matrices are denoted by A, B, C Matrix is a way of organizing the

    data in order to rows and columns. it can be w