数学建模b题程序.docVIP

附件1 %屋顶上的各月的每小时太阳辐照总量 %一月份 beta=atan(3/16); q=40.1.*(pi/180);%fei=q k=(0.5*(1+cos(beta))); Io=15.3; rho=0.2; B=xlsread(book2);It1=0; for n=1:31; s=23.45.*sin((2*pi*(284+n))/365).*(pi/180);%delta gama=1+0.033.*cos(2*pi*n/365); I=B(24*(n-1)+1:24*n,1); Id=B(24*(n-1)+1:24*n,2); In=B(24*(n-1)+1:24*n,3); for i=1:24 t=i-1; ome=15.*(t-12).*(pi/180); ap=sin(q)*sin(s)+cos(q)*cos(s)*cos(ome)+(pi/2);%sin(alphas) Ib=In(1,i)*ap; if I(1,i)==0; f=1; else f=sqrt(Ib/I(1,i)); end A=Ib/Io; Rb=(cos(q-beta)*cos(s)*cos(ome)+sin(q-beta)*sin(s))./(sin(q)*sin(s)+cos(q)*cos(s)*cos(ome)); % It1=It1+(Ib+Id(1,i).*A)*Rb+Id(1,i).*(1-A)*k*(1+f*sin(0.5*beta))^3+0.2*c; It1=It1+Ib*Rb+Id(1,i).*k+I(1,i)*k; end end It1 %二月份 beta=atan(3/16); q=40.1.*(pi/180);%fei=q k=(0.5*(1+cos(beta))); Io=20.3; rho=0.2; B=xlsread(book2);It2=0; for n=32:59; s=23.45.*sin((2*pi*(284+n))/365).*(pi/180);%delta gama=1+0.033.*cos(2*pi*n/365); I=B(24*(n-1)+1:24*n,1); Id=B(24*(n-1)+1:24*n,2); In=B(24*(n-1)+1:24*n,3); for i=1:24 t=i-1; ome=15.*(t-12).*(pi/180); ap=sin(q)*sin(s)+cos(q)*cos(s)*cos(ome)+(pi/2);%sin(alphas) Ib=In(1,i)*ap; if I(1,i)==0; f=1; else f=sqrt(Ib/I(1,i)); end A=Ib/Io; Rb=(cos(q-beta)*cos(s)*cos(ome)+sin(q-beta)*sin(s))./(sin(q)*sin(s)+cos(q)*cos(s)*cos(ome)); % It2=It2+(Ib+Id(1,i).*A)*Rb+Id(1,i).*(1-A)*k*(1+f*sin(0.5*beta))^3+0.2*I(1,i)*k; It2=It2+Ib*Rb+Id(1,i).*k+I(1,i)*k; end end It2 %三月份 beta=atan(3/16); q=40.1.*(pi/180);%fei=q k=(0.5*(1+cos(beta))); Io=27.4; rho=0.2; B=xlsread(book2);It3=0; for n=60:90; s=23.45.*sin((2*pi*(284+n))/365).*(pi/180);%delta gama=1+0.033.*cos(2*pi*n/365); I=B(24*(n-1)+1:24*n,1); Id=B(24*(n-1)+1:24*n,2); In=B(24*(n-1)+1:24*n,3); for i=1:24 t=i-1; ome=15.*(t-12).*(pi/180); ap=sin(q)*sin(s)+cos(q)*cos(s)*cos(ome)+(pi/2);%sin(alphas) Ib=In(1,i)*ap; if I(1,i)==0; f=1; else f=sqrt(Ib/I(1,i)); end A=Ib/Io; Rb=(cos(q-beta)*cos(s)*cos(ome)+s