Aula-Matlab

Prévia do material em texto

>> % comentario
>> 
>> % quadricular uma funçao (coordenadas)
>> 
>> %funçao de transferencia
>> 
>> clear all
>> num= [1 2];
>> den =[1 2 3 4];
>> A= [1 2 ; 2 3];
>> sys = tf(num,den)
 
Transfer function:
 s + 2
---------------------
s^3 + 2 s^2 + 3 s + 4
 
>> sys = tf(num,den,-1)
 
Transfer function:
 z + 2
---------------------
z^3 + 2 z^2 + 3 z + 4
 
Sampling time: unspecified
>> %tempo de amostragem nao especificado
>> [num,den]= tfdata (sys)
num = 
 [1x4 double]
den = 
 [1x4 double]
>> sys.num
ans = 
 [1x4 double]
>> %Exemplo
>> Ts=1; 
>> %periodo de amostragem
>> Z = [- 6; -5; 0]
Z =
 -6
 -5
 0
>> P = [-3 +2j; -3-2j ; -1];
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.
 
>> P = [-3 +2j; -3-2j ; -1; -2];
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.
 
>> P = [-3+2j; -3-2j ; -1; -2];
>> K=1;
>> sys= zpk ([-6;-5;0], [-3+2j;-3-2j;-2;-1],1)
 
Zero/pole/gain:
 s (s+6) (s+5)
---------------------------
(s+2) (s+1) (s^2 + 6s + 13)
 
>> sys = zpk (Z,P,K,Ts);
>> [num,den,Ts]= tfdata(sys,'v')
num =
 0 1 11 30 0
den =
 1 9 33 51 26
Ts =
 1
>> % declaraçao de variaveis para usar posteriormente, variaveis sem valor
>> syms n k w z
>> ztrans(z)
 
ans =
 
w/(w - 1)^2
 
>> q=ztrans(k)
 
q =
 
z/(z - 1)^2
 
>> r= iztrans(q)
 
r =
 
n
 
>> p = ztrans (2^(n+2))
 
p =
 
(4*z)/(z - 2)
 
>> s= iztrans (p)
 
s =
 
4*2^n
 
>> %Achar a fraçao parcial
>> [r p k]= residue (num,den)
r =
 -0.6500 - 1.3000i
 -0.6500 + 1.3000i
 4.8000 
 -2.5000 
p =
 -3.0000 + 2.0000i
 -3.0000 - 2.0000i
 -2.0000 
 -1.0000 
k =
 []
>> %Modelagem em Espaço de Estados