太阳能小屋设计模型

1、太阳能小屋的设计摘要本文根据天空散射辐射各向异性的 Hay 模型，计算出小屋顶部倾斜面上在不同时间 的太阳辐射量， 通过假定相邻整数时刻之间的辐射量满足线性的关系算出太阳对每个面 上每天和全年的辐射总量， 又根据太阳在不同时刻对小屋其它四个面的辐射量通过积分 计算得到每一种电池在小屋的五个面上一年的产生的总电量和每个时刻对应的功率值。 而问题一中要挑选出的电池与逆变器可以是一个双目标函数模型的优化解， 所以根据建 立的模型与题中所给出的要求，可相应的列出目标函数的约束条件，用软件计算即可求 出小屋每面所需要的电池和逆变器的型号和数量。具体结果如下 矚慫润厲钐瘗睞枥庑赖。型号东面南面西面北面顶面

2、电池板A2A312161643C8C101516SN11SN41逆SN61变SN71器SN111SN131SN141问题二主要是找出斜面的最佳倾角 0 ，使其满足全年接收的太阳光辐射总量达到最 大。倾角为 的光伏板接收的辐射量是一个关于变量 的函数，对倾角 求导并令导数 值为 0 查找出斜面辐射量的计算公式， 对 求导令其为 0得到一个函数的导数与倾角 的关系式，通过对该函数本身的函数表达式求导后的值应该相等的条件可以求解出最佳 的倾角 0 =37.6 度。再将求出的倾角值代回问题一所建立的模型中进行计算，得出在最 佳倾角的情况下，满足目标函数的最优解。 聞創沟燴鐺險爱氇谴净。对于问题三要设计

3、一个小屋，并对所设计小屋的外表面优化铺设光伏电池，我们可 以先把小屋的各种需要确定的指标全部假设出来，再根据给出要求列出不等式方程组， 将此方程组也代回到问题一的模型中进行计算， 得出小屋每面所需要的电池和逆变器的 型号和数量，即为所求问题的最优解。 残骛楼諍锩瀨濟溆塹籟。关键字： 双目标函数模型；太阳辐射；太阳能电池板；逆变器； Hay 模型1 问题重述近年来，大规模开发利用太阳能已经成为国际上应对气候变化、调整能源结构的重 要举措。我国有着丰富的太阳能资源和悠久的太阳能利用历史，近两年更是出现了蓬勃 的局面，太阳能发电，太阳能房屋等项目正在我国各地飞速发展。太阳能房屋的推广应 用对于节约常

4、规能源、减少环境污染、改善人们的生活水平有十分重要的意义，其中光 伏发电是太阳能房屋的核心技术。光伏发电要求在建筑物表面（屋顶及外墙）铺设光伏 电池，光伏电池组将太阳能转换成直流电，所产生的直流电经过逆变器转换成220V 交流电供家庭适用。不同种类的光伏电池每峰瓦的价格差别很大，且每峰瓦的实际发电效 率或发电量还受诸多因素的影响，如太阳辐射强度、光线入射角、环境、建筑物所处的 地理纬度、地区的气候与气象条件、安装部位及方式（贴附或架空）等。因此，在太阳 能小屋的设计中，研究光伏电池在小屋外表面的优化铺设是很重要的问题，科学合理的 铺设方法对太阳能小屋的设计起着关键作用。 酽锕极額閉镇桧猪訣锥。

5、附件 1-7 提供了相关信息。请参考附件提供的数据，对下列三个问题，分别给出小 屋外表面光伏电池的铺设方案，使小屋的全年太阳能光伏发电总量尽可能大，而单位发 电量的费用尽可能小， 并计算出小屋光伏电池 35 年寿命期内的发电总量、 经济效益（当 前民用电价按 0.5元/kWh 计算）及投资的回收年限。 在求解每个问题时， 都要求配有图 示，给出小屋各外表面电池组件铺设分组阵列图形及组件连接方式 （串、并联）示意图， 也要给出电池组件分组阵列容量及选配逆变器规格列表。 在同一表面采用两种或两种以 上类型的光伏电池组件时，同一型号的电池板可串联，而不同型号的电池板不可串联。 在不同表面上，即使是相

