Downloads.rar_MATLAB programs_lole_matlab LOLE

function COPT_LOLE(G,PR,A) format short g Lmax=input('Enter max Load in Load Model :'); Lmin=input('Enter min Load in Load Model :'); X=ff2n(G); InitiationMatrix=[zeros(1,2^G);zeros(1,2^G);ones(1,2^G);zeros(1,2^G)]; GeneratorCOPTMatrixTemp=InitiationMatrix'; for j=1:2^G for i=1:G if (X(j,i)==0) GeneratorCOPTMatrixTemp(j,1)=GeneratorCOPTMatrixTemp(j,1)+PR(i,1); GeneratorCOPTMatrixTemp(j,3)=GeneratorCOPTMatrixTemp(j,3)*A(i,1); else GeneratorCOPTMatrixTemp(j,2)=GeneratorCOPTMatrixTemp(j,2)+PR(i,1); GeneratorCOPTMatrixTemp(j,3)=GeneratorCOPTMatrixTemp(j,3)*(1-A(i,1)); end end end TemporaryMatrix=GeneratorCOPTMatrixTemp; for m=1:(2^G) for n=1:(2^G) if(GeneratorCOPTMatrixTemp(m,1)==GeneratorCOPTMatrixTemp(n,1)&& m~=n && n>m) GeneratorCOPTMatrixTemp(m,3)=GeneratorCOPTMatrixTemp(m,3)+GeneratorCOPTMatrixTemp(n,3); else end end end for m=1:2^G for n=1:2^G if(GeneratorCOPTMatrixTemp(m,1)==GeneratorCOPTMatrixTemp(n,1) && m<n && m~=n && GeneratorCOPTMatrixTemp(m,1)~=0) GeneratorCOPTMatrixTemp(n,:)=zeros; else end end end for m=1:1:((2^G)-1) for n=1:1:((2^G)-1) if (GeneratorCOPTMatrixTemp(n,1)<GeneratorCOPTMatrixTemp((n+1),1)) temp1=GeneratorCOPTMatrixTemp(n,1); temp2=GeneratorCOPTMatrixTemp(n,2); temp3=GeneratorCOPTMatrixTemp(n,3); GeneratorCOPTMatrixTemp(n,1)=GeneratorCOPTMatrixTemp((n+1),1); GeneratorCOPTMatrixTemp(n,2)=GeneratorCOPTMatrixTemp((n+1),2); GeneratorCOPTMatrixTemp(n,3)=GeneratorCOPTMatrixTemp((n+1),3); GeneratorCOPTMatrixTemp((n+1),1)=temp1; GeneratorCOPTMatrixTemp((n+1),2)=temp2; GeneratorCOPTMatrixTemp((n+1),3)=temp3; end end end GeneratorCOPTMatrix=GeneratorCOPTMatrixTemp; GeneratorCOPTMatrix(~any(GeneratorCOPTMatrixTemp,2),:)=[]; GeneratorCOPTMatrix; c=length(GeneratorCOPTMatrix(:,1)); suma=0; for i=c:-1:1 suma=suma+GeneratorCOPTMatrix(i,3); GeneratorCOPTMatrix(i,4)=suma; end l=length(GeneratorCOPTMatrix(:,1)); s=0; clc; fprintf('\n\nCalculated COPT for given Example as follows :\n\n'); fprintf('CAPACITY AVAILABLE\t\t CAPACITY UNAVAILABLE\t\t STATE PROBABILITY\t\t\t TIME (Ti) in Percentage\t\t Lost of Load Expected\n'); fprintf('=================================================================================================================================\n'); for i=1:c t=(Lmax-GeneratorCOPTMatrix(i,1))/(Lmax-Lmin)*100; if(t<0) t=0; end if(GeneratorCOPTMatrix(i,1)<Lmin) t=0; end lole=GeneratorCOPTMatrix(i,3)*t; if(lole<0) lole=0; end s=s+lole; fprintf('\t\t%d\t\t %d\t\t %10.8f\t\t %10.8f\t\t %10.8f\n',GeneratorCOPTMatrix(i,1),GeneratorCOPTMatrix(i,2),GeneratorCOPTMatrix(i,3),t,lole) end lo=s*365/100; fprintf('\n\nCalculated LOLE for above Problem ='); fprintf('%10.8f d/y\n\n',lo);