Discussion :

[foufi86]

bonjour
mon problème est que j'ai fait le code suivant pour calculer le différence entre deux courbes:

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
function b = fonct(lam)
global Uexp RFexp
%lecture des données expérimentales
fid =fopen('données expérimentales.m','r');
M = fscanf(fid,'%e',[1,1]);
xy=fscanf(fid,'%e%e',[2,M]);
xy1=xy';
Uexp=xy1(:,1);     %vesteur colonne U2- déplacement 
RFexp=xy1(:,2);    %vecteur colonne RF2- réaction
mef(lam);
rf4;
Fint=interp1(Uexp,RFexp,-c(:,3),'spline');
b=0;
for i=1:108
    X=Fint((i)-RFexp(i)).^2
    b=b+X
end
b

après compilation il me donne l'erreur suivante :

??? Attempted to access Fint(1.99984); index must be a positive integer or logical.

Error in ==> fonct at 16
    X=Fint((i)-RFexp(i)).^2

alors qu'est ce que je dois faire???
merci d'avance.

[Franck Dernoncourt]

Tu veux pas dire plutôt

Code :

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

X=(Fint(i)-RFexp(i)).^2

?

[green_fr]

C'est just que (i)-RFexp(i) n'est pas un nombre entier.
Tu es sûs que tu ne voulais pas écrire plutôt ça :

Code :

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

(Fint(i)-RFexp(i)).^2

[foufi86]

salut
oui merci beaucoup je c pas comment j'ai pas fait attention a ca
merciiiii