लाप्लास रूपान्तर

लाप्लास रूपान्तर (Laplace transform) एक प्रकार का समाकल रूपान्तर (integral transform) है। यह भौतिकी एवं इंजीनियरी के अनेकानेक क्षेत्रों में प्रयोग किया जाता है। उदाहरण के लिए परिपथ विश्लेषण में। इसको ${\displaystyle \displaystyle {\mathcal {L}}\left\{f(t)\right\}}$ से निरूपित करते हैं। यह एक रैखिक संक्रिया है जो वास्तविक अर्गुमेन्ट t (t ≥ 0) वाले फलन f(t) को समिश्र अर्गुमेन्ट वाले फलन F(s) में बदल देता है।

लाप्लास रूपान्तर, प्रसिद्ध गणितज्ञ खगोलविद पिएर सिमों लाप्लास के नाम पर रखा गया है। लाप्लास रूपान्तर का उपयोग अवकल समीकरण तथा समाकल समीकरण (इंटीग्रल इक्वेशन) हल करने में किया जाता है।

${\displaystyle F(s)={\mathcal {L}}\left\{f(t)\right\}=\int _{0}^{\infty }e^{-st}f(t)\,dt.}$

अनुबन्ध यह है कि उपरोक्त समाकलन का अस्तित्व हो। उपरोक्त प्रकार से परिभाषित लाप्लास रूपान्तर 'एकपक्षीय लाप्लास रूपान्तर' कहलाता है। लाप्लास रूपान्तर का द्विपक्षीय रूपान्तर निम्नलिखित प्रकार से पारिभाषित किया जाता है-

${\displaystyle F_{B}(s)={\mathcal {L}}\left\{f(t)\right\}=\int _{-\infty }^{\infty }e^{-st}f(t)\,dt.}$

${\displaystyle {\mathcal {L}}\left\{af(t)+bg(t)\right\}=a{\mathcal {L}}\left\{f(t)\right\}+b{\mathcal {L}}\left\{g(t)\right\}}$

