clc
clear all
close all



% Kettle
pd_k = makedist('Normal', 'mu' , 480 , 'sigma' , 50 );
Kettle_time(1)= random(pd_k);
Kettle_time(2)= random(pd_k);
Kettle_time(3)= random(pd_k);
pd_k_n = makedist('Normal', 'mu' , 2 , 'sigma' , 0.5 );
k_n= random(pd_k_n);
if k_n < 1.5
    n=1;
elseif 1.5<= k_n && k_n<2.5
    n=2;
else
    n=3;
end
   
Kettle_power=zeros(1,1440);
k_n_t=0;
k_t_1=1;

while k_n_t < n
    
for i=1:1440;
    if (Kettle_time(k_t_1) < i) && (i <= (Kettle_time(k_t_1)+2))
        Kettle_power(i)=2000;
    end
end
    k_n_t=k_n_t+1;
    k_t_1=k_t_1+1;
end 

% Washing Machine
pd_w_m = makedist('Normal', 'mu' , 1200 , 'sigma' , 100 );
Washing_machine_time= random(pd_w_m);
Washing_machine_power=zeros(1,1440);
for i_w_m=1:1440;
    if (Washing_machine_time < i_w_m) && (i_w_m <= (Washing_machine_time+18))
        Washing_machine_power(i_w_m)=1800;
    elseif (Washing_machine_time+18 < i_w_m) && (i_w_m <= (Washing_machine_time+36))
        Washing_machine_power(i_w_m)=220;
    elseif (Washing_machine_time+36 < i_w_m) && (i_w_m <= (Washing_machine_time+42))
        Washing_machine_power(i_w_m)=1800;
    elseif (Washing_machine_time+42 < i_w_m) && (i_w_m <= (Washing_machine_time+54))
        Washing_machine_power(i_w_m)=220;
    end
end

% Computer
pd_C = makedist('Normal', 'mu' , 1080 , 'sigma' , 50 );
Computer_time= random(pd_C);
pd_C_d = makedist('Normal', 'mu' , 60 , 'sigma' , 10 );
Computer_duration = random(pd_C_d);
Computer_power=zeros(1,1440);
for i=1:1440;
    if (Computer_time < i) && (i <= (Computer_time+Computer_duration))
        Computer_power(i)=200;
    else 
        Computer_power(i) = 5;
    end
end

Demand=Kettle_power + Washing_machine_power + Computer_power;
plot(Demand)