6、同型号的电池也不能进行串、并联连接。应注意分组连接方式 及逆变器的选配。 彈贸摄尔霁毙攬砖卤庑。问题一：根据大同市的气象数据，仅考虑贴附安装方式，选定光伏电池组件，对小 屋（见附件 2）的部分外表面进行铺设，并根据电池组件分组数量和容量，选配相应的 逆变器的容量和数量。 謀荞抟箧飆鐸怼类蒋薔。问题二：电池板的朝向与倾角均会影响到光伏电池的工作效率，请选择架空方式安 装光伏电池，重新考虑问题 1。问题三：根据附件 7 给出的小屋建筑要求，为大同市重新设计一个小屋，要求画出 小屋的外形图， 并对所设计小屋的外表面优化铺设光伏电池， 给出铺设及分组连接方式， 选配逆变器，计算相应结果。 厦礴恳蹒骈時

7、盡继價骚。1 问题分析对于问题一，将其分成三个部分进行求解。第一部分计算小屋顶部斜面受到的辐射强度， 因为附录四只给出了小屋东南西北四 个面以及水平面总辐射强度、水平面散射辐射强度、法向直射辐射强度每天在整数时刻 受到的总辐射强度值， 并没有直接给出小屋顶面的辐射强度值， 所以这里建立 Hay 模型 来对小屋顶面的辐射强度值进行计算； 第二部分计算出每种电池在一年内将太阳能转化 成电量的总量，因为太阳对每个面的辐射强度是关于时间的连续变化函数，而附件四中 给出的是整点时刻的离散值，所以需要对相邻整点之间的连续点的辐射量进行估算，这 里假定相邻整点之间的辐射量是呈线性变化的直线， 这样就便于求出

8、每个电池在每个墙 面上的功率及一年内将太阳能转化为总电量的值；第三部分则是建立双目标函数模型， 目标是求发电量的最大值与费用的最小值， 通过题目中所给的逆变器选择条件列出关于 目标函数的约束条件，进而求解出其的最优解。 茕桢广鳓鯡选块网羈泪。对于问题二，因为电池板东西南北四个方向架空所接收的太阳辐射量与贴附安装在 墙面上所接收的太阳辐射量变化不大，所以，只考虑小屋顶面的电板的架空情况，而架 空方向又可以分为朝向赤道方向和偏离赤道方向。 无论是讨论朝向赤道方向还是偏离赤 道方向，求解思路都可归纳如下： 鹅娅尽損鹌惨歷茏鴛賴。因为倾角为 的光伏板接收的辐射量是一个关于变量 的函数，查找出斜面辐射量

9、 的计算公式并对 求导令其为 0 得到一个函数的导数与倾角 的关系式，通过对该函数 本身的函数表达式求导后的值应该相等的条件可以求解出最佳的倾角。 籟丛妈羥为贍偾蛏练淨。对于问题三：按要求自己设计一个太阳能小屋，但是所设计小屋的外表面需要优化 铺设光伏电池，这样我们仍采用问题一中的第三部分的做法，使小屋的全年太阳能光伏 发电总量尽可能大，而单位发电量的费用尽可能小做为目标函数进行求解。 預頌圣鉉 儐歲龈讶骅籴。2 模型假设1. 小屋顶面朝北的斜面的辐射量相对朝南的斜面的辐射量很小，对经济的影响小 在计算过程中可忽略不计；2. 假设两整数点时间之间的太阳光对小屋每个面的辐射强度是与时间呈线性关系

