Mostrar código-fonte de Delta de Dirac

[[Ficheiro:Dirac distribution PDF.png|thumb|336x336px|Representação geométrica da função Delta de Dirac]]
Em [[matemática]], a função '''delta de Dirac''', também conhecida como função δ, é uma [[Teoria das distribuições|distribuição]] na [[reta real]], a qual vale infinito no ponto zero e é nula no restante da reta. A [[integral]] da '''função Delta de Dirac''' em toda reta é definida como tendo valor 1. Foi introduzida pelo físico teórico [[Paul Dirac]] em [[1930]] em seu livro ‘’''The Principles of Quantum Mechanics’''’.<ref>{{citar livro|nome = Paul|sobrenome = Dirac|título = Principles os Quantum Mechanics|ano = 1930|isbn = }}</ref> Seu
análogo no domínio discreto é o [[delta de Kronecker]],o qual vale 0 e 1.

Pode-se pensar no '''Delta de Dirac''' como um retângulo infinitamente estreito e infinitamente alto, com área igual à unidade. Em muitos casos, pode ser encarado como o [[limite]] de funções que tendem a estas condições.Além disso, se enfocarmos no contexto de processamento de sinais, ela é frequentemente interpretada como um impulso unitário.
[[Ficheiro:Dirac function approximation.gif|miniaturadaimagem|A função delta de Dirac como limite ( no sentido de distribuição ) da sequência da distribuição normal com centro em zero.
<math> \delta_{a}(x)=\tfrac {1}{\sqrt{\pi}a} \cdot e^{-\frac {x^2}{a^2}}</math>
]]
Matematicamente, o '''Delta de Dirac''' não pode ser caracterizado propriamente como uma [[Função (matemática)|função]], mas sim como um objeto matemático. Isso porque qualquer função que valha zero em todos os pontos exceto um, deve ter integral nula em toda a reta. Entretanto, quando em uma integral, ganha sentido matemático e, para a maioria dos propósitos, pode ser encarado e manipulado como uma função.

== Visão Geral ==
Apesar do nome, a '''função Delta de Dirac''' não é verdadeiramente uma função. '''''δ''(''x'')''' difere de uma função g(x)=0 apenas em zero, e portanto as duas são ''funções iguais quase sempre''. Segundo a teoria de integração de [[Henri Lebesgue|Lebesgue]], duas funções integráveis e iguais quase sempre, têm,necessariamente, integrais idênticas. Entretanto, a integral da função g(x) é zero, e a da '''função Delta de Dirac''' vale 1 por definição.

Desconsiderando o formalismo matemático, em termos práticos o Delta pode ser encarado como uma [[Função (matemática)|função]] na maioria dos casos. Geralmente é usado para modelar uma função impulso, isto é, modelar situações pontuais (como massas e cargas pontuais) ou instantâneas (adição instantânea de uma substância em uma [[reação química]] ou a dinâmica de um martelo batendo  fortemente em um sistema massa mola, por exemplo). 

