% CONTINGENCY RANKING BY JABER VALINEJAD
clear
clc
linedata=xlsread('linedata.xlsx');
busdata=xlsread('busdata.xlsx');
DCLOADFLOW
Pmax=P;
nl=linedata(:,1);
nr=linedata(:,2);
x=linedata(:,4);
pg=busdata(:,7)
nbus=length(busdata(:,1));
nline=length(linedata(:,1));
for i=1:nbus+nline
PI(i,1)=0;%PERFORMANCE INDEX
end
for i=1:nbus %DEPARTURE OF GENERATOR
if pg(i)>0
pg(i)=0;
DCLOADFLOW
for l=1:nbus
for k=1:nbus
if Pmax(l,k)~=0
PI(i,1)=((P(l,k)/Pmax(l,k))^10)+PI(i,1);
end
end
end
PI(i,1)=PI(i,1)/2;
end
end
for l=1:nline %DEPARTURE OF LINE
linedata=xlsread('linedata.xlsx');
linedata(l,:)=[];
DCLOADFLOW
for k=1:nbus
for h=1:nbus
if Pmax(k,h)~=0
PI(l+nbus,1)=((P(k,h)/Pmax(k,h))^10)+PI(l+nbus,1);
end
end
end
PI(i+nbus)=PI(i+nbus)/2;
end
disp(PI)
display('hghjj')
for i=1:nbus+nline+1
for j=1:nbus+nline+1
max(PI);
if max(PI)==0;break;end
if PI(j)==max(PI)
disp(max(PI))
disp(j)
PI(j)=0;
end
end
end