10、 的；3. 假设每一年的同一时刻大同市的太阳辐射量是相同的；4. 假设每个串联阵列，并联阵列都只能与一个逆变器相连3 符号说明E 表示小屋太阳能光伏发电一年的总量； eik,j表示第i月第 j日小屋的 k面一天的辐射总量；表示在前 年内小屋的太阳能光伏发电总量；Qi表示在前 年内小屋的光伏发电总量产生的经济效益；X i, j t 表示第 i 月第 j X 2 ti,j 表示第 i 月第 j X i, j t 表示第 i 月第 j X 4 ti,j 表示第 i 月第 j X 5 tX i,j t 表示第 i 月第 j日t 时刻小屋东面的辐射强度；日t 时刻小屋西面的辐射强度；日t 时刻小屋南面的

11、辐射强度，日t 时刻小屋北面的辐射强度，日t 时刻小屋屋顶倾斜面的辐射强度，Hb表示水平面的直接辐射量，HT表示法向直射辐射量，Hd表示水平面的散射量，表示水平面的总辐射量，H0表示大气外层上太阳总辐射量，表示倾斜面上的太阳辐射总量表示光电板的倾角，表示地球表面的反射率（一般取 0.2），Rb表示倾斜面与水平面的直射量之比k1表示小屋东面k2表示小屋西面k3表示小屋南面表示小屋北面k5表示小屋顶面k1表示小屋东墙的可用面积k2表示小屋西墙的可用面积k3表示小屋南墙的可用面积k4k5表示小屋北墙的可用面积表示小屋顶面的可用面积Sxj j表示电池 xj 的面积Syr表示电池 yr 的面积Vxj表示

12、电池 xj 的端电压Vyr表示电池 yr 的端电压Utj表示序号为 tj 的逆变器允许输入的电压范围Uhl表示序号为 hl 的逆变器允许输入的电压范围Pxj (ki)表示电池 xj 在 ki面的转化功率Pyr (ki)表示电池 yr 在 ki面的的转化功率Ptj额表示序号为 tj 的逆变器的配置容量（即额定功率）Phl 额表示序号为 hl 的逆变器的配置容量（即额定功率）4 模型建立与求解4.1 问题一4.1.1 斜面受到的辐射强度的计算由于光电板朝南（即朝向赤道方向） ，确定朝向赤道倾斜面上的辐射量通常采用Klein 所提出的计算方法，认为天空散射辐射量是均匀分布， 倾斜面上的太阳辐射总量

13、H是由直接太阳辐射量 H bt ，天空散射辐射量 H dt和地面反射辐射量 Hrt 三部分组成。而Hay 模型认为倾斜面上的天空散射辐射量 H dt是由太阳光盘的辐射量和其余天空穹顶均匀分布的散射辐射量两部分组成，可表达为 渗釤呛俨匀谔鱉调硯錦。Hdt Hd H0Rb 2 1HbH01 cos4.1)式中 H b ， H d分别为水平面上的直接辐射量和散射量，为斜面倾角 , Rb倾斜面与水平面直接辐射量之比。这样求斜面上的太阳辐射量的公式可改为： 铙誅卧泻噦圣骋贶頂廡。5 H 1i5,j t HbRb Hd HHb Rb 21 1HbH011 cos H 1 cos4.2)其中H0 24Hsc

14、 1 0.033cos360n 2 0cos cos sin 0 sin sin365 0 3604.3)cos cos sin s s sin sin180Rbcos cos sin 0 0 sin sin0 0 1804.4)所以 X i5,j t 的计算可转化为对 Hb 的计算。其中： 为当地的纬度， s, 0是倾斜面和水平面的日落时角， 23.45sin 2 283465 n 为太阳的赤纬角， n为一年中的第几天，s min 0 ,arccostan a tan， H sc为太阳系数取 1353W/m2 ， 0 为水平面的日落时角，可由 0 arccos tan tan 确定。因为太阳

15、在每刻的时角公式等于 t 15 t 12 ，倾斜面的日落时刻的时角为s min 0,arccos tan a tan时角 的计算公式为：t 15 t 12 度4.5)太阳高度角 的计算公式为：sin sin sin cos cos cos4.6)太阳方位角（倾斜面方位角） 的计算公式为sin cos4.7)sincos太阳光入射角 的计算公式为：cos si n coscos si n cos si ncos cossin si n cos cos cossi n（4.8） i n c ossi n当 0, 0 时，可得水平面上的太阳入射角为cos 0 sin sin cos cos cos