== História ==
[[Jean Baptiste Joseph Fourier|Joseph Fourier]] nos apresentou o que é chamado nos dias de hoje  de “Teorema de integral de Fourier”  na sua obra “ Théorie analytique de la chaleur” na forma: <ref name=Fourier>{{citar livro|título=The Analytical Theory of Heat |url=https://books.google.com/books?id=-N8EAAAAYAAJ&pg=PA408&dq=%22when+the+integrals+are+taken+between+infinite+limits%22+%22that+is+to+say,+that+we+have+the+equation%22&hl=en&sa=X&ei=rFe-T96cEIzKiQKcyfDtDQ&ved=0CD4Q6AEwAA#v=onepage&q=%22when%20the%20integrals%20are%20taken%20between%20infinite%20limits%22%20%22that%20is%20to%20say%2C%20that%20we%20have%20the%20equation%22&f=false |autor =JB Fourier |ano=1822 |página=408 |edição= English translation by Alexander Freeman, 1878 |publicado=The University Press}}, cf p 449 and pp 546-551.  The original French text can be found [https://books.google.com/books?id=TDQJAAAAIAAJ&pg=PA525&dq=%22c%27est-%C3%A0-dire+qu%27on+a+l%27%C3%A9quation%22&hl=en&sa=X&ei=SrC7T9yKBorYiALVnc2oDg&sqi=2&ved=0CEAQ6AEwAg#v=onepage&q=%22c%27est-%C3%A0-dire%20qu%27on%20a%20l%27%C3%A9quation%22&f=false '''here'''].</ref>
:[[Ficheiro:Fourier2.jpg|miniaturadaimagem|Jean-Baptiste Joseph Fourier]]<math>f(x)=\frac{1}{2\pi}\int_{-\infty}^\infty\ \ d\alpha f(\alpha) \ \int_{-\infty}^\infty dp\ \cos  (px-p\alpha)\ , </math>

A equação acima equivale a introdução da função- δ no seguinte formato: <ref name= Kawai>{{citar livro|título=Microlocal Analysis and Complex Fourier Analysis |editor1=Takahiro Kawai |editor2=Keiko Fujita |autor =Hikosaburo Komatsu |capítulo=Fourier's hyperfunctions and Heaviside's pseudodifferential operators |isbn=981-238-161-9 |ano=2002 |publicado=World Scientific |url=https://books.google.com/books?id=8GwKzEemrIcC&pg=PA200&dq=%22Fourier+introduced+the%22+%22+-function+much+earlier%22&hl=en&sa=X&ei=oJa6T5L2O6SriQKGloCUBw&ved=0CDQQ6AEwAA#v=onepage&q=%22Fourier%20introduced%20the%22%20%22%20-function%20much%20earlier%22&f=false |página=200 }}</ref>
:<math>\delta(x-\alpha)=\frac{1}{2\pi} \int_{-\infty}^\infty dp\ \cos  (px-p\alpha) \ . </math>
Anos mais tarde o matemático francês, [[Augustin-Louis Cauchy]] expressou o mesmo teorema no formato de funções exponenciais: <ref name= Myint-U>{{citar livro
|url=https://books.google.com/books?id=Zbz5_UvERIIC&pg=PA4&dq=%22It+was+the+work+of+Augustin+Cauchy%22&hl=en&sa=X&ei=RnW6T52LNovYiQLa9-mABw&ved=0CDgQ6AEwAA#v=onepage&q=%22It%20was%20the%20work%20of%20Augustin%20Cauchy%22&f=false  |autor1 =Tyn Myint-U. |autor2 =Lokenath Debnath |título=Linear Partial Differential Equations for Scientists And Engineers |isbn=0-8176-4393-1 |edição=4th  |ano=2007  |página=4 |publicado=Springer}}</ref><ref name=Debnath>{{citar livro|url=https://books.google.com/books?id=WbZcqdvCEfwC&pg=PA2&dq=%22It+was+the+work+of+Cauchy+that+contained%22&hl=en&sa=X&ei=Jym9T8L-NK6OigK-m_GYDg&ved=0CDQQ6AEwAA#v=onepage&q=%22It%20was%20the%20work%20of%20Cauchy%20that%20contained%22&f=false |título=Integral Transforms And Their Applications |autor1 =Lokenath Debnath |autor2 =Dambaru Bhatta |isbn=1-58488-575-0 |ano=2007 |edição=2nd |publicado=CRC Press |página=2}}</ref>
:<math>f(x)=\frac{1}{2\pi} \int_{-\infty} ^ \infty \ e^{ipx}\left(\int_{-\infty}^\infty e^{-ip\alpha }f(\alpha)\ d \alpha \right) \ dp. </math>
Cauchy percebeu que em alguns casos a ordem da integração teria um resultado significativo no  resultado final, conseguindo comprovar esta premissa por meio da  teoria das distribuições. <ref name=Grattan-Guinness>{{citar livro|título=Convolutions in French Mathematics, 1800–1840: From the Calculus and Mechanics to Mathematical Analysis and Mathematical Physics, Volume 2 |página=653 |url= https://books.google.com/books?id=_GgioErrbW8C&pg=PA653&dq=%22Further,+in+a+double+integral%22&hl=en&sa=X&ei=4gC9T7KVDvDRiALq-dTLDQ&ved=0CDgQ6AEwAA#v=onepage&q=%22Further%2C%20in%20a%20double%20integral%22&f=false |isbn=3-7643-2238-1 |ano=2009 |publicado=Birkhäuser |autor =Ivor Grattan-Guinness}}</ref>
A equação encontrada por Cauchy pode ser organizada para se assemelhar com a formula inicial encontrada por Fourier, representando a função- δ como:
:<math>\begin{align}
f(x)&=\frac{1}{2\pi} \int_{-\infty}^\infty e^{ipx}\left(\int_{-\infty}^\infty e^{-ip\alpha }f(\alpha)\ d \alpha \right) \ dp \\
&=\frac{1}{2\pi} \int_{-\infty}^\infty \left(\int_{-\infty}^\infty e^{ipx} e^{-ip\alpha } \ dp \right)f(\alpha)\ d \alpha =\int_{-\infty}^\infty \delta (x-\alpha) f(\alpha) \ d \alpha,
\end{align}</math>[[Ficheiro:Augustin-Louis Cauchy.jpg|miniaturadaimagem| Augustin-Louis Cauchy,1857. ]]
Na qual a função- δ e representada como:
:<math>\delta(x-\alpha)=\frac{1}{2\pi} \int_{-\infty}^\infty e^{ip(x-\alpha)}\ dp \ . </math>

A interpretação rigorosa da forma exponencial encontrada por Cauchy e como também as varias limitações sobre a função f necessárias para suas aplicações se estenderam por muitos séculos. Os problemas relatados pela interpretação clássica serão explicados a seguir: <ref name="Mitrović">{{citar livro|título=Fundamentals of Applied Functional Analysis: Distributions, Sobolev Spaces |autor1 =Dragiša Mitrović |autor2 =Darko Žubrinić |url=https://books.google.com/books?id=Od5BxTEN0VsC&pg=PA62&dq=%22greatest+drawback+of+the+classical+Fourier+transformation+is+a+rather+narrow+class+of+functions%22&hl=en&sa=X&ei=IKG6T_niFqWfiQLJoODdBg&ved=0CDQQ6AEwAA#v=onepage&q=%22greatest%20drawback%20of%20the%20classical%20Fourier%20transformation%20is%20a%20rather%20narrow%20class%20of%20functions%22&f=false |página=62 |isbn=0-582-24694-6 |ano=1998 |publicado=CRC Press}}</ref>
O maior inconveniente da classica transformacao de Fourier é a quantidade muito restrita de funções ( originais ) que podem ser efetivamente computadas. Ou seja, é necessário que essas funções decrescam suficientemente rápido para zero (próximo de infinito)  para assegurar que exista a integral de Fourier. Por exemplo, a transformação de Fourier para funções simples como as polinomiais não existe na forma clássica.

A extensão da transformação de Fourier clássica foi importante para que houvesse uma forma de aumentar o número de funções, que até então não eram muitas, para um número considerável de classes a serem transformadas. Tal ajudou a superar muitos obstáculos vivenciados por matemáticos da época.

Desenvolvimentos posteriores incluíram generalizações a  integral de Fourier, “inicialmente com a  teoria de Plancherel “pathbreaking L2-theory (1910)” , sucessivamente pelo trabalho de Wiener e  Bochner ( em meados de 1930) e culminando com a amalgamação para a teoria de distribuição de  L. Schwartz ( 1945 )”, <ref name=Kracht>{{citar livro|título=Topics in Mathematical Analysis: A Volume Dedicated to the Memory of A.L. Cauchy |url=https://books.google.com/books?id=xIsPrSiDlZIC&pg=PA553&dq=%22To+this+theory%22+%22and+even+more%22++%22that+one+was+able+to+generalize%22&hl=en&sa=X&ei=RJ66T-y7JOLjiAKuoeSUBw&ved=0CDQQ6AEwAA#v=onepage&q=%22To%20this%20theory%22%20%22and%20even%20more%22%20%20%22that%20one%20was%20able%20to%20generalize%22&f=false |autor1 =Manfred Kracht |autor2 =Erwin Kreyszig |página=553 |isbn=9971-5-0666-1 |editor=Themistocles M. Rassias |ano=1989 |publicado=World Scientific |capítulo=On singular integral operators and generalizations}}</ref>  guiando assim para o desenvolvimento da função delta de Dirac.

Em 1827, em um texto de Augustin Louis Cauchy,<ref>{{harvnb|Laugwitz|1989|p=230}}</ref> uma fórmula [[infinitesimal]] com altura infinita para uma função delta de impulso unitário  ( uma versão infinitesimal da distribuição de Cauchy) foi proposta.

Siméon Denis Poisson considerou esta questão apropriadamente associada com o estudo da propagacao de ondas como foi  também foi vista por [[Gustav Kirchhoff]]  posteriormente. Kirchhoff e [[Hermann von Helmholtz]]  também introduziram o impulse unitário como um limite gaussiano, que também corresponde ao ponto encontrado por Lord Kelvin como uma fonte de calor. No final do século dezenove, [[Oliver Heaviside]] usou a serie de Fourier para manipular o impulso unitário. <ref>A more complete historical account can be found in {{harvnb|van der Pol|Bremmer|1987|loc=§V.4}}.</ref> A função delta de Dirac como e conhecida nos dias de hoje foi introduzida como uma “notação conveniente” por Paul Dirac em seu livro The Principles of Quantum Mechanics,1930. Ele a chamou de “função delta”, pois foi analogamente encontrada pela função discreta de Kronecker delta.

== Definição<ref>{{citar livro|nome = Esequia|sobrenome = Sauter|título = Notas de Aula|ano = 2015|isbn = }}</ref> ==

Para definir matematicamente a função Delta de Dirac, é conveniente utilizar outra função singular, a chamada [[função de Heaviside]], também conhecida por função salto ou degrau. Considera-se o Delta como o limite de funções pulso unitário <math>f_\epsilon(t)</math>,  num curto intervalo de tempo <math>[-\epsilon,\epsilon]</math>, quando o parâmetro <math>\epsilon</math> tende a zero:

<math>f_\epsilon(t)=\frac{1}{2\epsilon}[u(t+\epsilon) - u(t-\epsilon)]=\begin{cases}
0 &t<-\epsilon\\
\frac{1}{2\epsilon} & -\epsilon<t<\epsilon\\
0 & t>\epsilon
\end{cases}</math>

Para representar um impulso num instante de tempo que não zero, realiza-se um deslocamento na equação anterior:

<math>f_\epsilon(t-a)=\frac{1}{2\epsilon}[u(t-(a-\epsilon)) - u(t-(a+\epsilon))]=\begin{cases}
0 &t<a-\epsilon\\
\frac{1}{2\epsilon} & a-\epsilon<t<a+\epsilon\\
0 & t>a+\epsilon
\end{cases}</math>

Observa-se que as áreas dos retângulos formados pelas funções pulso unitário <math>f_\epsilon(t)</math> valem 1 independentemente do valor de <math>\epsilon</math>. Ou seja:

<math>A=\int_{a-\epsilon}^{a+\epsilon} f_\epsilon(t-a)dt = \int_{a-\epsilon}^{a+\epsilon} \frac{1}{2\epsilon}dt = 1 </math>

Assim, considerando que quando <math> \epsilon \rightarrow 0 \quad f_\epsilon(t-a) \rightarrow \infty </math>, define-se a função Delta de Dirac como sendo

<math> \delta(t-a)=\lim_{\epsilon \to 0} f_\epsilon(t-a) </math>

A função Delta de Dirac pode também ser definida em termos de outras funções com propriedades análogas. Por exemplo:

<math> \delta(t)=\lim_{\epsilon \to 0} \frac{e^\frac{-t^2}{\epsilon^2}}{\epsilon\sqrt{\pi}} </math>

Resumidamente, pode se escrever:
: <math>\delta(x-x_0) = \begin{cases} \infty, & x = x_0 \\ 0, & x \ne x_0 \end{cases}</math>
: <math>\int_{-\infty}^{\infty} \delta(x) dx = 1</math>

É possível também encarar a distribuição Delta como um limite de [[aproximação da identidade|aproximações da identidade]].

== Propriedades<ref name=":0">Notas de Aula da Prof. Irene Strauch, UFRGS</ref> ==

=== Propriedade da Filtragem ===
Uma propriedade interessante da Delta de Dirac é a propriedade da filtragem:

: <math>\int_{-\infty}^{\infty} \delta(t-a)f(t)dt=f(a). </math>

Para chegar a este resultado utilizou-se o fato de que a <math>\delta(t-a) </math> é nula quando t≠a, então os limites de integração podem ser alterados para a-ε e a+ε. E como <math> f(t)</math> é contínua em t=a seu valor neste intervalo não será muito diferente de <math>f(a)</math>, pode-se dizer que aproximadamente,

<math>\int_{-\infty}^{\infty}\delta(t-a)f(t)dt=\int_{a-\epsilon}^{a+\epsilon}\delta(t-a)f(t)dt\approx f(a)\int_{a-\epsilon}^{a+\epsilon}\delta(t-a)dt=f(a)</math><ref name=":0" />

=== <math></math>Transformada de Laplace da Delta de Dirac ===
Partindo-se da definição da [[Transformada de Laplace]] e utilizando a Propriedade da Filtragem,

<math>\mathcal{L}(\delta(t-a))=\int_{0}^{\infty}\delta(t-a)e^{-st}dt=e^{-as}</math>.

=== Transformada de Fourier da Delta de Dirac ===
Assim como a Transformada de Laplace da função Delta de Dirac sua [[Transformada de Fourier]] também é obtida através da propriedade da filtragem,

<math>\mathcal{F}(\delta(t-a))=\int_{-\infty}^{\infty}\delta(t-a)e^{-i\omega t}dt=e^{-i\omega a}</math>.

== Aplicação em Física ==
Em [[física]], entre outras aplicações, ela é usada para representar [[densidade]]s de objetos pontuais (e.g., carga pontual) <ref name="physics.usu.riffe.3750.15">[[Utah State University]], ''Department of Physics'', ''Phys 3750'', ''Wave Phenomena'', ''Spring 2012'', ''The Dirac Delta Function'' [http://www.physics.usu.edu/riffe/3750/Lecture%2015.pdf <nowiki>[em linha]</nowiki>]</ref> e na normalização de operadores contínuos (e.g.,operador posição) da [[mecânica quântica]]. Além disso, podemos utilizá-la na [[mecânica clássica]] para representar um impulso sendo aplicado em um corpo.

=== Mecânica Estrutural   ===
A equação governante que descreve um sistema massa-mola sob uma ação de uma força (impulso) em t=0 é

<math>m{\operatorname{d^2}\!x\over\operatorname{d}\!t^2} +kx= I \delta (t)</math>

sendo m a massa, k a constante da mola e <math>x</math> a deflexão dela.

Além dessa equação, temos outras na área de mecânica estrutural, como a equação de Euler–Bernoulli, a qual descreve a deflexão de uma viga e é constituída por uma equação diferencial de quarta ordem.

<math>EI {d^4w \over dx^4}= q(x)</math>

onde EI é a rigidez a flexão, pois E é o módulo da elasticidade e I é o momento de inércia da seção transversal de uma viga entorno de um eixo. W é a deflexão e q(x) é a distribuição da carga. 

Se a viga apresenta uma força F em <math>x=x_0</math>, a distribuição da carga é dada pela seguinte equação:

<math>q(x)=F\delta(x-x_0).</math>

Porém ao considerarmos uma viga sob ação de duas forças opostas que distam d metros, elas geram um momento M=Fd na viga e ao aplicarmos o limite com d tendendo a zero temos a seguinte equação
[[Ficheiro:Bending.svg|miniaturadaimagem|Viga deformada devido a uma força uniforme]]
<math>q(x)=\lim_{d \to 0} (F\delta(x)-F\delta(x-d))</math>

<math>q(x)= \lim_{d \to 0} \bigg [ \left ( \frac{M}{d} \right )\delta(x) - \left ( \frac{M}{d} \right )\delta(x-d) \bigg]</math>

<math>q(x)= M\lim_{d \to 0} \left ( \frac{\delta(x)-\delta(x-d)}{d} \right )</math>

Observação: Consideramos na equação acima x=0

Essa equação além de ser aplicada no ramo da física, é também frequentemente utilizada na Engenharia Civil.<ref name=":4" />

== Aplicação em Engenharias ==

=== Engenharia Elétrica ===
Em [[Circuito elétrico|circuitos elétricos]] RLC, RL e RC, por exemplo, quando aplicamos uma fonte no circuito vemos uma saída diferente da entrada, essa saída atípica é composta pela resposta do circuito a excitação da fonte externa que depende das características da função de entrada e das características dos componentes do circuito como a capacitância, indutância e resistência. Como capacitores e indutores envolvem derivadas na relação corrente e tensão, teremos uma equação diferencial para o circuito. Esta equação pode ser resolvida usando a equação característica do circuito.

A fonte excitante desse circuito pode ser modelada por uma '''Delta de Dirac''', e o circuito responderá a esta excitação de uma forma muito particular como exemplificado adiante. Essa variação instantânea de tensão ou corrente pode ser pensada como uma chave abrindo ou fechando rapidamente, por exemplo, gerando um impulso de modelação instantânea. Este tempo de atuação é muito pequeno o que nos possibilita modelar o circuito usando as definições e propriedades da '''Delta de Dirac'''.

Exemplo: Circuito RLC em série sofrendo um chaveamento instantâneo. Seja <math>L</math> a indutância, <math>C</math> a capacitância, <math>R</math> a resistência e <math>\delta (0)</math>  delta de Dirac aplicada em t=0.

<math>L\frac{di}{dt}+Ri+\frac{q}{C}=\delta (0)</math> (Equação modeladora do circuito)

Deixando a equação em função da carga:

<math>Lq''+Rq'+\frac{q}{C}=\delta(0)</math>

Aplicando a Transformada de Laplace e considerando a carga e a corrente zero em t=0:

<math>LS^2Q(S)+SRQ(S)+\frac{Q(S)}{C}=1</math> (Equação no espaço das transformadas)

Observação: Interessante notar que caso a delta estivesse acompanhada de outra fonte, poderíamos usar a propriedade da filtragem listada acima quando calculada a transformada.

Para <math>R</math>=20Ω, <math>L</math>=1H e <math>C</math>=0.0001F, temos a tensão no [[capacitor]] modelada pela seguinte equação:

<math>V_c(t)=100.5e^{-10t}\sin(99.5t)</math>

Análise qualitativa da resposta:

Devido ao número de Euler estar elevado a um expoente negativo, quando o tempo tende ao infinito, a amplitude da tensão tende a 0. Esta resposta é uma característica de fenômenos em circuitos modelados pela delta de Dirac, pois o circuito recebe uma carga de energia no instante de aplicação que a longo prazo é drenada pela resistência do circuito que consome corrente, de acordo com o efeito Joule.

=== Engenharia Civil ===
Em [[Engenharia civil|Engenharia Civil]],  a função delta de Dirac é utilizada para modelagem de estruturas como, por exemplo, vigas. Nestes casos, a função delta de Dirac representa forças/momentos pontuais sendo aplicados na estrutura, como pode-se ver um exemplo abaixo:

Modelagem da deflexão em vigas sujeitas a cargas concentradas:

Considerando uma viga elástica horizontal de comprimento L sob a ação de forças verticais, coloca-se o eixo horizontal x com origem no extremo a esquerda da viga, logo, x=L é o outro extremo. Supõe-se que a viga está sujeita à uma carga ''W''(x) que provoca deflexão em cada ponto x <math>\in</math>[0, L]. Para pequenas deflexões pode-se aproximar a curvatura ''k''(x) pela variação instantânea de <math>\theta</math>(x), onde <math>\theta</math>(x) é o ângulo entre o eixo x e a tangente, ou seja,

<math>k(x) = d\theta(x)/dx</math>.

Como

<math>\operatorname{d}\!y(x)/\operatorname{d}\!x = tan(\theta(x))</math>

e para <math>\theta</math>(x) pequeno, tan(<math>\theta</math>(x)) pode ser considerada como <math>\theta</math>(x), portanto temos:

<math>\operatorname{d}\!y(x)/\operatorname{d}\!x = \theta(x)</math>

Derivando a equação anterior e substituindo, obtemos

<math>\operatorname{d^2}\!y(x)/\operatorname{d}\!x^2 = k(x)</math>.

Pela Lei de Hooke, k(x) = M(x)/EI, onde E é o módulo de Young, ''I'' é o momento de inércia da viga e ''M''(x) o momento fletor. Assim, substituindo na equação anterior:

<math>\operatorname{d^2}\!y(x)/\operatorname{d}\!x^2 = M(x)/EI </math>

A variação do momento de inércia ''M''(x) é a força de cisalamento ''V''(x):

<math>{d \over dx}M(x)=V(x) </math>

e a variação da força de cisalamento é a carga:

<math>{d \over dx}V(x)=W(x)</math>.

Logo,

<math>{d^2 \over dx^2}M(x)=W(x)</math>.

Considerando uma viga engastada, ou seja:

<math>y(0)=y'(0)=y(L)=y'(L)=0</math>

Neste exemplo, a carga estando concentrada na posição x=L/3 e tendo intensidade <math>P_0</math>, a expressão pode ser modelada por:

<math>W(x)=P_0\delta(x-L/3)</math>.<ref name=":4">{{citar livro|título=Transformada de Laplace|ultimo=AZEVEDO, SAUTER|primeiro=Fábio, Ezequia|editora=|ano=2016|local=|páginas=|}}</ref>

=== Engenharia Hídrica ===
Em [[Engenharia hídrica|Engenharia Hídrica]], a função delta de Dirac é utilizada na resolução da [[Equação de Helmholtz|Equação do Helmholtz]] para um sistema barragem-[[Lago artificial|albufeira]] <ref>{{Citar livro |nome=José |sobrenome=Simão Antunes do Carmo|título=Modelação em hidráulica fluvial e ambiente |local=Coimbra |editora=Imprensa da Universidade de Coimbra|ano=2004 |página=130|isbn=9728704283}}</ref>, porém sua maior utilização é na modelagem de águas subterrâneas, sendo na modelagem aquíferos confinados, onde a função delta Dirac representa a vazão bombeada por um poço em um ponto específico<ref>{{Citar livro |nome=Andrés|sobrenome=Sahuquillo|título=Modelos de uso conjunto de aguas superciales y subterraneas |local=España|editora=IGME |ano=2010|página=230|isbn=8478408525}}</ref> e também para calcular a fonte ou sorvedouro tridimensional de um aquífero (equação abaixo), que foi deduzida a partir da [[lei de Darcy]] e do princípio da conservação de massa em um volume elementar representativo (REV) de um aquífero:

<math>W=\sum_{i=1} Q_{i}\delta(x-x_{i})\delta(y-y_{i})\delta(z-z_{i})</math><ref>{{Citar periódico|ultimo=Roehrig|primeiro=J.|ultimo2=Wendland|primeiro2=E.|ultimo3=Martins|primeiro3=J. P.|ultimo4=Firmiano|primeiro4=A.|ultimo5=Roehrig|primeiro5=J.|ultimo6=Wendland|primeiro6=E.|ultimo7=Martins|primeiro7=J. P.|ultimo8=Firmiano|primeiro8=A.|data=2017-8|titulo=Geoprocessamento para a Solução Fraca do Transporte de Contaminantes Acoplado ao Fluxo de Água Subterrânea|url=http://www.scielo.br/scielo.php?script=sci_abstract&pid=S2179-84512017000200273&lng=en&nrm=iso&tlng=pt|jornal=TEMA (São Carlos)|lingua=pt|volume=18|numero=2|paginas=273–286|doi=10.5540/tema.2017.018.02.0273|issn=2179-8451}}</ref>

Sendo <math>Q_{i} \left [ \frac{L^3}{T} \right ]</math>a taxa de bombeamento ou injeção do poço  <math>i</math> em <math>(x_{i},y_{i},z_{i})</math>e <math>\delta</math> representa o armazenamento específico com unidade <math>\left [ \frac{1}{L} \right ]</math>.

== Aplicação em Estatística ==
Em [[estatística]], ela permite generalizar as fórmulas para variáveis aleatórias [[variável aleatória discreta|discretas]] e [[variável aleatória contínua|contínuas]], por exemplo:

O [[valor esperado]] de uma variável aleatória contínua é escrito como:
* <math>E[X] = \int x f(x) dx</math>

Por outro lado, o valor esperado de uma variável aleatória discreta é escrito como:
* <math>E[X] = \sum x_i p(x_i)</math>

O uso da '''Delta de Dirac''' permite unificar estas duas fórmulas, definindo-se a [[função densidade de probabilidade]] de uma variável aleatória discreta por:
* <math>f(x) = \sum p(x_i) \delta(x - x_i)</math>

<br />

== Aplicação no Metabolismo de uma Medicação ==
Com o delta de Direc nos permite descrever o metabolismo de uma medicação de acordo com uma taxa <math>(\tau)</math>(taxa que o organismo metaboliza o medicamento), o que é de extrema importância para determinar a dosagem ideal e também o intervalo de tempo entre as doses. A equação que descreve esse acontecimento é:

<math>m'=\left ( \frac{1}{\tau} \right )m(t)=x(t) </math> , para t>0

sendo x(t) a dosagem ao longo do tempo e m(t) a concentração do medicamento, lembrando que as doses entram no sangue de modo instantâneo e que a taxa de metabolização da droga é diretamente proporcional à concentração.

Considerando a concentração administrada instantaneamente com intervalos T, podemos escrever a equação que se refere a entrada x(t):

<math>x(t)= m_0 [\delta(t)+ \delta (t-T) + \delta(t-2T)+...] </math><ref>{{citar livro|título=TRANSFORMADA de LAPLACE em nove aulas|ultimo=Strauch|primeiro=Irene|editora=Departamento de Matemática Pura e Aplicada UFRGS|ano=2006|local=|páginas=5/10|}}</ref>

Como não há substância no organismo no primeiro momento da ingestão do medicamento, consideramos m(0)=0 e aplicamos a transformada de Laplace para calcular m(t):

<math>sM(s) +\left ( \frac{1}{\tau} \right )M(s)= m_0 (1+e^{-sT} + e^{-2sT} +...) = m_0 \sum_{n=0}^{\infin} (e^{-sT})^n   </math><ref name=":1">{{citar livro|título=Transformada de Laplace Livro Colaborativo|ultimo=Sauter,de Azevedo, Strauch|primeiro=Esequia,Fábio, Irene|editora=UFRGS|ano=2018|local=Porto Alegre|páginas=62-63|}}</ref>

Observação: <math>\mathcal{L}\{m(t)\}\ = sM(s) + m(0) </math>, porém como m(0)=0, esse termo não aparece na equação.

Manipulando a equação acima obtemos: 

<math>M(s) = \left ( \frac{m_0}{s+\left ( \frac{1}{\tau} \right )} \right ) \sum_{n=0}^{\infin} (e^{-sT})^n   </math><ref name=":1" />

Para encontrar a função m(t) basta aplicamos a inversa da transformada de Laplace usando a propriedade do deslocamento em s

<math>m(t)= m_0 ( e^\left ( \frac{-t}{\tau} \right )+ e^\left ( -\frac{t-T}{\tau} \right )u(t-T) + e^\left ( -\frac{t-2T}{\tau} \right )u(t-2T) + ...)   </math> <ref name=":1" /><br />
== Pente de Dirac ==
[[Ficheiro:Dirac comb.svg|miniaturadaimagem|Um Pente de Dirac é uma série infinita de funções delta de dirac espaçadas por intervalos de T.]]
Conhecido por ser um
trem de impulso uniforme encontrado por meio do delta de Dirac  ,também conhecido como  função shah, cria uma
função de amostragem, usualmente utilizada em processamento de sinais digitais e para análise de sinais de tempo discretos. O pente de Dirac é dado pela soma infinita, cujo limite pode ser entendido pelo sentido de distribuição,
:<math>\Delta(x) = \sum_{n=-\infty}^\infty \delta(x-n),</math>
que e a sequencia massas pontuais de cada uma das integrais.

Mediante a uma constante de normalizacao global, o pente de Dirac é igual a sua própria transformada de Fourier. Isso é muito significativo, pois se ''f'' for qualquer função de Schawartz, então a periodização de ''f'' será dada pela Convolução
:<math>(f*\Delta)(x) = \sum_{n=-\infty}^\infty f(x-n).</math>
Em particular,
:<math>(f*\Delta)^\wedge = \hat{f}\widehat{\Delta} = \hat{f}\Delta</math>
E precisamente como a “Poisson summation formula”. <ref>{{harvnb|Córdoba|1988}}; {{harvnb|Hörmander|1983|loc=§7.2}}</ref>

== Integral ==
Em certo sentido, pode-se dizer que a '''delta de Dirac''' é a derivada da [[função de passo Heaviside]], ou que a integral da '''delta de Dirac''' é a função de passo Heaviside:

:<math> H(x) = \int_{-\infty}^x { \delta(t)} \mathrm{d}t </math>

== Ver também ==

* [[Pente de Dirac]]
* [[Função de Heaviside]]

== Referências ==
<references />
{{Controlo de autoridade}}

{{DEFAULTSORT:Delta Dirac}}
[[Categoria:Funções generalizadas|Delta]]
[[Categoria:Paul Dirac]]