应用介绍

免费下载用于仿真PV + MPPT +部分阴影效果的matlab代码

function output=MPPT_PSO(input)

%%%----- PV_1 -----------------------

Tcell1 =input(1);% 25; %Ambient cell temperature (degrees celcius)

G1 = input(2);%1000%Solar Irradiation (W/m2)

Voc1 =43.2; % open circuit voltage at STC

Va1 = 0:0.1:Voc1;

Ipv1 = PVcharacteristics_func(Va1,G1,Tcell1);

factor1=Ipv1>=0;

Ipv1=factor1.*Ipv1;

Ppv1 = Va1.*Ipv1; 

%%%---------------

Pmax1 = max (Ppv1); %Finding the maximum power of the array

[row1,cumn1]=find(Ppv1<=Pmax1 & Ppv1>=Pmax1);

V_max1=Va1(row1,cumn1);

I_max1=Ipv1(row1,cumn1);

%%%----- PV_2 -----------------------

Tcell2 = input(3);%33; %Ambient cell temperature (degrees celcius)

G2 = input(4);%850;%Solar Irradiation (W/m2)

Voc2 = 43.2; %Open circuit voltage

Va2 = 0:0.1:Voc2;

Ipv2 = PVcharacteristics_func(Va2,G2,Tcell2);

factor2=Ipv2>=0;

Ipv2=factor2.*Ipv2;

Ppv2 = Va2.*Ipv2; %Calculating PV output power

%%%---------------

Pmax2 = max (Ppv2); %Finding the maximum power of the array

[row2,cumn2]=find(Ppv2<=Pmax2 & Ppv2>=Pmax2);

V_max2=Va2(row2,cumn2);

I_max2=Ipv2(row2,cumn2);

%%%----- PV_3 -----------------------

Tcell3 = input(5);%25; %Ambient cell temperature (degrees celcius)

G3 = input(6);%750%Solar Irradiation (W/m2)

Voc3 =43.2; %Open circuit voltage

Va3 = 0:0.1:Voc3;

Ipv3 = PVcharacteristics_func(Va3,G3,Tcell3);

factor3=Ipv3>=0;

Ipv3=factor3.*Ipv3;

Ppv3 = Va3.*Ipv3; %Calculating PV output power

%%%---------------

Pmax3 = max (Ppv3); %Finding the maximum power of the array

[row3,cumn3]=find(Ppv3<=Pmax3 & Ppv3>=Pmax3);

V_max3=Va3(row3,cumn3);

I_max3=Ipv3(row3,cumn3);

%%%%%------------ PV total calculation -------------------------------

PV_total=Va3.*(Ipv1+Ipv2+Ipv3);

Pmax_total = max (PV_total); %Finding the maximum power of the array

[row_total,cumn_total]=find(PV_total<=Pmax_total & PV_total>=Pmax_total);

V_max_total=Va3(row_total,cumn_total);

I_max_total=Ipv1(row_total,cumn_total)+Ipv2(row_total,cumn_total)+Ipv3(row_total,cumn_total);

%%---------------

output(1)=V_max_total(1);

output(2)=I_max_total(1);