光伏发电I-V、P-V特性matlab程序源码

光伏发电I-V、P-V特性matlab程序源码 新能源发电技术课程资料 电流I满足的方程式： 其中， 表示并联太阳能电池个数， 表示串联太阳能电池个数， 表示太阳能电池产生的原始电流， 表示等效的二极管的泄漏电流， 为光伏板的内部串联电阻， 为光伏板的内部并联电阻，A是理想系数，q为元电荷量，K为玻尔兹曼常数。 采用牛顿迭代法

clear all
Ta = 25;                    %设定温度
G = 2;                      %设定光照强度

A = 1.72;                   %太阳能电池板参数（BP生产商提供）；
q = 1.6*10^（-19）;
k = 1.380658*10^（-23）;
Eg = 1.1;
Ior = 19.9693*10^（-6）;
Iscr = 3.3;
ki = 0.0017;
ns = 40;
np = 4;
Rs = 5*10^（-5）;
Rp = 5*10^5;
Tr = 301.18;

Tc = Ta + 0.2*G + 273.18;         %计算电池表面温度
Isc = （Iscr + ki*（Tc-Tr））*G;      %计算短路电流
Is = Ior*（（Tc/Tr）^3）*exp（q*Eg*（1/Tr-1/Tc）/（k*A））;      %计算泄漏电流

V = 0;                  %设置V和I初值
I = 0;

U = 0:0.1:30;           %建立坐标轴
N = length（U）;
value = cell（1N ）;     %创建数组用于存放I值

for j = 1:N             %进入电压值循环   
   while（1）             %进入迭代循环计算I值
      FI = I-np*Isc + np*Is*（exp（（q*（V/ns+I*Rs/np））/（A*k*Tc））-1）+（V*np/ns+I*Rs）/Rp;   %f（I）函数
      dFI = 1+（np*Is*Rs*q）/（np*A*k*Tc）*exp（（q*（V/ns+I*Rs/np））/（A*k*Tc））+Rs/Rp;        %f（I）的导数
      I1 = I - FI/dFI;           %牛顿迭代公式
      if abs（（I1-I）/I1）>0.01     %判断是否满足停止迭代条件
          I = I1;