Application of residential demand response (DR) programs are currently realized up to a limited extent due to customers' difficulty in manually responding to the timedifferentiated prices. As a solution, this paper proposes an automatic home load management (HLM) framework to achieve the household minimum payment as well as meet the operational constraints to provide customer's comfort. The projected HLM method controls on/off statuses of responsive appliances and the charging/discharging periods of plugin hybrid electric vehicle (PHEV) and battery storage at home. This paper also studies the impacts of different timevarying tariffs, i.e., time of use (TOU), real time pricing (RTP), and inclining block rate (IBR), on the home load management (HLM). The study is effectuated in a smart home with electrical appliances, a PHEV, and a storage system. The simulation results are presented to demonstrate the effectiveness of the proposed HLM program. Peak of household load demand along with the customer payment costs are reported as the consequence of applying different pricings models in HLM.
1. Introduction
Smart Grid refers to the nextgeneration network that integrates information technology into the existing power grid to optimize energy efficiency through a twoway exchange of electricity/information between suppliers and customers in real time. In the context of energy efficiency, demand side management is one of the most important responsibilities of the future electric grid. To this end, smart grid operators design different programs to manage customers’ consumption patterns. Demand response (DR) programs play a major role in persuading customers to change their electric consumption patterns. The concept behind all the DR programs is to induce lower electricity use at times of high wholesale market prices or when system reliability is endangered.
DR offers a variety of financial and operational benefits for electricity customers and grid operators. Maintaining the supply of and demand for electricity in balance and preventing from high cost investment on generation and transmission are the most important benefits of DR application for grid operators
[1]
. The most important benefit of DR for customers is lower cost of electricity consumption. This would be resulted from adjusting load in response to the declared prices and decreasing market prices resulted from reduced supply costs. In addition, reductions in the probability of forced outages lead to imposing lower financial cost on the customer.
Since various timedifferentiated pricing models have yet been proposed, different pricebased DR programs can be offered to residential customers. Level of residential customers’ participation in the DR programs depends on the extent of comfort and benefits provided by the offered program. Several works verified that household customers would likely benefit from participating in DR programs if they properly respond to received signals.
The main obstacle of implementing pricebased DR programs is customers’ difficulty in manually responding to declared prices. Home load management (HLM) is the promise solution to alleviate this issue in the smart home. HLM system can automatically manage the electrical appliances consumption in response to the timevarying tariffs. Plugin hybrid electric vehicles (PHEVs) should be also incorporated in the HLM programs
[2]
. These vehicles are equipped with large batteries that are charged from the grid
[3]
. Analogously, the independent storage devices facilitate saving energy in the valleytime and sends energy back to the own consumption or to the grid in the peaktime of the load profile curve
[4]
. Vehicle to grid (V2G) describes an operation mode in which storage systems and electric vehicles are able to sell back electricity to the grid. This capability can effectively affect the HLM implementation results.
Several papers in the literature focus on controlling household appliances, such as washing machine, dryers, and other responsive plug loads
[2
5

6]
. Reference
[6]
presents a novel appliance commitment algorithm that schedules thermostatically controlled household loads based on the price and consumption forecasts. In this algorithm, the timevarying temperature range to reflect consumer choices on the appliances’ thermostat settings is specified as the user comfort. In
[7]
, an automatic energy consumption scheduler improvised in smart meters find the optimum scheduling based on the received price from the utility. The optimum scheduling is achieved with interactions among the users/customers and the utility company in the energy consumption game. The authors of
[8]
develop mathematical models for household appliances and propose the mathematical optimization models of residential energy consumption with the objectives of minimizing energy consumption, total cost of electricity and gas, emissions, peak load, and/or any combination of these objectives. These works study the issues of optimal appliance scheduling to minimize household energy consumption. To better decide on choosing and designing a proper pricebased DR program for customer and operator, it seems essential to investigate impacts of different pricebased DR programs on the HLM results. Neither of the previous works focuses on this challenge.
This paper aims to study HLM results, i.e., payment cost and peak of household load demand, incorporating various timedifferentiated tariffs. In doing so, at first, HLM is mathematically formulated as an optimizationbased problem to minimize the household payment cost of electric energy consumption during a day. Then, the result of HLM is studied under the implementation of different pricebased DR programs, i.e. real time pricing (RTP), time of use (TOU) pricing, peak time rebate (PTR) pricing, and the possible combinations of them. The on / off operation statuses of responsive loads including operation periods of household responsive appliances and charging/discharging periods of PHEVs and storage systems are outcomes of the proposed HLM problem. Vehicle to grid (V2G) capability is also considered in the problem. In order to provide customer comfort, the problem is subjected to practical constraints such as allowable interval for operation of appliances and charging PHEVs. Peak of household load demand along with the payment costs are broadly studied as the consequence of applying different timevarying tariffs. The results would explore the pros and cons of familiar pricedbased DR programs from the customer and operator points of view.
2. TimeVarying Tariffs
Economists have long advocated for electricity pricings that accurately reflect timevarying production costs. Electricity pricing usually peaks at certain predictable times of the day. In particular, if generation is constrained, prices can rise if power from other authorities or more costly generation is brought online. Billing customers by timevarying tariffs will encourage customers to adjust their consumption habits to be more responsive to market prices. Furthermore, regulatory and market design agency hope these price signals to delay the construction of additional generation or at least the purchase of energy from higher priced sources, thereby control the steady and rapid increase of electricity prices. In conclusion, timedifferentiated pricing models can potentially lead to economical and environmental advantages compared to the flat rates
[9]
. In this section, familiar timedifferentiated tariffs, i.e. real time pricing (RTP), time of use (TOU), and inclining block rate (IBR), are briefly described. Also, possible combinations of these tariffs are probed and their advantages are reported.
 2.1 TOU Tariff
