pso with wind 40uint

Using algorithm pso in economic dispatch 40 unit generation with wind power

%economic dispatch with considered valve point
clc
clear all
tic
pd=[70367110825884069480962897029776992410072101461022010072
    992497768554848086289776100729924962873327184];              % mw demand
pmin=[36 36 60 80 47 68 110 135 135 130 94 94 125 125 125 125 220
    220 242 242 254 254 254 254 254 254 10 10 10 47 60 60 60 90 90 90 25 25 25 242 1];
pmax=[114 114 120 190 97 140 300 300 300 300 375 375 500 500 500
    500 500 500 550 550 550 550 550 550 550 550 150 150 150 97 190 190 190 200 200 200 110 110 110 550 135];
s=70;
iter=100;
x（s41）=0;
kr=1;
d=1.12;
kp=1;
vi=5;
vr=15;
vo=45;
ll=（vr-vi）/vi;
%-----------------------------KC--------------------------------------%
v=[11.6 9 11.1 12.3 13.7 13.6;
    12.3 13.5 14.1 14.8 14.2 14.4;
    12.8 14 13.3 13.7 13.9 12.2;
    13.2 12.6 12 11.8 11.7 10.2;
    8.8 5.7 4.8 5.2 8.7 8.8;
    8.3 10.2 8.6 7.4 7.2 8.9;
    9.5 9.3 9.6 9.8 13.7 13.4;
    14.5 13.4 14.7 14.7 12.3 8.7;
    6.6 7.8 10.1 9.3 7.8 7;
    7.8 9.4 10 9.7 10.7 8;
    8.5 10.4 12.5 13.9 13.3 13.6;
    12.4 13.7 12.1 13.1 14.2 16.2;
    14.9 15.2 15.5 15.8 16.4 15.7;
    16.3 16.4 14.1 14.1 13.3 13.8;
    14.2 14.1 13.5 14 13.8 13;
    13.3 12.7 12.8 12.5 13.2 12.8;
    12.8 13 12.8 11.9 10.6 11;
    12.1 11.1 10 11 11 11.2;
    11 12 11 10.7 12.6 11.4;
    11.8 12.8 12.5 11.5 10.6 8;
    6.1 6.4 4.9 3.3 6.7 8.3;
    8 6.3 7.5 6.4 6 6.6;
    7 5.7 6.1 6.2 6.8 6.4;
    6.7 4.1 2.6 3.8 3 2.1];
%---------------------------------------------------------------------%
for i=1:24
   k=8;
   e=1;
   sum1=0;
   sum2=0;
   sum3=0;
   sum4=0;
   while abs（e）>=0.001
       sum1=0;
       sum2=0;
       sum3=0;
       sum4=0;
       for j=1:6
           sum1=sum1+（v（ij）^k）*log（v（ij））;
           sum2=sum2+log（v（ij））;
           sum3=sum3+v（ij）^k;
       end
       k1=（（sum1/sum3） -（sum2/6））^（-1）;
       e=k1-k;
       k=k1;
   end
   kg（1i）=k;
   for j=1:6
       sum4=sum4+v（ij）^k;
   end
   c（1i）=（sum4/6）^（1/k）;
end
%----------------------initialization process----------------//////
for tt=1:1:24
    for ttt=1:20
for i=1:1:s
    sum=0;
    flag=1;
    while flag==1
      k=floor（rand*41+1）;
        if k==42
          k=41;
        end
        for j=1:1:41
          if j~=k
              x（ij）=rand;
              while x（ij）<=（pmin（1j）/pmax（1j））
                  x（ij）=rand;
              end
              x（ij）=x（ij）*pmax（1j）;
              sum=sum+x（ij）;
          end
        end
      x（ik）=pd（1tt）-sum;
      if x（ik）<=（pmin（1k）/pmax（1k））
          flag=1;
          sum=0;
      else 
          flag=0;
      end
    end
end
for i=1:1:s
    for j=1:1:41
        pbest（ij）=x（ij）;
    end
end
cost=inf;
for i=1:1:s
    fp1=0.00690*x（i1）^2+6.73*x（i1）+94.705+abs（100*sin（0.084*（36-x（i1））））;
    fp2=0.00690*x（i2）^2+6.73*x（i2）+94.705+abs（100*sin（0.084*（36-x（i2））））;
    fp3=0.02028*x（i3）^2+

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      13652  2012-07-10 02:30  psowithwind40uint.m

----------- ---------  ---------- -----  ----

                13652                    1