14 lines
461 B
Matlab
14 lines
461 B
Matlab
function [ output_args ] = calIWr(Vr,Vi,sigmaP,sigmaQ,P,Q )
|
|
%计算电流实部的方差
|
|
dIrdP=(Vr./(Vr.^2+Vi.^2)).^2.*(sigmaP.^2);
|
|
dIrdQ=(Vi./(Vr.^2+Vi.^2)).^2.*(sigmaQ.^2);
|
|
dIrdVr=(P.*(Vr.^2+Vi.^2)-2*(P.*Vr+Q.*Vi).*Vr)./(Vr.^2+Vi.^2);
|
|
% dIrdVr=(dIrdVr).^2.*0.1;
|
|
dIrdVi=(Q.*(Vr.^2+Vi.^2)-2*(P.*Vr+Q.*Vi).*Vi)./(Vr.^2+Vi.^2);
|
|
% dIrdVi=(dIrdVi).^2.*0.1;
|
|
dIrdVr=zeros(size(dIrdVr));
|
|
dIrdVi=zeros(size(dIrdVi));
|
|
output_args=max(abs([dIrdP,dIrdQ,dIrdVr,dIrdVi]),[],2);
|
|
end
|
|
|