www.gusucode.com > 太阳能小屋的设计-数学建模【论文+程序】 > 2012高教社杯全国大学生数学建模竞赛太阳能小屋的设计【论文+程序】/文章+源码/B149源程序/M/qiangdu(!).m
clear;clc fushe=xlsread('cumcm.xls','sheet','C4:K8763'); w=[];dta=[]; for i=1:8760 w(i)=15*(fushe(i,1)-12)/57.3;%时角、弧度 n=ceil(i/24); dta(i)=23.45*sin(2*pi*(284+n)/365)/57.3;%赤纬角、弧度 sa(i)=sin(0.7)*sin(dta(i))+cos(0.7)*cos(dta(i))*cos(w(i));%太阳高度角、正弦 cA(i)=(sin(dta(i))-sa(i)*sin(0.7))/((sqrt(1-sa(i)^2))*cos(0.7));%太阳方位角、余弦 if w(i)<0 A(i)=acos(cA(i)); else A(i)=2*pi-acos(cA(i)); end sb(i)=sin(A(i)-0.5*pi)*sqrt(1-sa(i)^2); re(i)=fushe(i,7)-0.5*fushe(i,4)-fushe(i,5)*sb(i); end