Ajuste de modelos a dados experimentais - Determinação de parâmetros скачать в хорошем качестве

Ajuste de modelos a dados experimentais - Determinação de parâmetros

Acompanhe também conteúdos em:   / professormarcosfpmoreira     / cristianismo_na_industria   https://professormarcosfpm.wixsite.co...   / professormarcosfpmoreira      / @caminhando_com_jesus_e_maria   O Prof. Marcos Moreira possui graduação, mestrado, doutorado e 2 anos de pós-doutorado em Engenharia Química. É professor Associado do Departamento de Engenharia Química da Universidade Estadual do Oeste do Paraná (Unioeste, Toledo-PR, 24°43'29,77''S/53°45'11,84''O). http://buscatextual.cnpq.br/buscatext... É autor dos livros "Scilab aplicado à Engenharia Química (para graduação e pós-graduação)"; "Operações Unitárias da Engenharia Química utilizando o Excel/VBA"; "Balanços Globais de Massa na Engenharia Química (Uma abordagem atualizada e criativa)" e Matemática Aplicada à Engenharia Química (com Scilab). Para mais informações acesse: https://www.e-papers.com.br/product_a... Para entrar em contato escreva para professormarcosfpmoreira@gmail.com //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% //PROGRAMA 1 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc clear global x_exp y_exp alfa y0 //***************** // dados de entrada //***************** y0=0.75; x0=0; x_exp=[x0 5.2 16.2 26.8 38.1 50.2 63.5 80.6 99.7]; x_exp=x_exp'; y_exp=[y0 0.72 0.65 0.59 0.52 0.47 0.39 0.32 0.24]; y_exp=y_exp'; d=0.008; D=0.15; A_sai=%pi*d^2/4; A=%pi*D^2/4; g=9.81; //******** // modelo //******** function y_teorico=y_cal(x,C) y_teorico=(y0^0.5-(A_sai/(2*A))*C*(2*g)^0.5*x)^2; endfunction //********************** // função de otimização //********************** function F_otim=minha_funcao(C, x_exp, y_exp) F_otim=y_cal(x_exp,C)-y_exp endfunction //********************* // Definições de busca //********************* C_inicial=[1]; //************************* // Encontrando o parâmetro //************************* [fopt,C_otimo, gopt] = leastsq(list(minha_funcao,x_exp, y_exp),C_inicial) //******************** // imprimindo valores //******************** correlacao=correl(y_cal(x_exp,C_otimo),y_exp); printf('Coeficiente de descarga = %.2f
',C_otimo); printf('R^2 = %.6f
',correlacao^2); //********* // Gráfico //********* plot(x_exp,y_cal(x_exp,C_otimo),'b-',x_exp,y_exp,'ro'); plot(x_exp,y_cal(x_exp,1),'k-'); legend(["calculado", "experimental","C=1"],1); xtitle("AJUSTE MODELO","tempo (s)","Nível (m)"); //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% //PROGRAMA 2 //%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc clear global x_exp y_exp x0 y0 //***************** // dados de entrada //***************** y0=0.75; x0=0; x_exp=[x0 5.2 16.2 26.8 38.1 50.2 63.5 80.6 99.7]; x_exp=x_exp'; y_exp=[y0 0.72 0.65 0.59 0.52 0.47 0.39 0.32 0.24]; y_exp=y_exp'; d=0.008; D=0.15; A_sai=%pi*d^2/4; A=%pi*D^2/4; g=9.81; //************************** // Equação diferencial ***** //************************** function dy_dx=F(x,y); global y_exp dy_dx=-C*(A_sai/A)*(2*g)^0.5*y^0.5; endfunction //********************************** // resolvendo a equação diferencial //********************************** function y_calculados=y_calc(x,C) global x0 y0 tudo_junto=ode(y0,x0,x,F) y_calculados=[]; c=1; for i = 1:1:length(x_exp) y_calculados(c,1)= tudo_junto(1,i); c=c+1; end endfunction //******************************* // Função de otimização //******************************* function F_otim=minha_funcao(C_otimo, x_exp, y_exp) global y_exp x_exp F_otim=abs(y_calc(x_exp,C_otimo)-y_exp) endfunction //%***************** //iniciando a busca //%***************** C_inicial=[1]; // chute inicial [fopt,C_otimo, gopt] = leastsq(list(minha_funcao,x_exp, y_exp),C_inicial) //******************** // imprimindo valores //******************** printf('Coeficiente de descarga = %.2f
',C_otimo); correlacao=correl(y_calc(x_exp,C_otimo),y_exp); printf('R^2 = %.6f
',correlacao^2); //********* // Gráfico //********* plot(x_exp,y_calc(x_exp,C_otimo),'b-',x_exp,y_exp,'ro'); plot(x_exp,y_calc(x_exp,1),'k-'); legend(["calculado", "experimental","C=1"],1); xtitle("AJUSTE EQUAÇÃO DIFERENCIAL","tempo (s)","Nível (m)");