Experimento 3 MATLAB Analisis Frecuencia
Experimento 3 MATLAB Analisis Frecuencia
Experimento 3 MATLAB Analisis Frecuencia
Objetivos:
Alcance:
Experimento 3 1
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
Ejemplo 1.
ft = exp(a*t)*cos(w*t)
In Control System Toolbox™ software, you represent linear systems as model objects. Model
objects are specialized data containers that encapsulate model data and other attributes in
a structured way. Model objects allow you to manipulate linear systems as single entities
rather than keeping track of multiple data vectors, matrices, or cell arrays.
Experimento 3 2
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
The data encapsulated in your model object depends on the model type you use. For
example:
Other model attributes stored as model data include time units, names for the model inputs
or outputs, and time delays.
The following diagram illustrates the relationships between the types of model objects in
Control System Toolbox™, Robust Control Toolbox™, and System Identification Toolbox™
software. Model types that begin with id require System Identification Toolbox software.
Model types that begin with u require Robust Control Toolbox software. All other model types
are available with Control System Toolbox software.
Numeric Models
Experimento 3 3
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
Numeric LTI models are the basic numeric representation of linear systems or components
of linear systems. Use numeric LTI models for modeling dynamic components, such as
transfer functions or state-space models, whose coefficients are fixed, numeric values. You
can use numeric LTI models for linear analysis or control design tasks.
The following table summarizes the available types of numeric LTI models
Control System Toolbox software supports transfer functions that are continuous-time or discrete-
time, and SISO or MIMO. You can also have time delays in your transfer function representation.
Syntax
tf
sys = tf(num,den)
sys = tf(num,den,Ts)
sys = tf(M)
sys = tf(num,den,ltisys)
tfsys = tf(sys)
tfsys = tf(sys, 'measured')
tfsys = tf(sys, 'noise')
tfsys = tf(sys, 'augmented')
Experimento 3 4
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
Syntax
sys = zpk(z,p,k)
sys = zpk(z,p,k,Ts)
sys = zpk(M)
sys = zpk(z,p,k,ltisys)
s = zpk('s')
z = zpk('z',Ts)
zsys = zpk(sys)
zsys = zpk(sys, 'measured')
zsys = zpk(sys, 'noise')
zsys = zpk(sys, 'augmented'
pzplot: Pole-zero map of dynamic system model with plot customization options
Syntax:
h = pzplot(sys)
pzplot(sys1,sys2,...)
pzplot(AX,...)
pzplot(..., plotoptions)
Experimento 3 5
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
num = [1 0];
den = [1 3 2];
G = tf(num,den);
num and den are the numerator and denominator polynomial coefficients in descending powers of s. For
example, den = [1 3 2]represents the denominator polynomial s2 + 3s + 2.
G is a tf model object, which is a data container for representing transfer functions in polynomial form.
Tip Alternatively, you can specify the transfer function G(s) as an expression in s:
Create a transfer function model for the variable s.
s = tf('s');
1. Specify G(s) as a ratio of polynomials in s.
Experimento 3 6
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
:
Z = [0];
P = [-1-1i -1+1i -2];
K = 5;
G = zpk(Z,P,K);
Z and P are the zeros and poles (the roots of the numerator and denominator, respectively). K is the gain of
the factored form. For example,G(s) has a real pole at s = –2 and a pair of complex poles at s = –1 ± i. The
vector P = [-1-1i -1+1i -2] specifies these pole locations.
G is a zpk model object, which is a data container for representing transfer functions in zero-pole-gain
(factorized) form.
Ejemplo 2.
Para la función:
50 s
----------------------------------------------
s^4 + 6 s^3 + 50 s^2 + 150 s + 625
Encontrar:
a. Polos y ceros.
b. Diagrama de polos y ceros.
c. Expansión en fracciones parciales.
Experimento 3 7
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
Ejemplo 3.
Experimento 3 8
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
Ejemplo 4
𝟎, 𝟖
𝑮(𝒔) =
(𝟑𝒔 + 𝟏)(𝟏𝟎𝒔 + 𝟏)(𝟑𝟎𝒔 + 𝟏)
Experimento 3 9
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
Problema 1
Problema 2
1
𝐻(𝑠) =
𝑠4 + 2.6131𝑠 3 + 3.41412𝑠 2 + 2.6131𝑠 + 1
Problema 3
Experimento 3 10
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
6
x 10 Pole-Zero Map
3
-1
System: zin
Zero : -2.78e+05 - 2.57e+06i
Damping: 0.108
-2 Overshoot (%): 71.2
Frequency (rad/s): 2.58e+06
-3
-3 -2.5 -2 -1.5 -1 -0.5 0
-1 5
Real Axis (seconds ) x 10
Z(0) = 1,8 kΩ
Notas
Experimento 3 11
PONTIFICIA UNIVERSIDAD JAVERIANA
FACULTAD DE INGENIERIA - DEPARTAMENTO DE ELECTRONICA
4096 CIRCUITOS EN FRECUENCIA
Revisiones:
Experimento 3 12