Calculer delta t avec MATLAB. Depuis, en utilisant la méthode d'Euler, dP / dt peut être écrit comme (Pn + 1 - Pn) / delta t, décider d'abord sur cette valeur pour obtenir votre rapprochement. Tout ce que vous devez faire est de définir le temps de résolution, ainsi que le nombre d'intervalles de temps, et divisez-les pour obtenir delta t. Un exemple est la suivante:
t = 5
n = 50-
DeltaT = t / n-
Définissez une valeur pour tous les paramètres de l'équation. Si, par exemple, vous souhaitez mettre en œuvre Pn + 1 = Pn + x Delta T Pn, vous ne devez définir x dans MATLAB, puisque Delta T a été déjà définie à l'étape précédente:
r = .01-
Créer un tableau que vous allez utiliser pour calculer vos solutions. Remplissez-le avec 0 et définir une valeur de départ comme suit:
P = zéros (1, n + 1) -
P (1) = 2-
Utiliser un "pour" boucle pour obtenir la solution de l'Euler récurrente. Assurez-vous que vous utilisez également la "fin" commander pour éviter des erreurs de compilation.
pour i = 1: n
P (i + 1) = P (i) (1 + x deltat) -
fin
Tracer la solution afin de le visualiser. Cette étape est facultative, mais vous aurez une vue plus claire de votre rapprochement si vous le faites. Vous pouvez utiliser le "terrain" fonction comme suit:
parcelle (linspace (0, t, n + 1), p)