${\displaystyle {\mathcal {L}}\{f'(t)\}=s{\mathcal {L}}\{f(t)\}-f(0)}$

${\displaystyle {\mathcal {L}}\{f''(t)\}=s^{2}{\mathcal {L}}\{f(t)\}-sf(0)-f'(0)}$

${\displaystyle {\mathcal {L}}\left\{f^{(n)}(t)\right\}=s^{n}{\mathcal {L}}\{f(t)\}-s^{n-1}f(0)-\dots -f^{(n-1)}(0)}$ ${\displaystyle =s^{n}{\mathcal {L}}\{f(t)\}-\sum _{i=1}^{n}s^{n-i}f^{(i-1)}(0)}$

${\displaystyle {\mathcal {L}}\left\{\int _{0^{-}}^{t}f(\tau )d\tau \right\}={1 \over s}{\mathcal {L}}\{f\}}$

${\displaystyle {\mathcal {L}}\{tf(t)\}=-F'(s)}$

${\displaystyle {\mathcal {L}}\left\{e^{at}f(t)\right\}=F(s-a)}$

${\displaystyle {\mathcal {L}}\left\{f(t-a)u(t-a)\right\}=e^{-as}F(s)}$

${\displaystyle {\mathcal {L}}^{-1}\left\{e^{-as}F(s)\right\}=f(t-a)u(t-a)}$

टिप्पणी: ${\displaystyle u(t)}$ का अर्थ है यूनिट स्टेप फलन

${\displaystyle {\mathcal {L}}\{\,t^{n}f(t)\}=(-1)^{n}D_{s}^{n}[F(s)]}$

${\displaystyle {\mathcal {L}}\{f*g\}=F(s)G(s)}$

${\displaystyle {\mathcal {L}}\{f\}={1 \over 1-e^{-ps}}\int _{0}^{p}e^{-st}f(t)\,dt}$

${\displaystyle f(0^{+})=\lim _{s\to \infty }{sF(s)}}$

${\displaystyle f(\infty )=\lim _{s\to 0}{sF(s)}}$

फलन	समय डोमेन ${\displaystyle f(t)={\mathcal {L}}^{-1}\left\{F(s)\right\}}$	लाप्लास s-डोमेन ${\displaystyle F(s)={\mathcal {L}}\left\{f(t)\right\}}$	अभिसरण क्षेत्र (Region of convergence)	सन्दर्भ
यूनिट इम्पल्स	${\displaystyle \delta (t)\ }$	${\displaystyle 1}$	${\displaystyle \mathrm {all} \ s\,}$	inspection
delayed impulse	${\displaystyle \delta (t-\tau )\ }$	${\displaystyle e^{-\tau s}\ }$		time shift of unit impulse
unit step	${\displaystyle u(t)\ }$	${\displaystyle {1 \over s}}$	Re(s) > 0	integrate unit impulse
delayed unit step	${\displaystyle u(t-\tau )\ }$	${\displaystyle {e^{-\tau s} \over s}}$	Re(s) > 0	time shift of unit step
ramp	${\displaystyle t\cdot u(t)\ }$	${\displaystyle {\frac {1}{s^{2}}}}$	Re(s) > 0	integrate unit impulse twice
nth power (for integer n)	${\displaystyle t^{n}\cdot u(t)}$	${\displaystyle {n! \over s^{n+1}}}$	Re(s) > 0 (n > −1)	Integrate unit step n times
qth power (for complex q)	${\displaystyle t^{q}\cdot u(t)}$	${\displaystyle {\Gamma (q+1) \over s^{q+1}}}$	Re(s) > 0 Re(q) > −1	[1][2]
nth root	${\displaystyle {\sqrt[{n}]{t}}\cdot u(t)}$	${\displaystyle {\Gamma ({\frac {1}{n}}+1) \over s^{{\frac {1}{n}}+1}}}$	Re(s) > 0	Set q = 1/n above.
nth power with frequency shift	${\displaystyle t^{n}e^{-\alpha t}\cdot u(t)}$	${\displaystyle {\frac {n!}{(s+\alpha )^{n+1}}}}$	Re(s) > −α	Integrate unit step, apply frequency shift
delayed nth power with frequency shift	${\displaystyle (t-\tau )^{n}e^{-\alpha (t-\tau )}\cdot u(t-\tau )}$	${\displaystyle {\frac {n!\cdot e^{-\tau s}}{(s+\alpha )^{n+1}}}}$	Re(s) > −α	Integrate unit step, apply frequency shift, apply time shift
exponential decay	${\displaystyle e^{-\alpha t}\cdot u(t)\ }$	${\displaystyle {1 \over s+\alpha }}$	Re(s) > −α	Frequency shift of unit step
two-sided exponential decay	${\displaystyle e^{-\alpha |t|}\ }$	${\displaystyle {2\alpha \over \alpha ^{2}-s^{2}}}$	−α < Re(s) < α	Frequency shift of unit step
exponential approach	${\displaystyle (1-e^{-\alpha t})\cdot u(t)\ }$	${\displaystyle {\frac {\alpha }{s(s+\alpha )}}}$	Re(s) > 0	Unit step minus exponential decay
sine	${\displaystyle \sin(\omega t)\cdot u(t)\ }$	${\displaystyle {\omega \over s^{2}+\omega ^{2}}}$	Re(s) > 0	Bracewell 1978, पृष्ठ 227
cosine	${\displaystyle \cos(\omega t)\cdot u(t)\ }$	${\displaystyle {s \over s^{2}+\omega ^{2}}}$	Re(s) > 0	Bracewell 1978, पृष्ठ 227
hyperbolic sine	${\displaystyle \sinh(\alpha t)\cdot u(t)\ }$	${\displaystyle {\alpha \over s^{2}-\alpha ^{2}}}$	Re(s) > |α|	Williams 1973, पृष्ठ 88
hyperbolic cosine	${\displaystyle \cosh(\alpha t)\cdot u(t)\ }$	${\displaystyle {s \over s^{2}-\alpha ^{2}}}$	Re(s) > |α|	Williams 1973, पृष्ठ 88
exponentially decaying sine wave	${\displaystyle e^{-\alpha t}\sin(\omega t)\cdot u(t)\ }$	${\displaystyle {\omega \over (s+\alpha )^{2}+\omega ^{2}}}$	Re(s) > −α	Bracewell 1978, पृष्ठ 227
exponentially decaying cosine wave	${\displaystyle e^{-\alpha t}\cos(\omega t)\cdot u(t)\ }$	${\displaystyle {s+\alpha \over (s+\alpha )^{2}+\omega ^{2}}}$	Re(s) > −α	Bracewell 1978, पृष्ठ 227
natural logarithm	${\displaystyle \ln(t)\cdot u(t)}$	${\displaystyle -{1 \over s}\,\left[\ln(s)+\gamma \right]}$	Re(s) > 0	Williams 1973, पृष्ठ 88
Bessel function of the first kind, of order n	${\displaystyle J_{n}(\omega t)\cdot u(t)}$	${\displaystyle {\frac {\left({\sqrt {s^{2}+\omega ^{2}}}-s\right)^{n}}{\omega ^{n}{\sqrt {s^{2}+\omega ^{2}}}}}}$	Re(s) > 0 (n > −1)	Williams 1973, पृष्ठ 89
Error function	${\displaystyle \mathrm {erf} (t)\cdot u(t)}$	${\displaystyle {e^{s^{2}/4}\left(1-\operatorname {erf} \left(s/2\right)\right) \over s}}$	Re(s) > 0	Williams 1973, पृष्ठ 89
Explanatory notes: u(t) represents the Heaviside step function. ${\displaystyle \delta (t)\,}$ represents the Dirac delta function. Γ(z) represents the Gamma function. γ is the Euler–Mascheroni constant. t, एक वास्तविक संख्या है जो मोटे तौर पर "समय" को निरूपित करती है किन्तु t कोई भी स्वतंत्र राशि हो सकती है। s is the complex angular frequency, and Re(s) is its real part. α, β, τ, and ω are real numbers. n is an integer.

प्रतिलोम लाप्लास रूपान्तर (inverse Laplace transform) नीचे दिए गए समिश्र समाकल द्वारा निकाला जा सकता है। इस समाकल के कई नाम हैं, जैसे ब्रोमविच समाकल (Bromwich integra), फुर्ये-मेलिन समाकल (Fourier–Mellin integral) या मेलिन का प्रतिलोम सुत्र (Mellin's inverse formula):

${\displaystyle f(t)={\mathcal {L}}^{-1}\{F\}(t)={\frac {1}{2\pi i}}\lim _{T\to \infty }\int _{\gamma -iT}^{\gamma +iT}e^{st}F(s)\,ds,}$

जहाँ γ एक वास्तविक संख्या है ताकि समाकल का कन्टूर-पथ कन्वर्जेन्स के क्षेत्र F(s) के अन्दर हो। प्रतिलोम लाप्लास रूपान्तर निकालने का एक दूसरा सूत्र पोस्ट का प्रतिलोम सूत्र (Post's inversion formula) है।

नाभिकीय भौतिकी में जरेडियोसक्रिय क्षय को अभिव्यक्त करने वाला अवकल समीकरण नीचे दिया गया है। किसी नमूने में रेडियोसक्रिय परमाणुओं की संख्या N है तथा इसके क्षय की दर N के समानुपाती होती है। इसी को निम्नलिखित अवकल समीकरण द्वारा अभिव्यक्त किया जा सकता है-

${\displaystyle {\frac {dN}{dt}}=-\lambda N,}$

जहाँ λ, क्षय नियतांक (decay constant) है। इस समीकरण का हल लाप्लास रूपान्तर की सहायता से निकाला जा सकता है।

इस समीकरण को एक ही पक्ष (side) में ले जाकर लिखने पर,

${\displaystyle {\frac {dN}{dt}}+\lambda N=0.}$

अब हम इस समीकरण के दोनों पक्षों का लाप्लास रूपान्तर लेते हैं।

${\displaystyle \left(s{\tilde {N}}(s)-N_{0}\right)+\lambda {\tilde {N}}(s)=0,}$

जहाँ

${\displaystyle {\tilde {N}}(s)={\mathcal {L}}\{N(t)\}}$

तथा

${\displaystyle N_{0}=N(0).}$

इसका हल करने पर,

${\displaystyle {\tilde {N}}(s)={\frac {N_{0}}{s+\lambda }}.}$

अन्त में हम प्रतिलोम लाप्लास रूपान्तर लेते हें जिससे सामान्य हल प्राप्त होता है।

${\displaystyle {\begin{aligned}N(t)&={\mathcal {L}}^{-1}\{{\tilde {N}}(s)\}={\mathcal {L}}^{-1}\!\left\{{\frac {N_{0}}{s+\lambda }}\right\}\\&=\ N_{0}e^{-\lambda t},\end{aligned}}}$

क्षणिक परिपथों के विश्लेषण में प्रायः लाप्लास रूपान्तर का उपयोग किया जाता है। इसके लिए परिपथ के अवयवों को s-डोमेन में बदलकर काम को आगे बढ़ाते हैं। नीचे तुल्य परिपथ दिए गये हैं-

पार्श्व चित्र को देखें जिसमें दो लूप हैं। इनमें बहने वाली धारा ${\displaystyle i_{1}}$ तथा ${\displaystyle i_{2}}$ चित्र में दर्शायी गयी हैं। माना ${\displaystyle i_{1}}$ तथा ${\displaystyle i_{2}}$ के आरम्भिक मान शून्य हैं, अर्थात् ${\displaystyle i_{1}(0)=0}$ और ${\displaystyle i_{2}(0)=0}$। किरखॉफ के नियम के अनुसार,

${\displaystyle {di_{1}(t) \over dt}+5i_{1}(t)+40i(t)=110}$ (1)

${\displaystyle 2{di_{2}(t) \over dt}+10i_{2}(t)+40i(t)=110}$ (2)

चित्र से स्पष्ट है कि ${\displaystyle i(t)=i_{1}(t)+i_{2}(t)}$,

${\displaystyle {di_{1}(t) \over dt}+45i_{1}(t)+40i_{2}(t)=110}$

${\displaystyle {di_{2}(t) \over dt}+20i_{2}(t)+25i_{2}(t)=55}$ (समीकरण (2) को 2 से भाग देने पर)

इन पर लाप्लास रूपान्तर लगाने पर,

${\displaystyle sI_{1}(s)-i_{1}(0)+45I_{1}(s)+40I_{2}(s)={\frac {110}{s}}}$

${\displaystyle sI_{2}(s)-i_{2}(0)+20I_{1}(s)+25I_{2}(s)={\frac {55}{s}}}$

या,

${\displaystyle (s+45)I_{1}(s)+40I_{2}(s)={\frac {110}{s}}}$

${\displaystyle 20I_{1}(s)+(s+25)I_{2}(s)={\frac {55}{s}}}$

या,

${\displaystyle {\begin{bmatrix}(s+45)&40\\20&(s+25)\end{bmatrix}}{\begin{bmatrix}I_{1}(s)\\I_{2}(s)\end{bmatrix}}={\begin{bmatrix}110/s\\55/s\end{bmatrix}}}$

इसका हल निम्नलिखित है:

${\displaystyle {\begin{bmatrix}I_{1}(s)\\I_{2}(s)\end{bmatrix}}={\frac {1}{(s+25)(s+45)-800}}{\begin{bmatrix}(s+25)&-40\\-20&(s+45)\end{bmatrix}}{\begin{bmatrix}110/s\\55/s\end{bmatrix}}}$

अतः,

${\displaystyle I_{1}(s)={\frac {1}{s^{2}+70s+325}}{\biggl (}{\frac {110}{s}}(s+25)-{\frac {2200}{s}}{\biggr )}={\frac {1}{(s+5)(s+65)}}{\biggl (}110+{\frac {550}{s}}{\biggr )}}$

${\displaystyle I_{2}(s)={\frac {1}{s^{2}+70s+325}}{\biggl (}-{\frac {2200}{s}}+{\frac {55}{s}}(s+45){\biggr )}={\frac {1}{(s+5)(s+65)}}{\biggl (}55+{\frac {275}{s}}{\biggr )}}$

ध्यान दें कि ${\displaystyle I_{1}(s)=2I_{2}(s)}$ अतः हम केवल ${\displaystyle I_{2}(s)}$ की गणना ही करेंगे।

${\displaystyle I_{2}(s)={\frac {1}{(s+5)(s+65)}}{\biggl (}{\frac {55s+275}{s}}{\biggr )}={\frac {55}{(s+5)(s+65)}}{\biggl (}{\frac {s+5}{s}}{\biggr )}={\frac {55}{s(s+65)}}}$

इससे,

${\displaystyle i_{2}(t)={\frac {55}{65}}(1-e^{-65t})={\frac {11}{13}}(1-e^{-65t})}$

चूंकि ${\displaystyle i_{1}(t)=2i_{2}(t)}$, अतः

${\displaystyle i_{2}(t)={\frac {22}{13}}(1-e^{-65t})}$

1. ↑ Mathematical Handbook of Formulas and Tables (3rd edition), S. Lipschutz, M.R. Spiegel, J. Liu, Schuam's Outline Series, p.183, 2009, ISBN 978-0-07-154855-7 - provides the case for real q.
2. ↑ http://mathworld.wolfram.com/LaplaceTransform.html Archived 2013-01-30 at the वेबैक मशीन - Wolfram Mathword provides case for complex q

• फुर्ये रूपान्तर
• जेड रूपान्तर (Z-transform)
• अंतरण प्रकार्य (transfer function)