Representações diferentes do espaço de estado para o filtro de regressão automática e Kalman

8

Vejo que existem diferentes maneiras de escrever um modelo AR em uma representação no espaço de estados, para que possamos aplicar o filtro Kalman para estimar o sinal. Veja os exemplos 1, 2 e 3 aqui .

Gostaria de saber quais são as diferenças entre as diferentes representações do espaço de estados na estimativa pelo filtro de Kalman?

Obrigado!

Tim
fonte
Este é o lugar certo para isso, não para a Ciência da Computação . Se você não obteve respostas, tente atualizar a postagem mostrando seu esforço na semana passada - você já tentou pesquisar? Outra opção é adicionar uma recompensa ...
Lorem Ipsum
A discussão parece ser mais teórica do que aqui. O filtro Kalman é um método ideal de estimativa para um sistema dinâmico estocástico. Portanto, ele se encaixa perfeitamente na ciência computacional. Ainda não encontrei nada útil.
Tim
você já tentou colocar uma recompensa? Você apenas tem que ter mais atenção à sua pergunta, e existem maneiras de fazê-lo ...
Lorem Ipsum

Respostas:

8

Infelizmente, não sei muito sobre os filtros Kalman, mas acho que posso ajudá-lo com o material do espaço de estado.

No Exemplo 1, o modelo AR é exatamente sua boa e antiga definição recursiva de saída DSP:

yt=α+ϕ1yt1+ϕ2yt2+ηt

Nesse caso, escrevemos o modelo de espaço de estados com correspondência direta com a equação acima:

(ytyt1)=(ϕ1ϕ210)(yt1yt2)+(α0)+(10)ηt

Observe que, neste caso, os estados do sistema são valores atuais e anteriores da saída.

No segundo exemplo, você está separando seus estados dos seus valores de saída. Isso significa que os estados agora podem ser qualquer coisa, mesmo que ainda estejam diretamente mapeados nos valores de saída. Desta forma, obtemosc

yt=μ+ct

ct=ϕ1ct1+ϕ2ct2+ηt

E portanto

(ctct1)=(ϕ1ϕ210)(ct1ct2)+(10)ηt

Você também deve reconhecer isso como a representação no espaço de estados padrão de um sistema linear, porque as equações para evolução do estado e saída dependente do estado são duas equações diferentes . Essa separação é trivial no caso de um modelo de RA, mas essa última notação é como pensamos em todos os modelos lineares de espaço de estados em geral.

ϕ1ϕ2α

Você deve observar que dois sistemas lineares podem ser idênticos até uma mudança de base. Isso significa que podemos escolher uma base diferente para representar o mesmo sistema linear. Você pode se convencer de que foi exatamente isso que fizemos para passar do segundo para o terceiro exemplo. Particularmente, gostamos dessa transformação linear para transpor a matriz de transição de estados, para obtermos algum estado desconhecidos

yt=(10)αt

αt=(stst1)=(ϕ1ϕ210)(st1st2)+(α0)+(10)ηt

Agora podemos usar a mudança de base para descobrir o que esse estado deve ser com relação ao estado . E podemos calcular que sejasy

(stst1)=(ytϕ2yt1)

Essa forma (transposição da forma canônica de controlabilidade) é chamada de forma canônica de observabilidade porque, se podemos colocar um sistema nessa forma, podemos deduzir facilmente quais estados do sistema podem ser observados simplesmente observando a saída. Para obter uma descrição das formas canônicas, você pode ler este documento e, é claro, dar uma olhada na web. Observe que no documento os estados são invertidos, o que não altera nada na representação do sistema, simplesmente reordenando as linhas / colunas das matrizes.

Phonon
fonte
2

Em resumo, tudo depende do que você está tentando estimar, isto é, o que você sabe sobre o sinal e o que não sabe. O filtro Kalman tentará estimar o estado com base na sua definição do que é esse estado. O problema convencional é quando estamos tentando estimar os coeficientes de RA.

Vamos dar um exemplo de um modelo sem termo constante .AR(2)μ

yk=a1yk1+a2yk2+ηk

Para estimar o sistema acima, tudo que você precisa fazer é estimar os coeficientes de AR, e .a1a2

Configuração geral do espaço de estado do filtro Kalman:

xk=Fk1xk1+wk
yk=Hkxk+vk
wk=WGN(0,Qs) evk=WGN(0,Qo)

Nesse caso, precisamos estimar e . Portanto, é natural definir o estado como esses coeficientes. Neste exemplo, esses coeficientes são constantes ( ) e não há ruído nesses coeficientes -> .a1a2xk=[a1,a2]TFk=Fk1=Iwk=0Qs=0

Como tudo o que observamos é , eles se tornam as medidas para o nosso sistema. Uma vez que já definiu o que o vetor de estado é, para os nossos equações de medição para ser igual ao modelo AR dada, substituímos o nosso ruído de medição com e .ykvkηkHk=[yk1,yk2]

y k = H k x k + η k = [ y k - 1 y k - 2 ] [ a 1 a 2 ] + η k

xk=xk1=[a1a2]
yk=Hkxk+ηk=[yk1yk2][a1a2]+ηk

Agora, você pode usar o filtro Kalman para estimar seu estado e, consequentemente, seu sinal.

Nota: A única coisa estranha aqui é sua matriz depende das suas medidas . Algumas pessoas têm a idéia errada de que a Matriz de Ganhos de Kalman e Covariância Estatal é sempre independente da medição e que podem ser calculadas com antecedência. Este caso mostra claramente que não é esse o caso. Tanto a matriz de ganho de Kalman quanto a covariância estatal são estimadas com funções de , que neste caso é dependente da medição.y k H kHkykHk

ssk08
fonte
Discordo. Eu acho que você observability Estado compromisso com a inclusão de medição na matriz