intégration sur matlab avec trapz

[invitef90a29aa]

Bonjour,

j'ai fait un programme sur matlab qui utilise l'intégration, pour cela j'ai utilisé la méthode du trapèze, mais une erreur dont je connais pas la cause apparait:

??? Error using ==> permute
ORDER contains an invalid permutation index

Error in ==> trapz at 44
  y = permute(y,perm);

Error in ==> CB1wercor at 40
q=trapz(y,k);

voici mon programme:

function q=CB1wercor(t,y,n,teta,sigma,beta)

der=Mywernum(y,lamda1,lamda2,p,b)
dry=numder(der,y,h,n1)
Fty=Fwety(t,y,sigma,beta)

 y=0.01:0.01:t;
 k=(1-Fty).*dry;
q=trapz(y,k);

sachant que t,y,n,teta,sigma,beta sont des variables d'entrée que j'introduis en utilisant l'outil input; et les fonctions Mywernum, numder, Fwety sont des fonctions que j'ai programmé au préalable et qui ne génèrent aucune erreur puisque Matlab donne les résultats de ces fonctions et s'arrête à la ligne où se trouve trapz.

merci de me répondre
-----

[invite15928b85]

Bonsoir.

Cela n'a peut-être rien à voir avec l'origine du problème mais je trouve étrange que y, paramètre de la fonction, soit redéfini dans la fonction...

@+

[invitef90a29aa]

bonjour,
y n'est pas redéfinie dans la fonction, mais trapz nécéssite la notation:

y=0.01:0.01:t;
 k=(1-z).*dry;
q=trapz(y,k);

pour dire qu'on intégre la fonction k par rapport à y de 0.01 à t.

[invite15928b85]

Bonjour.

J'interprète les choses comme suit :

la fonction CB1wercor va :

de sa déclaration : function q=CB1wercor(t,y,n,teta,sigma,beta)

jusqu'à l'instruction : q=trapz(y,k);

L'argument formel y est utilisé dans l'appel d'autres fonctions puis redéfini par : y=0.01:0.01:t;

Je ferais deux choses:

1 : changer le nom de la variable y utilisée dans la fonction trapz. Il pourrait très bien y avoir interférence entre les deux usages, la redéfinition dans la fonction venant écraser quelque chose en mémoire (je pense à un comportement similaire à celui des pointeurs en C quand on commence à faire des bêtises ...)

2 : vérifier que les arguments de la fonction trapz sont de dimensions compatibles, avec le debugger par exemple.

Sorti de là, je n'ai pas d'autre idée car il n'y a aucun raison que la fonction trapz s'envoie en l'air comme cela si les arguments réels sont corrects.

Cordialement,

F

[A voir en vidéo sur Futura]

[invitef90a29aa]

bonjour,
dslée mais je n'ai pas compris comment je pourrais vérifier les dimensions des arguments de la fonction trapz avec le debugger
svp aidez moi

[invite0ab6bca5]

Bonsoir,
Il faut en effet vérifier que les deux arguments de trapz ont même dimension. Vous pourriez, juste au-dessus de la ligne "q = trapz(..." introduire deux lignes de sorties provisoires qui seraient, par exemple,
dimy=length(y)
dimk=length(k)

Bonne continuation