光伏电池输出特性曲线MATLAB代码

绘制不同光照强度下光伏电池输出特性曲线的MATLAB程序

clear
t=0:0.1:60;

S1=1000;                                         %一组不同的光照强度；
S2=1500;
S3=2000;
S4=2500;

T=25;                                            %设定温度；

Sref=1000;                                       %参考光照强度与参考温度；
Tref=25;

a=0.0001;b=0.0005;c=0.008;                     %补偿系数；

e=2.71828;

Vocref=44;Iscref=8.0;Vmref=35;Imref=7.5;   %太阳能电池板参数（BP生产商提供）；




Isc1=Iscref*S1/Sref*[1+a*（T-Tref）];               %工程用数学模型；
Voc1=Vocref*log（e+b*（S1-Sref））*[1-c*（T-Tref）];

Isc2=Iscref*S2/Sref*[1+a*（T-Tref）];
Voc2=Vocref*log（e+b*（S2-Sref））*[1-c*（T-Tref）];

Isc3=Iscref*S3/Sref*[1+a*（T-Tref）];
Voc3=Vocref*log（e+b*（S3-Sref））*[1-c*（T-Tref）];

Isc4=Iscref*S4/Sref*[1+a*（T-Tref）];
Voc4=Vocref*log（e+b*（S4-Sref））*[1-c*（T-Tref）];



Im1=Imref*S1/Sref*[1+a*（T-Tref）];
Vm1=Vmref*log（e+b*（S1-Sref））*[1-c*（T-Tref）];

Im2=Imref*S2/Sref*[1+a*（T-Tref）];
Vm2=Vmref*log（e+b*（S2-Sref））*[1-c*（T-Tref）];

Im3=Imref*S3/Sref*[1+a*（T-Tref）];
Vm3=Vmref*log（e+b*（S3-Sref））*[1-c*（T-Tref）];

Im4=Imref*S4/Sref*[1+a*（T-Tref）];
Vm4=Vmref*log（e+b*（S4-Sref））*[1-c*（T-Tref）];


C21=（Vm1/Voc1-1）/log（1-Im1/Isc1）;
C11=（1-Im1/Isc1）*exp（-Vm1/Voc1/C21）;

C22=（Vm2/Voc2-1）/log（1-Im2/Isc2）;
C12=（1-Im2/Isc2）*exp（-Vm2/Voc2/C22）;

C23=（Vm3/Voc3-1）/log（1-Im3/Isc3）;
C13=（1-Im3/Isc3）*exp（-Vm3/Voc3/C23）;

C24=（Vm4/Voc4-1）/log（1-Im4/Isc4）;
C14=（1-Im4/Isc4）*exp（-Vm4/Voc4/C24）;


i1=Isc1-Isc1*C11*[exp（t/C21/Voc1）-1];
i2=Isc2-Isc2*C12*[exp（t/C22/Voc2）-1];
i3=Isc3-Isc3*C13*[exp（t/C23/Voc3）-1];
i4=Isc4-Isc4*C14*[exp（t/C24/Voc4）-1];



k1=i1+abs（i1）;
k2=i2+abs（i2）;
k3=i3+abs（i3）;
k4=i4+abs（i4）;


plot（tk1tk2tk3tk4）;                    %plot（tk1tk2tk3tk4）;



xlabel（‘输出电压（v）‘‘FontSize‘16‘FontWeight‘‘bold‘）; 
ylabel（‘输出电流（A）‘‘FontSize‘16‘FontWeight‘‘bold‘）; 
title（‘光伏电池特性随光照的变化‘‘FontSize‘16‘FontWeight‘‘bold‘）; 
legend（‘1000W/㎡‘‘1500W/㎡‘‘2000W/㎡‘‘2500W/㎡‘‘Location‘ ‘East‘）
grid;  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2041  2014-04-25 18:03  PVguangzhao.m

----------- ---------  ---------- -----  ----

                 2041                    1