In recent years, the TOU energy price is researched by many scholars and also implemented by lots of utilities
[10]
. It is recognized that carrying out TOU energy price has significant effects on the load shifting
[11]
. In the TOU tariff, electricity price changes in definite levels during hours of the day. In this paper, a threelevel (onpeak, midpeak, offpeak) TOU tariff is utilized. Threelevel TOU pricing can be formulated as below:
where,
t
is the hourly time index,
R
(
t
) is the TOU electricity tariff at hour
t
,
λ
_{1}
,
λ
_{2}
, and
λ
_{3}
are respectively tariffs at offpeak periods (
T_{1}
), midpeak periods (
T_{2}
), and onpeak periods (
T_{3}
) during a day. Obviously,
T
_{1}
∪
T
_{2}
∪
T
_{3}
= 24 hrs and
λ
_{1}
≤
λ
_{2}
≤
λ
_{3}
. Transferring energy consumption from onpeak tariff periods to lower tariff ones leads to peak load shaving and valley filling which are also desirable from distribution system operator point of view.
 2.2 RTP
In RTP, electricity prices may change as often as hourly. This type of pricing better presents the market price variation. This pricing mechanism is based on the fact that the energy system’s operation cost in the peak load time is much higher than that in the valley time
[12]
. The above principle insures that the energy price is higher when the consumption is higher. With RTP prices, the price of electricity is time dependent. When demand is low, less expensive sources of electricity are used. When demand rises, more expensive forms of electricity generator are called upon, making prices higher.
 2.3 IBR
Many utilities such as Pacific Gas & Electric (PGE), San Diego Gas & Electric, and the Southern California Edison companies have used IBR pricing for years
[13]
. An IBR tariff,
R
(
E_{t}
) , is mathematically presented as:
where,
E_{t}
is the total amount of consumed energy at period
t
, and
γ
is the predetermined threshold of IBR pricing which is time independent. Tariffs of the energy consumption lower and higher than
γ
are respectively
α
and
β
where
β
is obviously greater than
α
.
Applying IBR can lead to load balancing and reducing peak load
[14]
. In the IBR pricing, energy consumption more than a predetermined threshold would impose a penalty cost to the customer
[15]
. This penalty is such that the amount of consumed energy more than the threshold should be paid by a higher tariff than that of below the threshold.
 2.4 Combination of tariffs
In this subsection, possible combinations of introduced tariffs are proposed. It seems that combining these nonflat rate pricings can collect their advantages together.
Combining IBR and three level TOU tariffs can be mathematically represented by:
where, {
α
_{1}
,
β
_{1}
} , {
α
_{2}
,
β
_{2}
} , and {
α
_{3}
,
β
_{3}
} are, respectively, blocked rates of offpeak, midpeak, and onpeak periods. It is clear that
α
_{3}
≥
α
_{2}
≥
α
_{1}
and
β
_{3}
≥
β
_{2}
≥
β
_{1}
. In (2),
α
and
β
were constant. Nevertheless, in (3), calculating the electrical energy price by
α_{i}
and
β_{i}
depend on not only the total energy consumption at periods, but also the time of day. According to (1) and (3), It would be rational to have
α_{i}
less than
γ_{i}
and
β_{i}
higher than
γ_{i}
.
IBR and RTP can be similarly combined. At each period, two rates for lower and higher energy consumption than predetermined threshold would be likely defined. This can be mathematically represented as bellow:
The cost of energy consumption at period t less than and beyond threshold
γ
are respectively calculated by tariffs
α_{t}
and
β_{t}
.
3. Proposed HLM Method
Home appliances are divided into two categories, responsive and nonresponsive. The consumption of responsive appliances unlike the nonresponsive ones can be reduced and shift to offpeak tariff periods. Here, HLM is defined to schedule the operation periods of responsive appliances while satisfying operation constraints as well as constraints determined by the customers. Also, other developing technologies, i.e. PHEVs and storage systems, are incorporated and V2G capability is taken into account. Needless to emphasize, the HLM is solved from customer viewpoint. The HLM determines the optimal operation time vector for responsive appliances. Also, HLM determines the PHEV and storage system charging/ discharging periods.
The tariff is determined by the operator and the user’s total electricity payment within the upcoming scheduling horizon is minimized a
where, PC is the payment cost function,
Tr_{h}
is the tariff at period h,
E_{h}
is energy consumption of appliances at period
h
,
is the charged / discharged energy of the battery at period
h
,
is the charged/discharged energy of PHEV at period
h
, and 𝐇 is the set of periods in which the responsive appliances are scheduled. Here, 𝐇 is one day ahead. Since each scheduling period is assumed to be equal to ten minutes in this paper, the total number of periods would be 24 × 6 = 144 periods. So,
H
= {1,2,3,...,144}.
E_{h}
is the summation of responsive and nonresponsive appliances’ consumption which is formulated as:
where,
FE_{h}
is the energy consumption of nonresponsive appliances at period
h
,
E_{j}
is the energy consumption of appliance
j
at each period [kWh/period], and
I_{jh}
is status binary indicator of appliance
j
at period
h
, where 1 means appliance
j
is on at period
h
.
E_{j}
is the known parameter and
I_{ jh}
is determined by solving the problem.
HLM should provide customer’s comfort. So, possible customer constraints should be included in the optimization process. The customer is allowed to determine an allowable time interval for each responsive appliance operation, i.e. both the beginning and the end periods of allowable interval for operation. For example, the customer wants to wash the dishes after 10 a.m. and before 4 p.m.. So, he sets the beginning of time interval (
b_{j}
) to 10 a.m. and the end of time interval (
e_{j}
) to 4 p.m.. The operation time of responsive appliance j should lie within determined allowable time interval of that appliance. This is certified in (7):
where,
U_{j}
is the up time needed for proper operation of the appliance
j
. Obviously,
b_{j}
is lower than
e_{j}
and
e_{j}

