bonjour,
je travaille actuellement sur des problèmes de champ de vitesse au sein d'un fluide, et je me sert de matlab pour l'exploitation des resultats. je n'arrive pas à tracer des courbes de niveau, que ce soit avec les écarts types, ou le champ moyen de vitesse. j'ai essayé avec les commandes: contour(x,y,nm), contourf(x,y,nm), etc....rien à faire.
pourriez vous m'aider SVP. Voici le programme; celui-ci fait appel à des fichiers resultat contenant des tableaux de dimension (1363x4),contenant, x,y,u,v. avec u et v les composantes de vitesses.merci d'avance:

clear all
close all
basename='E:\lit\LitFluidise_0 1\run_01\PIVgml.00-99_R01\i000';
extname='3i.res';
u=0;
v=0;
g=0;
n=0;
a=1363;
for i=0:9
for j=0:8

c=strcat(basename,int2str(i),i nt2str(j),extname);
f=load (c);
u=u+f(:,3);
v=v+f(:,4);
n=n+(u.*u)+(v.*v);

g=g+1;

end
end
x=f(:,1);
y=f(:,2);
um=double(u)/g;
vm=double(v)/g;
nm=double(n)/g;


%champ de vecteurs
%quiver(x,y,um,vm)
u=0;
v=0;
n=0;
for i=0:9
for j=0:8

c=strcat(basename,int2str(i),i nt2str(j),extname);
f=load (c);
e1=f(:,3)-um;
u=u+(e1.^2);
v=v+(f(:,4)-vm).^2;
n=n+(((u.*u)+(v.*v))-nm).^2;

g=g+1;

end
end
e_cu=sqrt(double(u)/g);
e_cv=sqrt(double(v)/g);
e_cn=sqrt(double(n)/g);

%vorticité

w=0;
g=0;
dx=1.2E+00
dy=1.8E+00

for i=0:9
for j=0:8

c=strcat(basename,int2str(i),i nt2str(j),extname);
f=load(c);
x=f(:,1);
y=f(:,2);
u1=f(:,3);
v1=f(:,4);
s=size(f);
T1=zeros(47,28);
T2=zeros(47,28);
h=1;
%for b=1:28
%for b=1:5


T1=reshape(u,47,29);
T2=reshape(v,47,29);
%h=h+47;
%end
%for o=1:46
%for m=1:27
for o=1:5
for m=1:5
%r1=r1+;
%r2=r2+;

w=w+(T1(o+1,m)-T1(o,m))./dy-(T2(o,m+1)-T2(o,m))./dx;
g=g+1;
end
end

end


end


wm=w./g;
wm=wm./2