【正文】 she=A(:,5)。Sanshe(:,2)=Sanshe(:,1).*(1+cosd(sita))/2。%求斜面散射輻射強(qiáng)度Zhishe(:,2)=Zhishe(:,1).*Arefa(:,3)。%求斜面直射輻射強(qiáng)度%%%%%%%%%%%%%%%%%%%%對法線方向修正%%%%%%%%%%%%%%%%%%%%%%%%%Zhishe_n=Zhishe。for i=1:n if(Zhishe(i,2)0) Zhishe_n(i,2)=0。 endend%%%%%%%%%%%%%%%%%%%%%%求出相對日照時(shí)數(shù)%%%%%%%%%%%%%%%%%%%%%%Jing_light=Zhishe_n(:,2)+Sanshe(:,2)。%求出晶硅日照時(shí)數(shù)for i=1:n if(Jing_light(i,1)=80amp。amp。Jing_light(i,1)=200) Jing_light(i,2)=*Jing_light(i,1)。 elseif(Jing_light(i,1)=200) Jing_light(i,2)=Jing_light(i,1)。 endendjing_day=sum(Jing_light)/1000。Fei_light=Zhishe_n(:,2)+Sanshe(:,2)。%求出非晶硅日照時(shí)數(shù)for i=1:n if(Fei_light(i,1)=30) Fei_light(i,2)=Fei_light(i,1)。 endendfei_day=sum(Fei_light)/1000。附件二：數(shù)據(jù)預(yù)處理分析的部分結(jié)果編號(hào)太陽時(shí)時(shí)角赤緯角太陽高度角太陽方位角斜面上的太陽光線入射角（度）01231074874975075175216751676167733823383338450505051505250535054691769186919866286638664866586668667附件三：轉(zhuǎn)換后經(jīng)緯度及標(biāo)準(zhǔn)日照的數(shù)據(jù)處理%導(dǎo)入數(shù)據(jù)%球坐標(biāo)與地面坐標(biāo)相互轉(zhuǎn)換%計(jì)算太陽時(shí)、時(shí)角、赤緯角%計(jì)算相應(yīng)經(jīng)度緯度下的標(biāo)準(zhǔn)日照量，取最大，得到標(biāo)準(zhǔn)日照量及對應(yīng)的經(jīng)緯度%由經(jīng)緯度得到斜角和轉(zhuǎn)換方位角clear。clc。A=importdata(39。39。)。n=length(A)。j_n=50:150。 %經(jīng)度搜索范圍w_n=50:50。 %緯度搜索范圍j_nshu=length(j_n)。 %經(jīng)度的個(gè)數(shù)w_nshu=length(w_n)。 %緯度的個(gè)數(shù)jing_b=zeros(j_nshu,w_nshu)。fei_b=zeros(j_nshu,w_nshu)。j_a=zeros(j_nshu,w_nshu)。%緯度w_a=zeros(j_nshu,w_nshu)。%經(jīng)度j_a0=。%大同緯度w_a0=。%大同經(jīng)度%%%%%%%%%%%%%%%%%%%%%地面坐標(biāo)系%%%%%%%%%%%%%%%%%%%%%% di=[sind(j_a0)*cosd(w_a0) sind(j_a0)*sind(w_a0) cosd(j_a0)。 sind(w_a0) cosd(w_a0) 0。 cosd(j_a0)*cosd(w_a0) cosd(j_a0)*sind(w_a0) sind(j_a0)]。%%%%%%%%%%%%%%%%%%j_n_i=1。for i=j_n %經(jīng)度 %%%%%%%%%%%%%%%%%求太陽時(shí)%%%%%%%%%%%%%%%%%% taiyangshi(:,1)=A(:,2)。 taiyangshi(:,2)=A(:,1)+(i105)/15。 %%%%%%%%%%%%%%%%%%%%求時(shí)角%%%%%%%%%%%%%%%%%%%%%%% shijiao=15*(taiyangshi(:,2)12)。 %%%%%%%%%%%%%%%%%%%%赤緯角%%%%%%%%%%%%%%%%%%%%%%% chiweijiao=*sin(2*pi*(284+floor(taiyangshi(:,1)/24)+1)/365)。 %%%%%%%%%%%%%%%%%%%%%% w_n_j=1。 for j=w_n %緯度 %太陽高度角的sin值 Temp=[cosd(j)*cosd(i),cosd(j)*sind(i),sind(j)]。% Qiu=[]。 Qiu=Temp*(di)39。%球坐標(biāo)中對應(yīng)的坐標(biāo) j_a(j_n_i,w_n_j)=asind(Qiu(3))。 cosw_a=Qiu(1)/cosd(j_a(j_n_i,w_n_j))。 sinw_a=Qiu(2)/cosd(j_a(j_n_i,w_n_j))。 w_a(j_n_i,w_n_j)=acosd(cosw_a)。 if sinw_a0 w_a(j_n_i,w_n_j)=w_a(j_n_i,w_n_j)。 end gaodujiao=sind(j).*sind(chiweijiao)+cosd(j).*cosd(chiweijiao).*cosd(shijiao)。%太陽高度角的sin值 %改變法向光強(qiáng)矩陣 guang=A(:,5)。 %%%%%%%%%%%%%%%%%計(jì)算斜面上的投影光強(qiáng)%%%%%%%%%%%%%%%%% guang1=gaodujiao.*guang。 for k=1:n if(guang1(k)0) guang1(k)=0。 end end %%%%%%%%%%%%%%%%%%%%%%得散射光強(qiáng)%%%%%%%%%%%%%%%%% guang22=A(:,4)。 guang2=guang22.*(1+cosd(90j_a(j_n_i,w_n_j)))/2。 %%%%%%%%%%%%%%%%%%%%%%%求總光強(qiáng)%%%%%%%%%%%%%%%%%%%%% zongguang=guang1+guang2。 %%%%%%%%%%%%%%%%%%%%%%%%晶硅%%%%%%%%%%%%%%%%%%%%% jing=zeros(n,1)。 for k=1:8760 if(zongguang(k)=80amp。amp。zongguang(k)=200)。 jing(k)=zongguang(k)*。 elseif(zongguang(k)200) jing(k)=zongguang(k)。 end end %將其轉(zhuǎn)化為標(biāo)準(zhǔn)日照時(shí)長 xiaoshishu1=jing/1000。 jing_b(j_n_i,w_n_j)=sum(xiaoshishu1)。 %%%%%%%%%%%%%%%%%%%%%%%%非晶硅%%%%%%%%%%%%%%%%%%% fei=zeros(n,1)。 for k=1:n if(zongguang(k)=30) fei(k)=zongguang(k)。 end end %將其轉(zhuǎn)化為標(biāo)準(zhǔn)日照時(shí)長 xiaoshishu2=fei/1000。 fei_b(j_n_i,w_n_j)=sum(xiaoshishu2)。 w_n_j=w_n_j+1。 %%%%%%%%%%%%%%%% end j_n_i=j_n_i+1。end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%t1=0。t2=0。for i=1:j_nshu for j=1:w_nshu if(jing_b(i,j)t1) temp_i=j_a(i,j)。%晶硅緯度 temp_j=w_a(i,j)。 %晶硅經(jīng)度 t1=jing_b(i,j)。 %晶硅標(biāo)準(zhǔn)日照 end if(fei_b(i,j)t2) temp_i2=j_a(i,j)。%非晶硅緯度 temp_j2=w_a(i,j)。 %非晶硅經(jīng)度 t2=fei_b(i,j)。 %非晶硅標(biāo)準(zhǔn)日照 end endend33