16、则倾斜面与水平面上接受到的直射辐射分别为HT Hn cos ,Hb Hncos 0H n为垂直于太阳光线的倾斜面的辐射强度，由上式可得：Hbcos 0cos4.9）由式（ 5）可以计算出倾斜面上日出，日落的时刻为rs s ss sTirs, j12 1 15s ,Tis,sj12 1 15s .4.1.2 电池将太阳能转化成电量的计算 因为太阳光的辐射强度是随时间变化的连续函数，而附件 4 中只给出不同面所受太 阳光辐射量的一些离散值， 由此不妨假设已知在 t0 时刻时小屋的某一面的辐射强度值为 a，在t0 1时刻的辐射强度值为 b时用 Xik,j t 表示太阳光某日在小屋的 k （ k为小屋

17、的 东西南北顶面五面中的一个面） 的辐射强度，当时间 t t0,t0 1 这段时间区域时认为太 阳光在小屋的 k 面的辐射值满足线性关系式 擁締凤袜备訊顎轮烂蔷。X ik,j tb a t a t0 1 bt0（ 4.10）则在第i月第 j日小屋的k面一天的辐射总量 eik,j 为k i, X rsTi,j rskiXmmssTj,i rsTi,mdt,kjiX ssTi,js,jsTi,经过简单的化简可得,jss Ti jb2aTiss,jTissj, 2 cT issj , Tss c Tis,sj222m 11（ 4.12）其中 T irs,j , Tis,sj 分别表示不大于日出时刻

18、T irs, j 和日落时刻 T is,sj 最近的那个整数时刻，c a t0 1 bt0 ，则小屋一年受太阳辐射的总量可表示为Eeik,j i 1, 2, ,j5; 1, 2,或 , 30. 31 (4.13)又根据附件给出的数据及建模的要求可以提出以下两个假设：X ik,jt 80 对于型号为 A 的单晶硅电池，若辐照强度在 80,200 w / m 2 时，A 电池的转换效率为5%120 对于薄膜电池 c1,c2，若辐照强度在 200,1000 w / m2 时， c1, c2电池的转换效率为X ik,j t 21001%，8其中 为取整符号。光伏电池组件启动发电时其表面所应接受到的最低

19、辐射量限值， 单晶硅和多晶硅电池启动发电的表面总辐射量大于等于 80W / m2 、薄膜电池表面总辐射量大于等于 30W / m2 ，可得不同类型的电池一天在 ki(i 1,2,3) 面将太阳能转化成电量的表达式如下贓熱俣阃歲匱阊邺镓騷。电池 Am ：tl 1WAkm1 i,j l tl Am Xik,1j tAmdt4.14)A 5%0，X k1 t 0,80 w/ m2 i,jX ik,j t 80120Am,X ik,j1,X k1 t 80,200 w/ m2 m 1,2, ,6 i,jt 200w / m2电池 Bm ：tl 1WBkm2 i,j l tl BmXik,j2 t B

20、mdt(4.15)0，X k2 t 0,80 w /m2Bmk22 m 1, 2, , 7BmBm,X ik,j2 t80w/m2电池 Cm ：tl 1WCkm3 i,j Cm Xik,j3 t C mdt(4.16)l tli,j0，X ik,j3 t 0, 30X ik,j3 t 2100CmCm81%, X ik,j3 t 200 m 1,2Cm ,当Xik,j3 t 30,2000，Xik,j3 t 0, 30Cm,Xik,j3 t 30, 200m 3, , 11其中k1, k2,k3 东西南北顶面 ，k1,k2,k3可能相等也可能不等； Am， Bm， Cm 分别表 示产品型号为

