Prévia do material em texto
Ajuste de curvas Método dos mínimos quadrados Jorge Díaz Calle Noção • Dados os valores experimentais de grandezas (considerando que contêm erros) pretende-se encontrar uma curva que “represente” os dados fornecidos e outros relacionados, de tal forma que seja minimizado o erro total cometido no sentido da média quadrática. Tabela com valores experimentais 0.11 2.9 8.7 0.7 5.6 8.4 5.3 0.3 2.2 0.1 Tensão 0.2 8.1 6.1 4.1 2.1 0.1 8.0 6.0 4.0 2.0 eIntensidad 10 9 8 7 6 5 4 3 2 1 ItemApós realizar 10 medições com um tensiômetro para diferentes intensidades de corrente, obtemos a Tabela de dados: Se geramos um gráfico dos dados, observa-se: Aparenta ser um comportamento linear com algumas oscilações. Tentando desenhar uma reta que unindo os dados: O conjunto de segmentos de reta, podem ser ajusta- dos com uma reta, pois os “picos” podem ser me- dições com erro. Ajustando os dados com uma reta • Observamos os dados. Então, a tensão será ajustada com uma reta cuja equação é: • Então os dados satisfazem o sistema: bxay y x Tensão eIntensidad ba ba ba ba bxay ii )0.2(0.11 )6.0(0.3 )4.0(2.2 )2.0(0.1 O sistema matricial é 110210 0.11 2.9 8.7 0.7 5.6 8.4 5.3 0.3 2.2 0.1 10.2 18.1 16.1 14.1 12.1 10.1 18.0 16.0 14.0 12.0 )0.2(0.11 )8.1(2.9 )6.1(8.7 )4.1(0.7 )2.1(5.6 )0.1(8.4 )8.0(5.3 )6.0(0.3 )4.0(2.2 )2.0(0.1 xx ii B b a A b a ba ba ba ba ba ba ba ba ba ba bxay O produto com a transposta • Lembrar que o produto com a transposta pode projetar em uma menor ou maior dimensão do sistema: ou • O primeiro caso permite reduzir a um sistema de duas equações com duas incógnitas, assim: 22210102 xx t x PAA 1010102210 x t xx PAA Sistema projetado a dimensão 2 0.11 2.9 8.7 0.7 5.6 8.4 5.3 0.3 2.2 0.1 1111 28.14.02.0 10.2 18.1 16.1 14.1 12.1 10.1 18.0 16.0 14.0 12.0 1111 28.14.02.0 b a 210102210102 x t xx t x BAXAA Tabela melhorada .......56.11 0.22 88.0 2.0 00.4 24.3 56.2 96.1 44.1 00.1 64.0 36.0 16.0 04.0 0.11 2.9 8.7 0.7 5.6 8.4 5.3 0.3 2.2 0.1 0.2 8.1 6.1 4.1 2.1 0.1 8.0 6.0 4.0 2.0 10 9 8 7 6 5 4 3 2 1 10 1 2 i iiiii yxxyxi 0.56 12.79 100.11 0.114.15 b a Resolvemos o sistema por Gauss Jordan 6.511.1 52.1703.3 6.511.1 12.790.114.15 25 6 10 55 292 01 6.511.1 55 292 01 25 6 55 292 xy 0.561011 12.790.114.15 :é estendidamatrix A Podemos projetar outros valores Sistema geral projetado a dimensão 2 n n nn n n nn y y y y y y y xxxx b a x x x x x x x xxxx 1 5 4 3 2 1 121 1 5 4 3 2 1 121 1111 1 1 1 1 1 1 1 1111 210102210102 x t xx t x BAXAA O sistema em dimensão 2 • Observe o sistema obtido • Isso significa que podemos resolver se completamos a tabela inicial com os produtos necessários e seus somatórios. n i i n i ii n i n i i n i i n i i y yx b a x xx 1 1 11 11 2 1 Caso um comportamento como: • Observar os dados e ao levar a um desenho: 13/1 7/1 25.0 4.0 0.1 5.2 0.5 0.3 0.2 5.1 0.1 8.0 6 5 4 3 2 1 ii yxitem Quais curvas posso utilizar? • Reta: • Parábola: • Polinômio: • Hipérbola: • Função exponencial: bxay * cxbxay ** 2 01 1 1 ... axaxaxay n n n n axbzaxCyeCyCey axax lnlnlnlnln bax y 1 Ajuste de Dados em Hipérbole • Uma hipérbole é representada por • Criamos uma nova variável • Assim, pensamos em um ajuste similar ao de uma reta, corrigindo a tabela. bax y 1 . 1 bax y z 0.13 0.7 0.4 5.2 0.1 4.0 13/1 7/1 25.0 4.0 0.1 5.2 0.5 0.3 0.2 5.1 0.1 8.0 6 5 4 3 2 1 item iii zyx Calculando a hipérbole 9.27 07.99 63.13 3.1389.41 b a 07.99 0.65 0.21 0.8 75.3 0.1 32.0 89.41 0.25 0.9 0.4 25.2 0.1 64.0 9.27 0.13 0.7 0.4 5.2 0.1 4.0 3.13 0.5 0.3 0.2 5.1 0.1 8.0 6 5 4 3 2 1 item 2 iiiii zxxzx .2.)3( xz . 23 1 x y Comparando 23 1 x y 53.0 19 1 2)7(3 1 )7( temos0.7 Para yx Ajuste quadrático • Para ajustar com um polinômio quadrático: • Exemplo: cxbxay 2 5.4 0.1 9.0 0.1 6.1 10.4 50.0 60.0 00.1 00.2 472.22 00.16 0625.5 00.1 4096.0 887.12 00.8 375.3 00.1 512.0 89.7 00.4 25.2 00.1 64.0 15.4 25.0 40.0 00.1 50.2 3.5 0.2 5.1 0.1 8.0 4 3 2 1 Item 2432 iiiiiiiii yxyxxxxyx Forma geral do sistema matricial 13 3 2 1 2 2 3 2 3 2 2 2 1 2 1 2 2 3 2 33 2 2 22 1 2 11 2 1 1 1 1 1 nxnx n i nn ii nnn iii iii B c b a A y y y y y c b a xx xx xx xx xx cxbxay cxbxay cxbxay cxbxay cxbxay cxbxay n i ni ni nn ii ni ni y y y y y xxxxx xxxxx c b a xx xx xx xx xx xxxxx xxxxx 3 2 1 321 222 3 2 2 2 1 2 2 3 2 3 2 2 2 1 2 1 321 222 3 2 2 2 1 11111 1 1 1 1 1 11111 n i i n i ii n i ii n i i n i i n i i n i i n i i n i i n i i n i i y yx yx c b a nxx xxx xxx 1 1 1 2 11 2 11 2 1 3 1 2 1 3 1 4 Para o exemplo 15.4 1.4 5.4 43.589.7 3.589.7887.12 89.7887.12472.22 c b a 924.7254.9725.2 2 xxy Para o caso de ajuste exponencial • Uma função exponencial pode ser expressada • Para ajustar os dados a uma curva exponencial observamos que podemos criar uma nova variável para transformar o ajuste em reta. axeCy axCeCeCyz axax lnlnlnlnln Baxz axBzCB é ajuste o portanto , : temos,ln Fazendo Dados para o ajuste • Dados iniciais 0.20 0.8 0.4 5.2 0.2 0.1 0.5 0.3 0.2 5.1 0.1 0.0 6 5 4 3 2 1 item ii yx 916.05.2ln 693.02ln 995.2 079.2 386.1 916.0 693.0 0.0 0.20 0.8 0.4 5.2 0.2 0.1 0.5 0.3 0.2 5.1 0.1 0.0 6 5 4 3 2 1 item iii zyx Ajuste em reta 916.05.2ln 693.02ln 051.2625.41069.85.12 995.2 079.2 386.1 916.0 693.0 000.0 000.5 000.3 000.2 500.1 000.1000.0 995.2 079.2 386.1 916.0 693.0 000.0 000.5 000.3 000.2 500.1 000.1 000.0 6 5 4 3 2 1 item 2 iiiii zxxzx n i i n i ii n i i n i i n i i z zx B a nx xx 1 1 1 11 2 069.8 051.26 65.12 5.1225.41 :resolvemos Logo, B a Ajuste exponencial • Resolvendo: • Substituindo: 079.0ln079.0e6076.0 eCCBa x x xax ey ey eeeCy 6076.0 6076.0079.0 6076.0079.0 0822.1 ou Aplicação • Uma pesquisadora em química, após um processo térmico sobre um alimento, tabelou a presença de uma componente em quatro tempos diferentes. Após uma análise dos dados conjetura que uma curva quadrática ajusta bem os mesmos. Determine a curva quadrática que melhor ajusta os dados: 23.3 08.0 25.0 40.0 50.2 3.9 0.5 0.2 5.1 8.0 4 3 2 1 Item ii yt Aplicação • Utilizando o ajuste realizado, a pesquisadora informa a presença da componente na terceira e quarta hora do processo térmico em uma tabela modificada. Quais os valores informados? É um bom ajuste? Existe outro ajuste melhor? 50.5 00.2 00.1 90.0 60.1 50.3 40.0 50.0 60.0 00.2 4721.646 000.625 000.16 0625.5 4096.0 887.136 000.125 000.8 375.3 512.0 89.31 00.25 00.4 25.2 64.0 23.3 08.0 25.0 4.0 5.2 3.9 0.5 0.2 5.1 8.0 4 3 2 1 Item 2432 iiiiiiiii ytyttttyt Ajuste quadrático • Ajustando com a curva: • Resolvemos 23.3 5.3 5.5 43.989.31 3.989.31887.136 89.31887.1364721.646 1 1 1 2 11 2 11 2 1 3 1 2 1 3 1 4 c b a y yt yt c b a ntt ttt ttt n i i n i ii n i ii n i i n i i n i i n i i n i i n i i n i i n i i cxbxay 2 Ajuste quadrático • A curva quadrática é • Logo os valores informados (Não físico – Errado) 7211.4336.3482.0 2 xxy 9109.0)0.4( 9489.0)0.3( y y Melhor ajuste seria hipérbole • Sugere-se uma hipérbole • Completando a tabela 89.31 00.25 00.4 25.2 64.0 57.74 50.62 00.8 75.3 32.0 4.19 5.12 00.4 50.2 40.0 08.0 25.0 40.0 50.2 3.9 0.5 0.2 5.1 8.0 4 3 2 1 Item 2 iiiiii tztzyt nmt y 1 . 1 nmt y z Hipérbole de ajuste • Resolvendo temos 82.187.2 1 t y 10352.0)0.4( 14728.0)0.3( y y