Los valores de la señal que 'perderemos' entre las instancias de muestreo durante el muestreo de señales de banda limitada

9

Según el teorema de muestreo de Nyquist-Shannon, cualquier señal de tiempo continuo con un ancho de banda menor que la frecuencia de Nyquist (con la frecuencia de muestreo), que se muestrea a la frecuencia de muestreo puede reconstruirse perfectamente mediante interpolación sinc (es decir la fórmula de interpolación de Whittaker-Shannon).BfN=fs/2fsfs

Supongamos que muestreamos una señal de tiempo continua desconocida, limitada en magnitud, con tiempo de muestreo constante en instancias de muestra , ( ), sin jitter de muestreo o cuantización. la restricción de que , con .T=1/fskTkZB=αfN0α1

Lo que me gustaría averiguar es el siguiente: En muestra del instante , quiero determinar para cada peor de los casos fraccional 'exceso' de cualquier señal de tiempo continuo entre las muestras y , que podría haber tenido. Es decir, cuánto fue mayor la señal de tiempo continuo que los valores muestreados (absolutos) más altos en los instantes de muestra y . La señal continua, o reconstrucción (¡ya que la interpolación sinc es perfecta!), Que hemos 'perdido' por muestreo.kαk1kkk1

Ejemplo: Configuramos y asumimos una señal de tiempo discreta [1,0,1,0,1,1,0,1,0,1] (observe el doble 1 cerca del medio, y esta señal incluso tiene ?). Su reconstrucción sinc (línea azul) de las muestras (impulsos negro) tiene el siguiente aspecto (I han trazado los sincs pertenecientes a cada muestra en gris): El 'sobreimpulso' entre las muestras y , se encuentra oα=1α=1Reconstrucción de Sinc para secuenciak=0k=10.770%. Así que perdimos un pico de valor 1.7 en nuestro tiempo continuo limitado de banda original, o la señal 'reconstruido perfectamente limitado de banda'. Si hubiera puesto 3 o más 1 consecutivos, el sobreimpulso habría sido menor (al final, el fenómeno de Gibbs es mucho más pequeño). Por lo tanto, 2 muestras continuas consecutivas como esta es el "peor de los casos".

Extender la señal en ambas direcciones hará que crezca el sobreimpulso: ingrese la descripción de la imagen aquí que muestra un sobreimpulso relativo de a un valor de casi 2.1.1.1

Para cualquier longitud de secuencia de , este 'sobreimpulso' crecerá indefinidamente, , que va a cuando . Esto se debe a que cada muestra de los sincs creará una 'interferencia' constructiva, y la suma de (las contribuciones de todas las envolventes de unidades sinc) para no converge.2mo(m)o(m)ln(m)m1/πnn

Esto (creo) similar a lo siguiente: si muestrea constantemente un valor 0, también podría reconstruir una señal de tiempo continuo con amplitud infinita que solo se muestrea en los nodos a valores de 0, por ejemplo, . Esto me dice lo mismo: que si permito que una señal esté en la frecuencia de Nyquist, el peor sobreimpulso que podría 'perder' es infinito.sinπfst

Ahora podemos afirmar que . Y podemos razonar que (el muestreo de una señal constante de la que sabe que está limitada en la banda tiene una reconstrucción constante única).o(m)|α=1=o(m)|α=0=0

¿Qué pasa si ?α<1

Si ahora asumimos que hacemos esta misma interpolación sinc, pero sabemos con certeza α<1, me gusta α=0.5. Entonces, (dice mi instinto) este efecto debería disminuir e incluso debería permanecer finito (cuandom) !. Dado que para cualquier señal de pared de ladrillo limitada al ancho de bandaαfs/2, obtenemos una respuesta de impulso de filtro de h(t)sinc(tkTαT)(¿Correcto?). Por lo tanto, las transiciones de señal no pueden ser tan rápidas como para el ejemplo del tren de impulsos cambiante anterior y, por lo tanto, las contribuciones de cada función sinc durante la reconstrucción no pueden crear interferencia constructiva infinita.

Mi problema: no sé cómo proceder desde aquí; cómo formar una 'prueba' del exceso de caso más desfavorable que podría haber encontrado entre 2 muestras consecutivas, sabiendo queα<1, para anyseñal (no necesariamente estas unidades de impulso-tren como ejemplos). Un valor dado paraα me da una pendiente dh(t)dt del núcleo de convolución limitante de banda h(t), lo que debería decirme algo sobre cuántas muestras consecutivas deben ser diferentes, pero no veo los pasos a seguir para llegar a una conclusión genérica.

