quarta-feira, 24 de setembro de 2014

Introdução a Programação C - Aula de 24/09

Senhores alunos:


Segue o link para acessar o material da nossa aula de hoje:
  • Classes de Armazenamento;
  • O Pré-processador C e as macros.

Até a noite!

terça-feira, 23 de setembro de 2014

Programação de Computadores I - Guia de referência rápida para o Scilab

Senhores alunos:

Estou disponibilizando uma primeira versão do "Cartão de referência do Scilab".

Peço que leiam e verifiquem possíveis erros.

Boa semana!


Keywords: Scilab, reference card.

sexta-feira, 19 de setembro de 2014

Programação de Computadores I - Aula de 22/09

Senhores alunos:


As notas relativas a nossa próxima aula podem ser lidas no link abaixo.

Incluo neste post exercícios que podem ser utilizados em uma avaliação escrita.

Até a próxima postagem!

Keywords: select, case, C programming language, Fortran, Scilab.

segunda-feira, 15 de setembro de 2014

Lista de exercícios sobre transitórios de segunda ordem

Senhores Alunos:



Disponibilizo, a lista de exercícios para preparação para a avaliação dos transitórios de segunda ordem. Apenas  o primeiro exercício não foi selecionado do livro texto.

A propósito, no primeiro exercício, o primeiro gráfico representa a corrente no indutor em um circuito RL série, medida em ampères. O segundo gráfico representa a tensão no capacitor de um circuito RLC série, medida em volts.

Keywords:  transient response, second-order, parallel RLC circuit, series RLC circuits, exercises.

Respostas:

1. R = 5 Ohms
    L = 0,32 H
    C = 3,7 mF
    Valores aproximados.

8.12 a) C > 6 mF  b) C = 6 mF c) C < 6 mF
8.13 R > 120 Ohms
8.14 v = -1,33 (e^-0,5t - e^-0,2t)
8.15 L = 25H, C= 0,2 mF e R = 750 Ohms
8.22 C = 0,05 uF
8.23 C = 40 mF
8.24 i = e^-10t ( 4 cos 19,36 t + 1,05 sen 19,36 t)
8.30 C = 150 nF
8.31 VL(0+) = 80
8.36 v = 35 e^-0,8t (15 sen 0,6t + 20 sen 0,6 t)
        i = 5 e^-0,8t sen 0,6t
8.43 R = 8 ohms e C= 2,075 mF
8.47 v = 200 t e^-10t
8.53 i'' + (5/4) i' + 100 i = 600

Transitórios RLC de segunda ordem - Script com Scilab

Senhores Alunos:

Segue, abaixo, o script Scilab completo para o cálculo das constantes da solução das equações diferenciais e o traçado dos gráfico da corrente e da tensão dos transitórios RLC série e paralelo de segunda ordem.


Bom uso!

Keywords: transient response, second-order, parallel RLC circuit, series RLC circuit, Scilab.

-----------------------------------------------------------------------

clc;
clear ;
xdel(winsid());
disp('Avaliação da Corrente no transitorio RLC de 2a ordem');
Tipo = input('Informe <1> para circuito serie <2> para paralelo: ')
disp('Entrada de Dados');
R = input('Valor de R =');
L = input('Valor de L =');
C = input('Valor de C =');
I0 = input('Valor de I0 =');
V0 = input('Valor de V0 =');
w0 = 1/(sqrt(L*C));
if Tipo==1
    Vs = input('Valor de Vs =');
    dI0 = -(-Vs+R*I0+V0)/L;
    dV0 = I0/C
    alpha = R/(2*L);
else 
    Is = input('Valor de Is =');
    dV0 = -(-Is+I0+V0/R)/C;
    dI0 = V0/L;
    alpha = 1/(2*R*C);
