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

टाइम डोमेन और फ्रेक्वेन्सी डोमेन

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

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

परिभाषा

$F(s)=ℒ\{f(t)\}=\underset{0}{\overset{\mathrm{\infty }}{\int }}{e}^{-st}f(t)dt.$

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

$F_B(s)=ℒ\{f(t)\}=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{-st}f(t)dt.$

गुण

रैखिकता

$ℒ\{af(t)+bg(t)\}=aℒ\{f(t)\}+bℒ\{g(t)\}$

अवकलन

$ℒ\{f^{\prime }(t)\}=sℒ\{f(t)\}-f(0)$

$ℒ\{f^{″}(t)\}=s^2ℒ\{f(t)\}-sf(0)-f^{\prime }(0)$

$ℒ\{f^{(n)}(t)\}=s^nℒ\{f(t)\}-s^{n-1}f(0)-\dots -f^{(n-1)}(0)$ $=s^nℒ\{f(t)\}-\sum _{i=1}^n{s}^{n-i}{f}^{(i-1)}(0)$

समाकलन

$ℒ\{\underset{0^{-}}{\overset{t}{\int }}f(\tau )d\tau \}=\frac{1}{s}ℒ\{f\}$

द्वैत

$ℒ\{tf(t)\}=-F^{\prime }(s)$

आवृत्ति विस्थापन

$ℒ\{e^{at}f(t)\}=F(s-a)$

समय विस्थापन

$ℒ\{f(t-a)u(t-a)\}=e^{-as}F(s)$

${ℒ}^{-1}\{e^{-as}F(s)\}=f(t-a)u(t-a)$

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

समय के n-घात से गुणा

$ℒ\{t^{n}f(t)\}=(-1{)}^n{D}_s^n[F(s)]$

संवलन (कॉन्वोलुशन)

$ℒ\{f*g\}=F(s)G(s)$

p आवर्तकाल वाले एक आवर्ती फलन का लाप्लास रूपान्तर

$ℒ\{f\}=\frac{1}{1-e^{-ps}}\underset{0}{\overset{p}{\int }}{e}^{-st}f(t)dt$

प्रारम्भिक मान प्रमेय

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

अन्तिम मान प्रमेय

$f(\mathrm{\infty })=\underset{s\to 0}{\lim }sF(s)$

प्रमुख फलनों के लाप्लास रूपान्तर

फलन	समय डोमेन $f(t)={ℒ}^{-1}\{F(s)\}$	लाप्लास s-डोमेन $F(s)=ℒ\{f(t)\}$	अभिसरण क्षेत्र (Region of convergence)	सन्दर्भ
यूनिट इम्पल्स	$\delta (t)$	$1$	$\mathrm{a}\mathrm{l}\mathrm{l}s$	inspection
delayed impulse	$\delta (t-\tau )$	$e^{-\tau s}$		time shift of unit impulse
unit step	$u(t)$	$\frac{1}{s}$	Re(s) > 0	integrate unit impulse
delayed unit step	$u(t-\tau )$	$\frac{e^{-\tau s}}{s}$	Re(s) > 0	time shift of unit step
ramp	$t\cdot u(t)$	$\frac{1}{s^2}$	Re(s) > 0	integrate unit impulse twice
nth power (for integer n)	$t^n\cdot u(t)$	$\frac{n!}{s^{n+1}}$	Re(s) > 0 (n > −1)	Integrate unit step n times
qth power (for complex q)	$t^q\cdot u(t)$	$\frac{\mathrm{\Gamma }(q+1)}{s^{q+1}}$	Re(s) > 0 Re(q) > −1
nth root	$\sqrt[n]{t}\cdot u(t)$	$\frac{\mathrm{\Gamma }(\frac{1}{n}+1)}{s^{\frac{1}{n}+1}}$	Re(s) > 0	Set q = 1/n above.
nth power with frequency shift	$t^n{e}^{-\alpha t}\cdot u(t)$	$\frac{n!}{(s+\alpha {)}^{n+1}}$	Re(s) > −α	Integrate unit step, apply frequency shift
delayed nth power with frequency shift	$(t-\tau {)}^n{e}^{-\alpha (t-\tau )}\cdot u(t-\tau )$	$\frac{n!\cdot e^{-\tau s}}{(s+\alpha {)}^{n+1}}$	Re(s) > −α	Integrate unit step, apply frequency shift, apply time shift
exponential decay	$e^{-\alpha t}\cdot u(t)$	$\frac{1}{s+\alpha }$	Re(s) > −α	Frequency shift of unit step
two-sided exponential decay	$e^{-\alpha |t|}$	$\frac{2\alpha }{{\alpha }^2-s^2}$	−α < Re(s) < α	Frequency shift of unit step
exponential approach	$(1-e^{-\alpha t})\cdot u(t)$	$\frac{\alpha }{s(s+\alpha )}$	Re(s) > 0	Unit step minus exponential decay
sine	$\sin (\omega t)\cdot u(t)$	$\frac{\omega }{s^2+{\omega }^2}$	Re(s) > 0	साँचा:Harvnb
cosine	$\cos (\omega t)\cdot u(t)$	$\frac{s}{s^2+{\omega }^2}$	Re(s) > 0	साँचा:Harvnb
hyperbolic sine	$\sinh (\alpha t)\cdot u(t)$	$\frac{\alpha }{s^2-{\alpha }^2}$	Re(s) > |α|	साँचा:Harvnb
hyperbolic cosine	$\cosh (\alpha t)\cdot u(t)$	$\frac{s}{s^2-{\alpha }^2}$	Re(s) > |α|	साँचा:Harvnb
exponentially decaying sine wave	$e^{-\alpha t}\sin (\omega t)\cdot u(t)$	$\frac{\omega }{(s+\alpha {)}^2+{\omega }^2}$	Re(s) > −α	साँचा:Harvnb
exponentially decaying cosine wave	$e^{-\alpha t}\cos (\omega t)\cdot u(t)$	$\frac{s+\alpha }{(s+\alpha {)}^2+{\omega }^2}$	Re(s) > −α	साँचा:Harvnb
natural logarithm	$\ln (t)\cdot u(t)$	$-\frac{1}{s}[\ln (s)+\gamma ]$	Re(s) > 0	साँचा:Harvnb
Bessel function of the first kind, of order n	$J_n(\omega t)\cdot u(t)$	$\frac{{(\sqrt{s^2+{\omega }^2}-s)}^n}{{\omega }^n\sqrt{s^2+{\omega }^2}}$	Re(s) > 0 (n > −1)	साँचा:Harvnb
Error function	$\mathrm{e}\mathrm{r}\mathrm{f}(t)\cdot u(t)$	$\frac{e^{s^2/4}(1-\mathrm{erf}(s/2))}{s}$	Re(s) > 0	साँचा:Harvnb
Explanatory notes: u(t) represents the Heaviside step function. $\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):

