Le versioni di matlab sono la 7.0.019901 r14 sul portatile e la 7.4.0287 r2007a sul fisso.
Per il programma ti posso postare una parte, perchè è molto lungo e composto da svariati m files:
Codice: Seleziona tutto
Mx=21;%x grid dimension
My=17;%y grid dimension
day=3;%total number of day
midday=6;
step=100;
stepf=50;%step fo the recostruction
N=day*step;
Nf=day*stepf;%total number of step for the recostruction
dt=0.01;
dtf=0.02;%dt for the recostruction function
R=2;%patch radius
cmask=0.08;%concentration mask for "fuzzy" velocities
%releasepatch coordinate
x0=10;
y0=11;
sig=5.0;% source intensity
lam=0.0; %feedback parameter
fcf=13.29;
fc=13.29;%new convertion factor
%fc=10;%old Leonid convertion factor
u=zeros(My,Mx);
v=zeros(My,Mx);
u_m=zeros(My,Mx);
v_m=zeros(My,Mx);
u_r=zeros(My,Mx);
v_r=zeros(My,Mx);
u_est_fuzzy=zeros(My,Mx);
v_est_fuzzy=zeros(My,Mx);
u_fuzzy=zeros(My,Mx);
v_fuzzy=zeros(My,Mx);
ufuzzy=zeros(My,Mx);
vfuzzy=zeros(My,Mx);
diffu=zeros(My,Mx);
diffv=zeros(My,Mx);
u_est_regression=zeros(My,Mx);
v_est_regression=zeros(My,Mx);
c0=zeros(My,Mx);
X=(1:Mx);
Y=(My:1);
a=[1 0.95 0.9 0.85 0.8 0.75 0.7 0.65 0.6 0.55 0.5 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1];
%a=[1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1];
%----------------------initial tracer distribution---------------------%
for i=1:My
y=i;
for j=1:Mx
x=j;
c0(i,j)= exp(-(x-x0)^2/R^2-(y-y0)^2/R^2);
if c0(i,j)<0.0
c0(i,j)=0.0;
else if c0(i,j)>0.0
c0(i,j)=c0(i,j);
end
end
end
end
figure()
contour(X,Y,c0,a)
caxis([0.1 1])
colorbar
title('Initial tracer distribution')
print -depsc Initial_tracer_distribution.eps
grid
%--------------------------------------------------------------%
comunque adesso ho notato che l'errore non me lo da solo per il comando contourf, ma, commentando quella parte di codice anche per il quiver successivo, e per tutti i contour o contourf successivi ??
Sembra quasi che siano i grafici i problemi

>:( >:( >:(
grazie mille per il tuo aiuto, se hai bisogno di maggiori info fammi sapere
grazie a presto
(b2b)
"Not everything that counts can be counted, and not everything that can be counted counts."