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

太阳能小屋外表面电池组件优化铺设的设计PAGEPAGE26太阳能小屋的设计摘要本文探讨了太阳能小屋外表面电池组件优化铺设的问题，并在此基础上给出了太阳能小屋的改进方案。针对问题一，本文首先计算了屋顶的太阳能辐射总量（见表1），在考虑小屋全年光伏电池发电总量最大、单位发电量费用最小的基础上建立了贴附式电池组件安装收益规划模型，对小屋屋顶和南面进行了部分铺设，得到了小屋表面电池板铺设方案（见表3，图2，图4）和逆变器的选择方案（见表6，图3，图5），并计算出了在当前方案下小屋的发电总量为1052097.37kw·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/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由于所有光伏组件在0～10年效率按100%，10～25年按照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（R2012a39）软件求得小屋房顶及各面的太阳辐射总量如下：表1小屋房顶及各面太阳辐射总量墙面东南西北顶辐射量578700104340087280024320037679005.3.2安装方案的确定根据屋顶及房屋各面的辐射总量以及目标函数：运用matlab（R2012a39）软件可求得小屋屋顶及房屋各面安装电池板后收益如下：表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.37kw·h，屋顶总发电量为1019279kw·h，1052097.37kw·h，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（R2012a39）求得根据建立模型运用matlab（R2012a39）软件算出架空后各面安装电池板收益如下：表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.37kw·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;clearall;closeall;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;%重置计数whilei<=8760%附件月或年结束位置ifb1(i)>=fushe1sum1(1)=sum1(1)+b1(i);num1(1)=num1(1)+1;endifb2(i)>=fushe1sum1(2)=sum1(2)+b2(i);num1(2)=num1(2)+1;endifb3(i)>=fushe1sum1(3)=sum1(3)+b3(i);num1(3)=num1(3)+1;endifb4(i)>=fushe1sum1(4)=sum1(4)+b4(i);num1(4)=num1(4)+1;endifb5(i)>=fushe1sum1(5)=sum1(5)+b5(i);num1(5)=num1(5)+1;endi=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=1580195615801651165019561650195614821640195619561668130013211414140014003106156159208181645;H=80899180899299199199199199299299299210001100711111411001100355180355355355712;N=16.8416.6418.7016.514.9815.1116.2116.3915.9814.815.9815.214.996.996.176.355.846.493.633.633.663.664.134.27;P=14.914.914.914.914.914.912.512.512.512.512.512.512.54.84.84.84.84.84.84.84.84.84.84.8;M=2153252002702452952653202102402802952501005810090100448121250;POWER_OUT=266.919418.9063253.55339.09323.7234396.7488341.5691408.7622279.888312.174373.184386.507331.0125168.3695.942163.35145.4041649.3458.8218.6928.03524.0396.25;LONG=71007100101001010010100;WIDTH=38003800320032006511.53;AREA=2423000026980000192353094444000060870453;FUSHE1=524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83581235.26581235.26581235.26581235.26581235.26581235.26581235.26581235.26581235.26581235.26581235.26;FUSHE2=795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47872801.59872801.59872801.59872801.59872801.59872801.59872801.59872801.59872801.59872801.59872801.59;FUSHE3=1008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161044044.061044044.061044044.061044044.061044044.061044044.061044044.061044044.061044044.061044044.061044044.06;FUSHE4=132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71243181.05243181.05243181.05243181.05243181.05243181.05243181.05243181.05243181.05243181.05243181.05;FUSHE5=375710037571003757100375710037571003757100375710037571003757100375710037571003757100375710037679003767900376790037679003767900376790037679003767900376790037679003767900;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=1580195615801651165019561650195614821640195619561668130013211414140014003106156159208181645;H=80899180899299199199199199299299299210001100711111411001100355180355355355712;N=16.8416.6418.7016.514.9815.1116.2116.3915.9814.815.9815.214.996.996.176.355.846.493.633.633.663.664.134.27;P=14.914.914.914.914.914.912.512.512.512.512.512.512.54.84.84.84.84.84.84.84.84.84.84.8;M=2153252002702452952653202102402802952501005810090100448121250;POWER_OUT=266.919418.9063253.55339.09323.7234396.7488341.5691408.7622279.888312.174373.184386.507331.0125168.3695.942163.35145.4041649.3458.8218.6928.03524.0396.25;LONG=71007100101001010010100;WIDTH=38003800320032006511.53;AREA=2423000026980000192353094444000060870453;FUSHE1=524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83524961.83581235.26581235.26581235.26581235.26581235.26581235.26581235.26581235.26581235.26581235.26581235.26;FUSHE2=795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47795135.47872801.59872801.59872801.59872801.59872801.59872801.59872801.59872801.59872801.59872801.59872801.59;FUSHE3=1008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161008343.161044044.061044044.061044044.061044044.061044044.061044044.061044044.061044044.061044044.061044044.061044044.06;FUSHE4=132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71132000.71243181.05243181.05243181.05243181.05243181.05243181.05243181.05243181.05243181.05243181.05243181.05;FUSHE5=375710037571003757100375710037571003757100375710037571003757100375710037571003757100375710037679003767900376790037679003767900376790037679003767900376790037679003767900;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/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/0.9*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基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系统基于单片机的控制系统在PLC虚拟教学实验中的应用研究基于单片机系统的网络通信研究与应用基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究基于单片机的模糊控制器在工业电阻炉上的应用研究基于双单片机冲床数控系统的研究与开发基于Cygnal单片机的μC/OS-Ⅱ的研究基于单片机的一体化智能差示扫描量热仪系统研究HYPERLINK"/