Line 1:
Line 1:
{{calculus}}
{{seealso|Integration by substitution|Chain rule}}
Often a [[partial differential equation]] can be reduced to a simpler form with a known solution by a suitable '''change of variables'''.
Often a [[partial differential equation]] can be reduced to a simpler form with a known solution by a suitable '''change of variables'''.
Often a partial differential equation can be reduced to a simpler form with a known solution by a suitable change of variables .
Change of variable for integral equations is discussed in Integration by substitution .
The article below discusses change of variable for PDEs in two ways:
By example
By giving the theory of the method
Technique explained by example
For example the following simplified form of the Black–Scholes PDE
∂
V
∂
t
+
1
2
S
2
∂
2
V
∂
S
2
+
S
∂
V
∂
S
−
V
=
0.
{\displaystyle {\frac {\partial V}{\partial t}}+{\frac {1}{2}}S^{2}{\frac {\partial ^{2}V}{\partial S^{2}}}+S{\frac {\partial V}{\partial S}}-V=0.}
is reducible to the Heat equation
∂
u
∂
τ
=
∂
2
u
∂
x
2
{\displaystyle {\frac {\partial u}{\partial \tau }}={\frac {\partial ^{2}u}{\partial x^{2}}}}
by the change of variables[ 1] :
V
(
S
,
t
)
=
v
(
x
(
S
)
,
τ
(
t
)
)
{\displaystyle V(S,t)=v(x(S),\tau (t))}
x
(
S
)
=
ln
(
S
)
{\displaystyle x(S)=\ln(S)}
τ
(
t
)
=
1
2
(
T
−
t
)
{\displaystyle \tau (t)={\frac {1}{2}}(T-t)}
v
(
x
,
τ
)
=
exp
(
−
(
1
/
2
)
x
−
(
9
/
4
)
τ
)
u
(
x
,
τ
)
{\displaystyle v(x,\tau )=\exp(-(1/2)x-(9/4)\tau )u(x,\tau )}
in these steps:
Replace
V
(
S
,
t
)
{\displaystyle V(S,t)}
by
v
(
x
(
S
)
,
τ
(
t
)
)
{\displaystyle v(x(S),\tau (t))}
and apply the chain rule to get
1
2
(
−
2
v
(
s
,
τ
)
+
2
∂
τ
∂
t
∂
v
∂
τ
+
S
(
(
2
∂
x
∂
S
+
S
∂
2
x
∂
S
2
)
∂
v
∂
x
+
S
(
∂
x
∂
S
)
2
∂
2
v
∂
x
2
)
=
0
{\displaystyle {\frac {1}{2}}(-2v(s,\tau )+2{\frac {\partial \tau }{\partial t}}{\frac {\partial v}{\partial \tau }}+S\left(\left(2{\frac {\partial x}{\partial S}}+S{\frac {\partial ^{2}x}{\partial S^{2}}}\right){\frac {\partial v}{\partial x}}+S\left({\frac {\partial x}{\partial S}}\right)^{2}{\frac {\partial ^{2}v}{\partial x^{2}}}\right)=0}
Replace
x
(
S
)
{\displaystyle x(S)}
and
τ
(
t
)
{\displaystyle \tau (t)}
by
ln
(
S
)
{\displaystyle \ln(S)}
and
1
2
(
T
−
t
)
{\displaystyle {\frac {1}{2}}(T-t)}
to get
1
2
(
−
2
v
(
ln
(
S
)
,
1
2
(
T
−
t
)
)
−
∂
v
(
ln
(
S
)
,
1
2
(
T
−
t
)
)
∂
τ
+
∂
v
(
ln
(
S
)
,
1
2
(
T
−
t
)
)
∂
x
+
∂
2
v
(
ln
(
S
)
,
1
2
(
T
−
t
)
)
∂
x
{\displaystyle {\frac {1}{2}}(-2v(\ln(S),{\frac {1}{2}}(T-t))-{\frac {\partial v(\ln(S),{\frac {1}{2}}(T-t))}{\partial \tau }}+{\frac {\partial v(\ln(S),{\frac {1}{2}}(T-t))}{\partial x}}+{\frac {\partial ^{2}v(\ln(S),{\frac {1}{2}}(T-t))}{\partial x}}}
Replace
ln
(
S
)
{\displaystyle \ln(S)}
and
1
2
(
T
−
t
)
{\displaystyle {\frac {1}{2}}(T-t)}
by
x
(
S
)
{\displaystyle x(S)}
and
τ
(
t
)
{\displaystyle \tau (t)}
and divide both sides by
1
2
{\displaystyle {\frac {1}{2}}}
to get
−
2
v
−
∂
v
∂
τ
+
∂
v
∂
x
+
∂
2
v
∂
x
2
=
0
{\displaystyle -2v-{\frac {\partial v}{\partial \tau }}+{\frac {\partial v}{\partial x}}+{\frac {\partial ^{2}v}{\partial x^{2}}}=0}
Replace
v
(
x
,
τ
)
{\displaystyle v(x,\tau )}
by
exp
(
−
(
1
/
2
)
x
−
(
9
/
4
)
τ
)
u
(
x
,
τ
)
{\displaystyle \exp(-(1/2)x-(9/4)\tau )u(x,\tau )}
and divide through by
−
exp
(
−
(
1
/
2
)
x
−
(
9
/
4
)
τ
)
u
(
x
,
τ
)
{\displaystyle -\exp(-(1/2)x-(9/4)\tau )u(x,\tau )}
to yield the heat equation.
Advice on the application of change of variable to PDEs is given by mathematician J. Michael Steele [ 2] :
"There is nothing particularly difficult about changing variables and transforming one equation to another, but there is an element of tedium and complexity that slows us down. There is no universal remedy for this molasses effect, but the calculations do seem to go more quickly if one follows a well-defined plan. If we know that
V
(
S
,
t
)
{\displaystyle V(S,t)}
satisfies an equation (like the Black-Scholes equation) we are guaranteed that we can make good use of the equation in the derivation of the equation for a new function
v
(
x
,
t
)
{\displaystyle v(x,t)}
defined in terms of the old if we write the old V as a function of the new v and write the new t and x as functions of the old t and S . This order of things puts everything in the direct line of fire of the chain rule; the partial derivatives
∂
V
∂
t
{\displaystyle {\frac {\partial V}{\partial t}}}
,
∂
V
∂
S
{\displaystyle {\frac {\partial V}{\partial S}}}
and
∂
2
V
∂
S
2
{\displaystyle {\frac {\partial ^{2}V}{\partial S^{2}}}}
are easy to compute and at the end, the original equation stands ready for immediate use."
Technique in general
Suppose that we have a function
u
(
x
,
t
)
{\displaystyle u(x,t)}
and a change of variables
x
1
,
x
2
{\displaystyle x_{1},x_{2}}
such that there exist functions
a
(
x
,
t
)
,
b
(
x
,
t
)
{\displaystyle a(x,t),b(x,t)}
such that
x
1
=
a
(
x
,
t
)
{\displaystyle x_{1}=a(x,t)}
x
2
=
b
(
x
,
t
)
{\displaystyle x_{2}=b(x,t)}
and functions
e
(
x
1
,
x
2
)
,
f
(
x
1
,
x
2
)
{\displaystyle e(x_{1},x_{2}),f(x_{1},x_{2})}
such that
x
=
e
(
x
1
,
x
2
)
{\displaystyle x=e(x_{1},x_{2})}
t
=
f
(
x
1
,
x
2
)
{\displaystyle t=f(x_{1},x_{2})}
and furthermore such that
x
1
=
a
(
e
(
x
1
,
x
2
)
,
f
(
x
1
,
x
2
)
)
{\displaystyle x_{1}=a(e(x_{1},x_{2}),f(x_{1},x_{2}))}
x
2
=
b
(
e
(
x
1
,
x
2
)
,
f
(
x
1
,
x
2
)
)
{\displaystyle x_{2}=b(e(x_{1},x_{2}),f(x_{1},x_{2}))}
and
x
=
e
(
a
(
x
,
t
)
,
b
(
x
,
t
)
)
{\displaystyle x=e(a(x,t),b(x,t))}
t
=
f
(
a
(
x
,
t
)
,
b
(
x
,
t
)
)
{\displaystyle t=f(a(x,t),b(x,t))}
Suppose
L
{\displaystyle {\mathcal {L}}}
is a differential operator such that
L
u
(
x
,
t
)
=
0
{\displaystyle {\mathcal {L}}u(x,t)=0}
Then it is also the case that
L
v
(
x
1
,
x
2
)
=
0
{\displaystyle {\mathcal {L}}v(x_{1},x_{2})=0}
where
v
(
x
1
,
x
2
)
=
u
(
e
(
x
1
,
x
2
)
,
f
(
x
1
,
x
2
)
)
{\displaystyle v(x_{1},x_{2})=u(e(x_{1},x_{2}),f(x_{1},x_{2}))}
and we operate as follows to go from
L
u
(
x
,
t
)
=
0
{\displaystyle {\mathcal {L}}u(x,t)=0}
to
L
v
(
x
1
,
x
2
)
=
0
{\displaystyle {\mathcal {L}}v(x_{1},x_{2})=0}
:
Apply the chain rule to
L
v
(
x
1
(
x
,
t
)
,
x
2
(
x
,
t
)
)
=
0
{\displaystyle {\mathcal {L}}v(x_{1}(x,t),x_{2}(x,t))=0}
and expand out giving equation
e
1
{\displaystyle e_{1}}
.
Substitute
a
(
x
,
t
)
{\displaystyle a(x,t)}
for
x
1
(
x
,
t
)
{\displaystyle x_{1}(x,t)}
and
b
(
x
,
t
)
{\displaystyle b(x,t)}
for
x
2
(
x
,
t
)
{\displaystyle x_{2}(x,t)}
in
e
1
{\displaystyle e_{1}}
and expand out giving equation
e
2
{\displaystyle e_{2}}
.
Replace occurrences of
x
{\displaystyle x}
by
e
(
x
1
,
x
2
)
{\displaystyle e(x_{1},x_{2})}
and
t
{\displaystyle t}
by
f
(
x
1
,
x
2
)
{\displaystyle f(x_{1},x_{2})}
to yield
L
v
(
x
1
,
x
2
)
=
0
{\displaystyle {\mathcal {L}}v(x_{1},x_{2})=0}
, which will be free of
x
{\displaystyle x}
and
t
{\displaystyle t}
.
References