Fórum de Matemática
DÚVIDAS? Nós respondemos!

Um Fórum em Português dedicado à Matemática
Data/Hora: 12 jul 2025, 15:48

Os Horários são TMG [ DST ]




Fazer Nova Pergunta Responder a este Tópico  [ 2 mensagens ] 
Autor Mensagem
 Título da Pergunta: Minimo quadrado polinomial
MensagemEnviado: 14 jan 2014, 16:55 
Offline

Registado: 14 jan 2014, 16:25
Mensagens: 1
Localização: Brazil
Agradeceu: 0 vez(es)
Foi agradecido: 0 vez(es)
Eu estou tentando encontrar o "a" e "b" que minimizam essa equação

yhat(k+1) = ay(k) + bu(k)

Usando o método dos mínimos quadrados tem-se essa equação

\(\sum_{k=1}^{N} \left [ y(k+1) - ay(k) -bu(k) \right ]^2\)

Para encontrar os valores de a e b tem-se as seguintes equações

\(\frac{\partial }{\partial a}J(a,b) = -2\sum_{k=1}^{N} y(k)\left [ y(k+1) - ay(k) -bu(k) \right ] = 0\)

\(\frac{\partial }{\partial b}J(a,b) = -2\sum_{k=1}^{N} u(k)\left [ y(k+1) - ay(k) -bu(k) \right ] = 0\)

Eu sei que a resposta correta é essa de baixo, mas não estou conseguindo manipular as equações para chegar a isso

\(a = \frac{\sum_{k=1}^{N}u^2(k) \sum_{k=1}^{N}y(k)y(k+1) - \sum_{k=1}^{N}u(k)y(k) \sum_{k=1}^{N}u(k)y(k+1)}{\sum_{k=1}^{N}y^2(k) \sum_{k=1}^{N}u^2(k) - (\sum_{k=1}^{N}u(k)y(k))^2}\)

\(b = \frac{\sum_{k=1}^{N}y^2(k) \sum_{k=1}^{N}u(k)y(k+1) - \sum_{k=1}^{N}u(k)y(k) \sum_{k=1}^{N}y(k)y(k+1)}{\sum_{k=1}^{N}y^2(k) \sum_{k=1}^{N}u^2(k) - (\sum_{k=1}^{N}u(k)y(k))^2}\)


Topo
 Perfil  
 
 Título da Pergunta: Re: Minimo quadrado polinomial
MensagemEnviado: 14 jan 2014, 18:57 
Offline

Registado: 05 jan 2011, 12:35
Mensagens: 2235
Localização: Lisboa
Agradeceu: 683 vezes
Foi agradecido: 346 vezes
partindo do princípio que tem os seus cálculos corretos pode igualar os somatórios

\(-2\sum_{k=1}^{N} y(k)\left [ y(k+1) - ay(k) -bu(k) \right ]=-2\sum_{k=1}^{N} u(k)\left [ y(k+1) - ay(k) -bu(k) \right ]\)

\(\sum_{k=1}^{N} y(k)\left [ y(k+1) - ay(k) -bu(k) \right ]=\sum_{k=1}^{N} u(k)\left [ y(k+1) - ay(k) -bu(k) \right ]\)

\(\sum_{k=1}^{N} y(k)\left [ y(k+1) - ay(k) -bu(k) \right ]-\sum_{k=1}^{N} u(k)\left [ y(k+1) - ay(k) -bu(k) \right ]=0\)

\(\sum_{k=1}^{N} \left(y(k)\left [ y(k+1) - ay(k) -bu(k) \right ]-u(k)\left [ y(k+1) - ay(k) -bu(k) \right ]\right)=0\)

\(\sum_{k=1}^{N} \left(y(k) y(k+1) - ay^2(k) -b.y(k)u(k) -u(k) y(k+1) + a .u(k)y(k) -b.u^2(k) \right)=0\)

\(\sum_{k=1}^{N} \left(y(k) y(k+1) -b.y(k)u(k) -u(k) y(k+1) -b.u^2(k) + a(u(k)y(k) - y^2(k)) \right)=0\)

\(\sum_{k=1}^{N} \left(y(k) y(k+1) -b.y(k)u(k) -u(k) y(k+1) -b.u^2(k)\right) + a\sum_{k=1}^{N}(u(k)y(k) - y^2(k)) =0\)

avance, agora é fácil colocar o \(a\) em evidência

para o \(b\) é semelhante

_________________
João Pimentel Ferreira
 
Partilhe dúvidas e resultados, ajude a comunidade com a sua pergunta!
Não lhe dês o peixe, ensina-o a pescar (provérbio chinês)
Fortalecemos a quem ajudamos pouco, mas prejudicamos se ajudarmos muito (pensamento budista)


Topo
 Perfil  
 
Mostrar mensagens anteriores:  Ordenar por  
Fazer Nova Pergunta Responder a este Tópico  [ 2 mensagens ] 

Os Horários são TMG [ DST ]


Quem está ligado:

Utilizadores a ver este Fórum: Nenhum utilizador registado e 43 visitantes


Criar perguntas: Proibído
Responder a perguntas: Proibído
Editar Mensagens: Proibído
Apagar Mensagens: Proibído
Enviar anexos: Proibído

Pesquisar por:
Ir para: