function I = solar(u)
%#eml
I=((1367*(1+0.034*cos(2*pi*172/365))*0.75^(1/(sin(30.6167*pi/180)*sin(23.45*pi/180)+...
cos(30.6167*pi/180)*cos(23.45*pi/180)*cos(u*pi/180))))*((sin(30.6167*pi/180)*cos(30*pi/180)-sin(30*pi/180)*cos(30.6167*pi/180))*sin(23.45*pi/180)+...
(cos(30.6167*pi/180)*cos(30*pi/180)+sin(30*pi/180)*sin(30.6167*pi/180))*cos(23.45*pi/180)*cos(u*pi/180))+(0.5*1367*(1+0.034*cos(2*pi*172/365))*((1-0.75^(1/(sin(30.6167*pi/180)*sin(23.45*pi/180)+cos(30.6167*pi/180)*cos(23.45*pi/180)*cos(u*pi/180))))/(1-1.4*log(0.75)))*(sin(30.6167*pi/180)*sin(23.45*pi/180)+...
cos(30.6167*pi/180)*cos(23.45*pi/180)*cos(u*pi/180)))*0.933+0.2*((1367*(1+0.034*cos(2*pi*172/365))*0.75^(1/(sin(30.6167*pi/180)*sin(23.45*pi/180)+...
cos(30.6167*pi/180)*cos(23.45*pi/180)*cos(u*pi/180))))*(sin(30.6167*pi/180)*sin(23.45*pi/180)+cos(30.6167*pi/180)*cos(23.45*pi/180)*cos(u*pi/180))+...
0.5*1367*(1+0.034*cos(2*pi*172/365))*((1-0.75^(1/(sin(30.6167*pi/180)*sin(23.45*pi/180)+...
cos(30.6167*pi/180)*cos(23.45*pi/180)*cos(u*pi/180))))/(1-1.4*log(0.75)))*(sin(30.6167*pi/180)*sin(23.45*pi/180)+...
cos(30.6167*pi/180)*cos(23.45*pi/180)*cos(u*pi/180)))*0.067)*0.682;
for u=-180:15:180
end
那位高手懂,能给我指点一下呢,公式这么长但是应该没有问题。