附录2：
function cx1
%======================导入数据=================== A=load('E:\\A.txt');%A种电池价格 B=load('E:\\B.txt');%B种电池价格 C=load('E:\\C.txt');%C种电池价格 NBQ=load('E:\\NBQ.txt');%逆变器价格
%======================南顶=================== AND=[0 0 0 0 0 0]';%A种电池使用数目 BND=[0 24 0 0 0 0 0]';%B种电池使用数目 CND=[0 0 0 0 0 0 0 0 0 0 0]';%C种电池使用数目
NBQND=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0]';%逆变器使用数目 %======================北墙=================== ABQ=[0 0 0 0 0 0]';%A种电池使用数目 BBQ=[0 0 0 0 0 0 0]';%B种电池使用数目
CBQ=[0 0 0 0 0 0 0 0 86 0 0]';%C种电池使用数目
NBQBQ=[0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0]';% 逆变器使用数目 %======================南墙=================== ANQ=[0 0 1 0 0 0]';%A种电池使用数目 BNQ=[0 0 0 0 0 0 0]';%B种电池使用数目
CNQ=[0 0 0 0 0 0 0 0 0 0 0]';%C种电池使用数目
NBQNQ=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]';% 逆变器使用数目 %======================西墙=================== AXQ=[1 1 1 1 0 0]';%A种电池使用数目 BXQ=[0 0 0 0 0 0 0]';%B种电池使用数目
CXQ=[0 0 0 0 0 0 0 0 0 0 0]';%C种电池使用数目
NBQXQ=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]';% 逆变器使用数目 %======================东墙=================== ADQ=[1 1 1 1 0 0]';%A种电池使用数目 BDQ=[0 0 0 0 0 0 0]';%B种电池使用数目
CDQ=[0 0 0 0 0 0 0 0 0 0 0]';%C种电池使用数目
NBQDQ=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]';% 逆变器使用数目 %======================北顶=================== ABD=[0 0 0 0 0 0]';%A种电池使用数目 BBD=[0 0 0 0 0 0 0]';%B种电池使用数目
CBD=[11 0 0 0 0 0 0 0 0 0 0]';%C种电池使用数目
NBQBD=[1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]';% 逆变器使用数目
A2=load('E:\\A2.txt');%A种电池效率 B2=load('E:\\B2.txt');%B种电池效率 C2=load('E:\\C2.txt');%C种电池效率 NBQ2=load('E:\\NBQ2.txt');%逆变器效率
AS=load('E:\\AS.txt');%A种电池面积
- 19 -
BS=load('E:\\BS.txt');%B种电池面积 CS=load('E:\\CS.txt');%C种电池面积
x=1:8760;
y1=load('E:\\y1jiu.txt');%屋顶单位辐射强度南 y2=load('E:\\y2.txt');%北墙单位辐射强度 y3=load('E:\\y3.txt');%南墙单位辐射强度 y4=load('E:\\y4.txt');%西墙单位辐射强度 y5=load('E:\\y5.txt');%东墙单位辐射强度 y6=load('E:\\y6.txt');%屋顶单位辐射强度北
%======================计算所需成本=================== disp('所需成本')
F1=AND'*A+BND'*B+CND'*C+NBQND'*NBQ F2=ABQ'*A+BBQ'*B+CBQ'*C+NBQBQ'*NBQ; F3=ANQ'*A+BNQ'*B+CNQ'*C+NBQNQ'*NBQ; F4=AXQ'*A+BXQ'*B+CXQ'*C+NBQXQ'*NBQ; F5=ADQ'*A+BDQ'*B+CDQ'*C+NBQDQ'*NBQ; F6=ABD'*A+BBD'*B+CBD'*C+NBQBD'*NBQ;
%======================计算有效辐射强度================= q=1;
while q<8761 if y1(q)<80 y1(q)=0; else end
if y2(q)<80 y2(q)=0; else end
if y3(q)<80 y3(q)=0; else end
if y4(q)<80 y4(q)=0; else end
if y5(q)<80 y5(q)=0;
- 20 -
else end
if y6(q)<80 y6(q)=0; else end q=q+1; end
%======================计算输出功率===================== q=1;
while q<8761 if y1(q)<200 y12(q)=0; else
y12(q)=y1(q); end
if y2(q)<200 y22(q)=0; else
y22(q)=y2(q); end
if y3(q)<200 y32(q)=0; else
y32(q)=y3(q); end
if y4(q)<200 y42(q)=0; else
y42(q)=y4(q); end
if y5(q)<200 y52(q)=0; else
y52(q)=y5(q); end
if y6(q)<200 y62(q)=0; else
- 21 -
y62(q)=y6(q); end q=q+1; end
M1=(trapz(x,y1)*(B2'*(BS.*BND)+C2'*(CS.*CND))+trapz(x,y12)*A2'*(AS.*AND));%屋顶南输出功率
M2=(trapz(x,y2)*(B2'*(BS.*BBQ)+C2'*(CS.*CBQ))+trapz(x,y22)*A2'*(AS.*ABQ));%北墙输出功率
M3=(trapz(x,y3)*(B2'*(BS.*BNQ)+C2'*(CS.*CNQ))+trapz(x,y32)*A2'*(AS.*ANQ));%南墙输出功率
M4=(trapz(x,y4)*(B2'*(BS.*BXQ)+C2'*(CS.*CXQ))+trapz(x,y42)*A2'*(AS.*AXQ));%西墙输出功率
M5=(trapz(x,y5)*(B2'*(BS.*BDQ)+C2'*(CS.*CDQ))+trapz(x,y52)*A2'*(AS.*ADQ));%东墙输出功率
M6=(trapz(x,y6)*(B2'*(BS.*BBD)+C2'*(CS.*CBD))+trapz(x,y62)*A2'*(AS.*ABD));%屋顶北输出功率
%======================计算35年实际发电量================== disp('35年实际发电量')
N1=M1*NBQ2'*NBQND*(10+15*0.9+10*0.8)/1000%屋顶南实际发电量 N2=M2*NBQ2'*NBQBQ*(10+15*0.9+10*0.8)/1000;%北墙实际发电量 N3=M3*NBQ2'*NBQNQ*(10+15*0.9+10*0.8)/1000;%南墙实际发电量 N4=M4*NBQ2'*NBQXQ*(10+15*0.9+10*0.8)/1000;%西墙实际发电量 N5=M5*NBQ2'*NBQDQ*(10+15*0.9+10*0.8)/1000;%东墙实际发电量 N6=M6*NBQ2'*NBQBD*(10+15*0.9+10*0.8)/1000;%屋顶北实际发电量 %======================计算经济效益==================== disp('经济效益') P1=N1*0.5-F1 P2=N2*0.5-F2; P3=N3*0.5-F3; P4=N4*0.5-F4; P5=N5*0.5-F5; P6=N6*0.5-F6;
%======================计算回收年限==================== disp('回收年限')
Q0=N1*0.5/(10+15*0.9+10*0.8); Q1=Q0; k=1;
while F1-Q1>0 if k<=10
Q1=Q1+Q0*0.1; else if k<=25 Q1=Q1+Q0*0.9*0.1; else if k<=35
- 22 -
Q1=Q1+Q0*0.8*0.1; else
disp('error') end end end k=k+0.1; end k end
附录3：
function cx2
%======================导入数据=================== SS=load('E:\\SS.txt');%水平面散射辐射强度 NZ=load('E:\\NZ.txt');%南向直射辐射强度 FZ=load('E:\\FZ.txt');%法向直射辐射强度 SZ=load('E:\\SZ.txt');%水平面总辐射强度
%======================选择有效值=================== q=1;
while q<8761 if SS(q)<200 SS(q)=0; else end
if NZ(q)<200 NZ(q)=0; else end
if FZ(q)<200 FZ(q)=0; else end
if SZ(q)<200 SZ(q)=0; else end q=q+1; end
- 23 -