You can find these codes online. For example yarpiz.com provides you with the basic of these codes and you can implement it for your problem statement. There is a video tutorial also explaining how the code works. Hope this helps you. 

hi ,

you can use this two m file i attach to this message

if you want more meta_heuristic algorithm message me.

best regards

Thank you very much for your cooperation.

I reformulate my problem. In the attached file you will find the simulation that I am doing. The problem is to have an undesired result. Since the variables come out of their operation interval.

for example:

****************************************************************************

****************************************************************************

%T=300;

%k1 = physconst('Boltzmann');

%K2=(T*k1)/1.6e-19;

%i=Iph-I0(exp((v+Rs*i)/n*K2)-1)-(v+Rs*i)/Rsh

%objf=1/24*sqrt(sum(i-[x1-x2*(exp((v+x3*i)/x4*K2)-1)-(v+x3*i)/x5]))  such that

%x1 =Iph in [0,1]

%x2 = I0 in [0,1]

%x3 = Rs in [0,1]

%x4 = n in [1,2]

%x5 = Rsh in [20 100]

**************************************************************

***************************************************************

Desired interval:

X1 = Iph in [0,1]

X2 = I0 in [0,1]

X3 = Rs in [0,1]

X4 = n in [1,2]

X5 = Rsh in [20,100]

** 1st case:

In rows 123 and 124 (in pso.m) the limitation of the active variables

* obtained result:

 The variables remain in the desired work area but the result is required. (See FIG. 1)

2nd case:

The limitation of the variables (in pso.m) is deactivated.

* obtained result:

The curve compute converges towards the experimental curve but the variables diverge from their operating domain. (Figure 2).

X1 = 0.042289610901086 false but in [0 1]

X2 = -1.703735888318927 false and diverge from [0 1]

X3 = -17.115210839470510 false and diverge from [0 1]

X4 = 2.990782015135789 false and diverge from [1 2]

X5 = 14.277914413465945 false and diverge from [20,100]

you can use PSO directly for your task,  if you have a license for GADS Toolbox  !!!!

thank you for your reply, but unfortunately I do not have a lisence gads toolbox

for solar parameter estimation of PSO, can any one having the matlab code for that.