b_{j}
is equal to or greater than
U_{j}
. For more clarification, assume that allowable operation time for washing machine is 9 a.m. and 11a.m. and required operation time to wash the clothes properly is 30 minutes. Having the scheduling period equal to 10 minutes,
U_{j}
=3,
b_{j}
=9 × 6+1=55, and
e_{j}
=11 × 6=66.
Also, appliances should operate uninterruptedly and the operation periods should be consecutive. This is mathematically represented by
where,
y_{jh}
is the start up indicator. Equation (9) is to relate the appliance status (
I_{jh}
), startup (
y_{jh}
), and shutdown (
z_{jh}
) indicators at period
h
for appliance
j
.
Since an appliance at a given period cannot be simultaneously started up and shut down, we have (10).
Other constraints related to PHEV and storage system, and their operation conditions are described in the following.
Two modes are defined for the battery:
B_{p}
indicates the charging mode and
B_{n}
signifies the discharging mode.
E_{Bp}
and
E_{Bn}
represent the chemical energy consumed or produced in the battery. Actual electrical energy is concluded with multiplying or dividing the chemical energy by AC/DC and DC/AC efficiency. Same as battery storage, two modes are defined for PHEV:
P_{p}
and
P_{n}
.
and
are respectively discharging and charging binary indicators of battery during period
h
.
and
are respectively discharging and charging binary indicators of PHEV during period
h
.
Therefore, the charged/discharged energy of battery and PHEV at each period incorporated in (6) are calculated by:
where,
η^{B}
,
η^{P}
,
η '^{B}
and
η '^{P}
are respectively charging efficiency of the battery, charging efficiency of PHEV, discharging efficiency of the battery, and discharging efficiency of the PHEV.
In addition, the batteries cannot be simultaneously in charging and discharging modes. Hence,
and,
For a PHEV, the outofhome interval, [
g_{p}
,
c_{p}
], is determined by the customer. So,
and
are zero in this interval.
The discharging value at each period should be less than the available battery charge. This is verified for battery storage in (15):
where,
is the initial charge of the battery storage. At the right hand side of (15),
E_{Bn}
appears with a minus sign to show the opposite direction of the discharging energy flow.
This constraint is slightly different for PHEV before and after the outofhome interval:
Where
is the initial charge of the PHEV battery. From (16), it is deduced that before the outofhome interval, the charge level of the PHEV battery is calculated similar to the storage system. But after arrival period of PHEV, its energy consumption during [
g_{p}
,
c_{p}
],
E_{p}
, is also considered for calculating the PHEV charge level.
Also, the capacity limit of the battery,
E
^{max B}
, is satisfied in (17).
The capacity limit of the PHEV battery,
E
^{max P}
, before and after outofhome interval is met in (18).
When the discharged energy is more than the home own consumption, the energy could be sold to the utility with the same tariff of that time, if profitable.
The other point is that PHEV charge level before
g^{P}
should be equal or greater than outofhome energy consumption. This is certified in (19).
Solving HLM determined the statuses of all appliances and other facilities for a dayahead period. Obviously, the results directly depend on the hourly tariffs. It seems that different pricings lead to not only different payment for the customer but also different home load curve, which is more important from operator viewpoint.
4. Numerical Study
This section analyzes the HLM formulation through different case studies to probe the impacts of different defined pricebased DR programs on the HLM results. The payment cost and load demand are reported in each case. The problem is solved by using a MIP solver, i.e. CPLEX 11.0.0 under GAMS environment. Computation times are not reported in detail as they all were less than 10 seconds on a PC equipped with a 3.2 GHz Pentium 4 processor and 2 GB of RAM.
We derive a variety of cases comprising
Case I: HLM considering TOU tariff; Case II: HLM considering IBR tariff; Case III: HLM considering combination of TOU and IBR tariffs;
As the concept behind the RTP and TOU pricing are the same, the results of applying RTP and combination of RTP and IBR are similar to applying TOU pricing and combination of TOU and IBR. Therefore, to prevent superfluous, RTP implementation is omitted from the case studies.
Nonresponsive appliances load demand is assumed as
Fig. 1
.
The nonresponsive load curve
Table 1
outlines the energy consumption (
E_{j}
), the number of operation periods (
U_{j}
), the beginning (
b_{j}
), and the end time (
e_{j}
) of allowable operation interval for each responsive appliance.
Responsive appliances data
Responsive appliances data
Battery storage and PHEV battery capacity is assumed 10 kWh and the efficiencies are equal to 0.88. Returning the energy back to the grid (V2G) capability is possible with the tariff of that period.
Threelevel TOU tariff, presented in
Fig. 2
, is obtained from
[16]
as the input of the HLM in Case I.
Threelevel tariffs for a day
In Case II, IBR pricing is applied to HLM algorithm. Generally, the distribution system company (DisCo) determines the tariff so that the DisCo’s cost plus a reasonable profit are covered. Consider a daily IBR pricing; if the energy consumption of the consumer during a definite period is less than a defined threshold, the customer will pay based on
α
[$/kWh] and if the energy consumption goes beyond the threshold, the amount of consumed energy below and beyond the threshold are calculated based on
α
and
β
, respectively. Due to proven positive impacts of IBR on the distribution system operation, the DisCo should specify
α
and
β
so that it convinces customers to participate in the IBR pricing. Based on the respected reviewer’s comment,
α
can be equal to the flat rate tariff to cover DisCo’s cost as well as its profit. However, in this way, no customer is willing to participate in IBR pricing; since, the payment cost in the case of applying IBR would be definitely more than the case of flat rate tariff application. Therefore,
α
should be less than flat rate tariff to encourage customers to participate in this program. Meanwhile,
β
should be higher than flat rate tariff. This creates incentives for end users to distribute their load at different times of the day in order to avoid paying for electricity at higher rates. Accordingly, the cost and a reasonable level of profit are guaranteed for the DisCo as well.
Based on the above discussion,
α
and
β
cannot be determined independently. They are directly related to each other as well as the threshold
γ
, and the flat rate tariff offered by DisCo.
The lower level of the IBR tariff at each period,
α
, is assumed to be the ¢ 9.866. The price of the higher level of IBR,
β
, is considered to be ¢ 14.898. The threshold
γ
for each 10minute period is assumed to be 1 kWh.
 4.1 Case I
