¿Existe una distorsión no lineal limitada en la banda?

Entonces, si genera una onda cuadrada simplemente cambiando una señal entre dos valores, en los límites de la muestra, produce una serie infinita de armónicos, que alias y producen tonos por debajo de su fundamental, que es muy audible. La solución es la Síntesis de banda limitada , ya sea usando síntesis aditiva o pasos de banda limitada para producir formas de onda que son las mismas que si hubiera limitado la banda de la onda cuadrada matemática ideal antes de muestrearla:

http://flic.kr/p/83JMjT

Pero me di cuenta de que si aplica una amplificación grande a una onda sinusoidal digital y luego la recorta digitalmente, producirá la misma forma de onda cuadrada, sin las ondas del fenómeno Gibbs. Entonces, también está produciendo productos de distorsión con alias, ¿verdad? Entonces, ¿ cualquier distorsión no lineal en el dominio digital que produzca armónicos fuera de los límites de Nyquist producirá productos de distorsión con alias? (Editar: he hecho algunas pruebas y he confirmado que esta parte es cierta).

¿Existe una distorsión de banda limitada para simular (en el dominio digital) los efectos de la distorsión (en el dominio analógico) antes de la limitación de banda y el muestreo? Si es así, como lo haces? Si busco "distorsión ilimitada de banda", encuentro algunas referencias a los polinomios de Chebyshev, pero no sé cómo usarlos o si solo funcionan para ondas sinusoidales o qué:

Este instrumento no intenta generar una distorsión de banda limitada. Los interesados ​​en la distorsión de banda limitada deben investigar el uso de polinomios de Chebyshev para generar el efecto. Distorsión de la tangente hiperbólica

 

"Chebyshev polinomio" - dar forma a las funciones con la importante propiedad de que son intrínsecamente de banda limitada, es decir, que no introducen armónicos espectrales espurias debido a la superposición etc. Wave Shaper

La aplicación de una función no lineal siempre introducirá armónicos, y la combinación de funciones no lineales con versiones muestreadas de señales continuas agrega la arruga que notas arriba (donde los armónicos de alta frecuencia tienen un alias a bajas frecuencias).

Se me ocurren algunas formas de proceder:

1. Puede usar un factor de sobremuestreo lo suficientemente alto como para capturar los armónicos adicionales (hasta una precisión arbitraria, por ejemplo, su piso de ruido),
2. Puede usar una función de recorte "más suave" (consulte, por ejemplo, aquí ) que tiene armónicos que se extinguen antes que el cortapelos duro. Esto es más fácil de modelar, pero introduce su propia distorsión a bajas frecuencias.
3. Sobre la base del enfoque que sugirió anteriormente, interpole su señal muestreada (por ejemplo, utilizando un interpolador Lagrange o Chebyshev) para construir un modelo de tiempo continuo. Luego, aplique el cortapelos duro y pase bajo en un dominio de tiempo continuo simulado. Muestra el resultado.

Puedes combinar (1) y (2). El tercer enfoque es complejo, pero le brinda el mejor control sobre la cantidad de distorsión que debe admitir, y probablemente escalará mejor a requisitos de muy alta fidelidad.

$f(x)$

$$f(x)=\sum_{n=0}^\infty \left[ \frac{f^{(n)}(0)}{n!} x^n\right]$$

$f(x)$$x=g(t)$$g(t)$$x^n$$g(t)^n$$n$$n$$n$veces la de tu señal. Para completar la imagen, debe calcular la amplitud asociada con cada término y decidir cuántos términos en la suma son relevantes.

(Con un poco de reflexión, también podría usar este formulario directamente para aproximar la no linealidad filtrada. Eso requeriría una buena representación en serie para el clipper).

Algunos enfoques para la distorsión no lineal sin alias (en orden creciente de dificultad):

1. $\frac{f_s}{2N}$$f$$f^{N+1}$

2. $N$$2N$

3. $f:\mathbb{R}\to\mathbb{R}$$f:\mathbb{R}^N\to\mathbb{R}^M$$N$$M>N$