21、Am ， Bm ， Cm的电池面积；所有光伏组件在 010年效率按 100%，1025年按照 90%折算，25年后按 80%折算。对( 4.14)，(4.15)，(4.16)式分别求和可得 到电池的一年在该面上转化的总电量。 坛摶乡囂忏蒌鍥铃氈淚。电池 Am一年在 k1面上转化的总电量为 WAkm1WAkm1 i,j ，电池Bm一年在 k2面上转化i,jk k k k的总电量 WBkm2WBkm2i,j ，电池Cm 一年在k3面上转化的总电量WCkm3WCkm3i, j 。i ,ji ,j4.1.3 建立双目标函数模型 根据所给附表及所知的数据和假设，可以按以下方法求解，具体思路如图4-1 中

22、所给出的流程图。假设在小屋 ki面有电池的串联阵列 nki个，并列阵列 mki个， j 1,2,.,nki ，第j个 串联阵列中的电池符号为 xj(ki)共有 aj(ki)个，l 1,2,.,mk ，第l个并联阵列中由 Tl(ki)个串联电路组成， 对于 r1,2.,Tl (ki ) ，第 l个并联阵列中 r个串联阵列的电池符号为 yr (ki) ，有br (ki )个。构造目标函数为 蜡變黲癟報伥铉锚鈰赘。nkinki Tl (ki )总电量： maxaj(ki)Wxj(ki)br (ki ) Wyr(ki) ，(4.17)j 1 l 1 r 118费用：Z ming1PGAmHAmg2PG

23、BmHBmg3HCmFf gf ，( 4.18)f1图 4.1 双目标函数建模过程其中， HAm,HBm,HCm 分别表示 A类单晶硅电池在五个面上的个数和， B类多晶硅电池 在五个面上的个数和， C类薄膜电池在五个面上的个数和， Ff 表示序号为 f 的逆变器 的个数， g1,g2,g3分别表示 A类单晶硅电池的价格， B类多晶硅电池的价格， C类薄膜 电池价格，g f 表示序号为 f 的逆变器的价格， PG Am ，PG Bm ，PG Cm 产品型号为 Am， Bm ，Cm 的电池的峰值功率 （组建功率）。目标函数的其约束条件为式 （4.19）（4.25）。 買鲷鴯譖昙膚遙闫撷凄。，分别为

24、电池总数约束条件nki 5 5 mki Tl （ki ）HAm HBm HCmaj（ki ）br （ki ），（4.19）j 1 i 1 i 1 l 1 r 1面积约束条件nk1nk1 Tl (k1)aj (k1) Sxjbr(k1) Syrk1j1l 1 r1nk2nk2 Tl (k2)aj (k2) Sxjbr (k2) Syrk2j1l 1 r1nk3nk3 Tl (k3)aj(k3) Sxjbr(k3) Syrk3j 1l 1 r14.20)nk4nk4 Tl (k4)aj(k4) Sxjbr(k4) Syrk4j1l 1 r1nk5nk1 Tl (k5)aj (k5 ) Sxjbr

25、(k5) Syrk5j 1l 1 r1ki表示小屋的 ki 面的面积，k1 表示小屋东面 , k 2为小屋西面 ,k3 表示小屋南面 ,k4 表示小屋北面, k5 表示小屋顶面 . Sx ， Syr为电池 xj(ki)和 yr(ki)的面积。綾镝鯛駕櫬鹕踪韦辚糴。根据附录中给出的要求我们可以得到串并联阵列与逆变器相连的约束条件， 对于串联有aj(ki) Pxj(ki) ptj额，i 1,2,.5 ； j 1,2,. nki( 4.21)aj(ki) Vxj(ki) Utj (ki)；对于并联有Tl (ki )Pyr(ki) br(ki) phl额 ，r=1br(ki) Vyr Uhlr1,r2

26、 1,2,.,Tl(ki) ，l 1,2,.mki ， yr1 iyr2 i 10% ，i max Vyr1(ki),Vyr2 (ki)4.22)4.23)4.24)4.25)且tj,hl 1,2,.,18 ， xj(ki),yr(ki) A1, A2,., A6 , B1,.B7 ,C1,.,C11 ，即为所有电池的集合 元素， Pxj (ki) ， Pyr (ki )表示电池 xj (ki)和yr (ki)在小屋 ki面上的转化功率且Pxj(ki)xjX(ik,ij)tSxj ， Pyr (ki)yrX (ik,ij)tSyr，xj 和 yr 示电池 xj(ki)和 yr (ki )的转换

