TRABALHO 02 Termodin mica
ACADÊMICOS:
Iago Braga Theodoro – RA 139209
Rayssa Gomes Guimarães – RA 139289
Roger Carvalho Kazama – RA 149392
Trabalho 02
Criar uma rede neural para dados tabulados de termodinâmica, achando suas entalpias e entropias equivalentes.
Tabela 01 – Dados termodinâmicos para a Água
Para a obtenção do código, as entradas foram definidas como sendo dados de temperatura e pressão, relacionando-as às saídas solicitadas (entalpias de líquido saturado, durante evaporação e vapor saturado e entropias de liquido e vapor saturado). Desta forma, obteve-se o seguinte código:
NN = input('Digite o numero de neuronios:') format long entrada = [ 50 0.1235 ; 75 0.3858; 100 1.014 ; 150 4.758; 200 15.54; 220 23.18; 250 39.73; 280 64.12; 300 85.81; 340 145.9]'; saida = [209.33 2382.7 2592.1 0.7038 8.0763; 313.93 2321.4 2635.3 1.0155 7.6824; 419.04 2257.0 2676.1 1.3069 7.3549; 632.20 2114.3 2746.5 1.8418 6.8379 ; 852.45 1940.7 2793.2 2.3309 6.4323; 943.62 1858.5 2802.1 2.5178 6.2861;1085.4 1716.2 2801.5 2.7927 6.0730;1236.0 1543.6 2779.6 3.0668 5.8571;1344.0 1404.9 2749.0 3.2534 5.7045;1594.2 1027.9 2622.0 3.6594 5.3357]'; net = newfit(entrada,saida,NN); net = train(net,entrada,saida);
T = input('Digite a Temperatura:');
P = input('Digite a Pressao:');
Y = [T P]'; fprintf(' Hliq\n Hevap\n Hvap\n Sliq\n Svap\n')
Tabela = sim(net,Y)
Assim, para uma rede neural com apenas 1 neurônio, obteve-se as seguintes saídas:
Digite o numero de neuronios:1
NN = 1
Digite a Temperatura:110
Digite a Pressao:1.433 Hliq Hevap Hvap Sliq Svap
Tabela = 1.0e+03 * 0.454570972513818 2.218287422999792 2.672905240622185 0.001368250329394 0.007371968717726
Analisando a tabela temos a seguinte relação entre resultado obtido pela rede neural e o especificado na tabela-base:
Variável
Rede Neural
Tabela-base
|Erro| hliq 454.57
461.30
6.73 hevap 2218.29
2230.2
11.91 hvap 2672.91
2691.5
18.59 sliq 1.3683
1.4185
0.0502 svap 7.372