太阳板光照强度与板倾斜角和方位角关系的MAtlab程序

%%%%%%%%%%%%%%%%每天的光照强度%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%以正午为例%%%%%%%%%%%%%%%%%%%
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）=