% fai=input('fai:')
% B=input('B:')
%
% Qyear0=0;
% Qyear1=0;
% s4=0;
% s2=0;
% s5=0;
% N=365;
% for day=1:N
%
% %大气层外的太阳辐射强度
% I=1367*(1+0.033*cos(2*pi*day/365));
% %I=2;
% % 赤纬角
% delta=23.45*sin(2*pi*(284+day)/365);
% %日出日落时角
% w=(acos(-tan((fai-B)*pi/180)*tan(delta*pi/180)))*180/pi;%%%%%%%%%%%%%%%%%
% %日出时间,小时
% Tchu=12-w/15;
% %日出时间,小时
% Tluo=12+w/15;
%
% Qday0=0;
% Qday1=0;
% s1=0 ;
% s3=0;
% n=1000;
% h=(Tluo-Tchu)/n;
%
%
% for i= 1:n; %%%%%%%%%%
% time =Tchu+(i-1)* h;%%%%%%%%%%%%%%%%%%%%%%
%
% %time=12;
% %入射角
% w0=15*time-180;
% % w0=0;
% rushejiao=acos(sin((fai-B)*pi/180)*sin(delta*pi/180)+cos((fai-B)*pi/180)*cos(delta*pi/180)*cos(w0*pi/180))*180/pi;
%
% %高度角
%
% gaodujiao=asin(sin(fai*pi/180)*sin(delta*pi/180)+cos(fai*pi/180)*cos(delta*pi/180)*cos(w0*pi/180))*180/pi;
% % rs=asin(cos(delta*pi/180)*sin(w*pi/180)/cos(gaodujiao*pi/180)); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%xinjia
% % Y=1.64*sin(B*pi/180)*cos(rs*pi/180)/tan(gaodujiao*pi/180); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% %太阳直射透明度系数
% AM0=((1229+(614*sin(gaodujiao*pi/180))^2)^0.5)-614*sin(gaodujiao*pi/180);
% AMh=AM0*0.9996;
%
% % ph0=((288-0.0065*H)/288)^5.256;
% % AMh=AM0*(ph0);
% % %
%
% tb=0.56*((2.71828^(-0.56*AMh))+(2.71828^(-0.0095*AMh)));
% td=0.271-0.294*tb; %散射光穿透系数
%
% t1=I*tb*cos(rushejiao/180*pi); %瞬时光直射强度
% t2=I*td*(cos(B*pi/180))^2*sin(gaodujiao*pi/180);%瞬时光散射强度,倾角
% %补偿
% t3=t1+t2;%瞬时光照强度
%
% % Im=0.008769*t3*(1+0.00255*(20-25));%瞬时最大功率点处电流%%%
% % Um=31.95*(1-0.00285*(20-25))*log(2.71828+0.55*(t3/1000-1));%瞬时最大功率点处电压
% % Pm1=Im*Um*0.001*180;%瞬时一组串(含9块光伏板)的发电功率
% % Pm2=Im*Um*0.001*200;
% %
% % Qday0=Qday0+Pm1*h;
% % Qday1=Qday1+Pm2*h;
%
% s1=s1+t1*h;
% s3=s3+t2*h;
% end
%
% s2=s2+s1; %一年光直射总辐射量
% s4=s4+s3; %一年光散射总辐射量
% s5=s2+s4; %一年光总辐射量
% % Qyear0=Qyear0+Qday0; %上限
% % Qyear1=Qyear1+Qday1; %上限
% % if N==2
% % f1=abs([s5]);
% % %g2=abs([Qyear3]);
% % else
% % f1=abs([f1 s5]);
% % %g2=abs([g2 Qyear3]);
% % end
% end
Asun=(9.5376*10^-6)*(0^2)-(5.3803*10^-6)*+0.00050652;%损失量