4. Restringir el diseño algebraico basado : en el elemento anterior, ha visto que la distorsión no lineal antialiasing conduce a filtros no lineales. Por supuesto, no todos los filtros no lineales están libres de alias, pero algunos pueden estarlo. Entonces, la pregunta obvia es un criterio para hacer que dicho filtro sea estrictamente libre de alias y cómo diseñarlo. Como resultado, una afirmación equivalente a estar libre de alias es que el filtro no lineal conmuta con traducciones de submuestra. Por lo tanto, debe asegurarse de que no haga una diferencia si traduce primero y luego filtra, o filtra primero y luego traduce. Esta condición conduce a restricciones de diseño muy estrictas.para filtros no lineales, pero depende de cómo realice la traducción de la señal. Por ejemplo, la traducción ideal requeriría infinitos coeficientes para el filtro no lineal. Por lo tanto, debe aproximar la traducción de la señal a un orden finito para obtener un filtro no lineal finito. La ausencia de alias se escala con la aproximación que usa, pero tiene un muy buen control sobre ella. Después de haber trabajado con las matemáticas de este enfoque, puede diseñar cualquier función de transferencia no lineal (no solo suave) como un modelo digital casi ideal en forma de filtro no lineal. No puedo bosquejar los detalles aquí, pero tal vez puedas encontrar algo de inspiración en esta descripción.

$$T_n(x) = \mathrm{cos}(n \, \mathrm{arccos}(x)).$$

$T_n(x)$

$$T_n(\mathrm{cos}(kx)) = \mathrm{cos}(n\,\mathrm{arccos}(\mathrm{cos}(kx))) = \mathrm{cos}(nkx). \tag{1}$$

Los propios polinomios se pueden generar fácilmente utilizando la siguiente relación de recurrencia :

$$T_0(x) = 1 \\ T_1(x) = x \\ T_n(x) = 2x\,T_{n-1}(x) - T_{n-2}(x).$$

Aquí están los primeros:

$$\eqalign{ T_0(x) &= 1 \\ T_1(x) &= x \\ T_2(x) = 2x(x - 1) &= 2x^2 - 1 \\ T_3(x) = 2x(2x^2 - 1) - x &= 4x^3 - 3x \\ T_4(x) = 2x(4x^3 - 3x) - (2x^2 - 1) &= 8x^4 - 8x^2 + 1 \\ \ldots }$$

$(1)$$T_2$$\mathrm{cos}(x)$

$$\eqalign{2\,\mathrm{cos}^2(x) - 1 &= 2\left(\frac{e^{ix} + e^{-ix}}{2}\right)^2 - 1 \\ &= \frac{2}{4}(e^{i2x} + 2e^{ix}e^{-ix} + e^{-i2x}) - 1 \\ &= \left(\frac{e^{i2x} + e^{-i2x}}{2}\right) + \frac{2}{2} - 1 \\ &= \mathrm{cos}(2x). }$$

Al calcular una serie de Chebyshev

$$f(x) = \sum_{n = 0}^{\infty}{a_n\,T_n(x)}$$

$n$$f(x)$

@ robert-bristow-johnson explica esto muy claramente en comp.dsp :

tienes que sobremuestrear hasta cierto punto. si representa la no linealidad (sin memoria, supongo) como un polinomio de orden finito (que se aproxima a la curva que está tratando de implementar), entonces cualquiera que sea el orden del polinomio es el mismo factor de sobremuestreo necesario y no se producirán alias. luego filtro de paso bajo (a esa frecuencia sobremuestreada) para deshacerse de todos los componentes de frecuencia más altos que su Nyquist original, luego disminuya la muestra y no tendrá alias.

En otras palabras, si su no linealidad es un polinomio, la frecuencia más alta que puede producir la distorsión será la frecuencia más alta en su señal multiplicada por el orden N del polinomio. (La no linealidad polinómica está multiplicando la señal por sí misma N veces, por lo que su espectro se enreda consigo mismo y se extiende por la misma relación).

Entonces conoce la frecuencia máxima (ya sea Nyquist o algún límite inferior para su aplicación), y conoce el orden del polinomio, por lo que puede sobremuestrear lo suficiente para evitar el aliasing, hacer la distorsión y luego el filtro de paso bajo y la disminución de la muestra.

De hecho, puede reducir la tasa de sobremuestreo permitiendo que ocurra algún aliasing, siempre que esté contenido en la banda que se eliminará antes de la disminución de muestreo:

Otro pequeño truco es que no necesita preocuparse por los alias que se pliegan hacia el área que LPF eliminará. por lo tanto, un polinomio de quinto orden solo necesita tener una relación de sobremuestreo de 3. esos 2 armónicos superiores pueden tener alias, pero no volverán a la banda base. al reducir la resolución, filtra los armónicos con alias. así que creo que la regla dura y rápida es

relación de sobremuestreo = (orden polinomial + 1) / 2