太阳能小屋的设计

太阳能小屋的设计
摘要
本文探讨了太阳能小屋外表面电池组件优化铺设的问题，并在此基础上给出了太阳能小屋的改进方案。

针对问题一，本文首先计算了屋顶的太阳能辐射总量（见表1），在考虑小屋全年光伏电池发电总量最大、单位发电量费用最小的基础上建立了贴附式电池组件安装收益规划模型，对小屋屋顶和南面进行了部分铺设，得到了小屋表面电池板铺设方案（见表3，图2，图4）和逆变器的选择方案（见表6，图3，图5），并计算出了在当前方案下小屋的发电总量为1052097.37 kw·h，经济总收益为315451元，回收成本年限为７年。

针对问题二，本文首先分析了太阳辐射量与倾斜面的方位角和倾角的关系，并建立了相应的模型，通过计算得到了小屋顶面和南墙电池组件的最佳倾斜角为37.3°（相对于水平面），方位角为０°。

在问题一模型的基础上建立了改进的架空式电池组件安装收益规划模型，进而得到了小屋表面电池组件的铺设方案（见表8，图8）和逆变器的选用方案（见表9，图9），此时小屋的发电总量为1168997，经济总收益372263元，回收年限为7年。

针对问题三，本文根据小屋的建筑要求及问题一、二中的结论，制定了小屋设计原则，重新构建了小屋（见图10至图14），并给出了改进后小屋的电池板铺设方案和逆变器的选择方案（图15至图18）。

最后，本文对所建模型的优点和缺点进行了客观的评价，认为本文研究的结果在实际应用中有很好的参考价值。

关键词：多目标规划，最佳倾斜角，太阳辐射总量，收益规划模型
1.问题重述
在设计太阳能小屋时，需在建筑物外表面（屋顶及外墙）铺设光伏电池，光伏电池组件所产生的直流电需要经过逆变器转换成220V交流电才能供家庭使用，并将剩余电量输入电网。

不同种类的光伏电池每峰瓦的价格差别很大，且每峰瓦的实际发电效率或发电量还受诸多因素的影响，如太阳辐射强度、光线入射角、环境、建筑物所处的地理纬度、地区的气候与气象条件、安装部位及方式（贴附或架空）等。

因此，在太阳能小屋的设计中，研究光伏电池在小屋外表面的优化铺设是很重要的问题。

附件1-7提供了相关信息。

请参考附件提供的数据，对下列三个问题，分别给出小屋外表面光伏电池的铺设方案，使小屋的全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小，并计算出小屋光伏电池35年寿命期内的发电总量、经济效益（当前民用电价按0.5元/kWh计算）及投资的回收年限。

在求解每个问题时，都要求配有图示，给出小屋各外表面电池组件铺设分组阵列图形及组件连接方式（串、并联）示意图，也要给出电池组件分组阵列容量及选配逆变器规格列表。

在同一表面采用两种或两种以上类型的光伏电池组件时，同一型号的电池板可串联，而不同型号的电池板不可串联。

在不同表面上，即使是相同型号的电池也不能进行串、并联连接。

应注意分组连接方式及逆变器的选配。

问题1：请根据山西省大同市的气象数据，仅考虑贴附安装方式，选定光伏电池组件，对小屋（见附件2）的部分外表面进行铺设，并根据电池组件分组数量和容量，选配相应的逆变器的容量和数量。

问题2：电池板的朝向与倾角均会影响到光伏电池的工作效率，请选择架空方式安装光伏电池，重新考虑问题1。

问题3：根据附件7给出的小屋建筑要求，请为大同市重新设计一个小屋，要求画出小屋的外形图，并对所设计小屋的外表面优化铺设光伏电池，给出铺设及分组连接方式，选配逆变器，计算相应结果。

2.模型假设
1、假设单晶硅和多晶硅电池启动发电的表面总辐射量≥80W/m
2、薄膜电池表面总辐射量≥30W/m2；
2、假设只允许同型号的光伏组件进行串联，不同平面上的任何电池板不能相联；
3、假设并联的光伏组件端电压相差不应超过10%；
4、假设本题所研究时间段内没有阴雨天气；
5、假设小屋光伏电池在35年寿命期内不会有任何损坏。

