太阳电池matlab模拟program1.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * 第 5章 太陽電池MATLAB模擬 Program 1: 太陽能電池之輸出電壓、電流、功率響應 Program 2: The output response of buck converter with load variation Program 3: Solar cell with Buck converter Program 1 太陽的日照強度=1 kW/m2時，且溫度25度C ， 當t=0~3000, 負載R=R1=6, 當t=30000~7000, 負載 R=R2=7, 當t=7000~10000, 負載R=R3=8。 畫出太陽能電池之: 輸出電壓響應圖(Fig.1)， 輸出電流響應圖(Fig.2) ， 輸出功率響應圖(Fig.3) 。 %太陽的日照強度=1 kW/m2)時，且溫度25度C ，(t=0~3000, 負載R=R1=6), (t=30000~7000, R=R2=7), (t=7000~10000, R=R3=8) %太陽能電池之輸出電壓響應圖(Fig.1)，輸出電流響應圖(Fig.2), 輸出功率響應圖(Fig.3) clear; clc; %參數調整% np = 5; %並聯個數 ns = 10; %串聯個數 k = 1.38*(10^-23); %波茲曼常數 q = 1.63*(10^-19); %一個電子所含之電荷量 T = 298; %太陽能表面溫度 Tr = 298; %參考溫度 A =1.8; %太陽能板理想因數 Irr = 0.3*(10^-19); %逆向飽和電流 Eg = 1.1; %半導體能量間隙所需之能量 Isc = 0.66; %太陽能板在參考溫度和1kW/m^2的日照條件下的短路電流值 Ki = 1.33; %太陽能板短路電流的溫度係數 Si = 1.0; %日到強度 1 dv=0.01; t1=3000; t2=7000; MATLAM Program (Photovoltaicwithload1) R1=6; R2=7; R3=8; V(1) =18; Iph = (Isc+Ki*(T-Tr)/1000)*Si; %公式(3) Isat = Irr*((T/Tr)^3)*exp(q*Eg*((1/Tr)-(1/T))/(k*A)); %公式(2) for t=1:10000 if tt1 R=R1; else if tt2 R=R2; else R=R3; end end I(t) = np*Iph-np*Isat*exp(((q*V(t))/(k*A*T*ns))-1); %公式(1) IL(t)=V(t)/R; if IL(t)I(t) V(t+1)=V(t)+dv; else if IL(t)I(t) V(t+1)=V(t)-dv; else V(t+1)=V(t); end end P(t) = V(t)*I(t); %公式(4) end %電壓響應圖% t = 1:10000; figure(1) plot(t,V(t),-b) legend(V(t),1) xlabel(t) ylabel(V (volt)) hold on %電流響應圖% figure(2) plot(t,I(t),-r,t,IL(t),-b) legend(I(t),IL(t),1) xlabel(t) ylabel(I(amp)) hold on %功率響應圖% figure(3) plot(t,P(t),-r) legend(P(t),1) xlabel(t) ylabel(P(W)) hold on Figure 1 Figure 2 Figure 3 Program 2