求解matlab高手-经管之家官网！

下面是我总matlab的一段程序上截取的一段：
for ii=1:5000;
if ii==1;
runs=1;simulation=0; markup=0; nop=250;
params=[0.0146 0.0208 0 0 0.3392 0.0181 0.0605 0 0.4432 2.0093];
if column==2;params=[0.0113 0 0.9964 0 0.4269 0 0.0608 0 0.4129 2.2209]; end;
if column==3;params=[ 0 0.0107 0 1.0101 0 0.1667 0.1290 0.0164 0.2163 3.4214]; end;
if column==4;params=[ 0 0 4.8440 0 0 0 0.0410 0 0.1710 2.0000]; nop=1; end;
if column==5;params=[ 0 0 0 0 0 0 0.0794 0.0108 0.1000 2.0129]; end;
if columnr==2; params=[0.0208 0.0028 0.4612 0.3598 0.2862 0.0103 0.0687 0 0.4899 2.1290]; enum=2; end;
if columnr==3; params=[0.0208 0.0174 0 0 0.2980 0 0.0613 0 0.4980 2.0000]; markup=1; end;
if columnr==4; params=[0.0089 0.0134 0.6157 0.1986 0.3031 0 0.0815 0 0.3931 2.1478]; nop=25; end;
if columnr==5; params=[0.0132 0 2.0562 0.0699 0.4698 0 0.0632 0 0.2481 2.1077]; nop=1; end;
if columnr==6; params=[0.0208 0.0200 0.2501 1.0392 0.4530 0.0201 0.1138 0.0162 0.3391 2.0000]; period=1; end;
%This initializes the parameters and criterion function. Value is the value function that spits out the criterion given a set of parameters
criterion=value(params(1,1),params(1,2),params(1,3),params(1,4),params(1,5),params(1,6),params(1,7),params(1,8),params(1,9),params(1,10));
params1=params; simcoeff1=simcoeff; criterion1=criterion; guess=0;scalar=0.25; v3=v2;

%This is the simulation annealing loop (for ii=2 onwards)
else;
jj=ii^(1/2); %Scaling parameter for cooling rate
randn('state',sum(100*clock)); z=[randn(1,2)/(100+jj),randn(1,2)/(2+0.02*jj),randn(1,2)/(10+0.1*jj),randn(1,2)/(20+0.2*jj),randn(1,1)/(10+0.1*jj),randn(1,1)/(2+0.02*jj)];
if crash==0; params=min(max(params+z+guess,[0 0 0 0 0 0 0 0 0.01 2]),upcap); end; %If not crashed take another guess, if crashed use best guess to reset
if column==2; params(2)=0; params(4)=0; params(6)=0; elseif column==3; params(1)=0; params(3)=0; params(5)=0; elseif column==4; params(1)=0; params(2)=0; params(5)=0; params(6)=0; elseif column==5; params(1)=0; params(2)=0;params(3)=0; params(4)=0; params(5)=0; params(6)=0; end;
criterion=value(params(1,1),params(1,2),params(1,3),params(1,4),params(1,5),params(1,6),params(1,7),params(1,8),params(1,9),params(1,10));
end;
%This shows the results after each loop
disp(' '); disp({'loop',ii,'column',column,'columnr',columnr,'nop',nop,'period',period,'gnum',gnum,'V2(n/2,n/2,1,1)',v2(round(gnum/2),round(hnum/2),1,1),'V2(n/2,n/2,2,2)',v2(round(gnum/2),round(hnum/2),2,2),'gstart',exp(gstart),'gfinish',exp(gfinish),get(0,'diaryfile')});
disp({ 'FCK ','FCL ','lamk','laml','irk ','irl ','mud ','mup ','var ','wagl'}); disp([params;params1]);
disp([simcoeff;simcoeff1;datacoeff]);
disp({'criterion', criterion,'criterion1',criterion1});
disp(datestr(now));
runs=2;
if isnan(criterion)==1; params=params1; criterion=criterion1; else; %if crashed return to baseline guess.
if criterion<criterion1; guess=0.33*(params-params1); params1=params; criterion1=criterion; simcoeff1=simcoeff; v3=v2;else guess=scalar*(params1-params); params=params1; end; end;
end;
这里面的for ,ifend 太多，把我弄混了，能否赐教，将对应的用不同的颜色标注。