Seno cardinal

Sinc(x) normalizada (azul) frente a la sinc desnormalizada (rojo) con la misma escala: x = −6π a 6π.

En matemática, el seno cardinal es una función especial denotada por $\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{c}(x)$; tiene dos definiciones, la «normalizada» y la «desnormalizada», que se definen de la siguiente forma:

1. En procesamiento digital de señales y teoría de la información, la función sinc normalizada comúnmente se define como:

   ${\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{c}}_N(x)=\frac{\mathrm{sen}(\pi x)}{\pi x}$

2. En matemática, la histórica función sinc desnormalizada, está definida por:

   $\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{c}(x)=\frac{\mathrm{sen}(x)}{x}$

En ambos casos el valor de la función tiene una singularidad evitable en cero, que generalmente se redefine específicamente como igual a 1. El seno cardinal es analítico en todo el dominio de los números reales, excepto para el valor $x=0$

La función «desnormalizada» es idéntica a la «normalizada» excepto por el factor de escala faltante en el argumento. La función sinc corresponde a la transformada de Fourier de un pulso rectangular, y la transformada inversa de Fourier de un espectro rectangular es una sinc.

Reseña histórica

Sinc es una contracción del nombre latino completo de la función sinus cardinalis («seno cardinal»). La notación fue introducida por el matemático e ingeniero británico Philip Woodward en su artículo de 1952 Information theory and inverse probability in telecommunication (Teoría de la Información y probabilidad inversa en las telecomunicaciones), en el que afirmó que la función «se reproduce con tanta frecuencia en el análisis de Fourier y en sus aplicaciones, que parece merecer alguna notación propia». Esta notación también aparece utilizada en su libro de 1953 "Probability and Information Theory, with Applications to Radar" (Probabilidad y Teoría de la Información, con aplicaciones al Radar).

Con anterioridad, ya en 1915, el también matemático británico Edmund Whittaker (1873 - 1956) había utilizado esta función aplicada a procesos de muestreo, aunque no le dio nombre ni notación específicos. La función permite resolver el problema de:

Determinar una función que pasa por los puntos $(a+kw,f(a+kw))$, donde $k$ es un número entero, y $w$ es un número complejo, obteniendo una interpolación tan suave como sea posible, sin singularidades y con rápidas oscilaciones de los valores tabulares dados de $f(x)$.

y toma la forma:

$C(x)=\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}f(a+kw)\frac{\mathrm{sen}(\frac{\pi }{w})(x-a-kw)}{(\frac{\pi }{w})(x-a-kw)}$

Propiedades

Los máximos y mínimos locales (pequeños puntos blancos) de la función sinc sin normalizar (curva de color rojo), coinciden con sus intersecciones con la función coseno (curva de color azul).

Parte real de sinc compleja Plantilla:Math.

Parte imaginaria de sinc compleja Plantilla:Math.

Valor absoluto
Plantilla:Math.

Los ceros (cortes con el eje horizontal) de la función sinc sin normalizar están en múltiplos enteros de π no nulos, mientras que los ceros de la función sinc normalizada se localizan en números enteros distintos de cero.

Los máximos y mínimos locales de la función sinc sin normalizar se corresponden con sus intersecciones con la función coseno. Es decir, Plantilla:Math para todos los puntos Plantilla:Mvar donde la derivada de Plantilla:Math es cero y por lo tanto se alcanza un valor extremo local.

Una buena aproximación de la coordenada Plantilla:Mvar del valor extremo enésimo Plantilla:Math, con Plantilla:Math positivo, es la coordenada

$x_n\approx (n+\frac{1}{2})\pi -\frac{1}{(n+\frac{1}{2})\pi },$

donde Plantilla:Math impar lleva a un mínimo local y Plantilla:Math par a un máximo local. Además de valores extremos en Plantilla:Math, la curva tiene un máximo absoluto en Plantilla:Math y debido a su simetría respecto al eje Plantilla:Math, también lo es para los valores de Plantilla:Math negativos Plantilla:Math.

La función sinc normalizada posee una representación simple como el productorio

${\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{c}}_N(x)=\frac{\mathrm{sen}(\pi x)}{\pi x}=\prod _{n=1}^{\mathrm{\infty }}(1-\frac{x^2}{n^2})$

y se relaciona con la función gamma Plantilla:Math través fórmula de reflexión de Euler,

${\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{c}}_N(x)=\frac{\mathrm{sen}(\pi x)}{\pi x}=\frac{1}{\mathrm{\Gamma }(1+x)\mathrm{\Gamma }(1-x)}.$

Euler descubrió que

$\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{c}(x)=\frac{\mathrm{sen}(x)}{x}=\prod _{n=1}^{\mathrm{\infty }}\cos \left(\frac{x}{2^n}\right).$

La transformada de Fourier de la función sinc normalizada (a la frecuencia ordinaria) es Plantilla:Math,

$\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{\mathrm{sinc}}_N(t)e^{i2\pi ft}dt=\mathrm{rect}(f),$

