O que há de tão bom nos solucionadores sem derivativos para SDEs?

8

Estou tentando me familiarizar com as SDEs e tenho lido alguns artigos de revisão sobre o assunto. Eles deixam a impressão de que uma grande quantidade de trabalho foi colocada em solucionadores livres de derivativos. Para meu entendimento, isso significa que, para um DDE como

dX=f(X)dt+g(X)dW,
as derivadas de f e g não são necessários para o método (me corrija se eu estiver errado).

Eu posso entender que essa propriedade é útil em algumas aplicações em que a derivada é difícil ou computacionalmente inviável de obter ou não existe. No entanto, eu não esperaria que esses problemas fossem muito relevantes na aplicação.

Isso sugere para mim que pelo menos um dos seguintes se aplica:

  • Há alguma outra vantagem relevante para os solucionadores sem derivativos que estou perdendo.

  • Problemas em que são necessários solucionadores sem derivativos (devido à razão acima) são mais relevantes do que eu acho.

  • A demanda por solucionadores livres de derivativos é menor que a “oferta”, ou seja, a atenção dada a eles por aqueles que desenvolvem solucionadores.

Qual é?

Wrzlprmft
fonte

Respostas:

7

Eu posso entender que essa propriedade é útil em algumas aplicações em que a derivada é difícil ou computacionalmente inviável de obter ou não existe. No entanto, eu não esperaria que esses problemas fossem muito relevantes na aplicação.

Se uma solução analítica para a derivada não for conhecida, é muito cara e propensa a erros. O cálculo das jacobianas é de entradas, mas as técnicas de diferenciação numérica precisarão fazer várias chamadas de função por entrada. Então, para acertar, as técnicas de diferenciação numérica precisam ser divididas por um pequeno número ao calcular a derivada, o que causa muitos problemas numéricos.n2

Com as ferramentas de autodiferenciação, esse custo é reduzido, mas ainda pode ser significativo. Portanto, quando os jacobianos analíticos não são prescritos, geralmente é bom ficar longe de métodos que exigem derivados.

No entanto, eu não esperaria que esses problemas fossem muito relevantes na aplicação.

Para a maioria das coisas, como SPDEs não lineares ou grandes sistemas de SDEs (milhares de anos) provenientes da biologia, escrever o jacobiano pode ser quase impossível e propenso a erros. Eu diria que é o contrário: esperar que um jacobiano analítico seja fornecido não é uma boa idéia.

Existem outras vantagens também. Os métodos Runge-Kutta são métodos livres de derivativos e podem fazer muita otimização de coeficiente.

A demanda por solucionadores livres de derivativos é menor que a “oferta”, ou seja, a atenção dada a eles por aqueles que desenvolvem solucionadores.

Esse não é o caso. No DifferentialEquations.jl, os métodos livres de derivativos foram implementados antes dos métodos KPS Stochastic Taylor Series porque, para a maioria dos usuários, levará a facilidade de uso e maior desempenho. Dito isto, no campo das equações diferenciais, sempre é possível encontrar um exemplo contrário, caso contrário, planejo implementar alguns métodos que explicitamente usam derivadas. No entanto, tenho certeza de que a maioria dos usuários provavelmente usará apenas os métodos livres de derivativos, porque a carga cognitiva no final é muito menor.

Chris Rackauckas
fonte
4

Não sou especialista em equações diferenciais especificamente estocásticas, mas suponho que minha resposta ainda terá algum valor.

  1. A computação da derivada pode ser desafiadora, como você mencionou na sua pergunta. No entanto, isso seria ainda mais pronunciado em um caso multidimensional, pois seria necessário calcular matrizes jacobianas ( entradas). Portanto, os solucionadores sem derivativos sofrerão com a maldição da dimensionalidade. A situação fica ainda pior quando derivativos de ordem superior são necessários para um esquema.n2
  2. A computação de um derivado por si só geralmente amplifica o ruído numérico. Portanto, por exemplo, se a função subjacente ( ou ) não for analítica, o erro na derivada poderá distorcer completamente a solução.fg
Anton Menshov
fonte
Obrigado pela sua resposta. Em relação ao ponto 1: Sim, você teria que calcular o jacobiano, mas na maioria dos casos, posso pensar que seria escasso, se não zero (pelo menos para ). Em relação ao ponto 2: Que tipo de aplicativo possui ou não analítico ? gfg
Wrzlprmft 30/03
Exigir que o usuário forneça um padrão de escarsidade para obter um desempenho decente é uma boa maneira de a maioria dos usuários ter um desempenho ruim. A maioria dos usuários deseja apenas um "solucionador automático de SDE" para a maioria dos problemas, como o que Dormand-Prince fornece, portanto, precisar dessa quantidade de entrada no "método mais básico" é um downgrade na usabilidade.
Chris Rackauckas
Sim, o jacobiano pode ser escasso. A questão é quão esparsa é, quão fácil é determinar o padrão de escarsidade e quantas avaliações de funções seriam necessárias para calculá-lo em comparação com um método sem derivadas. Em relação ao "ruído numérico". Também surgirá com funções analíticas, não tão severas (mas desafiadoras em termos computacionais o suficiente para examinar esquemas livres de derivativos). Funções não analíticas? Difícil responder quais são úteis no mundo estocástico da DE. Como especialista em equações integrais, sempre uso a função de Green como exemplo.
Anton Menshov