Informe Laboratorio 1 Control

download Informe Laboratorio 1 Control

of 11

Transcript of Informe Laboratorio 1 Control

INFORME LABORATORIO 1Caractersticas bsicas en el espacio del trabajo de MATLABOperaciones bsicas en matlab: OperacinSmboloExpresin en matlab

Suma+a+b

Resta-a-b

Multiplicacin*a*b

Divisin/a/b

Potencia^a^b

Matlab no toma en cuenta espacios. Si queremos que matlab evalu la lnea pero que no escriba respuesta, basta escribir punto y coma (;)al final de la sentencia.Variables:VariableDefinicinValor

AnsVariable usada por defecto para almacenar el ultimo resultado ???

PiRazn de una circunferencia a su dimetro3,1416

EpsNumero mas pequeo 2.2204e-016

InfInfinitoInf

NanNot a number magnitud no numricaNan

i y ji=j=0+1.000i

RealminEl numero real positivo mas pequeo que es utilizable2.2251e-308

RealmaxEl numero positivo mas grande que es utilizable1.7977e+308

Comandos:ComandoDefinicin

[X1 X2]=sizeAObtenemos el tamao de la fila y columna

diag(A)Se obtienen elementos de su diagonal

X3=length(A)Es equivalente a max(size(A))

X4=det(A)Se obtiene el determinante de A

A2=ones(N,M)Se obtiene una matriz de N por M de unos

A3=zeros(N,M)Se obtiene una matriz de N por M de ceros

Triu(A)Extrae la parte triangular superior

Tril(A)Extrae la parte triangular inferior

Inv(A)Se obtiene la inversa de A

ASe obtiene la transpuesta de la matriz A

Graficos:En 2D: plot, fplot, ezplot, subplot,stem,stairsEn 3D:plot3, stem3, surf, mesh, contourEjemplo:En 2D:>> t=linspace (-pi, pi, 30);>> y1=sin(2*t);>> y2=cos(5*t);>> plot(y1)

En 3D:>> t=linspace (-pi,pi,1000);y1=sin(2*t);y2=cos(5*t);plot3(t,y1,y2)

Mostrar ms de una grfica en una ventana:

>> figure>> subplot (2, 1,1);plot(t,y1)>> subplot(2,1,2);plot(y1,t)

Solucin a sistemas de ecuaciones:

p= [1 0 -8 9 -2]p = 1 0 -8 9 -2>> q= [1 0 4 -1]q = 1 0 4 -1>> r1=roots (p) (encontra las raices del polinomio p(x))r1 = -3.3028 2.0000 1.0000 0.3028>> r2=roots (q) (encontra las raices del polinomio q(x))r2 = -0.1231 + 2.0113i -0.1231 - 2.0113i 0.2463 Polyval(p,5) (encontra el valor de la ecuacin cuando x=5)ans = 468polyder(p) (encontrar la derivada del polinomio p(x))ans = 4 0 -16 9EJERCICIOS: cree las variables >> t=[2 4 6 8 10 12 14]t = 2 4 6 8 10 12 14>> y= [1 5 7; 2 5 pi]y = 1.0000 5.0000 7.0000 2.0000 5.0000 3.1416 - Extraiga la primera fila en F1 y la tercera columna en C3 de y.>> F1=y(1,:)F1 = 1 5 7>> C3=y(:,3)C3 = 7.0000 3.1416 - cree la matriz ypeq con la primera y tercera columna de y.>> a=y (:,1)a = 1 2>> ypeq=[a C3]ypeq = 1.0000 7.0000 2.0000 3.1416-borre la segunda fila de y.>> y=[1 5 7 ; 2 5 pi]y = 1.0000 5.0000 7.0000 2.0000 5.0000 3.1416m = 0 0 0 2.0000 5.0000 3.1416>> y-mans = 1 5 7 0 0 0