donde la función rectangular es 1 para argumentos entre -Plantilla:Sfrac y Plantilla:Sfrac, y cero en caso contrario. Esto corresponde al hecho de que el filtro sinc es el filtro de paso bajo ideal (lo que se denomina en inglés un brick-wall, es decir, un filtro electrónico idealizado, que presenta plena transmisión en la banda de paso, y atenuación completa en la banda restringida, con transiciones bruscas, que se conoce coloquialmente en su traducción literal como "filtro de muro de ladrillo").

Esta integral de Fourier, incluyendo el caso especial

$\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\frac{\mathrm{sen}(\pi x)}{\pi x}dx=\mathrm{rect}(0)=1$

es una integral impropia (según la definición de integral de Dirichlet), y no una integral de Lebesgue convergente, como

$\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\left|\frac{\mathrm{sen}(\pi x)}{\pi x}\right|dx=+\mathrm{\infty }.$

De la anterior integral de Fourier se deducen las expresiones siguientes:

$\underset{-\mathrm{\infty }}{\overset{+\mathrm{\infty }}{\int }}{\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{c}}_N(x)dx=1$

$\underset{-\mathrm{\infty }}{\overset{+\mathrm{\infty }}{\int }}\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{c}(x)dx=\pi$

La función sinc normalizada tiene propiedades que la hacen ideal en relación con la interpolación de funciones muestreadas con ancho de banda limitado:

• Se trata de una función de interpolación, es decir, Plantilla:Math, y Plantilla:Math para números enteros Plantilla:Math distintos de cero.
• Las funciones Plantilla:Math con (Plantilla:Math entero) forman una base ortonormal para funciones de ancho de banda limitado en el espacio funcional Plantilla:Math, con la frecuencia angular más alta Plantilla:Math (es decir, la frecuencia de ciclo más alta Plantilla:Math).

Otras propiedades de las dos funciones sinc incluyen:

• La función sinc no normalizada es la función de Bessel esférica de primera clase de orden cero, Plantilla:Math. La función sinc normalizada es Plantilla:Math.
• $\underset{0}{\overset{x}{\int }}\frac{\mathrm{sen}(\theta )}{\theta }d\theta =\mathrm{Si}(x)$

donde Plantilla:Math es la función integral senoidal.

• Plantilla:Math (no normalizada) es una de las dos soluciones linealmente independientes de la ecuación diferencial ordinaria lineal

$x\frac{d^{2}y}{d{x}^2}+2\frac{dy}{dx}+{\lambda }^{2}xy=0.$

La otra es Plantilla:Math, que no está limitada en Plantilla:Math, a diferencia de su contraparte, la función sinc.

• $\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\frac{{\mathrm{sen}}^2(\theta )}{{\theta }^2}d\theta =\pi \to \underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{\mathrm{sinc}}^2(x)dx=1,$

donde se hace referencia a la función sinc normalizada.

• $\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\frac{{\mathrm{sen}}^3(\theta )}{{\theta }^3}d\theta =\frac{3\pi }{4}$
• $\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\frac{{\mathrm{sen}}^4(\theta )}{{\theta }^4}d\theta =\frac{2\pi }{3}.$

Relación con la distribución delta de Dirac

La función sinc normalizada se puede utilizar como una función delta naciente, de acuerdo con la siguiente convergencia débil,

$\underset{a\to 0}{\lim }\frac{\mathrm{sen}\left(\frac{\pi x}{a}\right)}{\pi x}=\underset{a\to 0}{\lim }\frac{1}{a}{\text{sinc}}_N\left(\frac{x}{a}\right)=\delta (x).$

Este no es un límite ordinario, puesto que desde el lado izquierdo no converge. Más bien, significa que

$\underset{a\to 0}{\lim }\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\frac{1}{a}{\text{sinc}}_N\left(\frac{x}{a}\right)\phi (x)dx=\phi (0),$

para cualquier función continuamente diferenciable Plantilla:Math con soporte compacto.

En la expresión anterior, cuando Plantilla:Math, el número de oscilaciones por unidad de longitud del seno cardinal se aproxima a infinito. Sin embargo, la expresión siempre oscila dentro del intervalo Plantilla:Math, independientemente del valor de Plantilla:Mvar.

Esto complica la imagen informal de Plantilla:Math como cero para todos los valores de Plantilla:Mvar excepto para el punto Plantilla:Math, e ilustra el problema de la idea de la función delta como una función más que como una distribución. Una situación similar se aparece en el fenómeno de Gibbs.

Series

Todos los sumatorios en esta sección se refieren a la función sinc sin normalizar.

El sumatorio de Plantilla:Math sobre Plantilla:Mvar para los números enteros de 1 a Plantilla:Math es igual a Plantilla:Math.

$\sum _{n=1}^{\mathrm{\infty }}\mathrm{sinc}(n)=\mathrm{sinc}(1)+\mathrm{sinc}(2)+\mathrm{sinc}(3)+\mathrm{sinc}(4)+\cdots =\frac{\pi -1}{2}$

