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;
??