3.问题分析
3.1问题一分析
本题只考虑贴附安装方式，由题目可知在选择安装太阳能电池板时，要考虑如何安装才能使得小屋的全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小，因此在解决本题时我们需要列出只安装电池板时，发电量最大而费用最少的目标函数，并将发电量最大转化为收益最大，从而将多目标函数转化为单目标函数，建立优化模型解决。

通过运用lingo11.0软件计算出只安装电池板
的收益后，可将安装后亏损的面剔除，对剩余的面加上逆变器进行分析最终解决安装方案。

3.2问题二分析
众所周知，地球的“自转”形成昼夜的变化，地球的“公转”形成四季的变化。

地球的“自转”与“公转”的运动使得房屋的不同面在不同的时间分别受到相应的辐射量，而本题就是希望最大效率的吸收太阳辐射量。

由于在光伏发电系统的研究中，往往利用每小时平均光强和每小时平均温度等气象数据计算光伏阵列每小时的发电量。

由于电池板的朝向与倾角均会影响到光伏电池的工作效率，因此本题在问题一的基础上采用架空的方式安装光伏电池板。

光伏系统分为固定式光伏系统和安装了跟踪装置的太阳能光伏发电系统，由于跟踪装置比较复杂，初始成本和维护成本较高，目前太阳能光伏发电系统中使用跟踪装置的相对较少。

因此本文主要研究固定光伏电池系统。

由于本文需要全年负载符合均匀的光伏电池系统，因此确定光伏电池板的倾斜角是获得安装方案的重要环节。

本文预利用题目所给山西大同全年的辐射量来计算电池板的倾斜角。

运用所求倾斜角求得本斜面上的太阳辐射量。

最终根据太阳辐射量求得所需要的电池板类型及数量，进而确定逆变器型号以及安装方案。

3.3问题三分析
问题三要求根据相关要求为大同市重新设计一个小屋，并对小屋进行光伏电池组铺设。

考虑小屋外形的设计，由问题二分析可知电池板的朝向与倾角均会影响到光伏电池组的问题三要求根据相关要求为大同市重新设计一个小屋，并对小屋进行光伏电池组铺设。

工作效率，那么在重新设计小屋时，应尽可能的将小屋的外形设计成有利于提高光伏电池组工作效率的样式。

考虑小屋大小的设计，首先，在满足所给小屋设计要求的前提下，应该让小屋的尺寸尽可能大，这样能更有利于安装更多的光伏电池板；其次，小屋墙面尺寸的设计应有利于光伏电池板的铺设，这里我们应从墙面铺设紧密性的程度出发设计小屋墙面。

最后，在将小屋重新设计的基础上，对小屋进行优化铺设。

4.符号说明与名词解释
5.问题一模型的建立与求解
5.1 模型准备——屋顶接收辐射总量的计算
图1 太阳角关系示意图
由于屋顶有倾斜角，因此将屋顶作为斜面计算总的辐射量。

斜面的总辐射量计算如下：
朝向赤道斜面上的太阳辐射量，通常采用Klein 的计算方法，倾斜面上所接受到的太阳辐射总量T H 由直接辐射量bt H 、天空散射辐射量dt H 以及地面反射辐射量t H τ组成，即：
T bt dt t H H H H τ=++ （1）
bt H 与水平面上的直接辐射量'()bt b H H S =之间又如下关系：
b
b bt R H H ⋅= （2）
对于朝向赤道的倾斜面，b R 可以由下式确定：
()()cos cos sin sin sin 180
cos cos sin sin sin 180
ST ST b S S R π
φβδωωφβδπ
φδωωφδ
-⋅+
-⋅=
⋅⋅+
⋅ （3）
其中：ϕ是当地维度，β是倾斜角，δ是太阳赤纬角，b R 为倾斜面与水平面上直接辐射量之比。

水平面上的日落时角为：
[]1cos tan tan S ωφδ-=-⋅ (4)
斜面上的日落时角为： ()1
min cos tan tan s
ST
ωωφβδ-⎧⎪=⎨--⎡⎤⎪⎣⎦⎩ (5)
通过对附件1数据的分析，并借助matlab 软件，首先统计数据并运用方差
法计算得每月的典型日期，然后对12个月的典型日期，运用方差法同样可以得到当年的典型日期，以此作为n .
太阳赤纬角为：
2(284)23.45sin 365n πδ+⎛⎫
= ⎪⎝⎭
(6)
其中：n 为日期序号，n 由方差法确定。

