how to write hourly varying load consumption data in matlab code for ieee 33 bus system to get current values and losses

i am working on IEEE 33 bus system and i have static load data so that i want use dynamic data like hourly load consumption data to get current values and losses for proposed test system,hence i request please guide me to writing the code in matlab

[Saleh Ba-swaimi]

Daer sameer mohammad, if you still need help please let me knew.

[Lonah Segera]

Dear Saleh Ba-swaimi. Could you please help me, I am working on the IEEE 15 bus and I need to represent the bus load data in hourly format but I don't know how to do it.

Answer 1

[Thomson Mtonga]

Dear sameer mohammad,

For a start, try running the attached file "UntitledElpa". It plots the load variations as are given in:

clc;

clear all;

format long

Industrial = [6:1:18];

Residential = [2,3,4,5,19,20,21,22,23,24,25];

Commercial = [26,27,28,29,30,31,32,33];

mpc = case33Loss202;

mpc.bus(2,3)

LoadVariInd = [0.56;0.54;0.52;0.50;0.55;0.58;0.68;0.80;0.90;0.98;1.00;0.94;0.95;0.96;0.90;0.83;0.78;0.72;0.71;0.70;0.69;0.67;0.65;0.60];

LoadVariRes = [0.69;0.65;0.62;0.56;0.58;0.61;0.64;0.76;0.9;0.95;0.98;1;0.99;0.99;1;0.96;0.96;0.94;0.93;0.92;0.91;0.88;0.84;0.72];

LoadVariCom = [0.2;0.19;0.18;0.18;0.2;0.22;0.25;0.4;0.65;0.86;0.9;0.92;0.89;0.92;0.94;0.96;1;0.88;0.76;0.73;0.65;0.5;0.28;0.22];

for i=1:length(Industrial(1,:))

AA(i,:)=(mpc.bus(Industrial(1,i),3))*1e3;

BB(i,:)=(mpc.bus(Industrial(1,i),4))*1e3;

end

for i=1:length(Residential(1,:))

AA1(i,:)=(mpc.bus(Residential(1,i),3))*1e3;

BB1(i,:)=(mpc.bus(Residential(1,i),4))*1e3;

end

for i=1:length(Commercial(1,:))

AA2(i,:)=(mpc.bus(Commercial(1,i),3))*1e3;

BB2(i,:)=(mpc.bus(Commercial(1,i),4))*1e3;

end

P = sum(AA);

Q = sum(BB);

S = (P^2+Q^2)^0.5;

P1 = sum(AA1);

Q1 = sum(BB1);

S1 = (P1^2+Q1^2)^0.5;

P2 = sum(AA2);

Q2 = sum(BB2);

S2 = (P2^2+Q2^2)^0.5;

for j=1:length(LoadVariInd(:,1))

Demand(j,:)=LoadVariInd(j,:)*S;

Demand1(j,:)=LoadVariRes(j,:)*S1;

Demand2(j,:)=LoadVariCom(j,:)*S2;

end

E=max(Demand);

E1=max(Demand1);

E2=max(Demand2);

Demand=Demand/E;

Demand1=Demand1/E1;

Demand2=Demand2/E2;

graph=plot([1:length(LoadVariInd(:,1))],Demand, '-oy',[1:length(LoadVariInd(:,1))],Demand1, '-og',[1:length(LoadVariInd(:,1))],Demand2, '-ob');

legend('Industrial Load Profile', 'Residential Load Profile', 'Commercial Load Profile');

set(graph,'LineWidth',1.3);

xlabel('Time (h)');

ylabel('Demand (p.u.)');

xlim([1 24]);

ylim([0 1]);

grid on;

Optimal Placement and Sizing of D-STATCOM in Radial and Meshed Distribution Networks Using a Discrete-Continuous Version of the Genetic Algorithm By Cristian Mateo Castiblanco-Pérez, David Esteban Toro-Rodríguez, Oscar Danilo Montoya and Diego Armando Giral-Ramírez

[Thomson Mtonga]

You can get the "case33Loss202" file from the following link:

https://www.mathworks.com/matlabcentral/fileexchange/129209-optimal-shunt-capacitors-placement-and-sizing-using-mvo?s_tid=srchtitle

[Daravann MEL]

Dear Mr. Thomson

I am sory in advance If I call you Last name.

Thank you for your sharing.

However, I have some question that is related to your code for the following point?

1, How did you know or select bus for residential, commercial, and industrial?

2, How did you know load variation (load factor) for those 3 load type?

[sevgi]

Dear Mr. Thomson

How can I perform a 24-hour power flow analysis (bus voltages, line currents, etc.) of a 33-bus radial distribution system?