Discussion :

[thouraya_ad]

salut
j'ai écrit ce programme sur matlab,le problème c'est que dans la dernière étape il ne reconnait pas les x(i,j) et les y(i,j) bien qu'il les calcule au début(il dit que l'erreur est liée à inteline et feval,j'ai pas compris)

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
a=2
n=5
for j=1:2*n
       for i=1:2*n
           if i<=n      x(i,j)=(-a*(1+n-i))
		   else x(i,j)=(a*(1+n-i))
           end
        end  
    end
for i=1:2*n
       for j=1:2*n
           if j<=n      y(i,j)=(-a*(1+n-j))
		   else y(i,j)=(a*(1+n-j))
           end
        end  
    end
 
 
for i=1:2*n
    for j=1:2*n
        A(i,j)=dblquad('(x(i,j)-u)+(y(i,j)-v)',-5,5,-5,5)
    end
end

merci

[tug83]

explication de l'erreur:
http://www.mathworks.com/support/sol...lution=1-18I4D

[Jerome Briot]

Bonjour,

la fonction MESHGRID sera très utile ici. Si n devient grand, il faudra aussi penser à préallouer la mémoire pour la matrice de résultat A.

Code :

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

1
2
3
4
5
6
7
8
9
10
 
a=2;
n=5;
[y,x]=meshgrid([-a*n:a:0 -a:-a:-a*(n-1)]);
%A=zeros(2*n);
for i=1:2*n
    for j=1:2*n
        A(i,j)=...
    end
end

[rostomus]

Bonjour,
Voila une solution:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
a=2;
n=5;
[y,x]=meshgrid([-a*n:a:0 -a:-a:-a*(n-1)]);
A=zeros(2*n);
for i=1:2*n
    for j=1:2*n
        expression=sprintf('%d%s%d%s',x(i,j),'-u+',y(i,j),'-v');
        A(i,j)=dblquad(expression,-5,5,-5,5);
    end
end

[Jerome Briot]

Code :

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

expression=sprintf('%d%s%d%s',x(i,j),'-u+',y(i,j),'-v');

Ou :

Code :

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

expression=sprintf('%d-u+%d-v',x(i,j),y(i,j));

[rostomus]

Effectivement,

merci pour la remarque, Dutmatlab.

[thouraya_ad]

merci beaucoup à tous ce qui m'ont répondu ,ça marche;