Articulo de referencia

Ondícula de Cauchy

En matemáticas, las ondículas de Cauchy son una familia de ondículas continuas, utilizadas en la transformada continua de ondículas . Definición La ondícula de Cauchy de orden p...

En matemáticas, las ondículas de Cauchy son una familia de ondículas continuas, utilizadas en la transformada continua de ondículas .

Definición

La ondícula de Cauchy de ordenpag{\displaystyle p}se define como:

ψpag(t)=Γ(pag+1)2π(jt+j)pag+1{\displaystyle \psi _{p}(t)={\frac {\Gamma (p+1)}{2\pi }}\left({\frac {j}{t+j}}\right)^{p+1}}

dóndepag>0{\displaystyle p>0}yj=1{\displaystyle j={\sqrt {-1}}} Por lo tanto, su transformada de Fourier se define como

ψpag^(ξ)=ξpagmiξI[ξ0]{\displaystyle {\hat {\psi _{p}}}(\xi )=\xi ^{p}e^{-\xi }I_{[\xi \geq 0]}}.

A veces se define como una función con su transformada de Fourier [ 1 ].

ψpag^(ξ)=ρ(ξ)ξpagmiξI[ξ0]{\displaystyle {\hat {\psi _{p}}}(\xi )=\rho (\xi )\xi ^{p}e^{-\xi }I_{[\xi \geq 0]}}

dóndeρ(ξ)L(R){\displaystyle \rho (\xi )\in L^{\infty }(\mathbb {R} )}yρ(ξ)=ρ(aξ){\displaystyle \rho (\xi )=\rho (a\xi )}paraξR{\displaystyle \xi \in \mathbb {R} }casi en todas partes yρ(ξ)0{\displaystyle \rho (\xi )\neq 0}a pesar deξR{\displaystyle \xi \in \mathbb {R} }.

Además, solía definirse como [ 2 ]

ψpag(t)=(jt+j)pag+1{\displaystyle \psi _{p}(t)=({\frac {j}{t+j}})^{p+1}}

en investigaciones previas de la ondícula de Cauchy. Si definimos la ondícula de Cauchy de esta manera, podemos observar que la transformada de Fourier de la ondícula de Cauchy

ψpag^(ξ)dξ=02πΓ(pag+1)ξpagmiξdξ=2π{\displaystyle \int _{-\infty }^{\infty }{\hat {\psi _{p}}}(\xi )\,d\xi =\int _{0}^{\infty }{\frac {2\pi }{\Gamma (p+1)}}\xi ^{p}e^{-\xi }\,d\xi =2\pi }

Además, podemos ver que el máximo de la transformada de Fourier de la ondícula de Cauchy de ordenpag{\displaystyle p}sucedió enξ=pag{\displaystyle \xi =p}y la transformada de Fourier de la ondícula de Cauchy es positiva solo enξ>0{\displaystyle \xi >0}, significa que: (1) cuandopag{\displaystyle p}es bajo entonces la convolución de la ondícula de Cauchy es un filtro de paso bajo, y cuandopag{\displaystyle p}es alta la convolución de la ondícula de Cauchy es un filtro de paso alto. Dado que la transformada de ondícula es igual a la convolución a la ondícula madre y la convolución a la ondícula madre es igual a la multiplicación entre la transformada de Fourier de la ondícula madre y la función por el teorema de convolución . Y, (2) el diseño de la transformada de ondícula de Cauchy se considera con el análisis de la señal analítica.

Dado que la señal analítica es biyectiva a la señal real y solo hay frecuencia positiva en la señal analítica (la señal real tiene frecuencia conjugada entre positiva y negativa), es decir,

FT{incógnita}(ξ)¯=FT{incógnita}(ξ){\displaystyle {\overline {FT\{x\}(-\xi )}}=FT\{x\}(\xi )}

dóndeincógnita(t){\displaystyle x(t)}es una señal real (incógnita(t)R{\displaystyle x(t)\in \mathbb {R} }, para todostR{\displaystyle t\in \mathbb {R} }) Y la biyección entre la señal analítica y la señal real es que

incógnita+(t)=incógnita(t)+jincógnitaH(t){\displaystyle x_{+}(t)=x(t)+jx_{H}(t)}

incógnita(t)=Rmi{incógnita+(t)}{\displaystyle x(t)=Re\{x_{+}(t)\}}

