Programação de Computadores I - Aula de 01/12 - Scilab: Arquivos e funções

Senhores alunos:

Encontram-se disponíveis neste link as notas referentes a nossa próxima aula.

Clique aqui para conhecer a próxima tarefa.

Até amanhã!

Sistemas de Controle - Scripts dos controladores

Senhores alunos:

Abaixo apresento os scripts Scilab para projeto de controladores PID pelo Método Analítico e dos controladores em avanço ou atraso de fase pelo Lugar das raízes.

Até amanhã.

Script PID:

s=%s;
ng = input("Entre com o numerador da FTMA da planta - polinômio em s: ")
dg = input("Entre com o denominador da FTMA da planta - polinômio em s:")
wg = input("Entre com a frequencia de cruzamento de ganho: ")
pm = input("Entre com a margem de fase: ")
ki = input("Entre com Ki: ")
ngc=horner(ng,%i*wg); dgc=horner(dg,%i*wg); g=ngc/dgc;
thetar=(pm-180)*3.141592/180;
ejtheta=cos(thetar)+%i*sin(thetar);
sol=(ejtheta/g)+%i*(ki/wg)
x=imag(sol);
r=real(sol);
kp=r;
kd=x/wg;
if ki~=0
    dk=s;
    nk=kd*s^2+ kp*s+ ki;
else dk=1; 
     nk=kd*s + kp;
end;
GC=syslin('c',nk,dk);

print(%io(2),GC,"Compensador");

Script Avanço ou atraso - Root Locus:

s=%s;
z = input("Entre com o amortecimento desejado : ")
wn = input("Entre com a frequencia natural desejada:")
ng = input("Entre com o numerador de G- polinomio em s: ")
dg = input("Entre com denominador de G - polinomio em s: ")
H = input("Entre com H: ")
KC=input("Entre com o ganho KC: ")
s1 = -z*wn + %i*wn*sqrt(1-z^2);
ms=abs(s1);
ts=atan(imag(s1),real(s1));
ns1=horner(ng*H,s1); ds1=horner(dg,s1); g=ns1/ds1;
mg=abs(g);
tg=atan(imag(g),real(g));
tauz=(sin(ts)-KC*mg*sin(tg-ts))/(KC*mg*ms*sin(tg));
taup=-(KC*mg*sin(ts)+sin(tg+ts))/(ms*sin(tg));
nk=KC*(tauz*s+1);
dk=taup*s+1;
GC=syslin('c',nk,dk);
print(%io(2),s1,"s1");
print(%io(2),ms,"modulo de s1");
print(%io(2),ts*180/%pi,"fase de s1");
print(%io(2),mg,"modulo de G");
print(%io(2),tg*180/%pi,"fase de G");
print(%io(2),GC,"Compensador")
print(%io(2),tauz,"tauz");
print(%io(2),taup,"tauz");
print(%io(2),GC,"Compensador");
G=syslin('c',ng,dg);
G1=G*GC;
FTMF = G1/.H;
t=linspace(0,0.005,1000);
y=csim('step',t,FTMF);
plot2d(t,2.5*y);
xgrid();

Introdução a Programação em Linguagem C - Aula de 19/11

Senhores alunos:

Cliquem no link para acessar o material da aula de amanhã.

- Arquivos.

Até a quarta.

Programação de Computadores I - Aula de 24/11 - Arquivos

Senhores alunos:

Encontram-se disponíveis neste link as notas referentes a nossa próxima aula.

Conforme tratado com os presentes no início da aula passada, nossa próxima aula inicia às 7h e 30min.
Boa leitura!

Circuitos Elétricos - FIltros Passivos: Gabarito das listas de exercícios.

Senhores alunos:

Respostas dos exercícios propostos sobre filtros passivos do livro texto:

14.47 - fc = 796 kHz

14.48 -                                    1
            |G| = ------------------------------------------------
                     sqrt [ (1 - w^2 * LC) ^2   + (wL/R)^2 ]

           Se w = 0 --> G = 1
           Se  w = 1 / sqrt(LC) --> 0 <G < 1
           Se w = infinito --> G = 0

           Logo o Filtro é Passa Baixas

14.49 - 0,2 rad/s , -14 dB e -84,3 graus

14.50 - Filtro PA e fc = 318,2 Hz

14.51 - R = 1200 ohms (Valor comercial)

14.52 - R = 27k (valor comercial)

14.53 - R = 18k , L = 2,9 H e Q = 10,5

14.54 - Considerando a configuração RLC série:
             B = 0,5 rad/s
             Com C= 10 mF --> L = 1H
             R = 0.5 ohms

14.55 - 1,56 kHz< f < 1,62 kHz , Q = 25

14.56 - Considerar na demonstração a configuração RLC serie

14.57-  (a) wo = 1 rad/s e B = 3 rad/s
           (b) wo = 1 rad/s e B = 3 rad/s   

14.58 -  fo = 79,6 kHz
            fci = 78 kHz
            fcs = 81,2 kHz
            Q = 250

14.59 - B - 2400 rad/s e wo = 15811 rad/s

Respostas dos exercícios propostos sobre filtros passivos - 6 questões:

1) Para R1 = 100k:
 C1 = 10 nF, C2 = 2nF, , R2 = 10k,

2) Q = 0,29

3) f0 = 2 kHz
    Q = 0,67
Para L = 10mH --> C= 0,68uF e R =180 (Comerciais)

4) Correção - Ruido de 60 Hz!

   FPA : R = 3,9 k
           C = 68 nF
  Redução - 1,32dB

5) FPB: fc = 7860,9 Hz
           C = 1,2 nF
           R = 16,87 k -> 18k

6) Vo =                            1                                     
    Vi     ( 1 + R2 + wc1 )   +   j ( w  -   wc1 )
                    R1    wc2             wc2      w    

Programação de Computadores I - Aula de 17/11

Senhores alunos:

Encontram-se disponíveis neste link as notas referentes a nossa próxima aula.

Boa leitura!

Introdução a Programação I - Aula de 12/11

Senhores alunos:

No link abaixo, é possível obter as notas de nossa aula de hoje.

Estruturas.

Até a noite.

Sistema de Controle - Gabarito da Segunda Parte da Primeira Avaliação

Senhores Alunos

O gabarito das questões teóricas encontra-se abaixo. O Gabarito das questões práticas neste link.

1) Tempo de subida é o tempo que a resposta a uma entrada em degrau leva para excursionar de 10 a 90% do seu valor final.

2) Quando o amortecimento é razoavelmente maior do que 1.

3) a) O denominador apresenta o termo s^2.
    b) Degrau e rampa.

Boa Semana

Prog. Computadores I - Solução da Tarefa I no Scilab

Senhores alunos:

Em função das tarefas entregues não estarem bem adaptadas ao ambiente Scilab, apresento, no link abaixo, a solução da tarefa.

Script Scilab - Tarefa I.

Boa semana!

Programação de Computadores I - Aula de 10/11

Senhores alunos:

Encontram-se no link abaixo, as notas de nossa próximas aula.

Funções recursivas.

Bom início de de semana!