$f(t)={ℒ}^{-1}\{F\}(t)=\frac{1}{2\pi i}\underset{T\to \mathrm{\infty }}{\lim }\underset{\gamma -iT}{\overset{\gamma +iT}{\int }}{e}^{st}F(s)ds,$

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

उपयोग

उदाहरण १: किसी अवकल समीकरण का हल निकालना

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

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

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

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

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

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

$(s\tilde{N}(s)-N_0)+\lambda \tilde{N}(s)=0,$

जहाँ

$\tilde{N}(s)=ℒ\{N(t)\}$

तथा

$N_0=N(0).$

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

$\tilde{N}(s)=\frac{N_0}{s+\lambda }.$

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

$\begin{array}{cc}N(t)& ={ℒ}^{-1}\{\tilde{N}(s)\}={ℒ}^{-1}\left\{\frac{N_0}{s+\lambda }\right\}\\ & =N_0{e}^{-\lambda t},\end{array}$

s-डोमेन में तुल्य परिपथ और प्रतिबाधाएँ

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

साँचा:Math-डोमेन तुल्य परिपथ

अनुप्रयोग : दो लूप वाले एक परिपथ की क्षणिक अनुक्रिया (रिस्पॉन्स)

दो लूप वाला एक परिपथ

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

$\frac{d{i}_1(t)}{dt}+5i_1(t)+40i(t)=110$ (1)

$2\frac{d{i}_2(t)}{dt}+10i_2(t)+40i(t)=110$ (2)

चित्र से स्पष्ट है कि $i(t)=i_1(t)+i_2(t)$,

$\frac{d{i}_1(t)}{dt}+45i_1(t)+40i_2(t)=110$

$\frac{d{i}_2(t)}{dt}+20i_2(t)+25i_2(t)=55$ (समीकरण (2) को 2 से भाग देने पर)

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

$s{I}_1(s)-i_1(0)+45I_1(s)+40I_2(s)=\frac{110}{s}$

$s{I}_2(s)-i_2(0)+20I_1(s)+25I_2(s)=\frac{55}{s}$

या,

$(s+45)I_1(s)+40I_2(s)=\frac{110}{s}$

$20I_1(s)+(s+25)I_2(s)=\frac{55}{s}$

या,

$\left[\begin{array}{cc}(s+45)& 40\\ 20& (s+25)\end{array}\right]\left[\begin{array}{c}{I}_1(s)\\ I_2(s)\end{array}\right]=\left[\begin{array}{c}110/s\\ 55/s\end{array}\right]$

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

$\left[\begin{array}{c}{I}_1(s)\\ I_2(s)\end{array}\right]=\frac{1}{(s+25)(s+45)-800}\left[\begin{array}{cc}(s+25)& -40\\ -20& (s+45)\end{array}\right]\left[\begin{array}{c}110/s\\ 55/s\end{array}\right]$

अतः,

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

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

ध्यान दें कि $I_1(s)=2I_2(s)$ अतः हम केवल $I_2(s)$ की गणना ही करेंगे।

$I_2(s)=\frac{1}{(s+5)(s+65)}\left(\frac{55s+275}{s}\right)=\frac{55}{(s+5)(s+65)}\left(\frac{s+5}{s}\right)=\frac{55}{s(s+65)}$

इससे,

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

चूंकि $i_1(t)=2i_2(t)$, अतः

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

इन्हें भी देखें

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