El sumatorio de los cuadrados es también igual a Plantilla:Math.

$\sum _{n=1}^{\mathrm{\infty }}{\mathrm{sinc}}^2(n)={\mathrm{sinc}}^2(1)+{\mathrm{sinc}}^2(2)+{\mathrm{sinc}}^2(3)+{\mathrm{sinc}}^2(4)+\cdots =\frac{\pi -1}{2}$

Cuando los signos de las adiciones se alternan y comienzan con el signo +, el sumatorio es igual a Plantilla:Sfrac.

$\sum _{n=1}^{\mathrm{\infty }}(-1{)}^{n+1}\mathrm{sinc}(n)=\mathrm{sinc}(1)-\mathrm{sinc}(2)+\mathrm{sinc}(3)-\mathrm{sinc}(4)+\cdots =\frac{1}{2}$

El sumatorio de los términos alternantes de los cuadrados y los cubos también es igual a Plantilla:Sfrac.

$\sum _{n=1}^{\mathrm{\infty }}(-1{)}^{n+1}{\mathrm{sinc}}^2(n)={\mathrm{sinc}}^2(1)-{\mathrm{sinc}}^2(2)+{\mathrm{sinc}}^2(3)-{\mathrm{sinc}}^2(4)+\cdots =\frac{1}{2}$

$\sum _{n=1}^{\mathrm{\infty }}(-1{)}^{n+1}{\mathrm{sinc}}^3(n)={\mathrm{sinc}}^3(1)-{\mathrm{sinc}}^3(2)+{\mathrm{sinc}}^3(3)-{\mathrm{sinc}}^3(4)+\cdots =\frac{1}{2}$

Expansión de las series

Para Plantilla:Math sin normalizar:

$\mathrm{sinc}(x)=\frac{\mathrm{sen}(x)}{x}=\sum _{n=0}^{\mathrm{\infty }}\frac{{(-x^2)}^n}{(2n+1)!}$

Dimensiones superiores

El producto de funciones sinc 1-D proporciona fácilmente una función multivariable sinc para un sistema de referencia cartesiano de rejilla cuadrada (gráfico de celosía): Plantilla:Math cuya transformada de Fourier es la función indicatriz de un cuadrado en el espacio de frecuencia (es decir, el definido por un muro de ladrillo en el espacio 2-D). El seno cardinal para un gráfico de celosía no cartesiano (por ejemplo, una rejilla hexagonal) es una función cuya transformada de Fourier es la función indicatriz de la zona de Brillouin de la celosía. Por ejemplo, la función sinc para la red hexagonal es una función cuya transformada de Fourier es la indicatriz de la unidad hexágonal en el espacio de la frecuencia. Para una retícula cartesiana esta función no se puede obtener por un simple producto tensorial. Sin embargo, la fórmula explícita de la función sinc para rejillas hexagonales, sistemas cúbicos, celosías cúbicas centradas en las caras y otras celosías de dimensiones superiores pueden ser deducidas explícitamente utilizando las propiedades geométricas de las zonas de Brillouin y su conexión con los zonotopos.

Por ejemplo, una rejilla hexagonal puede ser generada por el sistema generador de vectores (enteros)

$u_1=\left[\begin{array}{c}\frac{1}{2}\\ \frac{\sqrt{3}}{2}\end{array}\right]u_2=\left[\begin{array}{c}\frac{1}{2}\\ -\frac{\sqrt{3}}{2}\end{array}\right].$

siendo

${\xi }_1=\frac{2}{3}{u}_1,{\xi }_2=\frac{2}{3}{u}_2,{\xi }_3=-\frac{2}{3}(u_1+u_2),𝐱=\left[\begin{array}{c}x\\ y\end{array}\right],$

puede deducirse el seno cardinal para esta red hexagonal como:

$\begin{array}{cc}{\mathrm{sinc}}_{\mathrm{H}}(𝐱)=\frac{1}{3}(& \cos (\pi {\xi }_1\cdot 𝐱)\mathrm{sinc}({\xi }_2\cdot 𝐱)\mathrm{sinc}({\xi }_3\cdot 𝐱)+\\ & \cos (\pi {\xi }_2\cdot 𝐱)\mathrm{sinc}({\xi }_3\cdot 𝐱)\mathrm{sinc}({\xi }_1\cdot 𝐱)+\\ & \cos (\pi {\xi }_3\cdot 𝐱)\mathrm{sinc}({\xi }_1\cdot 𝐱)\mathrm{sinc}({\xi }_2\cdot 𝐱))\end{array}$

Esta construcción se puede utilizar para diseñar ventanas de Lanczos para rejillas multidimensionales generales.

Véase también

• Filtro antialiasing
• Filtro Sinc
• Remuestreo de Lanczos
• Fórmula de interpolación de Whittaker–Shannon
• Proyección de Winkel-Tripel (cartografía)
• Integral senoidal
• Funciones trigonométricas de matrices
• Integral de Borwein
• Integral de Dirichlet

Enlaces externos

• Plantilla:MathWorld