www.gusucode.com > 太阳能小屋的设计-数学建模【论文+程序】 > 2012高教社杯全国大学生数学建模竞赛太阳能小屋的设计【论文+程序】/文章+源码/B149源程序/M/sun_jisuan_day_1.m

    
%%%%%%%%%%%%%%%%每天的光照强度%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%以正午为例%%%%%%%%%%%%%%%%%%%
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(365,1);
I_d0=zeros(365,1);
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)=(sin(fai)*cos(alpha)-cos(fai)*sin(alpha)*cos(gama))*sin(deta(k))+(cos(fai)*cos(alpha)+sin(fai)*sin(alpha)*cos(gama))*cos(deta(k))*cos(w(k))+sin(alpha)*sin(gama)*cos(deta(k))*sin(w(k));
    C_theta0(k)=sin(fai)*sin(deta(k))+cos(fai)*cos(deta(k))*cos(w(k));
    I_b(k)=I_b0(k)*(C_theta(k)/C_theta0(k));   
    I_d(k)=I_d0(k)*(I_b0(k)*(C_theta(k)/C_theta0(k))/I0(k)+0.5*(1-I_b0(k)/I0(k))*(1+cos(alpha)));
    I(k)=I_b(k)+I_d(k);
end
w=w'*pi/180;%时角
   deta=deta';%赤纬角
   

BB=[31,28,31,30,31,30,31,31,30,31,30,31];
CC=[0,31,59,90,120,151,181,212,243,273,304,334];
I_m=zeros(12,1);
for i=1:12
    for j=1:BB(i)
        I_m(i)=I_m(i)+I(CC(i)+j);
    end
end
II_m(:,ii)=I_m;
end
II_m=II_m';
plot(alpha00,II_m);