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

Question

sameer mohammad 2022-3-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1682219-how-to-write-hourly-varying-load-consumption-data-in-matlab-code-for-ieee-33-bus-system-to-get-curre

评论： Daravann MEL 2023-9-15

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

2 个评论
显示无隐藏无

Saleh Ba-swaimi 2023-1-9

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

Lonah Segera 2023-4-4

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 2023-6-28

在 MATLAB Online 中打开

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