`bvpinit`

Constrói a malha inicial e o palpite para um solver de problema de valor de contorno.

Categoria: Mathematics > Numerical Integration and Differential Equations

Sintaxe

solinit = bvpinit(x, yinit)

Usa a malha x e um vetor de palpite constante yinit.

solinit = bvpinit(x, yinitfun)

Calcula o palpite em cada ponto da malha com um handle y = yinitfun(x).

solinit = bvpinit(x, yinit, params)

Fornece também valores iniciais para parâmetros desconhecidos.

Descrição

Problemas de contorno não lineares dependem do palpite inicial, então escolhê-lo bem faz parte da solução. Um palpite constante é o mais simples; um handle de função permite codificar o que se sabe da forma da resposta.

A malha só precisa ser razoável, pois o solver a refina. Mais pontos ajudam onde se espera que a solução varie rápido.

Argumentos de entrada

x: Vetor de pontos da malha; a primeira e a última entrada são os pontos de contorno a e b.
yinit: Vetor com um palpite por componente da solução, ou um handle de função.
params: Vetor de palpites para parâmetros desconhecidos.

Argumentos de saída

solinit: Estrutura com os campos x e y (e parameters), para passar a bvp4c ou bvp5c.

Exemplo

Palpite constante

solinit = bvpinit(linspace(0, 1, 5), [0; 1]);
solinit.x
solinit.y

Saída

ans = 1×5

         0    0.2500    0.5000    0.7500    1.0000
ans = 2×5

     0     0     0     0     0
     1     1     1     1     1

Veja também

bvp4c bvp5c bvpset bvpxtend