pdeval

Avalia uma solução de pdepe e sua derivada em pontos escolhidos.

Categoria: Mathematics > Numerical Integration and Differential Equations

Sintaxe

[uout, duoutdx] = pdeval(m, xmesh, ue, xout)

Descrição

pdepe retorna valores apenas nos pontos da malha. pdeval interpola com a mesma discretização que o solver usa, de modo que a derivada devolvida é consistente com o fluxo usado na equação.

Para ver a solução ao longo do tempo, chame-a uma vez para cada linha de sol.

Argumentos de entrada

m

Simetria, como dada a pdepe.

xmesh

A malha usada com pdepe.

ue

Linha da solução em um instante fixo, por exemplo sol(end, :, 1).

xout

Pontos em que avaliar, dentro da faixa da malha.

Argumentos de saída

uout

Valores interpolados da solução.

duoutdx

Derivada espacial interpolada.

Exemplo

Valor entre pontos da malha

m = 0;
x = linspace(0, 1, 11);
t = linspace(0, 0.1, 5);
sol = pdepe(m, @(x, t, u, dudx) deal(1, dudx, 0), @(x) sin(pi*x), @(xl, ul, xr, ur, t) deal(ul, 0, ur, 0), x, t);
u = pdeval(m, x, sol(end, :, 1), 0.55);
fprintf('%.3f\n', u)

0.366

Veja também

pdepe interp1 deval