Cálculo Numérico

Prévia do material em texto

UNIVERSIDADE DE CAXIAS DO SUL
CENTRO DE CIÊNCIAS EXATAS E DA TECNOLOGIA
ENGENHARIA QUÍMICA
DISCIPLINA: CÁLCULO NUMÉRICO
PROFESSORA: VÂNIA MARIA PINHEIRO SLAVIEIRO
ALUNA: ELIANE ANTUNES ANDRADE
TRABALHO 02
Número da chamada: 09
Exercícios: 01, 03, 05 e 09
Questão 01.
>> A=[1 2 0 4;0 1 1 2;0 0 1 4;1 2 3 4];
>>b=[5040; 2160; 1520; 7200];
>>[x] = SLGauss(A,b)
A = 1 2 0 4
 0 1 1 2
 0 0 1 4
 1 2 3 4
b = 5040
 2160
 1520
 7200
x = 2160
 1040
 720
 200
Resposta: O tempo de vida útil das ferramentas consta abaixo.
 F1= 2160 horas de vida útil;
 F2= 1040 horas de vida útil;
 F3= 720 horas de vida útil;
 F4= 200 horas de vida útil.
 
Questão 03.
>> %demanda mínima
>> x=[2 3 4 5]';
>> y=[16.4 15.2 14.9 16.0]';
>> [c,v] = IVander(x,y,2)
c = 0.0833
 -0.3000
 -1.2833
 19.5000
v = 16.4000
>> t=3.5:0.01:4;
>> min(VPol(c,t))
ans = 14.8632
Gráfico:
>> %demanda máxima
>> x=[13 14 15 16 17]';
>> y=[28.0 36.5 43.0 34.0 31.2]';
>> [c,v] = IVander(x,y,13)
c = 1.0e+04 *
 0.0001
 -0.0087
 0.1940
 -1.9064
 6.9942
v = 28.0000
>> t=14.5:0.01:15;
>> max(VPol(c,t))
ans = 43.1006
Gráfico:
Resposta: A demanda máxima é 43,1006 MW e a demanda mínima é 14,8632 MW. 
Questão 05.
>> x=[20 40]';
>>y=[15e-6 16.98e-6]';
>> [C,v]=ISpline3(x,y,25)
C = 1.0e-04 *
 0 0 0.0010 0.1500
v = 1.5495e-05
>> Re=1e+08; %Número de Reynolds
>>V=50; %Velocidade do ar
>>v=v; %Viscosidade encontrada no spline3
>> %x é comprimento da placa a ser encontrado
>> x=(Re*v)/V
x = 30.9900
Resposta: O comprimento da placa é 30,99 metros.
Questão 09.
>> x=[2:2:30]';
>>y=[9.7;8.1;6.6;5.1;4.4;3.7;2.8;2.4;2;1.6;1.4;1.1;0.85;0.69;0.6];
>> [a,b,B,f,Sqe]=AjusteExp(x,y)
a = 11.9131
b = -0.1002
B = -0.1002
 2.4776
f = 9.7505
 7.9804
 6.5317
 5.3460
 4.3755
 3.5812
 2.9311
 2.3990
 1.9635
 1.6071
 1.3153
 1.0766
 0.8811
 0.7212
 0.5903
Sqe = 0.1251
>> b1=(-1)/(b*(5e6))
b1 = 1.9968e-06
Resposta: A capacitância é 1,9968µF.