matlab error operands to the || and & & operators must be convertible to logical scalar values
Ältere Kommentare anzeigen
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 Kommentare
Walter Roberson
am 29 Mär. 2019
Which line of code is the error on? You do not have any && or || in the code you posted
Mariam Gasra
am 29 Mär. 2019
Antworten (2)
Andrei Bobrov
am 29 Mär. 2019
Bearbeitet: Andrei Bobrov
am 29 Mär. 2019
0 Stimmen
Use
P(P< DA.Pl) | (P > DA.Ph) = 0;
instead
if (P< DA.Pl) | (P > DA.Ph)
P=0;
end
6 Kommentare
Walter Roberson
am 29 Mär. 2019
() are in the wrong place Andrei
Mariam Gasra
am 29 Mär. 2019
Walter Roberson
am 29 Mär. 2019
DA.Pl and DA.Ph are each vectors. What do you want it to mean to test values against those vectors?
Is the logic to follow that each element of P must be between one of the pairs? So in your data1 that uses [700 800; 100 400; 50 300] then you want to check to see whether (P >= 700 and P <= 800) or (P >= 100 and P<=400) or (P>=50 and P<=300) ? That is, it is to be accepted if it falls within any of the intervals and otherwise rejected ?
Mariam Gasra
am 29 Mär. 2019
Walter Roberson
am 29 Mär. 2019
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.
Mariam Gasra
am 30 Mär. 2019
Andrei Bobrov
am 29 Mär. 2019
0 Stimmen
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 Kommentar
Mariam Gasra
am 29 Mär. 2019
Kategorien
Mehr zu Logical finden Sie in Hilfe-Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!