In this case, TOU tariff is the input of the problem and the HLM is solved by the proposed formulation. The results of three cases are summarized in
Table 2
.
Operating Periods of Responsive Appliances in Case I
Operating Periods of Responsive Appliances in Case I
Referring to
Table 2
, the responsive appliances are on in the lowest tariff periods in their allowable interval of operation.
The net load resulted from the HLM scheduling is depicted in
Fig. 3
. The net load is the summation of the appliance consumption and the PHEV and battery storage charging/discharging. Positive and negative loads represent the energy flow from and toward the grid, respectively. It can be concluded from this figure that the peak and average loads are 6.851 kW and 1.315 kW, respectively. The payment in this case is ¢ 144.481. This cost is expectedly less than the imposed cost in cases without battery storage and PHEV, because of more potential of returning energy back to appliances or the main grid.
Net load in Case I
 4.2 Case II
Considering IBR tariff, operation periods of responsive appliances resulted from HLM application is presented in
Table 2
.
The results show that to prevent more payment cost due to exceeding the threshold
γ
, operation periods of appliances shift to the valley periods of load demand. The resulted load demand shows that the load is limited to the determined threshold such that the peak load is 5.05 kW. In comparison with the previous case, this is more desirable from operator viewpoint. Also, the results certify that applying IBR prevents simultaneous charging of PHEV and battery storage. The payment in this case is ¢ 147.396 which is more than the payment cost in the previous case.
 4.3 Case III
