HT_loss_oct.m

%calculation of HT loss for month oct
  %    SL.NO.   FR       IR        K
data=[       ]  %fr=fina reading.ir=initial reading.k=conts
      n=length(data(:,1)) 
      fr=data(:,2);%intilisation%
      ir=data(:,3);
      k=data(:,4);
      for i=1:1:n
      DTC_consumption(i)=(fr(i)-ir(i))*k(i)
      end
      DTC_con_sum=sum(DTC_consumption)
      es2_reading=input('enter the  month final reading =')
      es1_reading=input('enter the month initial reading =')
      Difference_month=(es2_reading-es1_reading)
      constant_K=input('K=')
      consumption_on_feeder_month=Difference_month*constant_K
      HT_installation_Consumption=input('HT Consumption=')
      Consumption_Rural_areas=input('rural area consupmtion=')
      Total_Consumption=DTC_con_sum+HT_installation_Consumption+Consumption_Rural_areas
      Loss_units=consumption_on_feeder_month-Total_Consumption
      Percentage_loss=(Loss_units/consumption_on_feeder_month)*100
%Calculation of LT Loss%        
    %Sl.No.    Fr     Ir    Diff   k  Tc_con  sold_eng  Diff
data= [1      2703.2 2473.3 229.9 30 6897    6551      346
       2      9509.8 8569.0 940   30 28218   23468     4750]
   n=length(data(:,1))
   for i=1:1:n
       Tc_con=data(:,6)
       sold_eng=data(:,7)
       Perc_lt_loss=((Tc_con(i)-sold_eng(i))/Tc_con(i))*100
   end
DTC_consumption=sum(Tc_con)
Total_consumption_Connected_DTC=sum(sold_eng)
LT_Line_loss=(DTC_consumption-Total_consumption_Connected_DTC)
Percentage_Line_loss=(LT_Line_loss/DTC_consumption)*100
% Losses using KVA-KM Method%
     % Sl.No. KVA KM  
data=[
      ]
  n=length(data(:,1))
  kva_as=data(:,2)
  km_as=data(:,3)
  for i=1:1:n
      product(i)=kva_as(i)*km_as(i)
      product_sum=sum(product) 
  end
      kva_sum=data(n,2)
      KM_sum=sum(km_as)
      Peak_load_input=input('enter the peak load current=')
      Annual_energy_sent=input('enter the energy sent=')
      Resistance_value=input('enter conductor resistance=')
      Load_factor=Annual_energy_sent/(sqrt(3)*11*0.8*8760* Peak_load_input)
      Loss_load_factor=(0.2*Load_factor+0.8*Load_factor^2)
      Diversity_factor=kva_sum/(sqrt(3)*11* Peak_load_input)
      Load_density_factor=(kva_sum*KM_sum)/product_sum
      Annual_energy_loss=(0.105*kva_sum^2*Resistance_value* KM_sum*Loss_load_factor)/(2*Load_density_factor*Diversity_factor^2)
      
      %voltage regulation and losses for  distribution conductors%
      [Rab_con,Wea_con,Coy_con,Rab_res,Wea_res,Coy_res] =textread('mydata.dat',...
                                 '%5.5f %9.5f %8.4f %5.5f %6.5f %7.6f',1)
      Percentage_reguln=(product_sum*Rab_con)/(2.5*100)
      Losses_conductor=(0.105*Loss_load_factor*product_sum*kva_sum*Rab_res)/(2*Diversity_factor^2)