کدهای رایگان مرتبط با سیستم های خورشیدی
پروژه ۹۰۴: شبیه سازی رایگان در متلب
انرژی خورشیدی یکی از منابع انرژیهای تجدیدپذیر و از مهمترین آنها می باشد. میزان تابش انری خورشیدی در نقاط مختلف جهان متغیر بوده و در کمربند خورشیدی زمین بیشترین مقدار را داراست. کشور ایران نیز در نواحی پرتابش واقع است و مطالعات نشان می دهد که استفاده از تجهیزات خورشیدی در ایران مناسب بوده و میتواند بخشی از انرژی مورد نیاز کشور را تأمین نماید. در پروژه حاضر شبیه سازی ولتاژ و جریان و همچنین منحنی های مربوط به سیستم های خورشیدی در متلب انجام گردیده است.
در صورت سوال در مورد محصول می توانید از بخش تماس با ما (منوی بالا)، با شماره تلفن مورد نظر مکاتبه نمایید.
برنامه شماره ۱:
clc, clear all
disp( ‘ ‘)
disp( ‘=========================== ‘)
disp( ‘ Program for calculating PV module current under ‘)
disp( ‘ given Voltage, Irradiance and Temperature ‘)
disp( ‘============================= ‘)
% Defining input variables
Va = input(‘ Enter Module operating voltage (V):’);
B = input(‘ Enter Irradiance : ‘);
Tw=input( ‘ Enter Module temperature in deg C : ‘ );% Tac
Ns = input(‘ Enter number of series connected cells :’);
Np=input(‘ Enter number of parallel connected cells :’);1;
Vo1=input( ‘ Enter the Open circuit voltage, Voc, at reference temp(T1=25):’);
Vo2=input( ‘ Enter the Open circuit voltage, Voc, At temp T2 : ‘);
Isc_T1=input(‘Enter the short circuit Current per cell at reference temp T1:’);
Isc_T2 = input(‘ Enter the short circuit Current per cell at temp T2 : ‘);
t2 = input(‘ Enter the given temp t2 : ‘);
% Calculating cell voltage
Vc = Va/Ns;
K = 1.38e-23; % Boltzmann’s constant
q = 1.60e-19; % Electron charge
N = 1.2; % “diode quality” factor
Vg = 1.12;
Eg=1.12; % the band gap energy of the semiconductor
% Calculating open circuit voltage per cell at temperatures T1 and T2
Vo_T1 = Vo1 /Ns;
Vo_T2 = Vo2 /Ns;
T1 = 273 + 25;
T2 = 273 + t2;Tref=T1;
Tac = 273 + Tw; % array working temp.
% Calculating photon generated current @ given irradiance
Iph_T1 = Isc_T1 * B;
% Calculate short-circuit current for Tac
a = (Isc_T2 – Isc_T1)/Isc_T1 * 1/(T2 – T1);
Iph = Iph_T1 * (1 + a*(Tac – T1));
% Calculating the thermo voltage of cell at reference temp
Vt_T1 = K * T1 / q; % = A * kT/q
% Calculating reverse saturation current for given temperature
Ir_T1 = Isc_T1 / (exp(Vo_T1/(N*Vt_T1))-1);
Ir_T2 = Isc_T2 / (exp(Vo_T2/(N*Vt_T1))-1);
b = Eg * q/(N*K);
Ir = Ir_T1 * (Tac/T1).^(3/N) .* exp(-b.*(1./Tac – 1/T1));
% Calculating series resistance per cell
X2v = Ir_T1/(N*Vt_T1) * exp(Vo_T1/(N*Vt_T1));
dVdI_Voc = – 1.15/Ns / 2; % dV/dI at Voc per cell, from manufacturers graph
Rs = – dVdI_Voc – 1/X2v; % series resistance per cell
Vt_Ta = N*K*Tac / q; % = A * kT/q
Ia = zeros(size(Vc));
disp( ‘============================= ‘)
for j=1:6;
Ia = Ia – (Iph – Ia – Ir.*( exp((Vc+Ia.*Rs)./Vt_Ta) -1))./ (-1 – (Ir.*(exp((Vc+Ia.*Rs)./Vt_Ta) -1)).*Rs./Vt_Ta);
end
———————————
برنامه شماره ۲:
function Ia = pviv( Va,B,Tw )
% pviv
% Built-in function for calculating PV module current under ‘)
% given Voltage Va, Irradiance B, and Temperature Tw
% pviv function requires defining cell or module electric
% characteristics
clc
disp( ‘ ‘)
disp( ‘============================ ‘)
disp( ‘ Built-in function for calculating PV module current under ‘)
disp( ‘ given Voltage, Irradiance and Temperature ‘)
disp( ‘============================= ‘)
% Defining input variables
Ns = input(‘ Enter number of series connected cells :’);
Np=input(‘ Enter number of parallel connected cells :’);1;
Vo1=input(‘Enter the Open circuit voltage, Voc, at reference temp (T1=25) :’);
Vo2=input( ‘ Enter the Open circuit voltage, Voc, At temp T2 : ‘);
Isc_T1 = input(‘ Enter the short circuit Current per cell at reference temp T1 : ‘);
Isc_T2 = input(‘ Enter the short circuit Current per cell at temp T2 : ‘);
t2 = input(‘ Enter the given temp t2 : ‘);
% Calculating cell voltage
Vc = Va/Ns;
K = 1.38e-23; % Boltzman’s constant
q = 1.60e-19; % Electron charge
N = 1.2; % “diode quality” factor
Vg = 1.12;
Eg=1.12; % the band gap energy of the semiconductor
% Calculating open circuit voltage per cell at temperatures T1 and T2
Vo_T1 = Vo1 /Ns;
Vo_T2 = Vo2 /Ns;
T1 = 273 + 25;
T2 = 273 + t2;Tref=T1;
Tac = 273 + Tw; % array working temp.
% Calculating photon generated current @ given irradiance
Iph_T1 = Isc_T1 * B;
% Calculate short-circuit current for TaK
a = (Isc_T2 – Isc_T1)/Isc_T1 * 1/(T2 – T1);
Iph = Iph_T1 * (1 + a*(Tac – T1));
% Calculating the thermo voltage of cell at reference temp
Vt_T1 = K * T1 / q; % = A * kT/q
% Calculating reverse saturation current for given temperature
Ir_T1 = Isc_T1 / (exp(Vo_T1/(N*Vt_T1))-1);
Ir_T2 = Isc_T2 / (exp(Vo_T2/(N*Vt_T1))-1);
b = Eg * q/(N*K);
Ir = Ir_T1 * (Tac/T1).^(3/N) .* exp(-b.*(1./Tac – 1/T1));
% Calculating series resistance per cell
X2v = Ir_T1/(N*Vt_T1) * exp(Vo_T1/(N*Vt_T1));
dVdI_Voc = – 1.15/Ns / 2; % dV/dI at Voc per cell, from manufacturers graph
Rs = – dVdI_Voc – 1/X2v; % series resistance per cell
Vt_Ta = N*K*Tac / q; % = A * kT/q
Ia = zeros(size(Vc));
disp( ‘=========================== ‘)
for j=1:6;
Ia = Ia – (Iph – Ia – Ir.*( exp((Vc+Ia.*Rs)./Vt_Ta) -1))./ (-1 – (Ir.*(
exp((Vc+Ia.*Rs)./Vt_Ta) -1)).*Rs./Vt_Ta);
end
end
برنامه شماره ۳:
clear;B = 1;
figure, hold on
for TaC=0:20:80
Va = linspace (0, 48-TaC/8, 200);
Ia = msx60i(Va, B, TaC);
plot(Va, Ia)
end
title(‘ Photovoltaic Module I-V Curve’)
xlabel(‘Module Voltage (V)’)
ylabel(‘Module Current (A)’)
axis([0 50 0 5])
gtext(‘0C’)
gtext(’20C’)
gtext(’40C’)
gtext(’60C’)
gtext(’80C’)
figure, hold on
for TaC=0:20:80
Va = linspace (0, 48-TaC/8, 200);
Ia = msx60i(Va, B, TaC);
Pa= Va.*Ia;
plot(Va, Pa)
end
title(‘ Photovoltaic Module P-V Curve’), xlabel(‘Module Voltage(V)’),ylabel(‘Module Current (A)’),axis([0 25 0 80]), gtext(‘0C’),
gtext(’20C’), gtext(’40C’), gtext(’60C’), gtext(’80C’)
برای عضویت روی عکس زیر کلیک کنید : (آخرین اخبار مرتبط با مهندسی و سایر تکنولوژی ها)
یا آدرس لینک زیر را در تلگرام خود جستجو نمایید:
m_b_coll@
مجموعه: انرژی های نو