Discussion :

[alfredo_butticaz]

Dans le code ci-dessous, a un moment, en utilisant G(ceil(ti)), l'erreur reçue est la suivante: Undefined function 'G' for input arguments of type 'double'. Alors que ceil est censé retourner un entier ?

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
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; 
fun_1=@(t)integral(@(ti)G(ceil(ti)).*H(ceil(t-ti)),0.1,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.1,t-1,'ArrayValued',true);
I=500;
J=1000; 
K=2;
fun_3=@(t)I*integral(@(ti)min(fun_2(ceil(ti)),J).*exp(-(K+I).*(t-ti)),0.1,t-1);
t=1:1:5;
figure(1)
fplot(fun_1,t);
figure(2)
fplot(fun_2,t);
figure(3)
fplot(fun_3,t);

Crdlmnt,
Alfredo B

[dante22wasa]

Bonsoir alfredo_butticaz,

J'ai testé votre script sur la version R2013a de Matlab et il a très bien fonctionné jusqu'à la ligne 20, où vous cherchez à tracer votre graphique de la fonction 'fun_2'. J'obtiens l'erreur suivante :

Code :

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

1
2
3
4
5
6
7
Subscript indices must either be real positive integers or logicals.

Error in
@(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))

Error in integralCalc/iterateArrayValued (line 157)
                fxj = FUN(t(1)).*w(1);

Soit une erreur dans le calcul de la fonction anonyme 'func_2'. J'ai essayé d'obtenir l'erreur que vous avez indiquée. Comment avez-vous procédé pour obtenir : 'Undefined function 'G' précisément ? Si vous avez posé un 'breakpoint' au sein de la fonction 'integralCalc' pour étudier le comportement de cette fonction en mode DEBUG : alors le vecteur 'G' n'existe plus tant que vous êtes au sein de 'integralCalc'.

Si j'en reviens à l'erreur que j'ai identifiée en lançant votre script, il vous suffit d'utiliser 'ceil' pour les vecteurs D et F dans 'fun_2' :

Code :

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

1
2
fun_2=@(t)integral(@(ti)G(ceil(ti)).*B(ceil(ti)).*(C.*D(ceil(t)).^E)./F(ceil(t)).*...
    exp(-integral(@(x)(C.*D(ceil(x)).^E)./F(ceil(x)),ti,5)-K.*(t-ti)),0.1,t-1,'ArrayValued',true);

Vous obtiendrez la courbe suivante :


Malheureusement, vous aurez encore une erreur sur la dernière ligne de fplot(fun_3,t) : vous demandez à Matlab de calculer une intégrale avec une borne 'b' vecteur ligne de longueur 150. En changeant 'fun_3' par :

Code :

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

fun_3=@(t)I*integral(@(ti)min(fun_2(ceil(ti)),J).*exp(-(K+I).*(t-ti)),0.1,t-1,'ArrayValued',true);

l'intégrale se calcule. Vous obtenez ceci en figure 3:


J'espère avoir pu vous aider sur votre problème, n'hésitez pas à me donner plus de détails s'il vous manque quelque chose ou si je suis à côté de la plaque .

En vous souhaitant une bonne soirée,

Cordialement,

D.

Le script modifié au complet :

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