dóndeincógnita+(t){\displaystyle x_{+}(t)}es la señal analítica correspondiente de la señal realincógnita(t){\displaystyle x(t)}, yincógnitaH(t){\displaystyle x_{H}(t)}es la transformada de Hilbert deincógnita(t){\displaystyle x(t)}.

Unicidad de la reconstrucción

Problema de recuperación de fase

Un problema de recuperación de fase consiste en reconstruir una función compleja desconocida.F{\displaystyle f}a partir de un conjunto de mediciones lineales sin fase. Más precisamente, dejemosV{\displaystyle V}Sea un espacio vectorial , cuyos vectores son funciones complejas, endo{\displaystyle \mathbb {C} }y {Li}iI{\displaystyle \{L_{i}\}_{i\in I}}un conjunto de formas lineales deV{\displaystyle V}ado{\displaystyle \mathbb {C} }. Se nos da el conjunto de todos{|Li(F)|}iI{\displaystyle \{|L_{i}(f)|\}_{i\in I}}, por algún desconocidoFV{\displaystyle f\in V}y queremos determinarF{\displaystyle f}. Este problema puede estudiarse desde tres puntos de vista diferentes: [ 1 ] (1) EsF{\displaystyle f}determinado de forma única por{|Li(F)|}iI{\displaystyle \{|L_{i}(f)|\}_{i\in I}}(hasta una fase global)? (2) Si la respuesta a la pregunta anterior es positiva, ¿es la aplicación inversa?{|Li(F)|}iIF{\displaystyle \{|L_{i}(f)|\}_{i\in I}\implies f}¿Es “estable”? Por ejemplo, ¿es continuo ? ¿ Uniformemente Lipschitz ? (3) En la práctica, ¿existe un algoritmo eficiente que recupereF{\displaystyle f}de{|Li(F)|}iI{\displaystyle \{|L_{i}(f)|\}_{i\in I}}¿

El ejemplo más conocido de un problema de recuperación de fase es el caso en el queLi{\displaystyle L_{i}}representar los coeficientes de Fourier: por ejemplo:

Lnorte(F)=12πππF(t)mijnortetdt{\displaystyle L_{n}(f)={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(t)e^{-jnt}\,dt}, paranorteZ{\displaystyle n\in \mathbb {Z} },

dóndeF{\displaystyle f}es una función de valor complejo en[π,π]{\displaystyle [-\pi ,\pi ]} Entonces,F{\displaystyle f}puede reconstruirse porLnorte(F){\displaystyle L_{n}(f)}como

F(t)=norte=Lnorte(F)mijnortet{\displaystyle f(t)=\sum _{n=-\infty }^{\infty }L_{n}(f)e^{jnt}}.

y de hecho tenemos la identidad de Parseval

||F||2=norte=|Lnorte(F)|2{\displaystyle ||f||^{2}=\sum _{n=-\infty }^{\infty }|L_{n}(f)|^{2}}.

dónde||F||2=12πππ|F(t)|2dt{\displaystyle ||f||^{2}={\frac {1}{2\pi }}\int _{-\pi }^{\pi }|f(t)|^{2}\,dt}es decir, la norma definida enL2([π,π]){\displaystyle L^{2}([-\pi ,\pi ])}Por lo tanto , en este ejemplo, el conjunto de índicesI{\displaystyle I}es el número enteroZ{\displaystyle \mathbb {Z} }, el espacio vectorialV{\displaystyle V}esL2([π,π]){\displaystyle L^{2}([-\pi ,\pi ])}y la forma linealLnorte{\displaystyle L_{n}}es el coeficiente de Fourier . Además, el valor absoluto de los coeficientes de Fourier{|Lnorte(F)|}norteZ{\displaystyle \{|L_{n}(f)|\}_{n\in \mathbb {Z} }}solo puede determinar la norma deF{\displaystyle f}definido enL2([π,π]){\displaystyle L^{2}([-\pi ,\pi ])}.

Teorema de unicidad de la reconstrucción

En primer lugar, definimos la transformada wavelet de Cauchy como:

Wψpag[incógnita(t)](a,b)=1bincógnita(t)ψpag(tab)¯dt{\displaystyle W_{\psi _{p}}[x(t)](a,b)={\frac {1}{b}}\int _{-\infty }^{\infty }x(t){\overline {\psi _{p}({\frac {t-a}{b}})}}\,dt}.

Entonces, el teorema es: Teorema. [ 1 ] Para un fijopag>0{\displaystyle p>0}, si existen dos números diferentesb1,b2>0{\displaystyle b_{1},b_{2}>0}y la transformada wavelet de Cauchy definida como se indicó anteriormente. Entonces, si hay dos funciones de valor realF,gramoL2(R){\displaystyle f,g\in L^{2}(\mathbb {R} )}satisfecho

|Wψpag[F(t)](a,b1)|=|Wψpag[gramo(t)](a,b1)|{\displaystyle |W_{\psi _{p}}[f(t)](a,b_{1})|=|W_{\psi _{p}}[g(t)](a,b_{1})|},aR{\displaystyle \forall a\in \mathbb {R} }y

|Wψpag[F(t)](a,b2)|=|Wψpag[gramo(t)](a,b2)|{\displaystyle |W_{\psi _{p}}[f(t)](a,b_{2})|=|W_{\psi _{p}}[g(t)](a,b_{2})|},aR{\displaystyle \forall a\in \mathbb {R} },

entonces hay unαR{\displaystyle \alpha \in \mathbb {R} }de tal manera queF+(t)=mijαgramo+(t){\displaystyle f_{+}(t)=e^{j\alpha }g_{+}(t)}. F+(t)=mijαgramo+(t){\displaystyle f_{+}(t)=e^{j\alpha }g_{+}(t)}implica que

Rmi{F+(t)}=Rmi{mijαgramo+(t)}F(t)=porqueαgramo(t)pecadoαgramoH(t){\displaystyle Re\{f_{+}(t)\}=Re\{e^{j\alpha }g_{+}(t)\}\implies f(t)=\cos {\alpha }g(t)-\sin {\alpha }g_{H}(t)}y

Imetro{F+(t)}=Imetro{mijαgramo+(t)}FH(t)=pecadoαgramo(t)+porqueαgramoH(t){\displaystyle Im\{f_{+}(t)\}=Im\{e^{j\alpha }g_{+}(t)\}\implies f_{H}(t)=\sin {\alpha }g(t)+\cos {\alpha }g_{H}(t)}.

Por lo tanto, obtenemos la relación

F(t)=(porqueαpecadoαbroncearseα)gramo(t)broncearseαFH(t){\displaystyle f(t)=(\cos {\alpha }-\sin {\alpha }\tan {\alpha })g(t)-\tan {\alpha }f_{H}(t)}

yF(t),gramoH(t)spaganorte{FH(t),gramo(t)}=spaganorte{F(t),FH(t)}=spaganorte{gramo(t),gramoH(t)}{\displaystyle f(t),g_{H}(t)\in span\{f_{H}(t),g(t)\}=span\{f(t),f_{H}(t)\}=span\{g(t),g_{H}(t)\}}. Volviendo al problema de recuperación de fase, en el caso de la transformada wavelet de Cauchy, el conjunto de índicesI{\displaystyle I}esR×{b1,b2}{\displaystyle \mathbb {R} \times \{b_{1},b_{2}\}}conb1b2{\displaystyle b_{1}\neq b_{2}}yb1,b2>0{\displaystyle b_{1},b_{2}>0}, el espacio vectorialV{\displaystyle V}esL2(R){\displaystyle L^{2}(\mathbb {R} )}y la forma linealL(a,b){\displaystyle L_{(a,b)}}se define comoL(a,b)(F)=Wψpag[F(t)](a,b){\displaystyle L_{(a,b)}(f)=W_{\psi _{p}}[f(t)](a,b)}. Por eso,{|L(a,b)(F)|}a,bR×{b1,b2}{\displaystyle \{|L_{(a,b)}(f)|\}_{a,b\in \mathbb {R} \times \{b_{1},b_{2}\}}}determina el subespacio bidimensionalspaganorte{F,FH}{\displaystyle span\{f,f_{H}\}}enL2(R){\displaystyle L^{2}(\mathbb {R} )}.

Referencias

  1. 1 2 3 Mallat, Stéphane; Waldspurger, Irène (2015). "Recuperación de fase para la transformada wavelet de Cauchy" . Journal of Fourier Analysis and Applications . 21 (6): 1251– 1309. arXiv : 1404.1183 . Bibcode : 2015JFAA...21.1251M . doi : 10.1007/s00041-015-9403-4 .
  2. Argoul, Pierre; Le, Thien-phu (2003). "Indicadores instantáneos del comportamiento estructural basados ​​en el análisis continuo de ondículas de Cauchy" . Mechanical Systems and Signal Processing . 17 (1): 243– 250. Bibcode : 2003MSSP...17..243A . doi : 10.1006/mssp.2002.1557 .