%%%%%%%%%%%%%%%%每天的光照強度%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%以正午為例%%%%%%%%%%%%%%%%%%%
clear?all
clc
alpha00=1:1:90;
alpha0=alpha00*pi/180;
n0=length（alpha0）;
fai=40.1*pi/180;%緯度
gama=0*pi/180;%房屋方位角
Isc=1353;
AA=xlsread（‘工作表.xlsx‘）;
n=length（AA（:1））;
I_b01=AA（:4）-AA（:5）;
I_d01=AA（:5）;
I_b0=zeros（3651）;
I_d0=zeros（3651）;
for?i=1:365
for?j=1:24
????I_b0（i）=I_b0（i）+I_b01（（i-1）*24+j）;
????I_d0（i）=I_d0（i）+I_d01（（i-1）*24+j）;
end
end
for?ii=1:n0
????ii
%alpha=10*pi/180;%斜面偏角
alpha=alpha0（ii）;
for?i=1:365????

???deta（i）=pi/180*23.45*sin（（2*pi*（284+i））/365）;%
???w0（i）=acos（-tan（fai）*tan（deta（i）））;%水平面上日落時
???I0（i）=24/pi*Isc*（1+0.033*cos（360*i/365））*（cos（fai）*cos（deta（i））*sin（w0（i））+2*pi*w0（i）/360*sin（fai）*sin（deta（i）））;%
???%%%大氣外水平面太陽輻射強度
???w（i）=0;
???k=i;
????C_theta（k）=