Hay 模型认为倾斜面上天空散射辐射量是由太阳光盘的辐射量和其余天空穹顶均匀分布的散射辐射量两部分组成，可表达为： ()0
01
1cos 12d d dt d b H H H H H H R H H β⎡⎤⎛⎫--=++-⎢⎥ ⎪⎝⎭⎣⎦ (7)
其中：b H 和d H 分别为水平面上直接和散射辐射量；b R 为倾斜面与水平面上
直接辐射量之比，0H 为大气外层水平面上太阳辐射量；β为倾角。

大气外层的水平辐射量可由下式求出：
024
36010.033cos cos sin sin sin 365180SC S S n H I πφδωωφδπ⎡⎤⎛⎫
=
+⋅⋅+⋅⋅ ⎪⎢⎥⎣⎦⎝⎭
(8) 其中：sc I 为太阳常数。

由于本题是将光伏电池以贴附方式安装，电池板与屋顶的缝隙较小，因此地
面反射辐射量可以忽略不计
综上所述，屋顶上太阳辐射总量的表达式为：
()0
01
1cos 12d d T b b d b H H H H H H R H R H H β⎡⎤
⎛⎫--=+++-⎢
⎥ ⎪⎝⎭⎣⎦
(9) 5.2模型建立——电池组件安装收益规划模型
本题仅考虑贴附安装方式对小屋的部分表面进行铺设，由于要保证发电量最大、成本最低，因此建立一个多目标优化模型。

其目标为全年太阳能最大发电总量和建设成本最低。

5.2.1目标函数
由于所有光伏组件在0～10年效率按100%，10～25年按照90%折算，25年后按80%折算。

因此各组件的35年发电时间可折合成等效发电时间为：10+15*90%+10*80%=24.3年。

小屋屋顶及墙面上电池组件在等效时间内的总发电量可按下式计算： 24
1124.3()p i j i i f n H S η==⋅⋅⋅∑ （10）
其中，i n 为使用光伏电池的数目，j H 为屋顶及房屋各面所受的辐射强度，i S 为所需电池板的面积，i η为光伏电池的转换效率。

由于小屋光伏电池板在其寿命时间内不存在任何的损坏维修，因此电池板的花费成本为：
24
1()c i i i i f n P W ==⋅⋅∑ （11）
其中，i P 为电池每峰瓦的价格，i W 为电池组件功率。

由公式（1）、（2）可列出最大收益模型为：
24
24
3
11
1
max 24.3()100.5()j i i i i i i i f H n S n P W η-===⋅⋅⋅⋅⋅-⋅⋅∑∑ （12）
5.2.2约束条件
（1）对电池板铺设总面积的限制：
24
1
24
1
i i j i i i j i n S S n S S ==⎧⋅≤⎪⎪⎨⎪⋅≤⎪⎩∑∑
（2）电池板尺寸的限制：
24
1
24
1
i i j i i i j i n L L n D D ==⎧⋅≤⎪⎪⎨⎪⋅≤⎪⎩∑∑ （3）逆变器成本限制：
5
1
1
2
k p i SN f =<
∑ 5.2.3贴附式电池组件安装收益规划模型
综合上述目标函数及约束条件，可列出贴附式电池组件安装收益规划模型如下：
24
24
3
11
1
max 24.3()100.5()j i i i i i i i f H n S n P W η-===⋅⋅⋅⋅⋅⋅-⋅⋅∑∑
24
124
1
24
15
1
12i i j i i i j i i i j
i k p i n S S n L L n D D
SN f ====⎧⋅<⎪⎪⎪⋅≤⎪⎪⎨⎪⋅≤⎪⎪⎪<⎪⎩∑∑∑∑ (13)
5.3模型求解
5.3.1太阳辐射总量的计算
根据模型准备数据以及题目所给资料，运用matlab （R2012a 7.14.0.739）软件求得小屋房顶及各面的太阳辐射总量如下：
根据屋顶及房屋各面的辐射总量以及目标函数：
24
24
3
11
1
max 24.3()100.5()j i i i i i i i f H num S num P M η-===⋅⋅⋅⋅⋅⋅-⋅⋅∑∑
运用matlab （R2012a 7.14.0.739）软件可求得小屋屋顶及房屋各面安装电池板后收益如下：
虑在北面安装电池板。

