Osciladores Electronicos No Lineales
Osciladores Electronicos No Lineales
Osciladores Electronicos No Lineales
DOCUMENTO I. MEMORIA
Parte I. Memoria pág. 1 a 28 28 páginas
DOCUMENTO II. PRESUPUESTO
1. Mediciones pág. 3 a 3 1 páginas
2. Precios unitarios pág. 4 a 4 1 páginas
3. Sumas parciales pág. 5 a 5 1 páginas
4. Presupuesto general pág. 6 a 6 1 páginas
Autorizada la entrega del proyecto del alumno:
Fdo.: . . . . . . . . . . . . . . . . . . . . . . . . Fecha: . . . . . . / . . . . . . / . . . . . . . . .
Fdo.: . . . . . . . . . . . . . . . . . . . . . . . . Fecha: . . . . . . / . . . . . . / . . . . . . . . .
UNIVERSIDAD PONTIFICIA COMILLAS
ESCUELA TÉCNICA SUPERIOR DE INGENIERÍA (ICAI)
INGENIERO INDUSTRIAL
Para realizar este estudio se han empleado dos herramientas informáticas, MATLAB
y PSpice. Utilizando MATLAB se ha diseñado un modelo matemático que reproduce
con cierta fidelidad lo que ocurre en la realidad. Para diseñar este modelo se ha tenido
en cuenta la saturación de la etapa amplificadora, que es la que hace que los osciladores
reales no sean lineales. Además, para hacer el modelo aun más parecido a la realidad se
ha introducido el efecto que producen los codos de saturación en la característica del
amplificador. Se ha podido observar durante el desarrollo del modelo que la respuesta de
un oscilador de segundo orden depende principalmente de 3 parámetros que en el texto
hemos llamado µ, δ y m, que representan el efecto de las no linealidades, la magnitud
de estas y el efecto de los codos de saturación respectivamente.
Cuando se compraran los modelos no lineales con el lineal, se puede apreciar que el
espectro del modelo no lineal no muestran una única frecuencia, y que las trayectorias en
el plano de fases no son circulares. Para tratar de reducir el efecto de las no linealidades
debemos conseguir que µ 1 para que la señal obtenida se lo mas parecida a una
senoidal pura. En el caso del puente de Wien esto no es posible, ya que el valor de µ
viene fijado por los elementos que componen el circuito.
Por último, se han realizado ensayos en el laboratorio que muestran que los modelos
anteriormente citados son correctos. Para mejorar el modelo matemático del oscilador
sería necesario introducir ciertas modificaciones al circuito, ya que en la realidad
la saturación de la etapa amplificadora no es simétrica, lo que hace que aparezcan
componentes en frecuencia inesperadas cuando analizamos el espectro de la señal que
produce el oscilador. Este problema no es muy grave ya que existen métodos muy
sencillos y baratos para hacer que la característica de la saturación sea simétrica.
This project is a study about nonlinear oscillators. The method followed has
consisted in the analysis of the temporal response of these devices followed by a
comparison with the linear models found in textbooks. In fact, this research has focused
on the study of the Wien bridge, a second-order oscillator widely used in low-frequency
electronic applications.
To carry out this research two computational tools have been used: MATLAB and
PSpice. Using MATLAB a mathematical model has been designed which reproduces
with certain fidelity what happens in reality. To design this model the saturation of the
amplifier stage had to be taken into account, which makes real oscillators non-linear.
Furthermore, to make the model more accurate it was taken into consideration the effect
produced by the saturation elbows of the characteristic of the amplifier. During the
process it was observed that the model’s response of a second order oscillator depends
mainly on three parameters. In the text these are called µ, δ and m, representing the
nonlinear effect, the magnitude of these and the effect of saturation elbows respectively.
Once analyzed, the mathematical model was compared to a Wien bridge design in
PSpice. Comparing the spectrum and the phase planes, it was noted that there were
no big differences between the two of them, existing only one small difference in the
amplitudes of the signals, which was due to the fact that the model developed with
PSpice was much more accurate than the model in MATLAB. Even so, the mathematical
model shows with sufficient precision the behaviour of the Wien bridge.
When comparing non-linear models with linear ones, it can observed that the
spectrum of the nonlinear model does not show a single frequency and the phase plane
trajectories are not circular. To attempt to reduce the effect of these nonlinearities, µ
was set to be µ 1 so that the signal obtained would be as close as possible to a pure
sine signal. In the case of the Wien bridge this is not possible because the value of µ is
fixed by the elements of the circuit.
Finally, laboratory tests were carried out which show that the models developed are
correct. To improve the mathematical model of the oscillator, it would be necessary to
introduce certain changes to the circuit because in a real environment the saturation of
the amplifier stage is not symmetrical. This fact causes the apparition of unexpected
frequency components when analyzing the spectrum of the signal produced by the
oscillator. This problem is not very serious as there are very simple and inexpensive
methods to make the saturation characteristic symmetric.
1. Introducción 4
2. Oscilador Puente de Wien 6
2.1. Condiciones de oscilación . . . . . . . . . . . . . . . . . . . . . . . . 6
2.2. Puente de Wien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3. Oscilador de Lienard 9
3.1. Teorema de Lienard . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
3.2. Trayectorias en el plano de fase . . . . . . . . . . . . . . . . . . . . . 10
3.3. Ejemplo: Puente de Wien normalizado . . . . . . . . . . . . . . . . . 11
3.3.1. Normalización . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.3.2. Puente de Wien normalizado . . . . . . . . . . . . . . . . . . . . 13
4. Comparación teórica de los modelos lineal y no lineal 15
4.1. Método de resolución . . . . . . . . . . . . . . . . . . . . . . . . . . 15
4.2. Modelo lineal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
5. Comparación numérica de los modelos lineal y no lineal 20
5.1. Espectro . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
5.2. Trayectorias en el plano de fases . . . . . . . . . . . . . . . . . . . . 22
6. Comparación experimental de los modelos lineal y no lineal. Conclusiones. 24
6.1. Ensayos del puente de Wien . . . . . . . . . . . . . . . . . . . . . . 24
6.2. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
Bibliografía 28
Introducción
Wien que se usan mucho en aplicaciones de baja frecuencia que no requieren alta pureza
espectral o estabilidad.
ẍ + ω02 · x = 0 (1)
Vo A
= (3)
Vin 1−A·B
A·B =1 (4)
|A · B| = 1 (5)
Esto ocurre a una frecuencia, ω0 , conocida como frecuencia de oscilación. Hay que
tener en cuenta que en la práctica es imposible conseguir que A · B sea exactamente
1. En la práctica |A · B| > 1 de modo que la amplitud aumenta de forma exponencial
hasta que sature la etapa amplificadora. Por lo tanto, la señal de salida está limitada por
la saturación del amplificador.
Puesto que el oscilador debe ser capaz de generar la señal por sí mismo, es decir,
sin ninguna señal de entrada ni excitación externa, cabe preguntarse por que arranca la
oscilación. La respuesta se encuentra en el ruido eléctrico. En el sistema hay presentes
pequeñas señales armónicas de todas las frecuencias. Entre ellas, el oscilador amplifica
de manera exponencial aquella cuya frecuencia coincida con ω0 hasta que su amplitud
alcanza la de saturación. Este resultado es el que se observa a la salida como oscilación,
es decir, una señal sinusoidal con cierto grado de distorsión armónica.
R2 C2
+
R1 C1
Vin
-
A V0
Vo s·R·C
= (7)
Vin 1 + 3 · s · R · C + (s · C · R)2
1
Sustituyendo s = j · ω0 con ω0 = R·C y aplicando las condiciones (5) y (6)
obtenemos que la realimentación atenúa la oscilación en 1/3 por lo que la ganancia en
la etapa amplificadora debe ser al menos 3.
R3
A=1+ =3 (8)
R4
Ésta es la ganancia crítica. Para que el circuito comience a oscilar por sí mismo
deberá tener cierto grado de inestabilidad, esto significa que la ganancia deberá ser
ligeramente superior a 3. La amplitud de la oscilación estará limitada por la saturación
del operacional, que depende de su polarización, con la consiguiente distorsión armónica
que esto provoca.
Oscilador de Lienard
donde las funciones ’g’ y ’φ’ existen en el dominio de R, además ’g’ es una función
impar y ’φ’ una función par ( ẍ y ẋ son derivadas con respecto al tiempo).
Φ(x)
1
+
- 1 + 1 x
-δ
x
Figura 3. Función φ(x)
El teorema de Lienard dice que si φ(x) es par, φ(0) < 0 y φ(x0 ) > 0 para x0 > δ+1,
entonces existe una solución periódica aislada, es decir un ciclo estable. Por ejemplo,
φ(x) de la Figura 3 cumple esta condición.
ẋ = y − µ · F (x) (10)
ẏ = −x (11)
Siendo F (x)
Z x
F (x) = φ(s) · ds (12)
0
Además, cabe destacar que si se cumplen las condiciones de Lienard, las soluciones
dependen únicamente de dos parámetros, µ y δ. A continuación se muestran varias
soluciones del sistema formado por las ecuaciones (10) y (11), estas están representadas
en los planos de fase. Estas trayectorias han sido obtenidas con MATLAB.
Como se puede observas si se comparan las dos figuras anteriores se observa que
para un valor grande de µ las trayectorias en el plano de fase de la solución dejan de ser
circulares (Figura 5). Por el contrario en la Figura 4 se puede ver que para un valor de
µ = 0, 1 la respuesta en el plano de fases describe trayectorias casi circulares.
3.3.1. Normalización
La estructura clásica de un oscilador es un amplificador realimentado con una red
selectiva en frecuencia (Figura 6).
k·s
H(s) = (13)
s2 + γ · s + ω02
x
para |x| < 1
sat(x) = 1 para x > 1 (15)
−1 para x < −1
v0 = E · sat(x) (16)
A · vin
x= (17)
E
µ=T ·γ (21)
k·A
φ(x) = 1 − · sat0 (22)
γ
A·k
δ= −1 (23)
γ
1
k = ω0 = (24)
R·C
3
γ= (25)
R·C
Por lo tanto,
µ=3 (26)
A
δ= −1 (27)
3
Como se puede ver, el ciclo límite no es circular. Esto quiere decir que a pesar de
que el puente de Wien es un buen oscilador no nos sirve si queremos señales senoidales
de alta pureza espectral.
ẋ = y (28)
y = −x − µ · φ(x) · ẋ (29)
Como ya se ha mencionado anteriormente, los modelos lineales nos dicen que las
soluciones de los osciladores son señales senoidales puras. En la práctica esto no es así,
en verdad la solución es una suma de señales senoidales.
Si nos centramos en los osciladores electrónicos vemos que las no linealidades que
aparecen se pueden modelar con dos parámetros, µ y δ. El parámetro µ nos indica
en cuánto afectan los términos no lineales mientras que δ indica la magnitud de las
no linealidades (Figura 3). Si además los efectos no lineales son producidos por la
saturación de la etapa amplificadora, también podemos tener en cuenta los codos de
saturación obteniendo así un modelo más parecido a la realidad. Para ello introducimos
Como se puede observar en la figura 10, los efectos de las discontinuidades finitas
de la función φ(x) son apreciables. Estos efectos se hacen cada vez menores al aumentar
m, es decir al hacer más suaves la pendientes de la función φ(x) (Figura 9).
Como se puede ver en las figuras, cuando el valor de µ es muy pequeño (µ 1),
la diferencia con el modelo lineal no es muy apreciable, pero para valores grandes de
µ el efecto de las no linealidades es importante. Tambien es interesante comparar las
respuestas temporales que aparecen en las figuras 12 y 13. Se puede apreciar que si
µ 1 la respuesta ya no es senoidal (Figura 13).
5.1. Espectro
El oscilador puente de Wien tiene armónicos en múltiplos impares de la frecuencia
fundamental. Como se puede ver en la figura, en la práctica solo son representativos los
armónicos 1, 3 y 5, ya que el resto tienen una amplitud demasiado pequeña. Aun así, el
tercer y quinto armónico son también muy pequeños comparados con el fundamental.
3.000V
2.000V
1.000V
0V
0.305KHz 1.000KHz 2.000KHz 3.000KHz 4.000KHz 5.000KHz 6.000KHz 7.000KHz 8.000KHz
V(C1:2)
Frequency
Date: June 11, 2010 Page 1 Time: 08:38:59
Las diferencias que aparecen en los valores de la amplitud se deben a que el modelo
matemático de MATLAB es mucho menos preciso que el modelo de PSpice. Además
el parámetro m ha sido fijado mediante experimentación comparando los planos de
fase. A continuación se muestra el espectro del modelo lineal. Se puede observar, que a
diferencia del modelo no lineal, solo lo muestra un pulso en la frecuencia fundamental.
Además la amplitud de la oscilación es mucho menor.
2.0A
0A
-2.0A
-4.0A
-6.0V -5.0V -4.0V -3.0V -2.0V -1.0V -0.0V 1.0V 2.0V 3.0V 4.0V 5.0V 6.0V
-I(C1) * EXP(7)
V(C1:2)
Date: June 11, 2010 Page 1 Time: 08:43:17
Como se puede ver en la Figura 18, las soluciones representadas en el plano de fases
utilizando PSpice, no son trayectorias circulares, que como se ha comentado antes se
debe a la característica no lineal de la etapa amplificadora. Esta no linealidad hace que
aparezcan componentes en frecuencia múltiplos de la fundamental tal y como se puede
ver en la Figura 15.
0
-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7
-1
-2
-3
-4
-5
-6
-7
5,00
4,00
3,00
2,00
1,00
0,00
-0,003 -0,002 -0,001 0 0,001 0,002 0,003
-1,00
-2,00
-3,00
-4,00
-5,00
6.2. Conclusiones
Como se ha podido observar, la respuesta de un oscilador de segundo orden depende
principalmente de 3 parámetros que hemos denominado µ, δ y m. Estos parámetros
ayudan a modelar el efecto de las no linealidades y la magnitud de estas. Los efectos no
lineales que aparecen en este tipo de osciladores son debidos a que la característica de
la etapa amplificadora es no lineal, ya que satura al llegar a ciertos niveles de tensión.
Además se a podido observar en las simulaciones realizadas que las aproximaciones al
modelo lineal son aceptables si µ 1, ya que el efecto de las no linealidades se hace
muy pequeño. En la práctica, conseguir un valor tan pequeño de µ no es posible debido
a que cuando se construye un modelo físico del oscilador el valor de este parámetro
viene fijado por los componentes del circuito. En concreto, en el oscilador puente de
Wien el menor valor posible es de µ = 2.
Otra conclusión que se puede sacar del estudio realizado, es que el uso de este tipo
de osciladores está limitado a aplicaciones que no requieren una alta pureza espectral
de la señal porque en el espectro aparecen frecuencias adicionales a la fundamental. En
el caso particular del puente de Wien aparecen componentes de frecuencia en todos los
múltiplos impares, aunque como se ha demostrado la amplitud de estos armónicos es
pequeña en comparación con la de la frecuencia fundamental del circuito. Además, en
la práctica, la saturación de los operacionales no es simétrica, lo que agrava el problema,
ya que introduce aun más componentes en frecuencias no deseadas. Para solucionar
este problema habría que limitar la saturación de la etapa amplificadora por debajo de
los niveles de alimentación del amplificador operacional.
El oscilador puente de Wien es un sistema de segundo orden que utiliza una red
RC. Este tipo de redes se utilizan para frecuencias bajas ( por debajo de 1MHz). Si
queremos obtener frecuencias mayores deberemos utilizar una red de realimentación
LC, como la que utilizan los osciladores de Hartley, Colpitts (sistemas de tercer orden)
o el oscilador resonante (de segundo orden). En concreto, el oscilador de Colpitts es
utilizado cuando se requiere una señal con muy pocas componentes armónicas.
En futuros desarrollos sería interesante estudiar estos circuitos y averiguar de qué
parámetros adicionales dependen. El método de estudio sería parecido al desarrollado
en este proyecto porque a pesar de tener un orden mayor, los criterios de oscilación
son los mismos que los del puente de Wien. Por ejemplo, una diferencia a la hora de
estudiar la estabilidad sería que las trayectorias que sigue la respuesta hasta llegar a un
ciclo estable recorren el espacio en lugar de estar contenidas en un plano. Además sería
interesante realizar el estudio a frecuencias mayores ya que este tipo de circuitos tienen
un buen comportamiento a altas frecuencias.
[1] Oetiker, T., Partl, H., Hyna, I., Schlegl, E., The Not So Short Introduction to
LATEX 2ε, Diciembre 2009.
[2] Sedra, A., Smith. K., Microelectronic Circuits, Oxford Univ. Press. 1998.
[3] Joedan, D. W., Smith, P., Nonlinear Differential Equations, Oxford Univ. Press.
1987 2nd edition.
[4] Pérez, C., MATLAB y sus Aplicaciones en las Ciencias y la Ingeniería, PEARSON
EDUCACIÓN. S.A. Madrid 2002.
1. Mediciones 3
2. Precios unitarios 4
3. Sumas parciales 5
4. Presupuesto general 6
Mediciones
Precios unitarios
En la tabla 2 se definen los precios de las actividades. Estas se han dividido en tres
tipos: ingeniería, técnico de laboratorio y redacción y maquetación. A continuación se
define que partidas corresponden a cada tipo de actividad:
Ingeniería:
• Investigación previa.
• Programación y diseño del modelo matemático.
Técnico de laboratorio:
• Trabajo en el laboratorio.
Escritura:
• Redacción y maquetación.
Sumas parciales
Presupuesto general
Tabla 4. Presupuesto