Retinita
fuente
Discutimos tales secuencias patológicas en comp.dsp en 2002, tema: A Sampled Data Interpolation Poser, groups.google.com/d/msg/comp.dsp/EQ31d-2SS2o/wT5HXbjQpogJ y 2003, tema: Worst Case Signal for Reconstruction, groups.google.com/d/msg/comp.dsp/xwb9p3awrOg/zl20Wl2EiesJ
Olli Niemitalo el
Creo que hay un teorema que relaciona el ancho de banda de una función con un límite superior de su densidad promedio de cruces por cero. Ahora, para una función infinita en casi todas partes, esas muestras de valores finitos son tal vez como los cruces por cero son para una función de valores finitos: su densidad promedio tiene un límite superior.
Olli Niemitalo
Gracias, leeré las discusiones grupales de Google en detalle cuando tenga más tiempo (¿a dónde fue su respuesta con la figura?). Aún así, la respuesta de MBaz parece sugerir que existe una derivada absoluta máxima en el peor de los casos, que, sisup|x|es finito, será finito. Por lo tanto, para cualquier señal de banda limitada, no puede ir a un valor infinito. ¿Cómo se relaciona eso con lo que estás sugiriendo?
Retinite
Eliminé mi respuesta porque no tuve en cuenta que la secuencia discreta debe ser tal que sobreviva al filtrado de paso bajo intacto. Por lo tanto, es posible que tenga razón sobre lo que sucede en , y mi comentario anterior estaría de acuerdo con eso. podría encontrarse entre muestras, por lo que no dice mucho. α<1sup|x|
Olli Niemitalo el
Me pregunto si las matemáticas serían más simples para secuencias periódicas de longitud infinita con el período
Olli Niemitalo

Respuestas:

7

No tengo una respuesta real, pero tengo la sensación de que este resultado te ayudará: la desigualdad de Bernstein dice que, si la señal está limitada en banda a , entonces donde significa "límite superior mínimo".x(t)|f|B

|dx(t)dt|4πBsupτR|x(τ)|,tR
sup

Encontré esta desigualdad en el excelente libro de Amos Lapidoth (y gratuito en formato PDF) "A Foundation in Digital Communication". Se puede encontrar una prueba en MA Pinsky, "Introducción al análisis de Fourier y Wavelets".

MBaz
fuente
¡Gracias! Eso es realmente útil; razonamiento a partir de la señal horaria continua. Esto debería significar que al extrapolar linealmente 'hacia adelante' en la muestra y 'hacia atrás' en la muestra , obtendríamos un triángulo para el que sabemos que la señal de tiempo continuo debe estar debajo de ese triángulo. Si no fuera así, contendría contenido de mayor frecuencia. ¿No podríamos decir que max¿está realmente limitado por la frecuencia más alta ( ) (co) seno a la amplitud máxima que yo permita? (Tendré que leer la prueba en Pinsky si puedo encontrarla para ver cómo se relaciona)k1k|dx(t)/dt|αfN
Retinite
Todavía no puedo encontrar una prueba que entiendo, y no quiero gastar más de 100 USD para el libro de Pinksy solo para obtener 1 prueba. Mi instinto dice que podríamos estar seguros de que lugar de ( vs ), con el valor de señal máximo permitido. Encontré algunas pruebas generales aquí, pero no entiendo el uso de la norma con , y no estoy seguro de si la respuesta implica que sería una función ( aproximación) en el dominio de la frecuencia. |dx(t)dt|2πBAmax2π4πAmaxL1||g||1g(ξ)rect
Retinite
1
No importa, simplemente completando la explicación en la publicación mencionada anteriormente construí la prueba.
Retinite
5

Observaciones

He usado +1 y -1 en la secuencia en lugar de su 1 y 0. Con , la función continua limitada por banda en sus dos primeras figuras (con la modificación mencionada anteriormente) es:α=1fm(T)

(1)fm(T)=k=1mmsign(sinc(πkπ/2))sinc(πTπk),

dónde:

sinc(T)={sin(T)/Tif T01if T=0sign(x)={1if x<00if x=01if x>0.

fm(1/2) crece linealmente al logaritmo de :m

Crecimiento del pico en función de $ m $
Figura 1. trazada en función de El eje horizontal logarítmico linealiza el crecimiento como .fm(1/2)log2(m)m

Podemos simplificar con la ayuda de Wolfram Alpha :fm(1/2)

(2)fm(1/2)=2k=1m|sin(π(k0.5))π(k0.5)|={4/πif n=12(ψ(0)(1/2)ψ(0)(m+1/2))πotherwise,

donde es la función digamma . El término dominante de la serie de sobre es:ψ(0) (2)m=

2log(m)π,

que explica la linealización vista en la Fig. 1. Ahora podemos construir una versión normalizada de la función que hereda su limitación de banda pero no como :gm(T)fm(T)m

gm(T)=πfm(T)2log(m)

Como , parece aproximarse a una frecuencia sinusoide de Nyquist muestreada en sus ceros:mgm(T)

$ g_ {100000} (T) $
Figura 2. no .g100000(T)

El teorema de muestreo original de Nyquist-Shannon requiere que la frecuencia más alta esté por debajo de la mitad de la frecuencia de muestreo, por lo que parece que tenemos un caso límite que no está cubierto por él. Sin embargo, todavía están cubiertos arbitrariamente grandes finitos consecuencia, arbitrariamente finitos grandes .mfm(1/2)

Esquema de prueba

Aquí hay un resumen para una prueba de su declaración original: deje que el período de muestreo sea 1. Deje estar limitado a la frecuencia inferior , donde representa una frecuencia con un período de 2 y . Deje es finito para todos entero . Excluir el caso trivial para todo . Deje . De ello se desprende que para algunos . Ya sea:f(T)αππα<1f(T)Tf(T)=0Tg(T)=f(T)/supTf(T)g(T)0T

Caso 1. para algún entero . es finito para todos .g(T)0TsupTf(T)T

Caso 2. para todos número entero . es infinito por alguna . Hasta un factor de escala, está determinado por una fracción de sus ceros. Utilice uno más de los ceros restantes para hacer desaparecer la función: para todo . Esta es una contradicción, ya que anteriormente se determinó que para algunos . El caso 2 no puede ser cierto.g(T)=0TsupTf(T)Tg(T)αg(T)=0Tg(T)0T

De ello se desprende que el caso 1 es verdadera y es finito para todos .f(T)T

Sería bueno encontrar una prueba definitiva de que una parte de los ceros distribuidos uniformemente se puede utilizar para reconstruir la función dado su ancho de banda relativamente bajo en comparación con la densidad media de esos ceros. Supongo que si , el teorema de muestreo es suficiente para hacer que desaparezca. En la literatura, he encontrado algunas declaraciones de interés:α<1g(T)

La prueba de la parte 2 del teorema 4.1 mostró que una señal de banda limitada con con la propiedad de que la señal desapareció en los puntos debe desaparecer de forma idéntica.=πx=nZ

Jeffrey Rauch, " Serie de Fourier, integrales y muestreo de análisis complejo básico ".

Es bien sabido que puede tener más ceros, más o menos, que cos sin desaparecer de forma idéntica.gλt

BF Logan, Jr. " Información en los cruces cero de las señales de paso de banda ", Bell System Technical Journal, vol. 56, págs. 487-510, abril de 1977

La mayoría de los resultados sobre la especificación única de señales unidimensionales se basan en el hecho de que una función limitada de banda es completa (analítica en todas partes) y, por lo tanto, se especifica de manera única por sus ceros (reales y complejos) dentro de un factor constante y exponencial. Una función arbitraria de límite de banda se especifica de manera única por sus cruces por cero (reales) si se garantiza que todos sus ceros sean reales.
...
El trabajo adicional ha implicado la identificación de señales que se especifican de manera única por sus cruces por cero (reales) a pesar de que también contienen ceros complejos. Esto es posible si la velocidad de cruce por cero es, en cierto sentido, más alta que la velocidad de información o el ancho de banda de una señal.

SR Curtis, " Reconstrucción de señales multidimensionales de cruces por cero ", tesis, MIT, 1985.

Olli Niemitalo
fuente
Interesante. El enfoque derivativo máximo de la otra publicación demuestra dar una estimación del peor de los casos. Quiero abordar este problema nuevamente desde mi lado inicial (y el tuyo). Básicamente, podríamos decir que la señal consta de dos cosenos (uno corriendo hacia adelante y el otro hacia atrás) unidos en las muestras 0 y 1. Me parece que debería ser 'bastante fácil' rehacer su análisis para y etc. para crear una función o varias estimaciones de esta g (m) o f (1/2) que también dependen de . α=0.50.25α
Retinite
@Retinite quizás no sea tan fácil porque primero habría que asegurarse de que las muestras realmente codifiquen una función limitada en banda como se anuncia.
Olli Niemitalo el
Gracias por la prueba! Para me ocurrió: . Esto proporciona la secuencia [... 1 0 -1 0 1 1 0 -1 0 1 ...]. (¡¿esto es realmente BL a ?!) Esto no podría simplificarse (automáticamente) de modo que pudiera obtener una expansión en serie alrededor de . Lo que sí obtengo es una serie geométrica bastante obvia en para la que puedo verificar si la suma es finita (y lo es) y cuál podría ser ese valor. Pero este sigue siendo un método de fuerza bruta parcial. α=0.5fm(T)=k=1mmsign(cos(kπαπ4(sign(k12)+1))))sinc(kT)α=0.5mm
Retinite
Puede probar si la secuencia representa una función de banda limitada de media banda. Compare la "interpolación" dada por el núcleo sinc de banda completa y el núcleo sinc de media banda. Si en algún los dos no convergen como , la respuesta es no. (Comillas porque también puede probar en los puntos de muestra.)Tm
Olli Niemitalo
1

Considere la función de banda limitada con la transformada de Fourier que se puede recuperar perfectamente (¡desde la interpolación!) De sus muestras espaciadas a segundo de distancia, aunque las muestras incluyen solo la central pico y perder todos los otros máximos y mínimos locales de la función sinc. Retrase la función sinc en segundos para que la muestra pierda el pico central por completo, pero en su lugar obtenga muestras adyacentes con valores idénticos El exceso del máximo es, por lo tanto,sinc(t)rect(f)112

sinc(12)=sinπ/2π/2=2π.
12π. No tengo una prueba, pero sospecho que esto será el sobreimpulso máximo para el caso . Los valores más pequeños de darán muestras más cercanas al valor pico de y el sobreimpulso será correspondientemente más pequeño.α=1α1
Dilip Sarwate
fuente