2013-04-17 10:04:43 +08:00
|
|
|
|
function [ output_args ] = obj( this,x )
|
|
|
|
|
|
%OBJ Summary of this function goes here
|
|
|
|
|
|
% Detailed explanation goes here
|
|
|
|
|
|
PDi=this.PDi;
|
|
|
|
|
|
QDi=this.QDi;
|
|
|
|
|
|
PD=x(1:length(PDi));
|
|
|
|
|
|
QD=x(length(PDi)+1:length(PDi)+length(QDi));
|
|
|
|
|
|
wPD=this.wPD;
|
|
|
|
|
|
wQD=this.wQD;
|
|
|
|
|
|
PD0=this.PD0;
|
|
|
|
|
|
QD0=this.QD0;
|
2013-04-25 22:48:52 +08:00
|
|
|
|
t4=(wPD(PDi).*(PD-PD0(PDi))).^2;
|
|
|
|
|
|
% t4=wPD(PDi).*((PD-PD0(PDi))./PD0(PDi)).^2;
|
|
|
|
|
|
t5=(wQD(QDi).*(QD-QD0(QDi))).^2;
|
|
|
|
|
|
% t5=wQD(QDi).*((QD-QD0(QDi))./QD0(QDi)).^2;
|
2013-04-21 10:46:45 +08:00
|
|
|
|
Volt0=this.Volt0;
|
|
|
|
|
|
SEVolt=x(length(PDi)+length(QDi)+1:length(PDi)+length(QDi)+length(Volt0));
|
|
|
|
|
|
wVolt=this.wVolt;
|
2013-04-25 22:48:52 +08:00
|
|
|
|
t6=(wVolt.*(SEVolt-Volt0)).^2;
|
|
|
|
|
|
% t6=wVolt.*((SEVolt-Volt0)./Volt0).^2;
|
2013-04-21 10:46:45 +08:00
|
|
|
|
Balance=this.Balance;
|
|
|
|
|
|
t6(Balance)=0;%<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ƽ<EFBFBD><EFBFBD><EFBFBD>ڵ<EFBFBD><EFBFBD><EFBFBD>
|
2013-04-25 09:52:16 +08:00
|
|
|
|
% xx=t6(8);
|
|
|
|
|
|
% t6=xx;
|
|
|
|
|
|
% t6(9:end)=0;
|
2013-04-21 10:46:45 +08:00
|
|
|
|
% <EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|
|
|
|
|
SEVAngel=x(length(PDi)+length(QDi)+length(SEVolt)+1:end);
|
|
|
|
|
|
loadCurrent=LoadCurrent( SEVolt,SEVAngel,PD0(PDi),QD0(QDi),PDi,QDi );
|
|
|
|
|
|
wLoadCurrent=this.wLoadCurrent;
|
|
|
|
|
|
mLoadCurrent=this.mLoadCurrent;
|
2013-04-25 09:52:16 +08:00
|
|
|
|
% t7=(wLoadCurrent.*(loadCurrent-mLoadCurrent)).^2;
|
|
|
|
|
|
t7=wLoadCurrent.*((loadCurrent-mLoadCurrent)./mLoadCurrent).^2;
|
|
|
|
|
|
output_args=sum(t4)+sum(t5)+sum(t6)+sum(0);
|
2013-04-17 11:55:04 +08:00
|
|
|
|
output_args=full(output_args);
|
2013-04-17 10:04:43 +08:00
|
|
|
|
end
|
|
|
|
|
|
|
