Please read the scale_load() <http://www.pserc.cornell.edu//matpower/docs/ref/matpower5.1/scale_load.html> documentation carefully. To scale the load at a single bus (e.g. bus 3) you need to set the load_zone argument to select that bus …
nb = size(mpc(t).bus, 1); load_zone = zeros(nb, 1); load_zone(3) = 1; [mpc(t).bus, mpc(t).gen] = scale_load(0.5, mpc(t).bus, mpc(t).gen, load_zone); Ray > On Feb 19, 2016, at 7:59 AM, shruti seemala <[email protected]> wrote: > > sir > i am doing my project on demand response using time based program. > 1.i am taking an ieee 14 bus system. > 2.i apply demand response on one of the buses which have more demand. > 3.for that i require peak,off peak and valley period demand and price under > different periods. > 4.for 8hrs,4hrs and 10hrs the load variation is 0.5*peak > load,1*peakload,0.8*peakload respectively i assumed. > 5.i want to select 3rd bus and assume peakload as 41MW. > 6.i use this code for first 8hrs. > mpc0 = loadcase('case14'); > for t = 1:8 > mpc(t) = mpc0; > mpc(t).bus = scale_load(0.5, mpc(t).bus); > results(t) = runpf(mpc(t)); > end > i get some result.but i want the load to be change at third bus only as > mentioned w.r.t peakload > 7.for that i try this code > mpc0 = loadcase('case14'); > >> for t = 1:8 > mpc(t) = mpc0; > mpc(t).bus = scale_load(0.5, mpc(t).bus(4,4)); > results(t) = runpf(mpc(t)); > end > Attempted to access bus(:,3); index out of bounds because numel(bus)=1. > > Error in scale_load (line 168) > load_zone(bus(:, PD) ~= 0 | bus(:, QD) ~= 0) = 1; %% FIXED loads > it is showing like this. > i use mpc(t).bus(4,4) because it is the representation of 3rd bus load in the > bus data.is <http://data.is/> it correct approach. > 8.how can we set price for different periods. > 9.please suggest me
