Discussion :

[bailabar]

Aidez moi à résoudre numériquement ce système avec c++ ou Matlab.

[Invité]

Bonjour,
Résolution par différences finies (schéma explicite) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
tau = 0.01;
dt = 0.01;
k = 100;
b = 1;
t = [-tau:dt:10];
h(1) = 1;
h(2) = 1;
for i = 2:length(t)-1
    h(i+1) = dt*(-tanh(k*h(i-1))+b*cos(2*pi*t(i)))+h(i);
end
plot(t, h)