太阳能小屋外表面电池组件优化铺设的设计-毕业论文.doc

太阳能小屋外表面电池组件优化铺设的设计太阳能小屋的设计摘要本文探讨了太阳能小屋外表面电池组件优化铺设的问题，并在此基础上给出了太阳能小屋的改进方案。针对问题一，本文首先计算了屋顶的太阳能辐射总量（见表1），在考虑小屋全年光伏电池发电总量最大、单位发电量费用最小的基础上建立了贴附式电池组件安装收益规划模型，对小屋屋顶和南面进行了部分铺设，得到了小屋表面电池板铺设方案（见表3，图2，图4）和逆变器的选择方案（见表6，图3，图5），并计算出了在当前方案下小屋的发电总量为1052097.37 kwh，经济总收益为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/m2、薄膜电池表面总辐射量30W/m2 ；2、假设只允许同型号的光伏组件进行串联，不同平面上的任何电池板不能相联；3、假设并联的光伏组件端电压相差不应超过10%；4、假设本题所研究时间段内没有阴雨天气；5、假设小屋光伏电池在35年寿命期内不会有任何损坏。3.问题分析3.1问题一分析本题只考虑贴附安装方式，由题目可知在选择安装太阳能电池板时，要考虑如何安装才能使得小屋的全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小，因此在解决本题时我们需要列出只安装电池板时，发电量最大而费用最少的目标函数，并将发电量最大转化为收益最大，从而将多目标函数转化为单目标函数，建立优化模型解决。通过运用lingo11.0软件计算出只安装电池板的收益后，可将安装后亏损的面剔除，对剩余的面加上逆变器进行分析最终解决安装方案。3.2问题二分析众所周知，地球的“自转”形成昼夜的变化，地球的“公转”形成四季的变化。地球的“自转”与“公转”的运动使得房屋的不同面在不同的时间分别受到相应的辐射量，而本题就是希望最大效率的吸收太阳辐射量。由于在光伏发电系统的研究中，往往利用每小时平均光强和每小时平均温度等气象数据计算光伏阵列每小时的发电量。由于电池板的朝向与倾角均会影响到光伏电池的工作效率，因此本题在问题一的基础上采用架空的方式安装光伏电池板。光伏系统分为固定式光伏系统和安装了跟踪装置的太阳能光伏发电系统，由于跟踪装置比较复杂，初始成本和维护成本较高，目前太阳能光伏发电系统中使用跟踪装置的相对较少。因此本文主要研究固定光伏电池系统。由于本文需要全年负载符合均匀的光伏电池系统，因此确定光伏电池板的倾斜角是获得安装方案的重要环节。本文预利用题目所给山西大同全年的辐射量来计算电池板的倾斜角。运用所求倾斜角求得本斜面上的太阳辐射量。最终根据太阳辐射量求得所需要的电池板类型及数量，进而确定逆变器型号以及安装方案。3.3问题三分析问题三要求根据相关要求为大同市重新设计一个小屋，并对小屋进行光伏电池组铺设。考虑小屋外形的设计，由问题二分析可知电池板的朝向与倾角均会影响到光伏电池组的问题三要求根据相关要求为大同市重新设计一个小屋，并对小屋进行光伏电池组铺设。工作效率，那么在重新设计小屋时，应尽可能的将小屋的外形设计成有利于提高光伏电池组工作效率的样式。考虑小屋大小的设计，首先，在满足所给小屋设计要求的前提下，应该让小屋的尺寸尽可能大，这样能更有利于安装更多的光伏电池板；其次，小屋墙面尺寸的设计应有利于光伏电池板的铺设，这里我们应从墙面铺设紧密性的程度出发设计小屋墙面。最后，在将小屋重新设计的基础上，对小屋进行优化铺设。4.符号说明与名词解释符号说 明光照辐射强度，表示房屋的受辐射面表示所用光伏电池的个数，表示光伏电池的型号表示每种型号光伏电池的转换效率表示每种型号光伏电池的价格表示每种光伏电池的组件功率，分别表示每种型号光伏电池的长度，宽度，分别表示各个光伏电池的面积，小屋各面面积，分别表示小屋每个面的长度，宽度；表示房屋的各个面5.问题一模型的建立与求解5.1 模型准备屋顶接收辐射总量的计算东北南水平面倾斜面太阳图1 太阳角关系示意图由于屋顶有倾斜角，因此将屋顶作为斜面计算总的辐射量。斜面的总辐射量计算如下：朝向赤道斜面上的太阳辐射量，通常采用Klein的计算方法，倾斜面上所接受到的太阳辐射总量由直接辐射量、天空散射辐射量以及地面反射辐射量组成，即： （1）与水平面上的直接辐射量之间又如下关系： （2）对于朝向赤道的倾斜面，可以由下式确定： （3）其中：是当地维度，是倾斜角，是太阳赤纬角，为倾斜面与水平面上直接辐射量之比。水平面上的日落时角为： (4)斜面上的日落时角为： (5)通过对附件1数据的分析，并借助matlab软件，首先统计数据并运用方差法计算得每月的典型日期，然后对12个月的典型日期，运用方差法同样可以得到当年的典型日期，以此作为.太阳赤纬角为： (6)其中：为日期序号，由方差法确定。Hay模型认为倾斜面上天空散射辐射量是由太阳光盘的辐射量和其余天空穹顶均匀分布的散射辐射量两部分组成，可表达为： (7)其中：和分别为水平面上直接和散射辐射量；为倾斜面与水平面上直接辐射量之比，为大气外层水平面上太阳辐射量；为倾角。大气外层的水平辐射量可由下式求出： (8)其中：为太阳常数。由于本题是将光伏电池以贴附方式安装，电池板与屋顶的缝隙较小，因此地面反射辐射量可以忽略不计综上所述，屋顶上太阳辐射总量的表达式为： (9)5.2模型建立电池组件安装收益规划模型本题仅考虑贴附安装方式对小屋的部分表面进行铺设，由于要保证发电量最大、成本最低，因此建立一个多目标优化模型。其目标为全年太阳能最大发电总量和建设成本最低。5.2.1目标函数由于所有光伏组件在010年效率按100%，1025年按照90%折算，25年后按80%折算。因此各组件的35年发电时间可折合成等效发电时间为：10+15*90%+10*80%=24.3年。小屋屋顶及墙面上电池组件在等效时间内的总发电量可按下式计算： （10） 其中，为使用光伏电池的数目，为屋顶及房屋各面所受的辐射强度，为所需电池板的面积，为光伏电池的转换效率。由于小屋光伏电池板在其寿命时间内不存在任何的损坏维修，因此电池板的花费成本为： （11） 其中，为电池每峰瓦的价格，为电池组件功率。由公式（1）、（2）可列出最大收益模型为： （12）5.2.2约束条件（1）对电池板铺设总面积的限制：（2）电池板尺寸的限制：（3）逆变器成本限制：5.2.3贴附式电池组件安装收益规划模型综合上述目标函数及约束条件，可列出贴附式电池组件安装收益规划模型如下： (13)5.3模型求解 5.3.1太阳辐射总量的计算根据模型准备数据以及题目所给资料，运用matlab（R2012a 7.14.0.739）软件求得小屋房顶及各面的太阳辐射总量如下：表1 小屋房顶及各面太阳辐射总量墙面东南西北顶辐射量578700104340087280024320037679005.3.2安装方案的确定根据屋顶及房屋各面的辐射总量以及目标函数：运用matlab（R2012a 7.14.0.739）软件可求得小屋屋顶及房屋各面安装电池板后收益如下： 表2 房屋各面收益一览表屋面北面南面东面西面顶面收益（元）-2667.96810197.873732.14010718.53370647.9由表1中的数据可看出北面墙在不考虑逆变器的情况下已经亏损，因此不考虑在北面安装电池板。其余各面在不考虑逆变器的情况下所需的电池板数如下：表3 南、东、西、顶四面需要安装的最大电池板数目电池板型号南面东面西面顶面A300042B10001B20001B40001B50001C11116180C21000C31000C76331C80200C91010C100220注：表中未显示的电池型号均未使用根据表3中小屋屋顶及各面墙选择电池板的数目从而确定各面墙输入逆变器的功率，通过对附件5的数据处理得到各逆变器的标准输入功率如下：表4 各逆变器的标准输入功率逆变器SN1SN2SN3SN4SN5SN6SN7SN8SN9输入功率600120011522304350455203300561011110逆变器SN10SN11SN12SN13SN14SN15SN16SN17SN18输入功率2222011002200334455668338106482600026000根据表4可得到各面墙安装电池后产生功率如下：表5 各墙面安装电池板后产生功率屋面南面东面西面顶面墙产生功率（W）2192.7322805.663133.03512093.61考虑到经济以及安装的复杂程度，本文优先选择电量接近的1个逆变器。通过表3及表4的数据统计可得小屋屋顶及各面墙应选择的逆变器如下：表6 屋顶及各墙面选择的逆变器屋面南面东面西面顶面所选逆变器SN12SN7SN7SN17南面安装电池示意图如下：图2 南面墙电池安装示意图注：绿色矩形为C1单晶硅光伏电池组件，蓝色矩形为C9非晶硅薄膜光伏电池组件。南立面电池组件连接分组阵列及功率示意图如下：图3 南立面电池组件连接及功率示意图顶面电池组件安装图如下：图4 顶面电池组件安装示意图注：黄色矩形为A3单晶硅光伏电池组件。顶面电池组件连接分组阵列及功率示意图如下：图5 顶面电池组件连接分组阵列及功率示意图在以上方案下，计算得35年的南墙总发电量32818.37kwh ，屋顶总发电量为 1019279 kwh， 1052097.37 kwh，6.问题二模型的建立与求解6.1最佳倾角的求解1-4朝向赤道斜面上的太阳辐射量，通常采用Klein的计算方法，倾斜面上所接受到的太阳辐射总量由直接辐射量、天空散射辐射量以及地面反射辐射量组成，即：与水平面上的直接辐射量之间又如下关系：对于朝向赤道的倾斜面，可以由下式确定：其中：是当地维度，是倾斜角，是太阳赤纬角。水平面上的日落时角为：斜面上的日落时角为：太阳赤纬角为：其中：为日期序号，由方差法确定。Hay模型认为倾斜面上天空散射辐射量是由太阳光盘的辐射量和其余天空穹顶均匀分布的散射辐射量两部分组成，可表达为：其中：和分别为水平面上直接和散射辐射量；为倾斜面与水平面上直接辐射量之比，为大气外层水平面上太阳辐射量；为倾角。大气外层的水平辐射量可由下式求出： (14)其中，为太阳常数。地面反射辐射常量的表达式为： (15)其中：为地面反射率，一般情况下，倾斜面上太阳辐射总量的表达式为：(16)由上式直接求出最大太阳辐射倾斜角的数学表达式有些困难，在冬半年，太阳赤纬角为负值。如当地纬度为，通常总有，因此倾斜面上日落时角和水平面日落时角相等，这时可直接推导出最佳倾角的数学表达式。此时，水平面上的总辐射量为： (17)对于确定的地点，其太阳辐射量及地面反射率等均为常数，将对求导，并令，得： (18)当直接辐射量为0时，最佳倾角不等于0，这正是天空散射辐射量各向异性模型所需的结果。对于夏半年，最佳倾角往往接近于0，甚至为负值。6.2电池板排列间距的计算影子倍率法计算太阳高度角并计算间距5图6 电池阵列安装示意图一般在水平面垂直竖立的高为的木杆，其南北方向的角度为，方位角为，那么影子倍率可由下式表示： (19)其中： 影子倍率阵列高度影子长度太阳高度角太阳方位角由图6以及公式可以得出：由于本文研究的房屋是正南方向，因此太阳方位角为0.根据球面三角函数分析认为太阳高度与观测者的地理纬度、太阳赤纬和方位角有着一定的关系，他们之间的关系式为： (20)其中为太阳高度角，为时角，为当时的太阳赤纬，为当地的纬度(大同的纬度为)，根据最佳倾斜角的求解公式可知： 斜面上的日落时角为： (21)太阳赤纬角为： (22)由于本文研究的太阳方位角为0，则，太阳高度角公式变为： (23)得出计算太阳高度角的基本公式： (24)通过以上模型及附件给出的数据可以计算出电池板间的阵列间距，但是由于电池板的数量级规格较多，使用此模型计算大小不同电池板的阵列间距极为困难，且经过复杂计算后，准确度也有很大偏差。因此，本文使用阴影转换的方法，可以方便的计算出阵列间距及所需电池板数目。方法示意图如下：图7 阴影面及间距求法在本题中已经得到最佳倾斜角，通过倾斜角及电池板面积可以求得电池板在屋顶平面上产生的阴影，进而可以将电池板的面积等效为电池板在屋顶平面上产生的阴影的面积。此阴影面积可以等效为问题一中贴附电池板的面积，从而将求阵列间距问题转化为问题一中求贴附电池板的面积。运用此方法可建立模型求得在屋顶及房屋各面上所需的电池板数目。6.3架空式电池组件安装收益规划模型本题是在问题一的基础上求得最佳倾斜角建立优化模型。本题同样是多目标优化模型，其目标为全年太阳能最大发电总量和建设成本最低。根据第一问的贴附式电池组件安装收益规划模型可以改进的到架空式电池组件安装收益规划模型如下： (25)6.4模型求解根据前面公式(25)运用matlab（R2012a 7.14.0.739）求得根据建立模型运用matlab（R2012a 7.14.0.739）软件算出架空后各面安装电池板收益如下：表7 将电池组件架空安装后房顶及各面的收益屋面北面南面东面西面顶面收益（元）-2667.968120213732.14010718.53530975同问题一，经计算分析北墙亏损，因此北墙不考虑安装电池组件，根据房屋屋顶及除去背面的其他面吸收的辐射量，可确定在架空安装方式下选择的电池型号及数目如下：表8 选择的电池板数量及型号南面东面西面顶面A300055B10001B30002C11816180C73331C80200C91010C102220房顶及各面安装电池板后产生功率以及相应逆变器选择情况如下：表9 房顶及各面安装电池板后产生功率屋面南面东面西面顶面墙产生功率（W）2192.7322805.663133.03512093.61表10 房顶及各面选择的逆变器屋面南面东面西面顶面所选逆变器SN12SN5SN5SN17屋顶安装电池板的架构图如下：图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 kwh，经济总收益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=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;endsum1num1Lingo程序：question_1.lg4model:sets:!FUSHE*：各个面的辐射强度，POWER_OUT：年输出电量;xinxi/1.24/:L,H,P,N,num,M,FUSHE1,FUSHE2,FUSHE3,FUSHE4,FUSHE5,POWER_OUT;about/1.5/:AREA,LONG,WIDTH;endsetsdata: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-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*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);endLingo程序：question_2.lg4model:sets:!FUSHE*：各个面的辐射强度，POWER_OUT：年输出电量,SR:电池板的阴影面积;xinxi/1.24/:L,H,P,N,num,M,FUSHE1,FUSHE2,FUSHE3,FUSHE4,FUSHE5,POWER_OUT;about/1.5/:AREA,LONG,WIDTH;endsetsdata: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 111