MATLAB: Too much time taking for output

Question

while running this its taking too much time please help me,

thanks in advance

Answer 1

I might have messed your code up some, but you will get the idea. Your biggest problem was you were not preallocating. Several for-loops can be eliminated. If you want people to help you, you must format your code properly.

%Modeling of PV system using MATLAB
filename = 'book1.xls';
sheetname = 'data';
I_L=xlsread(filename, sheetname,'D2:D8761'); %Load Current (A)
T= xlsread(filename, sheetname, 'A2:A8761'); %Ambient Temp (°C )
S= xlsread(filename, sheetname, 'B2:B8761'); %Solar Radiation(W/m²)
NOCT=43.6;
Alpha=0.068; %Temperature Coefficient
IPV_M=6.89; %Current at MPP
T_cell=T+S*(NOCT-20)/800; %Cell Temp
I_PV=S/1000+Alpha*(T_cell-25);
I_PV(I_PV<0)=0;
%(Routine 2):Initialization
PV_eff=0.16; %PV module efficiency
Wire_eff=0.98; %Wire efficiency
INV_eff=0.95; %inverter efficiency
V_sys=230; %System voltage (V)
PSH=(mean(S)*12)/1000;
V_B=12; %Battery voltage
DCharge_eff=0.8; %Battery charging efficiency
DOD=0.8; %Allowed Depth of charge (%)
Alpha=0.068;
PV_WP=120;
%(Routine3):Intuitive method to find the initial N of PV modules needed
P_L=I_L*230;
A_PV=1.408*0.56;
E_L=sum(P_L)/length(I_L)*24;
N_PV = ceil(E_L/(PV_eff*INV_eff*Wire_eff*PSH*1000*A_PV));
N_PVmin=ceil(N_PV/3);
N_PVmax=5*N_PV;
nh=2;
C_battery=(E_L*nh)/(DOD*DCharge_eff);
C_batterymin=C_battery/3;
C_batterymax=5*C_battery;
%(Routine 4):Find LLP at each Battery Capacity and PV module trial
SOC=zeros(1,length(I_L));
I_Load=zeros(1,length(I_L));
I_Charge=zeros(1,length(I_L));
I_Discharge=zeros(1,length(I_L));
I_Deficit=zeros(1,length(I_L));
I_Damp=zeros(1,length(I_L));
I_Battery=zeros(1,length(I_L));
xxx=N_PVmin:N_PVmax;
yyy=C_batterymin:500:C_batterymax;
C_batteryf=zeros(length(xxx),length(yyy));
C_PVf=zeros(length(xxx),length(yyy));
LLP_calculated=zeros(length(xxx),length(yyy));
t=1;
Vmp=17.4;
NOCT=43.6;
K=0.8;
D=1e-5;
ns=6; %Number of cells per battery
SOC1=1; %Max. battery SOC=100% of total capacity
SOC2=SOC1;
SOC3=0.2; %Min. battery SOC=20% of total capacity
x=1;
I_net=I_PV-I_L;
for m=xxx %Number of PV modules
    y=1;
    for n=yyy %capacity of battery
        SOCmax=n; %Battery Capacity (Wh)
        SOCmin=SOCmax.*(1-DOD);
        w=0;
        for k=1:length(I_L) %Number of hours
            %(Routine 5):In Case of the battery is not empty
            if SOCmax>0
                %(Routine 5.1):In Case of I_PV=I_Load
                if I_net(k)==0
                    I_Load(k)=I_PV(k);
                    if k==1
                        SOC(k)=SOC1;
                    elseif w==0 %In case of battery is not discharged on the previous stage
                        SOC(k)=SOC1;
                    elseif w==1%In case of battery is discharged on the previous stage
                        SOC(k)=SOC(k-1);
                    end
                    %(Routine 5.2):In Case of I_PV>I_Load
                elseif I_net(k)>0
                    I_Load(k)=I_L(k);
                    if k==1
                        SOC(k)=SOC1;
                        I_Damp(k)=I_net(k);
                    elseif k>1
                        if w==0
                            SOC(k)=SOC1;
                            I_Damp(k)=I_net(k);
                        elseif w==1
                            if SOC(k-1)>=SOC1
                                SOC(k)=SOC(k-1);
                                I_Damp(k)=I_net(k);
                            elseif SOC(k-1)<SOC1
                                I_Charge(k)=I_net(k);
                                %%Charging mode
                                B=SOC2;
                                V1=(2+0.148*B)*ns;
                                R1=((0.758+(0.1309/(1.06-B)))*ns)/SOCmax;
                                %syms v;


                                ee=K*V1*I_net(k)-D*SOC2*SOCmax;%,v,0,t));
                                SOCx=SOC1+ee/SOCmax;
                                SOC2=SOCx;
                                SOC(k)=SOC(k-1)+abs(SOC1-SOC2);
                            end
                        end
                    end
                    %(Routine 5.3):In Case of I_PV<I_Load
                elseif I_net(k)<0
                    if w==0
                        I_Discharge(k)=I_L(k)-I_PV(k);
                        I_Battery(k)=I_Discharge(k);
                        I_Load(k)=I_PV(k)+I_Battery(k);
                        %%Discharging mode
                        B=SOC2;
                        V1=(1.926+0.124*B)*ns;
                        R1=(0.19+(0.1037/(B-0.14)))*(ns/SOCmax);
                        %syms v;
                        ee=K*V1*I_net(k)-D*SOC2*SOCmax;
                        SOCx=SOC1+ee/SOCmax;
                        SOC2=SOCx;
                        SOC(k)=SOC2;
                        w=w+1;
                    elseif w==1
                        if SOC(k-1)>SOC3
                            I_Discharge(k)=I_L(k)-I_PV(k);
                            I_Battery(k)=I_Discharge(k);
                            I_Load(k)=I_L(k);
                            B=SOC2;
                            V1=(1.926+0.124*B)*ns;
                            R1=(0.19+(0.1037/(B-0.14)))*(ns/SOCmax);
                            %syms v;
                            ee=K.*V1.*I_net(k)-D*SOC2*SOCmax;
                            SOCx=SOC1+ee/SOCmax;
                            SOC2=SOCx;
                            SOC(k)=SOC(k-1)-abs(SOC1-SOC2);
                        elseif SOC(k-1)<=SOC3
                            SOC2=SOC3;
                            SOC(k)=SOC3;
                            I_Deficit(k)=I_net(k);
                        end
                    end
                end
                    % In Case of the battery is empty
            else
                % In Case of I_PV=I_Load
                if I_net(n)==0
                    SOC(k)=SOC1;
                    I_Load(n)=I_PV(n);
                    I_Damp(n)=0;
                    I_Charge(n)=0;
                    I_Deficit(n)=0;
                    I_Discharge(n)=0;
                    I_Battery(n)=0;
                    %In Case of I_PV>I_Load
                elseif I_net(n)>0
                    SOC(k)=SOC1;
                    I_Load(n)=I_L(n);
                    I_Damp(n)=I_net(n);
                    I_Charge(n)=0;
                    I_Deficit(n)=0;
                    I_Discharge(n)=0;
                    I_Battery(n)=0;
                    %(Routine 6.3):In Case of I_PV<I_Load
                elseif I_net(n)<0
                    SOC(k)=SOC1;
                    I_Damp(n)=0;
                    I_Charge(n)=0;
                    I_Deficit(n)=I_L(n)-I_PV(n);
                    I_Discharge(n)=0;
                    I_Battery(n)=0;
                    I_Load(n)=I_PV(n);
                end
            end
        end
        E_Excess=sum(I_Damp.*230);
        SOC_per=SOC./SOCmax;
        C_batteryf(x,y)=n; %Battery Capacity for each hour and PV module
        C_PVf(x,y)=m; %Number of PV modules for each hour
        LLP_calculated(x,y)=abs(sum(I_Deficit)/sum(I_L));
        y=y+1;
    end
    x=x+1;