27、效率。利用软件求出产品型号为 Am，Bm ，Cm的电池在 ki 面上电池的个数 nAkim ，nBkim ， nCkim 在 年内小屋的太阳能光伏发电总量 驅踬髏彦浃绥譎饴憂锦106 5 k k 7 5 k k 11 5 k knkAim WAkmi 1nBkim WBkmi 1nCkim WCkmi 1 ，0 10m 1 i 1 m 1i 1 m 1i 16 5 k k 7 5 k k 11 5 k k0.9 1nkAimWAkmi1nBkim WBkmi1nCkimWCkmi1，1025m 1i 1 m 1i 1 m 1i 16 5 kiki7 5 kiki11 5 kiki0.8 3.5

28、nkAim WAkmi 1nBkim WBkmi 1nCkim WCkmi 1m 1 i 1 m 1i 1 m 1i 125 35令35 ，可得6 5 k k 7 5 k k 11 5 k k4.26)35 31.5nkAim WAkmi 1nBkim WBkmi 1nCkim WCkmi 1 ，m1i1 m 1i1 m 1i1其中WAkmi 1 ，WBkmi 1 ,WCkmi 1 表示型号为 Am，Bm ，Cm的电池使用的第一年在 k1面上转 化的总电能。从而，前 年内小屋的光伏发电总量产生的经济效益 猫虿驢绘燈鮒诛髅貺庑。4.27)Q 0 .5 。利用 matlab软件，可求得总费用为 3

29、00880.9，经 24 年时间能回收成本。各型号的电池 板及逆变器的数量如表 4.1 所示，计算得到的逆变器的型号及数量表 4.2 所示，电池板、 逆变器的选取及费用见表 4.3。 锹籁饗迳琐筆襖鸥娅薔。表 4.1 各面各型号电池板的数量表 / 个型号东面南面西面北面顶面总数A211A3126164377C81616C101515表 4.2 各面各型号逆变器的数量表 / 个型号东面南面西面 北面 顶面SN11SN41SN61SN71SN111SN131SN14111表 4.3 电池与逆变器的价格型号单价个数合计元A24842.514842.5电池板A3298077229460C838.416

30、614.4C1057.615864SN1290012900SN4690016900SN615000115000逆变器SN710200110200SN11450014500SN1310300110300SN1415300115300总计元300880.9依据太阳能电池的型号及数量，按照逆变器连接要求，分析计算得出太阳能电池与逆变器的连接方式，在各面的铺设连接方式如下示意图： 構氽頑黉碩饨荠龈话骛。顶面铺设连接示意图东面铺设连接示意图：南面铺设连接示意图：12134.2 问题二因为电池板的朝向与倾角均会影响到光伏电池的工作效率， 选择架空方式安装光伏 电池可以提高电池板的发电量，只考虑小屋屋顶的电

31、池板的朝向与倾角。 輒峄陽檉簖疖網儂號 泶。求导并取值1、朝向赤道倾斜面的上的最佳倾角计算 由（4.2）式，因为对于确定的地点，其太阳辐射量及地面反射率等均为常量，倾角 为 的光伏板接收的辐射量是一个关于变量 的函数，所以将 X i5,j t 对 为 0 可得 尧侧閆繭絳闕绚勵蜆贅。4.28)dRb sin Hd H0 HbHH 0d 2 H0 H d H b又因为cos cos sin s s Rbcos cos sin 0 0sin sin180sin sin180求导可得:ddRb Rbdsin cos sin s s cos sinsin sin180cos cos sin 0 04.

