Scribd
Upload
8 views2 pages

User and GPT Code

The MATLAB code calculates the total solar radiation on a flat plate collector based on user inputs for the day of the year, latitude, longitude, tilt angle, and azimuth angle. It uses various constants and formulas to compute solar angles and radiation values for each hour of the day, ultimately plotting the total radiation against the hour. The code includes necessary conversions and calculations to ensure accurate results.

Uploaded by

shamikhahmed7
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
8 views2 pages

User and GPT Code

The MATLAB code calculates the total solar radiation on a flat plate collector based on user inputs for the day of the year, latitude, longitude, tilt angle, and azimuth angle. It uses various constants and formulas to compute solar angles and radiation values for each hour of the day, ultimately plotting the total radiation against the hour. The code includes necessary conversions and calculations to ensure accurate results.

Uploaded by

shamikhahmed7
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

MATLAB Code - User's Original (Corrected)

clear; clc;
format short;

%% USER INPUT
n = input('Enter day of the year (1 to 365): ');
phi = input('Enter the latitude of the location (in degrees): ');
LL = input('Enter the longitude of the location (in degrees): ');
beta = input('Enter tilt angle of collector (in degrees): ');
gamma = input('Enter azimuth angle of collector (in degrees): ');

%% CONSTANTS & CONVERSIONS


rho = 0.2;
phi = deg2rad(phi);
beta = deg2rad(beta);
gamma = deg2rad(gamma);
delta = deg2rad(23.45 * sind((360/365)*(284 + n)));
B = deg2rad((360/365)*(n - 1));
E = 229.2 * (0.000075 + 0.001868*cos(B) - 0.032077*sin(B) ...
- 0.014615*cos(2*B) - 0.04089*sin(2*B));
GMT = ceil(LL/15);
SL = GMT * 15;
ST_LT = (E - 4 * (SL - LL)) / 60;

hours = 0:23;
omega = zeros(1, 24);
cos_theta = zeros(1, 24);
alpha_s = zeros(1, 24);
theta_z = zeros(1, 24);
m = zeros(1, 24);
A = zeros(1, 24);
k = zeros(1, 24);
IB = zeros(1, 24);
IBC = zeros(1, 24);
IDC = zeros(1, 24);
IRC = zeros(1, 24);
IC = zeros(1, 24);

for h = 1:24
ST = hours(h) + ST_LT;
omega(h) = deg2rad((ST - 12) * 15);
cos_theta(h) = sin(delta)*sin(phi)*cos(beta) ...
- sin(delta)*cos(phi)*sin(beta)*cos(gamma) ...
+ cos(delta)*cos(phi)*cos(beta)*cos(omega(h)) ...
+ cos(delta)*sin(phi)*sin(beta)*cos(gamma)*cos(omega(h));
cos_theta(h) = max(min(cos_theta(h), 1), -1);
alpha_s(h) = asin(cos(phi)*cos(delta)*cos(omega(h)) + sin(phi)*sin(delta));
theta_z(h) = acos(cos(phi)*cos(delta)*cos(omega(h)) + sin(phi)*sin(delta));
if alpha_s(h) > 0
m(h) = 1 / (sin(alpha_s(h)) + 0.50572 * (6.07995 + rad2deg(alpha_s(h)))^(-1.6364));
else
m(h) = Inf;
end
A(h) = 1160 + 75 * sin(deg2rad((360/365)*(n - 275)));
k(h) = 0.174 + 0.035 * sin(deg2rad((360/365)*(n - 100)));
IB(h) = A(h) * exp(-k(h) * m(h));
IBC(h) = IB(h) * cos_theta(h);
if IBC(h) < 0
IBC(h) = 0;
end
C = 0.095 + 0.04 * sin(deg2rad((360/365)*(n - 100)));
IDC(h) = C * IB(h) * (1 + cos(beta)) / 2;
IRC(h) = rho * IB(h) * (sin(alpha_s(h)) + C) * (1 - cos(beta)) / 2;
IC(h) = IBC(h) + IDC(h) + IRC(h);
end

figure;
plot(hours, IC, 'r-', 'LineWidth', 2);
xlabel('Hour of Day');
ylabel('Total Radiation on Collector (W/m^2)');
title(['Total Solar Radiation on Flat Plate Collector (Day ' num2str(n) ')']);
grid on;

You might also like