¿Cómo utilizar MATLAB para resolver la ecuación diferencial parcial del problema de valor límite de dos puntos?
Por ejemplo, si la condición inicial se cambia a x(2)=x'(2)=0 , entonces tspan(1) Debe ser 2.
Sin embargo, a menudo nos encontramos con problemas de valores límite en aplicaciones de ingeniería, lo que está más allá de las capacidades de la función ode**. Por supuesto, podemos escribir nuestra propia función para resolverlo.
Como disparar), pero eso está más allá de las capacidades de algunas personas. Afortunadamente, Matlab proporciona un solucionador de bvp.
solinit = bvpinit(x, yinit, params)
sol = bvpsolver(odefun, bcfun, solinit, options)
Porque el problema del valor límite puede tener many Una solución, para garantizar que la solución sea la que necesitamos, se debe utilizar la función bvpinit para estimar el valor inicial.
BVP solver: bvp4c y bvp5c se proporcionan en MATLAB, este último tiene mejor control de errores.
Parámetros de entrada:
x: el punto de la cuadrícula a calcular, equivalente al tspan de ode**
Yinit: valor estimado, que puede ser un Valor específico. También puede ser una función, similar a x0 de ode**.
Parámetros: Otros parámetros desconocidos también son conjeturas.
Un identificador de función que describe el problema de valor límite de una ecuación diferencial.
Bcfun: función de valor límite, generalmente valores bilaterales (los límites superior e inferior de X se consideran dos límites), pero también se admiten valores multilaterales (consulte la ayuda para obtener más detalles).
solinit: Grilla de inicialización generada por BVP init.
Opciones: Parámetros de optimización del solucionador BVP, que se pueden configurar a través de bvpset. Consulte la ayuda para conocer parámetros específicos.
Parámetros de salida:
La mayoría de ellos son los mismos que ode45.