Discussion :

[bissa1410]

Bonjour, j'ai un problème dans la programmation de la méthode de trapèze, j'ai l'intégrale suivante:
int(mu(s(w)))dw entre les bornes tau(i) et Tau avec Tau=0.05 et tau(i) on le calcul comme suit: tau(i+1)=tau(i)+h où h est le pas et est égal à 0.1 et i=1,50, f(s(w)) est une fonction de s qu'on peut la tirée du système de 3 équations qui contient l'une d'elle s'(t)=d(Sin-s(t))-mu(s(t))x1(t) où d et Sin sont deux constantes et x1 est la solution d'une autre équation différentielle du même système, ma piste est la suivante:
j'ai fait un petit programme d'Euler pour résoudre le système différentielle et après j'ai introduit la fonction trapz pour résoudre l'intégrale, en posant

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
h=0.1;
t=0:h:50;
tau=0.05;
Sin=1;
d=1;
K=0.94;
s(1)=2;
x1(1)=0.5;
x2(1)=1;
for i=1:length(t)-1,
mu(i)=s(i)/s(i)+K;
Le système differentielle avec la méthode d'Euler;
tau(i+1)=tau(i)+h;
w=linspace(tau(i),Tau,500)
y(i)=mu(s(i))
T(i)=trapz(w,y(i))

et l'erreur de matlab me dit qu'il y a problème de dimension????
Je ne sais pas est ce que c'est clair mais est ce possible de m'aider???
merci d'avance

[Jerome Briot]

Il nous faudrait le code complet et exact

De même que le message d'erreur exact et complet

[bissa1410]

Bonjour, voilà le code complet et exacte

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
close all , clear all, clc
h=0.1;		
t=0:h:50;	
 
k1  = 1;
k2  = 1;
K1  = 2;
K2  = .5;
Sin = 1;
m1  = 2;
m2  = 1;
Tau=0.05;
Dtilde1  =.01;
Dtilde2  =.03;
 
S(1)=2;	
X1(1)=0.5;
X2(1)=1;
 
tau(1)=0;
 
  for i=1:length(t)-1,
if      i<=50
             D(i) = 0.8;
 elseif i<=80
             D(i) = 0.05;
 elseif i<=150
             D(i) = 0.6;
 elseif i<=200
             D(i) = 0.2;
 elseif i<=300
             D(i) = 0.9;
 elseif i<=370
             D(i) = 0.3;
 elseif i<=400
             D(i) = 0.05;
 elseif i<=450
             D(i) = 0.7;
 elseif i<=500
             D(i) = 0.5;
 elseif i<=550
             D(i) = 0.8;
 elseif i<=580
             D(i) = 0.1;
 elseif i<=650
             D(i) = 0.6;
 elseif i<=700
             D(i) = 0.2;
 elseif i<=800
             D(i) = 0.9;
 elseif i<=870
             D(i) = 0.3;
 elseif i<=900
             D(i) = 0.05;
 elseif i<=950
             D(i) = 0.7;
else
             D(i) = 0.5;
end
 
D1(i)=D(i)+Dtilde1;
D2(i)=D(i)+Dtilde2;
 
    %le système
    mu1  = m1*(S(i)/(S(i)+K1));
    mu2  = m2*(S(i)/(S(i)+K2));
    f=((Sin-S(i))*D(i))-(k1*mu1*X1(i))-(k2*mu2*X2(i)); 
    g=(mu1-D1(i))*X1(i);
    l=(mu2-D2(i))*X2(i);
    S(i+1)=S(i)+h*f;		
    X1(i+1)=X1(i)+h*g;	
    X2(i+1)=X2(i)+h*l;
    %trapèze
    tau(i+1)=tau(i)+h;
 
    w=linspace(tau(i),Tau,500)
    y(i)=mu1(i)-D1(i);
    T(i)=trapz(w,y);
   end
 
S=S';
X1=X1';
X2=X2';
 
figure(1)
 plot(t,S','k','LineWidth',2); title('Concentration du substrat');xlabel('le temps t');ylabel('S');grid
figure(2)
 plot(t,X1','k','LineWidth',2); title('Concentration de l''espèce 1');xlabel('le temps t');ylabel('X1');grid
figure(3)
 plot(t,X2','k','LineWidth',2); title('Concentration de l''espèce 2');xlabel('le temps t');ylabel('X2');grid

[prgasp77]

De même que le message d'erreur exact et complet

Un membre prend le temps de se pencher sur ton problème. Prends le temps de lire son message. En entier. Oui, Toutes les DEUX phrases.