其余各面在不考虑逆变器的情况下所需的电池板数如下：
表
根据表3中小屋屋顶及各面墙选择电池板的数目从而确定各面墙输入逆变器的功率，通过对附件5的数据处理得到各逆变器的标准输入功率如下：
根据表4可得到各面墙安装电池后产生功率如下：
通过表3及表4的数据统计可得小屋屋顶及各面墙应选择的逆变器如下：
图2 南面墙电池安装示意图
注：绿色矩形为C1单晶硅光伏电池组件，蓝色矩形为C9非晶硅薄膜光伏电池组件。

图3 南立面电池组件连接及功率示意图顶面电池组件安装图如下：
图4 顶面电池组件安装示意图注：黄色矩形为A3单晶硅光伏电池组件。

图5 顶面电池组件连接分组阵列及功率示意图
在以上方案下，计算得35年的南墙总发电量32818.37kw ·h ，屋顶总发电量为 1019279 kw ·h ， 1052097.37 kw ·h ，
6.问题二模型的建立与求解
6.1最佳倾角的求解[1-4]
朝向赤道斜面上的太阳辐射量，通常采用Klein 的计算方法，倾斜面上所接受到的太阳辐射总量T H 由直接辐射量bt H 、天空散射辐射量dt H 以及地面反射辐射量t H τ组成，即：
t dt bt T H H H H τ++=
bt H 与水平面上的直接辐射量)('S H H b bt =之间又如下关系：
b b bt R H H ⋅=
对于朝向赤道的倾斜面，b R 可以由下式确定：
()()cos cos sin sin sin 180
cos cos sin sin sin 180
ST ST b S S R π
φβδωωφβδπ
φδωωφδ
-⋅+
-⋅=
⋅⋅+
⋅
其中：ϕ是当地维度，β是倾斜角，δ是太阳赤纬角。

水平面上的日落时角为：
[]1cos tan tan S ωφδ-=-⋅
斜面上的日落时角为：
()1
min cos tan tan s
ST
ωωφβδ-⎧⎪=⎨--⎡
⎤⎪⎣⎦⎩ 太阳赤纬角为：
2(284)23.45sin 365n πδ+⎛⎫
=
⎪⎝⎭
其中：n 为日期序号，n 由方差法确定。

Hay 模型认为倾斜面上天空散射辐射量是由太阳光盘的辐射量和其余天空穹顶均匀分布的散射辐射量两部分组成，可表达为：
()0
01
1cos 12d d dt d b H H H H H H R H H β⎡⎤
⎛⎫--=++-
⎢⎥ ⎪⎝⎭⎣⎦
其中：b H 和d H 分别为水平面上直接和散射辐射量；b R 为倾斜面与水平面上
直接辐射量之比，0H 为大气外层水平面上太阳辐射量；β为倾角。

大气外层的水平辐射量可由下式求出：
024
36010.033cos cos sin sin sin 365180SC S S n H I πφδωωφδπ⎡⎤⎛⎫
=
+⋅⋅+⋅⋅ ⎪⎢⎥⎣⎦⎝⎭
(14) 其中，sc I 为太阳常数。

地面反射辐射常量的表达式为：
()1
1cos 2
H ρβ⋅⋅+ (15)
其中：ρ为地面反射率，一般情况下2.0=ρ，倾斜面上太阳辐射总量的表
达式为：
()()0
01
1
1cos 11cos 22
d d T b b d b H H H H H H R H R H H H βρβ⎡⎤⎛⎫--=+++-
+⋅⋅-⎢⎥ ⎪⎝⎭⎣⎦(16) 由上式直接求出最大太阳辐射倾斜角的数学表达式有些困难，在冬半年，太
阳赤纬角为负值。

如当地纬度为ϕ，通常总有ϕβϕω<-)(，因此倾斜面上日落时角s ω和水平面日落时角st ω相等，这时可直接推导出最佳倾角的数学表达式。

此时，水平面上的总辐射量T H 为：
0001
()()(1cos )(1cos )22
d b d T b b b H H H H H R H H H H βρβ=+
+-⨯+- (17)
对于确定的地点，其太阳辐射量及地面反射率等均为常数，将T H 对ρ求导，
并令
0T
dH d β
=，得： ()2221tan tan 180arctan tan tan 11180b b b S S b b b S S H H H H H H H H H H H H πφωωβπφωωρ⎡⎤⎛⎫+-+ ⎪⎢⎥⎝⎭⎢
⎥=⨯⎛⎫⎛⎫⎢⎥
-+-+- ⎪ ⎪⎢⎥⎝⎭⎝⎭⎣⎦
(18) 当直接辐射量为0时，最佳倾角不等于0，这正是天空散射辐射量各向异性
模型所需的结果。

对于夏半年，最佳倾角往往接近于0，甚至为负值。

6.2电池板排列间距的计算
影子倍率法计算太阳高度角并计算间距
[5]
图6 电池阵列安装示意图
一般在水平面垂直竖立的高为L 的木杆，其南北方向的角度为h ，方位角为α，那么影子倍率R 可由下式表示：
cot cos s L
R h L
α==⋅ (19)
其中： R ——影子倍率 L ——阵列高度 s L ——影子长度
h ——太阳高度角
α——太阳方位角
由图6以及公式cot cos s
L R h L
α=
=⋅可以得出： cot cos s L L h α=⋅⋅
由于本文研究的房屋是正南方向，因此太阳方位角为0.
cot s L L h =⋅
根据球面三角函数分析认为太阳高度与观测者的地理纬度、太阳赤纬和方位角有着一定的关系，他们之间的关系式为：
sinh sin sin cos cos cos ,ϕδϕδα=⋅+⋅⋅ (20) 其中h 为太阳高度角，ω为时角，δ为当时的太阳赤纬，φ为当地的纬度(大同的纬度为o 1.40)，根据最佳倾斜角的求解公式可知：
斜面上的日落时角为：
()1min cos tan tan s
ωωφβδ-⎧⎪=⎨--⎡
⎤⎪⎣⎦⎩ (21)
太阳赤纬角为：
2(284)23.45sin 365n πδ+⎛⎫
= ⎪⎝⎭
(22)
由于本文研究的太阳方位角为0，则1=αcos ，太阳高度角公式变为：
sinh sin sin cos cos cos()sin[90()]φδφδφδφδ=⋅+⋅=-=-- (23)
得出计算太阳高度角的基本公式：
90()h φδ=-- (24)
通过以上模型及附件给出的数据可以计算出电池板间的阵列间距，但是由于电池板的数量级规格较多，使用此模型计算大小不同电池板的阵列间距极为困难，且经过复杂计算后，准确度也有很大偏差。

因此，本文使用阴影转换的方法，可以方便的计算出阵列间距及所需电池板数目。

方法示意图如下：
图7 阴影面及间距求法
在本题中已经得到最佳倾斜角β，通过倾斜角及电池板面积可以求得电池板
在屋顶平面上产生的阴影，进而可以将电池板的面积等效为电池板在屋顶平面上产生的阴影的面积。

此阴影面积可以等效为问题一中贴附电池板的面积，从而将求阵列间距问题转化为问题一中求贴附电池板的面积。

运用此方法可建立模型求得在屋顶及房屋各面上所需的电池板数目。

6.3架空式电池组件安装收益规划模型
本题是在问题一的基础上求得最佳倾斜角建立优化模型。

本题同样是多目标优化模型，其目标为全年太阳能最大发电总量和建设成本最低。

根据第一问的贴附式电池组件安装收益规划模型可以改进的到架空式电池组件安装收益规划模型如下：
24
24
3
11
1
max 24.3(/cos(arccos(0.98287)))100.5()
j i i i i i i i f H n S n P W βη-===⋅⋅⋅-⋅⋅-⋅⋅∑∑ 24
124
1
24
15
1
12i i j i i i j i i i j
i k p i n S S n L L n D D
SN f ====⎧⋅≤⎪⎪⎪⋅≤⎪⎪⎨⎪⋅≤⎪⎪⎪<⎪⎩∑∑∑∑ (25) 6.4模型求解
根据前面公式(25)运用matlab （R2012a 7.14.0.739）求得 337.=β 根据建立模型运用matlab （R2012a 7.14.0.739）软件算出架空后各面安装电池板收益如下：
屋顶及除去背面的其他面吸收的辐射量，可确定在架空安装方式下选择的电池型号及数目如下：
房顶及各面安装电池板后产生功率以及相应逆变器选择情况如下：
屋顶安装电池板的架构图如下：
图8 问题二屋顶架空排列图
屋顶安装逆变器后组件连接图如下：
图中：省略号为16个电池组
图9 屋顶安装逆变器后组件连接图
在以上铺设方式下计算可得发电总量为 1168997，经济总收益372263元，回收年限为7年。

7．问题三模型的分析与建立
7.1 小屋改进原则分析
1 整体最大设计原则
根据第三问的要求，重新设计的小屋在满足一定建筑要求的情况下应尽量增加受光面积，这里我们考虑在满足条件的前提下尽可能的增大房屋体积。

2 增进电池板功效设计原则
根据问题二的相关结论，可以将小屋的外形设计为有利于增加电池板功效形式，这里我们考虑将屋顶的倾角尽量的接近于37.3°。

3 墙面尺寸合适设计原则
考虑墙面大小和门窗大小设计时，应考虑到实际铺设光伏电池板时墙面的利用问题，这里我们设计的墙面大小应使得铺设电池板的墙面利用率最高。

7.2 小屋改进模型
根据分析设定的改进模型，我们重新设计的小屋模型如下：
图10 小屋整体透视图
图11小屋东视图图12小屋西视图
图13小屋南面视图
图14 小屋北面视图
7.3 改进小屋的铺设
屋顶：
图15 改进后小屋屋顶电池板安装图
图中：省略号为8个电池板
图16 改进后小屋屋顶安装逆变器连接图
图示颜色为B2类型的电池板
图17 改进后小屋南面墙电池板安装示意图
图18改进后小屋南面墙逆变器电池连接图
8.模型评价
8.1模型优点
对于问题一，本文将多目标优化模型合理转化为单目标优化模型，便于程序实现，小屋电池板的铺设方案通过此模型得到很好的优化。

在现行方案下发电总量1052097.37 kw·h，经济总收益315451元，回收成本年限为７年。

对于问题二，在通过对问题一改进后的模型，很好的完成了小屋架空式电池板的铺设，通过模型的对比和实际的统计的结果，此模型达到了很好的优化效果。

发电总量为 1168997，经济总收益372263元，回收年限为7年。

可见电池板架空放回比贴附方式铺设能效更高，这符合实际情况。

对于问题三，通过对一、二题的分析，对房屋接受阳光面积极限大的思想，建立了斜屋顶的房屋的方案。

通过计算，此模型明显提高了房屋的发电量和经济收益。

8.2模型缺点
对于问题三，由于给定电池板的大小和形状限制，导致房屋不能设计的接触阳光的面积尽可能大，导致发电量不能最大化。

参考文献
[1] 杨金焕，毛家俊，陈中华，不同方位倾斜面上太阳辐射量及最佳倾角的计算，上海交通大学学报，2002,37（7）：1032-1036.
[2] 唐润生，吕恩荣，集热器最佳倾角的选择，太阳能学报，1988，9 (4)：269-375.
[3] 杨金焕，于化从，葛亮，太阳能光伏发电应用技术[M]，北京：电子工业出版社，2009.
[4] 杨金焕，固定式光伏阵列最佳倾角的分析，太阳能学报， 1992（13）:86.
[5] 吴玉庭，任建勋，梁新刚，聚光条件下太阳电池热电特性分析，太阳能学报，2004（25）:337.
附录：
Matlab程序：
get_fushe.m
%根据限制求各月或年的辐射量
clc;clear all;close all;
sum1 = zeros(1,5); %辐射强度总数统计
num1 = zeros(1,5); %辐射强度相加个数
aver1 = zeros(1,5);
fushe1 = 30; %辐射强度限制
i=4; %附件月或年起始位置
%设置月或年的读取区域
range1 = strcat('H',num2str(i),':H',num2str(8763));%东
range2 = strcat('I',num2str(i),':I',num2str(8763));%南
range3 = strcat('J',num2str(i),':J',num2str(8763));%西
range4 = strcat('K',num2str(i),':K',num2str(8763));%北
range51 = strcat('E',num2str(i),':E',num2str(8763));%水平
range52 = strcat('F',num2str(i),':F',num2str(8763));%散射
range53 = strcat('G',num2str(i),':G',num2str(8763));%直射
%从附件读取月或年的数据
b1 = xlsread('cumcm2012B附件4_山西大同典型气象年逐时参数及各方向辐射强度.xls','逐时气象参数',range1);%东
b2 = xlsread('cumcm2012B附件4_山西大同典型气象年逐时参数及各方向辐射强度.xls','逐时气象参数',range2);%南
b3 = xlsread('cumcm2012B附件4_山西大同典型气象年逐时参数及各方向辐射强度.xls','逐时气象参数',range3);%西
b4 = xlsread('cumcm2012B附件4_山西大同典型气象年逐时参数及各方向辐射强度.xls','逐时气象参数',range4);%北
b51 = xlsread('cumcm2012B附件4_山西大同典型气象年逐时参数及各方向辐射强度.xls','逐时气象参数',range51);
b52 = xlsread('cumcm2012B附件4_山西大同典型气象年逐时参数及各方向辐射强度.xls','逐时气象参数',range52);
b53 = xlsread('cumcm2012B附件4_山西大同典型气象年逐时参数及各方向辐射强度.xls','逐时气象参数',range53);
b5 = b51+b52+b53/0.98295;%顶
i=1; %重置计数
while i<=8760 %附件月或年结束位置
if b1(i)>=fushe1
sum1(1)=sum1(1)+b1(i);
num1(1)=num1(1)+1;
end
if b2(i)>=fushe1
sum1(2)=sum1(2)+b2(i);
num1(2)=num1(2)+1;
end
if b3(i)>=fushe1
sum1(3)=sum1(3)+b3(i);
num1(3)=num1(3)+1;
end
if b4(i)>=fushe1
sum1(4)=sum1(4)+b4(i);
num1(4)=num1(4)+1;
end
if b5(i)>=fushe1
sum1(5)=sum1(5)+b5(i);
num1(5)=num1(5)+1;
end
i=i+1;
end
sum1
num1
Lingo程序：
question_1.lg4
model:
sets:
!FUSHE*：各个面的辐射强度，POWER_OUT：年输出电量;
xinxi/1..24/:L,H,P,N,num,M,FUSHE1,FUSHE2,FUSHE3,FUSHE4,FUSHE5,POW ER_OUT;
about/1..5/:AREA,LONG,WIDTH;
endsets
data:
L=1580 1956 1580 1651 1650 1956 1650 1956 1482 1640 1956 1956 1668 1300 1321 1414 1400 1400 310 615 615 920 818 1645;
H=808 991 808 992 991 991 991 991 992 992 992 992 1000 1100 711 1114 1100 1100 355 180 355 355 355 712;
N=16.84 16.64 18.70 16.5 14.98 15.11 16.21 16.39 15.98 14.8 15.98 15.2 14.99 6.99 6.17 6.35 5.84 6.49 3.63 3.63 3.66 3.66 4.13 4.27;
P=14.9 14.9 14.9 14.9 14.9 14.9 12.5 12.5 12.5 12.5 12.5 12.5 12.5 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 ;
M=215 325 200 270 245 295 265 320 210 240 280 295 250 100 58 100 90
100 4 4 8 12 12 50;
POWER_OUT=266.919 418.9063 253.55 339.09 323.7234 396.7488 341.5691 408.7622 279.888 312.174 373.184 386.507 331.0125 168.36 95.942 163.35 145.404 164 9.345 8.82 18.69 28.035 24.03 96.25;
LONG=7100 7100 10100 10100 10100;
WIDTH=3800 3800 3200 3200 6511.53;
AREA=24230000 26980000 19235309 44440000 60870453;
FUSHE1=524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26;
FUSHE2=795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59;
FUSHE3=1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06;
FUSHE4=132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05;
FUSHE5=3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3767900 3767900 3767900 3767900 3767900 3767900 3767900 3767900 3767900 3767900 3767900;
enddata
!i代表各个房屋的各个面，其中 1；东，2：西，3：南，4：北，5：订满。

FUSHE*和i的值相对应;
i=5;
!目标函数;
max=24.3*@sum(xinxi:(FUSHE5/1000000)*num*L*H*(N/100))*0.5/1000-@s um(xinxi:num*P*M);
!YEAR*@sum(xinxi:(FUSHE5/1000000)*num*L*H*(N/100))*0.5/1000=@sum( xinxi:num*P*M);!计算回本年限，YEAR为回本年数;
out_power=24.3*@sum(xinxi:(FUSHE5/1000000)*num*L*H*N/100)/1000;
Acost=@sum(xinxi:num*P*M);
@for(xinxi:@gin(num));
@sum(xinxi:num*L*H)>=AREA(i)-1000;
@sum(xinxi:num*L*H)<=AREA(i);
!in_power:变压器输入功率;
in_power=@sum(xinxi:num*POWER_OUT);
end
Lingo程序：
question_2.lg4
model:
sets:
!FUSHE*：各个面的辐射强度，POWER_OUT：年输出电量,SR:电池板的阴影面积;
xinxi/1..24/:L,H,P,N,num,M,FUSHE1,FUSHE2,FUSHE3,FUSHE4,FUSHE5,POW ER_OUT;
about/1..5/:AREA,LONG,WIDTH;
endsets
data:
L=1580 1956 1580 1651 1650 1956 1650 1956 1482 1640 1956 1956 1668 1300 1321 1414 1400 1400 310 615 615 920 818 1645;
H=808 991 808 992 991 991 991 991 992 992 992 992 1000 1100 711 1114 1100 1100 355 180 355 355 355 712;
N=16.84 16.64 18.70 16.5 14.98 15.11 16.21 16.39 15.98 14.8 15.98 15.2 14.99 6.99 6.17 6.35 5.84 6.49 3.63 3.63 3.66 3.66 4.13 4.27;
P=14.9 14.9 14.9 14.9 14.9 14.9 12.5 12.5 12.5 12.5 12.5 12.5 12.5 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 4.8 ;
M=215 325 200 270 245 295 265 320 210 240 280 295 250 100 58 100 90 100 4 4 8 12 12 50;
POWER_OUT=266.919 418.9063 253.55 339.09 323.7234 396.7488 341.5691 408.7622 279.888 312.174 373.184 386.507 331.0125 168.36 95.942 163.35 145.404 164 9.345 8.82 18.69 28.035 24.03 96.25;
LONG=7100 7100 10100 10100 10100;
WIDTH=3800 3800 3200 3200 6511.53;
AREA=24230000 26980000 19235309 44440000 60870453;
FUSHE1=524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 524961.83 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26 581235.26;
FUSHE2=795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47 795135.47
872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59 872801.59;
FUSHE3=1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1008343.16 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06 1044044.06;
FUSHE4=132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 132000.71 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05 243181.05;
FUSHE5=3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3757100 3767900 3767900 3767900 3767900 3767900 3767900 3767900 3767900 3767900 3767900 3767900;
enddata
!i代表各个房屋的各个面，其中 1；东，2：西，3：南，4：北，5：顶面。

FUSHE*和i的值相对应;
i=5;
!目标函数;
max=24.3*@sum(xinxi:(FUSHE5/1000000)*num*L*H/0.9*(N/100))*0.5/100 0-@sum(xinxi:num*P*M);
!YEAR*@sum(xinxi:(FUSHE5/1000000)*num*L*H*(N/100))*0.5/1000=@sum( xinxi:num*P*M);!计算回本年限，YEAR为回本年数;
out_power=24.3*@sum(xinxi:(FUSHE5/1000000)*num*L*H/0.9*N/100)/100 0;
Acost=@sum(xinxi:num*P*M);
@for(xinxi:@gin(num));
@sum(xinxi:num*L*H)>=AREA(i)-1000;
@sum(xinxi:num*L*H)<=AREA(i);
!in_power:变压器输入功率;
in_power=@sum(xinxi:num*POWER_OUT);
end。