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);