end
cc=1;
LLP_ff=zeros(1,size(LLP_calculated,1)*size(LLP_calculated,1));
C_PV_ff=zeros(1,size(LLP_calculated,1)*size(LLP_calculated,1));
C_battery_ff=zeros(1,size(LLP_calculated,1)*size(LLP_calculated,1));
for ii=1:size(LLP_calculated,1)
    for jj=1:size(LLP_calculated,1)
        if LLP_calculated(ii,jj)>=0.0095 && LLP_calculated(ii,jj)<=0.0105
            LLP_ff(cc)=LLP_calculated(ii,jj);
            C_PV_ff(cc)=C_PVf(ii,jj);
            C_battery_ff(cc)=C_batteryf(ii,jj);
            cc=cc+1;
        end
    end
end
%%Initialization of the cost Function
%PV
CC_PV=456; %Capital Cost for one PV
MC_PV=6.5; %Maintenance cost of one PV module per year
Ls=25; %The duration of operation of the system in years
L_PV=25; %The total lifetime period for PV array
%Battery
Ca_battery=1200; %Capacity for one battery (Wh)
CC_batwh=4.8; %Capital Cost for Wh
CC_bat=CC_batwh*Ca_battery; %Capital Cost for one battery
MC_bat=3.4; %Maintenance cost of one battery per year
L_bat=5; %The total lifetime period for battery
Y_bat=(Ls/L_bat)-1; %The expected numbers of the storage battery replacement during the system lifetime
B_rep=50; %Replacement cost for one battery
CC_B_rep=Y_bat*B_rep; %Cost of the storage battery replacement during the system lifetime
%Charge Controller
CC_cc=400; %Capital Cost for one Charge Controller
MC_cc=0; %Maintenance cost of one Charge Controller per year
L_cc=25; %The total lifetime period for Charge Controller
Y_cc=(Ls/L_cc)-1; %The expected numbers of the Charge
%Controller replacement during the system lifetime
N_cc=1; %Number of Charge Controllers during the system lifetime
%Inverter
CC_inv=800; %Capital Cost for one Inverter
MC_inv=0; %Maintenance cost of one Inverter per year
L_inv=25; %The total lifetime period for Inverter
Y_inv=(Ls/L_inv)-1; %The expected numbers of the Inverter replacement during the system lifetime
N_inv=1; %Number of inverters during the system lifetime
%Other Costs
%Circuit Breaker
N_CB=4; %Number of Circuit breaker
C_CB=25; %Cost for one circuit breaker
CC_CB=N_CB*C_CB; %Cost for all circuit breakers
%Support Structure
CC_SS=200;
%Civil work
CC_CW=400;
%Total cost for the Other Costs
CC_OC=CC_CB+CC_SS+CC_CW+CC_B_rep;
ir=0.035; %Real Interest Rate
fr=0.015; %Inflation Rate
ndr=((1+ir)/(1+fr))-1; %Net of discount?inflation rate
% Total life cycle cost calculation
LCCx=zeros(1,length(C_PV_ff));
CC_D=zeros(1,length(C_PV_ff));
LCC=zeros(1,length(C_PV_ff));
for kk=1:length(C_PV_ff)
    LCCx(kk)=(CC_OC/((((1+ndr)^Ls)-1)/((ndr*((1+ndr)^Ls)))))+((C_PV_ff(kk)*(CC_PV+Ls*MC_PV))/L_PV)+(((ceil(C_battery_ff(kk)/Ca_battery)*CC_bat*(1+Y_bat))+(MC_bat*(LsY_bat)))/L_bat)+(((N_cc*CC_cc*(1+Y_cc))+(MC_cc*(Ls-Y_cc)))/L_cc)+(((N_inv*CC_inv*(1+Y_inv))+(MC_inv*(Ls-Y_inv))/L_inv));
    CC_D(kk)=(C_PV_ff(kk)*(CC_PV))+((ceil(C_battery_ff( kk ) / Ca_battery ) * CC_bat * ( 1 + Y_rbat ) ) ) +(N_cc*CC_cc*(1+Y_cc))+(N_inv*CC_inv*(1+Y_inv))+CC_CB+CC_SS;
    LCC(kk) = LCCx(kk)-(((0.13*(CC_D(kk))))/((((1+ndr)^Ls)-1)/((ndr*((1+nd)^Ls)))));
end
%The optimum sizes of PV and battery combination
[MM,II]=min(LCC);
MM;
C_PV_best= C_PV_ff(II);
C_battery_best= C_battery_ff(II);
fprintf('Best PV Moduls Number is: %d\n',C_PV_best);
fprintf('Best Battery Capacity is: %d\n',C_battery_best);