Facultan de produccion y servicios

UNIVERSIDAD NACIONAL DE
SAN AGUSTIN

FACULTAN DE PRODUCCION Y SERVICIOS

ESCUELA PROFECIONAL DE
ING. ELECTRONICA

TEORIA DE CONTROL AUTOMATICO 1

LABORATORIO 03

RESPUESTA EN FRECUENCIADE LOS SCA

GRUPO 1:

PACHECO HUALLLA ELMER CUI: 20070953

AREQUIPA-PERÚ

2010

2010 – 02 – 24

Parte 1

%////////////////////////////////////////////////
%Sistema de Suspensión%////////////////////////////////////////////////
clear %borra las variables antiguas
%////////////////////////////////////////////////
%Parámetros
m=250; %masa suspensa (Kg)
k=10000; %rigidezdel resorte (N/m)
b=316.23; %amortecimiento (Ns/m)
%Matrices del sistema
A=[0 1;-k/m -b/m];
B=[-1;b/m];
C=[-k/m -b/m];
D=[b/m];
%Sistema en el espacio de estados
G=ss(A,B,C,D);
%Respuesta enfrecuencia
%Diagrama de bode del sistema (observe escalas de bode (G))
bode(G)
%La escala horizontal es la logarítmica, y muestra la frecuencia en rd/s
%La escala vertical del grafico superior eslineal, y muestra la ganancia
%en db
%La escala vertical del grafico inferior es lineal y muestra la fase en
%grados
%MATLAB asigna valores por defecto a nuestro grafico de bode, el Matlab%mostro que la amplitud de la ganancia está entre -0Db y +40Db, la fase esta
%entre -90 y 90 grados, y que la frecuencia está entre 10 y 1000 rad/s.
%Vamos a construir un vector w donde el logaritmo dela distancia entre un
%valor y otro es constante, comenzando en 10 elevado a (-3), terminando en
%10 elevado a (4), con 2000 puntos, y convencionando que la unidad es
%rad/s:
w=logspace(-3,4,2000);%Determinamos los vectores de magnitud y fase,
[mag,phase]=bode(G,w);
% mag y fase son arrays de 3 dimensiones, esto quiere decir que la función
% bode trabaja con sistemas multivariables (mimo)% Entonces, para nuestro caso, debemos solo capturar el vector
% correspondiente
mag=mag(1,:);
phase=phase(1,:);
%ploteando el bode de ganancia
subplot(2,1,1)
semilogx(w,20*log10(mag))...