%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