Discussion :

[alfredo_butticaz]

Je tente en ce moment de résoudre une équation complexe sur matlab, pour commencer je souhaiterais calculer une intégrale plus simple du style:

Code :

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

1
2
 
integral(@(x)x.^2+2.*integral(@(t)t,-1,x),-2,2)

Où est l'erreur?

Crdlmnt,
Alfredo B

[phryte]

Bonjour,

Je ferais comme cela:

Code :

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

1
2
3
syms x
A=int(@(t)t,-1,x)
integral(@(x) x.^2+2*(eval(A)),-2,2)

[alfredo_butticaz]

Cela fonctionne aussi, entre temps j'ai réussi en utilisant l'option 'ArrayValued',true

Un nouveau problème survient maintenant dans mon code, I(ceil(x)) dans l'intégrale engendre l'erreur suivante: Undefined function 'I' for input arguments of type 'double'. Alors que par définition ceil devrait me retourner un entier o_O

Merci d'avance,
Alfredo B.

[alfredo_butticaz]

Bonjour,

Je ferais comme cela:

Code :

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

1
2
3
syms x
A=int(@(t)t,-1,x)
integral(@(x) x.^2+2*(eval(A)),-2,2)

Bonsoir Phryte, tout d’abord merci pour votre contribution tout à l’heure,
ensuite, auriez-vous une idée de pourquoi l’erreur que j'ai est maintenant du type "Error using integral (line 85)
A and B must be floating-point scalars.*»

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
A=1:1:100; 
B=1:1:100; 
C=1; 
D=1:1:100; 
E=2; 
F=1:1:100;
G=1:1:100; 
H=0.1:0.1:10; 
I=500;
J=1000; 
K=2;
fun_1=@(t)integral(@(ti)G(ceil(ti)).*H(ceil(t-ti)),0,t-1);
fun_2=@(t)integral(@(ti)G(ceil(ti)).*B(ceil(ti)).*(C.*D(t).^E)./F(t).*...
    exp(-integral(@(x)(C.*D(ceil(x)).^E)./F(ceil(x)),ti,5)-K.*(t-ti)),0,t-1);
fun_3=@(t)I*integral(@(ti)min(fun_2(ceil(ti)),J).*exp(-(K+I).*(t-ti)),0.1,t-1);
t=2:1:5; figure(1) y = feval(fun_1,t); plot(t,y) figure(2) y = feval(fun_2,t); plot(t,y) figure(3) y = feval(fun_3,t); plot(t,y)

Crdlmnt,
Alfredo B