32、29)求得arctansin sin s scos sin180cos cos sin 0 00 0 180sin sinsin sin cos sin180.(4.30)其中Hd H0 HbHH 0 .其中 2 H0 Hd H b .2、偏离赤道斜面上的太阳辐射量计算朝向赤道的倾斜面上能接收到更多的太阳辐射条件，然而在有些情况下，可能受地 形、地物的限制，无法正对赤道安装，对于偏离赤道的倾斜面上太阳辐射的计算，普遍 采用 Klein 推导，后经修正的计算方法，其中此时的 Rb 可以用下式 4 计算 :识饒鎂錕缢灩筧嚌俨淒。14Tis,sj Tirs,j sin sin cos cos sin

33、 cosRdcos sinTis,sj sinTir,sj cos cos sin sin cos cosTis,sj cosTir,sj cos sin sinss rs2 cos cos sin 0 sin sin 0 180 .（. 4.31） 对（ 4.31）式求导并与（ 4.28）联立即可算出偏离赤道斜面上的最佳倾角。将所得到的 最佳倾角带回问题一中进行计算即可得到最优化的解。 凍鈹鋨劳臘锴痫婦胫籴。经计算，求的与水平面最佳倾斜角为 37.6 度。根据最佳倾斜角及第一问的计算结果 可知，当把顶面的太阳能电板的倾斜角调整为最佳倾斜角，其它铺设方案不变时，效果 最佳，图形如下 ;恥諤銪灭

34、萦欢煬鞏鹜錦。此时所需费用不变但是只需 22 年就能回本4.3 问题三在问题三中，需要重新设计一个小屋， 并对所设计小屋的外表面优化铺设光伏电池， 即要求所设计的小屋全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能 小。假设建筑平面体型长边为 x1，短边为 x2 ，最低净空高度为 h ，最高点高度为 H ，东 西南北墙的面积分别用 S1,S2,S3,S4表示，东西南北面上的窗户面积分别用 1, 2, 3, 4 表 示，顶面的倾斜角仍用 表示，根据附件 7 给出的小屋建筑要求有 鯊腎鑰诎褳鉀沩懼統庫。1 H hS1 S2 hx2H h x2x24.32)S3 x1h x1 H hsinS

35、4 x1hH hx2 Htanh 22 tan15且有 h 2.8,H 5.4, x1 x2 74,x1 15,x2 3121 2 0.35S1 S23 0.5 (4.33)S34 0.3S4 要求所设计的小屋全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小， 则可采用问题一中的双目标函数模型进行计算， 其中(4.20)式中的 硕癘鄴颃诌攆檸攜驤蔹。k1S11,k2S 2,2k3S 3,3k4S 4, 4 k5S 5 5 (4.34)根据要求，可设计以下图形：由 Matlab 计算可得到如下铺设方式： 南面的铺设方法：16顶面的铺设方法：西面铺设连接示意图：东面铺设连接示意图：花费

36、18203.5 元，大约 19 年回本。5 模型评价与改进5.1 模型评价1 模型优点：(1) .变量之间的联系紧凑，误差范围减少。(2) 建立的模型将问题分解为几个步骤来求，增大了能算出答案的概率2 模型缺点：17（1）文中涉及的量很多，在一定程度会使计算相对比较复杂，不易理解；（2）模型未考虑到门和窗的不可拆开性 （即门和窗的面积不能由分布在各处的小块面 积组成）；5.2 模型改进上述问题也可以看成一个矩形和圆形的装填问题，采用最大穴度算法来计算。首先 根据问题一中的模型先判断出每个面的需要铺的电池种类和个数， 假设计算后得到在 ki 面上需要铺设的电池为 R1,R2, , Rn ，每一块

37、电池 Rj的宽为 j，长为 l j ，让小屋的每个 墙面的右下角和直角坐标平面的原点重合， 墙面的短边为宽记为 Wi ，墙面的长边为宽记 为 Li 。现在的问题是要把这 n块电池尽可能多的放到 ki 墙上，也就是说，要使墙面的面 积利用率尽可能大，并满足如下限制条件： 阌擻輳嬪諫迁择楨秘騖。（1）任一放入容器的电池块的边必须和容器平行； （2）任意两块放入墙面的电池块没有嵌入。几个概念： 格局：在某一时刻，容器中已经互不嵌入地放了若干物体，还有若干物体在容器 外面等待放入，这种状态，称为一个格局。若容器中还没有放任何物体。此时的格局称 为初始格局；若全部物体都按放置规则放进了容器或容器中已放不

