1er Laboratorio
Teoria de Control Automtico 1
INTRODUCCION AL MATLAB
1. REPETIR Y EJERCITAR LOS SIGUIENTES COMANDOS EN MATLAB Definicin de una constante:
a= 1 b= [1 2]
Matriz constante:
A=[2 2 3 0 0 7 5 9 -1] A=[2 2 3; 0 0 7;5 9 -1]
Escribiendo nmeros complejos:
a=2+i b=-5-3*i
Podemos usar operaciones con anteriormente:
a=1;b=2;
matrices usando objetos definidos
Obsrvese que si colocamos punto y coma al final de la expresin, no es mostrado en la pantalla, lo que puede ser conveniente en unas situaciones. Expresin booleana:
A==1 A=[a+b pi 3 b^2 0 atan(a) 5 sin(b) -1]
Vector constante:
v=[1 2 3 4 5] v=1:5
Podemos formar matrices y vectores de zeros:
B=zeros( )
1er Laboratorio Matriz de ceros con 2 filas y 3 columnas:
B=zeros(2,3)
Teoria de Control Automtico 1
Matriz de zeros con las dimensiones de la matriz A:
A=[2 2 3;0 0 7;5 9 -1]; B=zeros(size(A))
Formando una matriz diagonal con los elementos de la diagonal principal de una matriz:
C=diag(diag(A))
OPERACIONES: De modo semejante, podemos formar matrices y vectores de unos: Matrices de unos con 2 filas y 3 columnas:
C=ones(2,3)
Matriz identidad:
A=diag(ones(1,3)) A=eye(3)
MATRICES DIAGONALES: Matriz diagonal con los elementos de diagonal principal yendo de 1 a 5:
D=diag(1:5)
Suma de matrices (recuerde las matrices deben tener la misma dimensin):
B=A+A
Sumar 1 a todos los elementos de una matriz:
C=B+1
Extrayendo los elementos de la diagonal principal:
A=[1 2 3 4 5 6 7 8 9] B=diag(A)
1er Laboratorio Multiplicacin de matrices:
A=[1 2 3;4 5 6;7 8 9] C=[1 2 0;0 0 1;0 2 3] D=A.*B
Teoria de Control Automtico 1 Rank (rango) de una matriz:
r=rank(A)
Matriz transpuesta:
B=A
Multiplicacin elemento a elemento:
A=[1 0 0;0 2 3;5 0 4] B=[2 0 0;0 2 2;0 0 3] C=A.*B
Inversa de una matriz:
A=[0 1;-2 -3] B=inv(A) A*B
Extraccin de la fila 2:
a=C(2,:)
Extraccin de columna 3:
b=C(:,3)
Determinante de una matriz:
d=det(A)
POLINOMIOS: Traza de una matriz:
A=[1 2 3;4 5 6;7 8 9] t=trace(A)
Polinomio p1 con races en 0 e -1:
v=[0 -1] p1=poly(v)
1er Laboratorio Polinomio p2 con coeficientes 1 e 2 e 1:
P2=poly([1 2 1])
Teoria de Control Automtico 1
CLCULO DE RACES:
P=roots(p1)
AUTOVALORES E AUTOVECTORES: Autovalores:
A=[0 1;-2 -3] r=eig(A) %r vector de autovalores
[V,D]=eig(A) % produce una matriz diagonal D de autovalores % y una matriz completa V cuyas columnas son % sus correspondientes autovectores. As (A*V=V*D)
1er Laboratorio Funciones:
Teoria de Control Automtico 1
function[y]=mifuncion(x) %definimos un archivo mifuncion.m if x<0 %donde estar nuestra 'y=f(x)' y=x^2 else y=sin(x*(pi/180)) end
Una vez definida la funcin podemos evaluarla en el punto x=30 (grados)
>> y=mifuncion(30) y = 0.5000
Otra funcin (mifuncin2): la funcin
en Matlab.
function[y]=mifuncion(x) x=input('ingrese su vector x=') y=x+x.^(2)+sin((2*pi)*x) figure(1) % abre nueva ventana de grfico plot(x,y,'r') % la curva es de color rojo figure(2) % abre nueva ventana de grfico plot(x,y,'b:') % la curva es de color azul y entrepunteada end
Cree un vector variando de -2 a 3 con paso de 0.5 Calculando la funcin:
y=mifuncion2(x) ingrese su vector x=-2:0.5:3 x = Columns 1 through 8 -2.0000 1.5000 -1.5000 -1.0000 -0.5000 0 0.5000 1.0000
Columns 9 through 11 2.0000 y = Columns 1 through 8 2.0000 3.7500 0.7500 0.0000 -0.2500 0 0.7500 2.0000 2.5000 3.0000
Columns 9 through 11 6.0000 8.7500 12.0000
1er Laboratorio
Teoria de Control Automtico 1
Ploteando el resultado: Creando nuevas ventanas grficas (figure(1), figure(2),,figure(n)):
figure(1) plot(x,y,`r)
12
% abre nueva ventana de grfico % la curva es de color rojo
10
-2 -2
-1.5
-1
-0.5
0.5
1.5
2.5
figure(2) plot(x,y,b:`)
% abre nueva ventana de grfico % la curva es de color azul y entrepunteada
12
10
-2 -2
-1.5
-1
-0.5
0.5
1.5
2.5
1er Laboratorio 2. EJERCICIOS Implemente en Matlab la siguiente funcin, luego plotee
Teoria de Control Automtico 1
Implementando la funcin en exm1.m
function y =exm1(x) if x==1 y = 'indeterminado'; elseif x==-1 y = 'indeterminado'; else y = 1./((x.^2)-1); end
Definiendo el vector x : x=-10:0.1:10 Definiendo y Graficando : y=exm1(x) : plot(x,y)
1er Laboratorio Implemente en Matlab la siguiente funcin, luego plotee
Teoria de Control Automtico 1
Implementando la funcin en exm2.m
function y =exm2(x1,x2) if x1>0 if x2>0 y = x1 + x2; end else y = sqrt((x1.^2)+(x2.^2)); end
Graficando mediante matrices en exm2_1.m :
1er Laboratorio Haga un m. file que ayude a encontrar el minimo de intervalo (0,2) La funcin esta implementada en exm3.m:
function y =exm3(x) y=(x.^3) -(2*x) - (5);
Teoria de Control Automtico 1 , dentro del
Para hallar el minimo se usa el archivo exm3_1.m
paso=input('ingrese paso') x=[0:paso:2]; y=exm3(x) minimo=min(y)
Construya una seal escalon unitario de 0 a 50 segundos, con step inicial en 25 seg. El paso deber ser de 0,5seg. Plotee el resultado. La funcin implementada esta en exm4.m:
t0 = [0:0.5:25]; u0 = zeros(size(t0)); t1 = [25:0.5:50]; u1 = ones(size(t1)); t = [t0 t1]; u = [u0 u1]; plot(t,u)
La grafica es:
1er Laboratorio Construya una seal peine de dirac, plotee el resultado. La funcin esta en exm5.m :
r=-5:1:5; p=zeros(size(r)); p(r>=-30)=1; stem(r,p); axis([-5 5 -1 2]);
Teoria de Control Automtico 1
-La grafica es :
3. CONCLUSIONES Se ha aprendido a utilizar los comandos bsicos de Matlab. Matlab es un programa que nos ayuda a hacer, calcular y graficar funciones en 2d o 3d, adems de hacer funciones especuales como el impulso unitario y funcin delta de dirac. Matlab tiene una forma fcil de almacenar archivos, archivos.m
10