Cree la matriz transpuesta y grafquela. Agregue las leyendas con los textos lineal,cuadrtica,cubica y pngale por ttulo grafico A>> A= [1 2 3 4 5 6;1 4 9 16 25 36;1 8 27 64 125 216]A = 1 2 3 4 5 6 1 4 9 16 25 36 1 8 27 64 125 216>> plot(A)>> plot(A)>> title('Grafico A')>> legend('lineal','cuadratica','cubica')

-ingrese los siguientes polinomios, multiplquelos y proceda a guardarlos en la variable num.

>> a=[1 0 -8 9 -2]a = 1 0 -8 9 -2>> b=[1 0 4 -1]b = 1 0 4 -1num=conv(a,b)num = 1 0 -4 8 -34 44 -17 2Ingrese los siguientes polinomios, multiplquelos y proceda a guardarlos en la variable den.

>> u= [1 2 7 -1]u = 1 2 7 -1>> v=[1 1 -3]v = 1 1 -3>> den=conv(u,v)den = 1 3 6 0 -22 3

Siendo num el numerador y den el denominador de una funcin de transferencia, genere esta funcin y gurdela en la variable datadata=tf([num],[den])Transfer function:s^7 - 4 s^5 + 8 s^4 - 34 s^3 + 44 s^2 - 17 s + 2------------------------------------------------s^5 + 3 s^4 + 6 s^3 - 22 s + 3-ingrese la siguiente funcin de transferencia y presntala en forma de ganancia, ceros y polos .>>num1=[4 7 0 0 0 3 1];>> den1=[3 0 0 0 1 0 5 8];>> s1=tf([num1],[den1])Transfer function:4 s^6 + 7 s^5 + 3 s + 1-----------------------3 s^7 + s^3 + 5 s + 8>> s2=zpk(s1)Zero/pole/gain 1.3333 (s+1.807) (s+0.3263) (s^2 - 1.245s + 0.6454)(s^2 + 0.8618s + 0.6569)-------------------------------------------------------(s+0.9625) (s^2 - 2.2s + 1.569) (s^2 + 1.599s + 1.279)(s^2 - 0.3616s + 1.38)

Reducir las siguientes funciones de transferenciaa)

>> G1=tf([3 5 -7],[1 7 -2 1]);>> G2=tf([2 -1],[1 2 6]);>> G3=tf([4],[1 4]);>> G4=series(G1,G2)Transfer function: 6 s^3 + 7 s^2 - 19 s + 7----------------------------------------s^5 + 9 s^4 + 18 s^3 + 39 s^2 - 10 s + 6>> G5=series(G3,G4)

Transfer function: 24 s^3 + 28 s^2 - 76 s + 28-----------------------------------------------------s^6 + 13 s^5 + 54 s^4 + 111 s^3 + 146 s^2 - 34 s + 24

b)

>> G1=tf([1 2],[1 7 1]);>> G2=tf([2],[1 2 1]);>> G3=tf([1],[1 1]);>> G4=G1/(G1*G2+1)Transfer function: s^5 + 11 s^4 + 34 s^3 + 41 s^2 + 19 s + 2---------------------------------------------------s^6 + 16 s^5 + 80 s^4 + 132 s^3 + 98 s^2 + 46 s + 5>> G5=G1+G4Transfer function: 2 s^7 + 36 s^6 + 224 s^5 + 582 s^4 + 702 s^3 + 418 s^2 + 130 s + 12--------------------------------------------------------------------------s^8 + 23 s^7 + 193 s^6 + 708 s^5 + 1102 s^4 + 864 s^3 + 425 s^2 + 81 s + 5

Calcule la solucin del siguiente sistema de ecuaciones:

>> A=[1 2 3;4 2 1;7 5 9]A = 1 2 3 4 2 1 7 5 9>> B=[402;521;638] B = 402 521 638>> C=A^-1C = -0.4815 0.1111 0.1481 1.0741 0.4444 -0.4074 -0.2222 -0.3333 0.2222>>X=C*BX = -41.1481 403.4074 -121.2222