38、下任何物体。此时的 格局称为终止格局。 氬嚕躑竄贸恳彈瀘颔澩。占角动作：在某一格局下，若即将放进去的物体与容器中已有物体包中的某两个 物体贴边则称放此物体的这种作法为一个占角动作 。釷鹆資贏車贖孙滅獅赘。 两物体间的距离：在两电池 i, j 中各取一点 i i, j j ，设这两点的欧氏 距离为 dij ，则两物体间的距离 d min i i , j j dij 。 占角动作的穴度将电池 R 按占角动作放进矩形容器之后， R与形成这个角以外的其它物体的距离 的最小值为 dmin ，则此占角动作的穴度 Cde 1 dmin ， q为一性价比常数。用引算法算法表示为Greedy packing p

39、rocedure（Con此） 过程中的参数 Con 表示当前格局。在格局 Con 之下，枚举所有的占角动作，并计算每一个占角动作的穴度；While 存在占角动作 do按占角动作选择原则选择唯一一个占角动作； 作这个占角动作，作完之后得到新的格局，仍记为 Con;在格局 Con 之下，枚举所有的占角动作，并计算每一个占角动作的穴度；End while返回格局 Con;End18参考文献1日本太阳能学会，太阳能利用新技术 M ，北京：科学出版社， 2009。2 高虹，张爱黎，新型能源技术与应用 M ，北京：国防工业出版社， 2007。3 韩菲，潘玉良，苏忠贤，固定式太阳能光伏板最佳倾角设计方法研究

40、 J,工程设计学 报 2009， 16:348352. 怂阐譜鯪迳導嘯畫長凉。4 杨金焕,毛家俊 ,陈中华,不同方位倾斜面上太阳辐射及最佳倾角的计算 J，上海交通大 学学报， 2002，36:10311036 谚辞調担鈧谄动禪泻類。5 申政,吕建,杨洪兴 ,蒋英,太阳辐射角接受面最佳倾角的计算与分析 J,天津城市建设学 院学报 2009，15:6164。 嘰觐詿缧铴嗫偽純铪锩。6 陈端兵,黄文奇 ,求解矩形和圆形装填问题的最大穴度算法 J，计算机工程与应用， 2007，43:13。熒绐譏钲鏌觶鷹緇機库。附录Matlab 程序：function fankuijieguo(weizhi)num =

41、 xlsread(辐射强度 .xls,逐时气象参数 ,weizhi);n1 total1 = geshu(num,0,30)n2 total2 = geshu(num,30,80)n3 total3 = geshu(num,80,200)n4 total4 = geshu(num,200,10000) function n total = geshu(num,shu1,shu2) n1 = find(shu1 = num & num shu2); n = length(n1);total = 0;for i = 1:1:ntotal = total + num(n1(i);end functi

42、on x = Xg() %求当地维度的正弦、余弦值weidu = 40.1; %weidu 表示当地的纬度，角度制 weidu = weidu/180*pi; %换成弧度制 sin_wd = sin(weidu);cos_wd = cos(weidu);%求斜面倾斜角度的正弦、余弦值a = atan(1200/6400);%a为斜面倾斜角度，弧度制cos_a = cos(a);19sin_a = sin(a);%求每天太阳赤纬角的正弦、余弦、正切值n = 1:365; %表示天数b = 23.45 .* sin(2*pi*(284+n)/365);%b 表示太阳赤纬角，角度制 鶼渍螻偉阅劍鲰腎邏蘞。b = b ./180*pi; %换成弧度制sin_b = sin(b); cos_b = cos(b); tan_b = tan(b);%求一天中每时刻所对应的时角 ,及对