matlab error operands to the || and &amp &amp operators must be convertible to logical scalar values

2 次查看(过去 30 天)
clc;
clear all
close all
%%%
unit=input('input the number of unit:');
D=input('input total load :');
dP=D;
Bdat1=input('transmission losses considered matrix:');
B=Bdat1;
dB=diag(B);
data1=input('input the data:n a b c min max');
DA=array2table(data1,'V',{'Unit' 'a' 'b' 'c' 'Pl' 'Ph'});
x=max(DA.b);
n=input('insert number of iteration:');
for i=1:n
while abs(dP)>0.00001
P=(x-DA.b)./(2*(DA.c+x*dB));
P=min(P,DA.Ph);
P=max(P,DA.Pl);
dP=D+P'*B*P-sum(P);
x=x+dP*2/(sum(1./DA.c));
end
end
if (P< DA.Pl) | (P > DA.Ph)
P=0;
end
C=DA.a+DA.b.*P+DA.c.*P.*P;
totalCost=sum(C);
display(P);
display(totalCost);
lamda=x;
display(lamda);
Loss=P'*B*P;
display(Loss);
display(C);
display(P);
if i insert the data :Unit=3
Total load=850
Bdat1=[0.00003 0 0;0 0.00009 0;0 0 0.00012]
data1=[1 605 7.92 0.001562 700 800; 2 310 7.785 0.00194 100 400;3 78 7.97 0.00482 50 300]
n=1
i get the value of P=
500.0000
255.4601
109.3481
but i want the the first value of P which is 500 become zero
in general if one P less than min or greater than max set this P=0
  2 个评论
Mariam Gasra
Mariam Gasra 2019-3-29
clc;
clear all
close all
%%%
unit=input('input the number of unit:');
D=input('input total load :');
dP=D;
Bdat1=input('transmission losses considered matrix:');
B=Bdat1;
dB=diag(B);
data1=input('input the data:n a b c min max');
DA=array2table(data1,'V',{'Unit' 'a' 'b' 'c' 'Pl' 'Ph'});
x=max(DA.b);
n=input('insert number of iteration:');
for i=1:n
while abs(dP)>0.00001
P=(x-DA.b)./(2*(DA.c+x*dB));
P=min(P,DA.Ph);
P=max(P,DA.Pl);
dP=D+P'*B*P-sum(P);
x=x+dP*2/(sum(1./DA.c));
end
end
if P < DA.Ph || P > DA.Pl
P = 0;
end
C=DA.a+DA.b.*P+DA.c.*P.*P;
totalCost=sum(C);
display(P);
display(totalCost);
lamda=x;
display(lamda);
Loss=P'*B*P;
display(Loss);
display(C);
display(P);

请先登录,再进行评论。

回答(2 个)

Andrei Bobrov
Andrei Bobrov 2019-3-29
编辑:Andrei Bobrov 2019-3-29
Use
P(P< DA.Pl) | (P > DA.Ph) = 0;
instead
if (P< DA.Pl) | (P > DA.Ph)
P=0;
end
  6 个评论
Walter Roberson
Walter Roberson 2019-3-29
Your P is calculated from DA.b and DA.c, both of which are vectors, so your P will be a vector.
It would look to me to make more sense to check to see whether P is within the corresponding range for that element, rather than it being within any of the ranges.
However, I think Andrei's solution of
P(P < DA.Pl | P > DA.Ph) = 0;
should be all that is needed.
You talk about the value of P being 500, but P is going to be a vector, and whether 500 is valid or not would depend upon the row.

请先登录,再进行评论。


Andrei Bobrov
Andrei Bobrov 2019-3-29
Thanks Walter! My typo.
Use
P((P< DA.Pl) | (P > DA.Ph)) = 0;
instead
if (P< DA.Pl) | (P > DA.Ph)
P=0;
end
  1 个评论
Mariam Gasra
Mariam Gasra 2019-3-29
if i use this code P((P< DA.Pl) | (P > DA.Ph)) = 0;
the value of P=500 but i want it equal zero
if i use P((P> DA.Pl) | (P < DA.Ph)) = 0;
i get all value of P equal zero

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Mathematics 的更多信息

产品


版本

R2015b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by