end     
printf('di(0)/dt = %.2f, alpha = %.2f, w0 = %.2f\n',dI0,alpha,w0); 
printf('dv(0)/dt = %.2f\n',dV0); 
if Tipo==1
    if alpha<w0
        disp('Caso subamortecido');
        wd = sqrt(w0^2-alpha^2);
        A1 = I0;
        A2 = (dI0 + alpha*A1)/wd;
        A3 = V0 - Vs;
        A4 = (dV0 + alpha*A3)/wd;
        printf('s1 = %.2f + j%.2f, s2 = %.2f - j%.2f\n',-alpha,wd,-alpha,wd);
        printf('wd = %.2f\n',wd);
        t = 0:0.001:5/alpha;
        i = exp(-alpha*t).*(A1*cos(wd*t) + A2*sin(wd*t));
        v = exp(-alpha*t).*(A3*cos(wd*t) + A4*sin(wd*t)) + Vs;
    elseif alpha>w0
        disp('Caso superamortecido');
        delta = sqrt(alpha^2-w0^2);
        s1 = -alpha + delta;
        s2 = -alpha - delta;
        A2 = (s1*I0 -dI0)/(2*delta);
        A1 = I0 - A2;
        A4 = (s1*(V0-Vs)-dV0)/(2*delta);
        A3 = V0 - Vs - A4;
        printf('s1 = %.2f, s2 = %.2f\n',s1,s2);
        tau = (1/s1)*(A1/A2); // pondera em relação as constantes
        t = 0:0.001:6*abs(tau);
        i = A1*exp(s1*t) + A2*exp(s2*t);
        v = A3*exp(s1*t) + A4*exp(s2*t) + Vs;
    else
        disp('Caso de amortecimento crítico');
        s = -alpha;
        A1 = I0;
        A2 = dI0 + alpha*A1;
        A3 = V0 - Vs;
        A4 = dV0 + alpha*A3;
        printf('s = %.2f\n',s);
        t = 0:0.001:5/alpha;
        i = (A1 + A2*t).*exp(s*t);
        v = (A3 + A4*t).*exp(s*t) + Vs;
    end
else  
    if alpha<w0
        disp('Caso subamortecido');
        wd = sqrt(w0^2-alpha^2);
        A3 = V0;
        A4 = (dV0 + alpha*A3)/wd;
        A1 = I0 - Is;
        A2 = (dI0 + alpha*A1)/wd;
        printf('s1 = %.2f + j%.2f, s2 = %.2f - j%.2f\n',-alpha,wd,-alpha,wd);
        printf('wd = %.2f\n',wd);
        t = 0:0.001:5/alpha;
        v = exp(-alpha*t).*(A3*cos(wd*t) + A4*sin(wd*t));
        i = exp(-alpha*t).*(A1*cos(wd*t) + A2*sin(wd*t)) + Is;
    elseif alpha>w0
        disp('Caso superamortecido');
        delta = sqrt(alpha^2-w0^2);
        s1 = -alpha + delta;
        s2 = -alpha - delta;
        A4 = (s1*V0 -dV0)/(2*delta);
        A3 = V0 - A4;
        A2 = (s1*(I0-Is)-dI0)/(2*delta);
        A1 = I0 - Is - A2;
        printf('s1 = %.2f, s2 = %.2f\n',s1,s2);
        tau = (1/s1)*(A1/A2);// pondera em relação as constantes
        t = 0:0.001:6*abs(tau);
        v = A3*exp(s1*t) + A4*exp(s2*t);
        i = A1*exp(s1*t) + A2*exp(s2*t) + Is;
    else
        disp('Caso de amortecimento crítico');
        s = -alpha;
        A3 = V0;
        A4 = dV0 + alpha*A3;
        A1 = I0 - Is;
        A2 = dI0 + alpha*A1;
        printf('s = %.2f\n',s);
        t = 0:0.001:5/alpha;
        v = (A3 + A4*t).*exp(s*t);
        i = (A1 + A2*t).*exp(s*t) + Is;
    end
end 
printf('Constantes ED corrente, A1 = %.2f, A2 = %.2f\n',A1,A2);
printf('Constantes ED tensao, A1 = %.2f, A2 = %.2f\n',A3,A4);
subplot(211);
plot(t,i)
xlabel('Tempo (s)');
ylabel('Corrente (A)');
xgrid;
subplot(212);
plot(t,v)
xlabel('Tempo (s)');
ylabel('Tensão (A)');
xgrid;

Gabarito exercícios de controle

Senhores alunos:

Neste link é possível baixar o gabarito da lista de exercícios disponibilizada para preparação para a primeira avaliação.

Até amanhã.




quinta-feira, 11 de setembro de 2014

quarta-feira, 10 de setembro de 2014

Programação de Computadores I - Aula de 15/09

Senhores alunos:


As notas relativas a nossa próxima aula podem ser lidas no link abaixo.

Estruturas de controle de fluxo por laço - parte ii.

Até a próxima postagem!

Keywords: flow control structures, C programming language, Fortran, Scilab.

quinta-feira, 4 de setembro de 2014

Programação de Computadores I - Aula de 08/09

Senhores alunos:


As notas relativas a nossa próxima aula podem ser lidas no link abaixo.

Estruturas de controle de fluxo por laço.

Até a próxima postagem!

Keywords: flow control structures, C programming language, Fortran, Scilab.

quarta-feira, 3 de setembro de 2014