function?[P?Fcost?Pl]=eld（dataBPd）
if?nargin?~=?3
?????error（‘Wrong?number?of?input?arguments‘）
end
ss=size（data）;
if?ss（2）~=5
????P=‘data?is?wrong?the?matrix?should?have?six?columns‘;
????Fcost=‘verify?your?data‘;
?????Pl=‘verify?your?data‘;
else
end
warning?off
n=length（data（:1））;
?????Aeq=ones（1n）;
?????a=data（:1）;
??????????b=data（:2）;
???????????????c=data（:3）;
????????????????????l=data（:4）;
?????????????????????????u=data（:5）;
???????????????????????????P=l;
?????????????????????????for?i=1:10
?????????????????????????????Pl=P‘*B*P;
?????????????????????????????Pd1=Pd+Pl;
?????????????????????????????ll=diag（1-2*B*P）;
?????????????????????????????A1=inv（ll）*a;
??????????????????????????????B1=inv（ll）*b;
??????????????????????????????H=2*diag（A1）;
??????????????????????????????P=quadprog（HB1[][]AeqPd1lu）;
?????????????????????????end
??????????????????????????????Fcost=sum（a.*P.*P+b.*P+c）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????953??2014-02-12?12:44??pso\eld.m
?????文件?????????172??2014-02-12?12:44??pso\forcecol.m
?????文件?????????181??2014-02-12?12:44??pso\forcerow.m
?????文件????????5788??2014-02-12?12:44??pso\goplotpso.m
?????文件???????28672??2014-02-12?12:44??pso\introduction.doc
?????文件??????126739??2014-02-12?12:44??pso\ll1411.pdf
?????文件????????4588??2014-02-12?12:44??pso\normmat.m
?????文件?????????922??2014-02-12?12:44??pso\psoeld.m
?????文件????????1032??2014-02-12?12:44??pso\psotest.m
?????文件???????22221??2014-02-12?12:44??pso\pso_Trelea_vectorized.m
?????文件?????????758??2014-02-12?12:44??pso\test.m
?????文件????????1530??2014-02-12?12:44??license.txt