This case combines IBR and TOU pricing as mentioned in Section 2. The lower level of the IBR tariff at period
t
,
α_{t}
, is assumed to be 0.9 of TOU tariff at that period. The price of the higher level of IBR,
β_{t}
, is considered to be 1.4 of
α_{t}
[5]
. The threshold
γ
is also assumed to be 1 kWh.
HLM leads to the presented results in
Table 2
. Although the resulted payment cost, ¢ 150.627, is more than the previous cases, the resulted peak load is 4.902 kW which is less than the previous cases. So, it seems that applying IBR is more desirable from operator facet. However, payment cost directly depends on the tariff rates which can be properly designed to reduce payment cost of customer in addition to the peak load in the case of applying IBR.
5. Conclusion
In this paper, an inhome optimizationbased algorithm is proposed to schedule demand side facilities and determine on / off status of responsive appliances and other household facilities. Different timevarying tariffs are employed in HLM to probe the payment cost resulted from solving the problem. The results show that in the case of applying TOU, highconsumption responsive appliances, storage system, and PHEVs might simultaneously operate/ charge at offpeak time and the solution probably creates a new peak during these periods. Application of IBR to TOU pricing prevents the negative impacts of the proposed method on the household load curve; since, IBR application causes appliances operation to be distributed during the allowable interval of their operation. Designing a proper tariff rate as an open research area is currently under our study.
BIO
Mohammad Rastegar (S’12) received the B.S. and M.Sc. degrees from Electrical Engineering at Sharif University of Technology, Tehran, Iran, in 2009 and 2011, respectively. He is currently pursuing the Ph.D. degree in the Electrical Engineering Department at Sharif University of Technology, Tehran. His research interests include power system reliability, demand response, smart home, and smart grid. He is Technical Manager at Zeytoon Foundation.
Mahmud FotuhiFiruzabad (SM’99) received the B.S. and M.S. degrees in electrical engineering from Sharif University of Technology and Tehran University, Tehran, Iran, in 1986 and 1989, respectively, and the M.S. and Ph.D. degrees in electrical engineering from the University of Saskatchewan, Saskatoon, SK, Canada, in 1993 and 1997, respectively. He joined the Department of Electrical Engineering at Sharif University of Technology in September 2002. Currently, he is a Professor and Head of the Department of Electrical Engineering, Sharif University of Technology. He is also an Honorary Professor in the Universiti Teknologi Mara (UiTM), Shah Alam, Malaysia. Dr. FotuhiFiruzabad is a member of the Center of Excellence in Power System Management and Control. He serves as an Editor of the IEEE TRANSACTIONS ON SMART GRID. Recently, he is a fellow member of IEEE PES.
Jaeseok Choi(S’88, M’91, SM’05) was born in Kyeongju, Korea, in 1958. He received the B.Sc., M.Sc., and Ph.D. degrees from Korea University, Seoul, in 1981, 1984, and 1990, respectively. He was a Postdoctoral at the University of Saskatchewan, Saskatoon, SK, Canada, in 1996. He was a visiting professor at Cornell University, Ithaca, NY, USA, in 2004. Since 1991, he has been on the faculty of Gyeongsang National University, Jinju, Korea, where he is a professor. He is also adjunct professor at IIT, IL, USA since 2007. His research interests include fuzzy applications, probabilistic production cost simulation, reliability evaluation, and outage cost assessment of power systems.
Albadi MH.
,
ElSaadany EF.
2007
“Demand response in electricity markets: an overview”
presented at IEEE Power Eng. Soc. Gen. Meeting
ClementNyns K.
,
Haesen E.
,
Driesen J.
2010
“The Impact of Charging PlugIn Hybrid Electric Vehicles on a Residential Distribution Grid”
IEEE Trans. Power Syst.
25
371 
380
DOI : 10.1109/TPWRS.2009.2036481
Vazquez S.
,
Lukic S. M.
,
Galvan E.
,
Franquelo L. G.
,
Carrasco J. M.
2010
“Energy storage systems for transport and grid applications”
IEEE Trans. Indus. Electronics
57
3881 
3895
DOI : 10.1109/TIE.2010.2076414
MohsenianRad A. H.
,
LeonGarcia A.
2010
“Optimal residential load control with price prediction in realtime electricity pricing environment,”
IEEE Trans. Smart Grid
1
120 
133
MohsenianRad A. H.
,
Wong V. W. S.
,
Jatskevich J.
,
Schober R.
,
LeonGarcia A.
2010
“Autonomous DemandSide Management Based on GameTheoretic Energy Consumption Scheduling for the Future Smart Grid,”
IEEE Trans. Smart Grid
1
(3)
320 
331
DOI : 10.1109/TSG.2010.2089069
Bozchalui M. Ch.
2012
“Optimal Operation of Residential Energy Hubs in Smart Grids,”
IEEE Trans. Smart Grid
3
(4)
1755 
1766
DOI : 10.1109/TSG.2012.2212032
2006
Benefits of Demand Response in Electricity Markets and Recommendations for Achieving Them. Report to the United States Congress
http://eetd.lbl.gov.
Time of use tariff rates for BGE
Available at: http://www.ieso.ca/imoweb/siteshared/tou_rates.asp
Gellings C.
,
Chamberlin J.
1988
Demand Side Management: Concepts and Methods
the Fairmont Press Inc
Lilburn
Tang Y.
,
Song H.
,
Hu F.
,
Zou Y.
2005
“Investigating on TOU pricing principles”
in Proc. IEEE PES Trans. Dist. Conf. Expos.
Dalian, China
Borenstein S.
2008
Equity effects of increasingblock electricity pricing Center for the Study of Energy Markets, Working Paper 180
MohsenianRad A. H.
,
Wong V.
,
Jatskevich J.
,
Schober R.
2010
“Optimal and autonomous incentivebased energy consumption scheduling algorithm for smart grid,”
presented at the IEEE PES Conf. Innov. Smart Grid Technol.
Gaithersburg, MD
Available at: http://www.ieso.ca/imoweb/siteshared/tou_rates.asp?sid=ic&Chart=SUMMER#chart
Jin Young Gyu
,
Choi TaeSeop
,
Park Sung Chan
,
Yoon Yong Tae
2013
“The Design of an Optimal Demand Response Controller Under Real Time Electricity Pricing”
Journal of Electrical Engineering & Technology
(http://dx.doi.org/10.5370/JEET.2013.8.3. 436)
8
(3)
436 
445
DOI : 10.5370/JEET.2013.8.3.436