jump to navigation

Introducing Integral Transforms October 14, 2009

Posted by Stephen Godfrey in Mathematics.
Tags: ,
trackback
Welcome to my first post, I have decided to start off this blog with a post about the first field of maths that I did some independent study in. This was back in 2005 when I was doing Honours, I looked into this area because I had just done a course where we looked at solving PDE’s with Fourier and Laplace Transforms (among other things) and thought the whole idea was really interesting.

As I spent the year studying I slowly found out that the current research in the field was rather different to what was doing in the classics, it appeared that I was born about 60-70 years too late. Most of the current research looks at new ways to invert Laplace transforms, deals with certain classes of distributions (also called Generalised functions), or plug a Hypergeometric function in the kernel. Not easy for work for an undergrad.

Now enough of this chit chat lets get into some maths!

The idea of using transformations in Mathematics is an old one. The idea is to change your problem into a simpler but equivalent problem, then change back to get the solution of the original problem. In this post I shall briefly explain why integrals transforms are useful when solving differential equations. The simplest answer is that one integral undoes one lot of differentiation, so if we integrate a ordinary differential equation we should only have an algebraic one.

The key point to remember is that we have mapped differentiation to something simpler like multiplication.

Integral transforms have been in wide use during the past two centuries as a tool to solve various problems in pure and applied mathematics. Many integral transforms were originally introduced to solve specific problems, but over the course of time have been found to be of use in the solution of other problems as well.

Definition 1
Let {\Omega \subseteq \mathbf{R}^n} and let {K: \Omega \times \Omega \rightarrow \mathbf{R}} Let

\displaystyle A=\left\{ f:\Omega \rightarrow \mathbf{R} {\Big |} \int_\Omega |f(x)K(x,y)| dx <\infty \right\}. \ \ \ \ \ (1)

Then the mapping
\displaystyle \mathcal{T} \{f(x):y\}=\int_\Omega f(x) K(x,y) dx, \ \ \ \ \ (2)

is an integral transform with domain {A}, kernel {K(x,y)} and transform variable {y}.

There are several key questions that should be asked about any integral transform {\mathcal{T}}. For this post and possibly later ones we shall only look at the following 4 key questions.

  1. What functions {f(x)} are in {A}, the domain of the operator {\mathcal{T}}?
  2. For every function {f \in A} is there a unique function {\mathcal{T} (f(x))}?
  3. Given a function {\mathcal{T} \{f(x)\}}, does there exist an operator that recovers {f(x)}? If so does the domain of this operator match exactly with {\mathcal{T} (A)} ?
  4. What differential operator does this transform diagonalize?

For practical purposes the most important question that of if we can undo the transformation. With a bit of Functional Analysis we can answer this question for every integral transform. Since integration is linear it is easy to see that {\mathcal{T}} is linear. That is, the operator {\mathcal{T}} has the following property,\displaystyle \mathcal{T} (\alpha f(x) +\beta g(x) )=\alpha \mathcal{T} (f(x)) +\beta \mathcal{T} (g(x)). \ \ \ \ \ (3)

The following lemma shows that if {\mathcal{T}} is linear and one-to-one, then {\mathcal{T}^{-1}} exists and it is also a linear operator.

Lemma 1
If a linear operator {L} from {X} into {Y} is one-to-one, then there exists an operator {M}, called the inverse of {L}, such that {ML=I_X} and {LM=I_Y}, where {I_X} and {I_Y} are identity operators on {X} and {Y}. The operator {M} is also linear. Often {M} is denoted by {L^{-1}}.

It is important to note that although the linear operator {\mathcal{T}} is an integral transform the lemma does not state how {\mathcal{T}^{-1}} is defined. In many cases {\mathcal{T}^{-1}} is also an integral transform, however this is not always the case. It is important to keep in mind that Lemma 1 is a statement about the existence of {\mathcal{T}^{-1}}. In many cases there are several different ways of inverting an integral transform.

Many integral transforms share similar properties. What I intend to do is a quick survey of the well know integral transforms to look at the similarities as motivation to look for a more general framework of Integral transforms.

In an attempt to give a more general theory of integral transforms we shall attack this problem in three different ways, though self-adjoin DE’s, determining the relationship between the kernel in {\mathcal{T}} and {\mathcal{T}^{-1}}, and by having the kernel of an integral transform being a Hypergeometric function.

Overall, as I mentioned above I will be interested in solving differential equations by integral transform methods. This is more out of personal interest as it is not the only place these transforms arise. The reason that integral transform methods are well suited to this problem is that many integral transforms convert certain differential operators into operators that act by multiplication.

As a simple example we consider the ordinary differential equation (ODE)

\displaystyle \frac{d^2y}{dt^2}-5\frac{dy}{dt}+6y=0, \ \ \ \ \ (4)

with the initial conditions {y(0)=0} and {y'(0)=1}. The Laplace transform of this ODE is
\displaystyle s^2 Y(s) -sy(0)-y'(0)-5(sY(s)-y(0))+6Y(s)=0, \ \ \ \ \ (5)

which in turn can be solved for {Y(s)} the Laplace transform of the unknown function {y(t)}, after some work you find that
\displaystyle Y(s)=\frac{1}{s^2-5s+6}. \ \ \ \ \ (6)

To find the solution to the differential equation we apply the inverse Laplace transform to {Y(s)}.
For Partial Differential Equations (PDE’s) we often do not reduce the equation to an algebraic problem. Often we reduce it to the problem of solving an ODE. Consider the inhomogeneous equation of telegraphy in the {x} variable,

\displaystyle \frac{\partial^2 u}{\partial x^2} -\frac{1}{c^2}\frac{\partial^2 u}{\partial t^2}-\frac{4\pi \sigma}{c^2}\frac{\partial u}{\partial t}=\frac{4\pi}{c} f(x). \ \ \ \ \ (7)

By applying the Fourier transform in the {x} variable it is transformed into
\displaystyle \frac{d^2 U}{dt^2}+4\pi \sigma \frac{dU}{dt}-4\pi^2\xi^2c^2U=-4\pi c F(\xi), \ \ \ \ \ (8)

which is a much simpler problem to solve. In fact the solution to equation (8) can be shown to equal,
\displaystyle U(\xi,t)=-4\pi c F(\xi) \int_0^t e^{-2\pi \sigma (t-\tau )} \frac{\sin (\sqrt{\xi^2c^2-4\pi^2\sigma^2}(t-\tau ))}{\sqrt{\xi^2c^2-4\pi^2\sigma^2}} d \tau . \ \ \ \ \ (9)

To find the solution to equation (7) we just apply the inverse Fourier transform to {U(\xi,t)} in the first variable.
Although I will mainly be interested in solving problems like the previous two examples, I will also look at solving certain types of integral equations.

Comments»

No comments yet — be the first.

Leave a comment