太阳能小屋设计模型

太阳能小屋的设计摘要本文根据天空散射辐射各向异性的Hay模型，计算出小屋顶部倾斜面上在不同时间的太阳辐射量，通过假定相邻整数时刻之间的辐射量满足线性的关系算出太阳对每个面上每天和全年的辐射总量，又根据太阳在不同时刻对小屋其它四个面的辐射量通过积分计算得到每一种电池在小屋的五个面上一年的产生的总电量和每个时刻对应的功率值。而问题一中要挑选出的电池与逆变器可以是一个双目标函数模型的优化解，所以根据建立的模型与题中所给出的要求，可相应的列出目标函数的约束条件，用软件计算即可求出小屋每面所需要的电池和逆变器的型号和数量。具体结果如下问题二主要是找出斜面的最佳倾角a，使其满足全年接收的太阳光辐射总量达到最0大。倾角为a的光伏板接收的辐射量是一个关于变量a的函数，对倾角a求导并令导数值为0查找出斜面辐射量的计算公式，对a求导令其为0得到一个函数的导数与倾角a的关系式，通过对该函数本身的函数表达式求导后的值应该相等的条件可以求解出最佳的倾角a0=37.6度。再将求出的倾角值代回问题一所建立的模型中进行计算，得出在最佳倾角的情况下，满足目标函数的最优解。对于问题三要设计一个小屋，并对所设计小屋的外表面优化铺设光伏电池，我们可以先把小屋的各种需要确定的指标全部假设出来，再根据给出要求列出不等式方程组，将此方程组也代回到问题一的模型中进行计算，得出小屋每面所需要的电池和逆变器的型号和数量，即为所求问题的最优解。关键字：双目标函数模型；太阳辐射；太阳能电池板；逆变器；Hay模型近年来，大规模开发利用太阳能已经成为国际上应对气候变化、调整能源结构的重要举措。我国有着丰富的太阳能资源和悠久的太阳能利用历史，近两年更是出现了蓬勃的局面，太阳能发电，太阳能房屋等项目正在我国各地飞速发展。太阳能房屋的推广应用对于节约常规能源、减少环境污染、改善人们的生活水平有十分重要的意义，其中光伏发电是太阳能房屋的核心技术。光伏发电要求在建筑物表面（屋顶及外墙）铺设光伏电池，光伏电池组将太阳能转换成直流电，所产生的直流电经过逆变器转换成220V交流电供家庭适用。不同种类的光伏电池每峰瓦的价格差别很大，且每峰瓦的实际发电效率或发电量还受诸多因素的影响，如太阳辐射强度、光线入射角、环境、建筑物所处的地理纬度、地区的气候与气象条件、安装部位及方式（贴附或架空）等。因此，在太阳能小屋的设计中，研究光伏电池在小屋外表面的优化铺设是很重要的问题，科学合理的铺设方法对太阳能小屋的设计起着关键作用。附件1-7提供了相关信息。请参考附件提供的数据，对下列三个问题，分别给出小屋外表面光伏电池的铺设方案，使小屋的全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小，并计算出小屋光伏电池35年寿命期内的发电总量、经济效益（当前民用电价按0.5元/kWh计算）及投资的回收年限。在求解每个问题时，都要求配有图示，给出小屋各外表面电池组件铺设分组阵列图形及组件连接方式（串、并联）示意图，也要给出电池组件分组阵列容量及选配逆变器规格列表。在同一表面采用两种或两种以上类型的光伏电池组件时，同一型号的电池板可串联，而不同型号的电池板不可串联。在不同表面上，即使是相同型号的电池也不能进行串、并联连接。应注意分组连接方式及逆变器的选配。问题一：根据大同市的气象数据，仅考虑贴附安装方式，选定光伏电池组件，对小屋（见附件2）的部分外表面进行铺设，并根据电池组件分组数量和容量，选配相应的逆变器的容量和数量。问题二：电池板的朝向与倾角均会影响到光伏电池的工作效率，请选择架空方式安装光伏电池，重新考虑问题1。问题三：根据附件7给出的小屋建筑要求，为大同市重新设计一个小屋，要求画出小屋的外形图，并对所设计小屋的外表面优化铺设光伏电池，给出铺设及分组连接方式，选配逆变器，计算相应结果。1问题分析对于问题一，将其分成三个部分进行求解。第一部分计算小屋顶部斜面受到的辐射强度，因为附录四只给出了小屋东南西北四个面以及水平面总辐射强度、水平面散射辐射强度、法向直射辐射强度每天在整数时刻受到的总辐射强度值，并没有直接给出小屋顶面的辐射强度值，所以这里建立Hay模型来对小屋顶面的辐射强度值进行计算；第二部分计算出每种电池在一年内将太阳能转化成电量的总量，因为太阳对每个面的辐射强度是关于时间的连续变化函数，而附件四中给出的是整点时刻的离散值，所以需要对相邻整点之间的连续点的辐射量进行估算，这里假定相邻整点之间的辐射量是呈线性变化的直线，这样就便于求出每个电池在每个墙面上的功率及一年内将太阳能转化为总电量的值；第三部分则是建立双目标函数模型，目标是求发电量的最大值与费用的最小值，通过题目中所给的逆变器选择条件列出关于目标函数的约束条件，进而求解出其的最优解。对于问题二，因为电池板东西南北四个方向架空所接收的太阳辐射量与贴附安装在墙面上所接收的太阳辐射量变化不大，所以，只考虑小屋顶面的电板的架空情况，而架空方向又可以分为朝向赤道方向和偏离赤道方向。无论是讨论朝向赤道方向还是偏离赤道方向，求解思路都可归纳如下：因为倾角为a的光伏板接收的辐射量是一个关于变量a的函数，查找出斜面辐射量的计算公式并对a求导令其为0得到一个函数的导数与倾角a的关系式，通过对该函数本身的函数表达式求导后的值应该相等的条件可以求解出最佳的倾角。对于问题三：按要求自己设计一个太阳能小屋，但是所设计小屋的外表面需要优化铺设光伏电池，这样我们仍采用问题一中的第三部分的做法，使小屋的全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小做为目标函数进行求解。2模型假设小屋顶面朝北的斜面的辐射量相对朝南的斜面的辐射量很小，对经济的影响小在计算过程中可忽略不计；假设两整数点时间之间的太阳光对小屋每个面的辐射强度是与时间呈线性关系的；假设每一年的同一时刻大同市的太阳辐射量是相同的；假设每个串联阵列，并联阵列都只能与一个逆变器相连3符号说明E表示小屋太阳能光伏发电一年的总量；#)、表示第i月第j日小屋的k面一天的辐射总量;G,j)①人表示在前人年内小屋的太阳能光伏发电总量；QiX(1)(t)«,j)X(2)(t)«,j)X(3)(t)«,j)X(4)(t)«,j)X(5)(t)k,j)表示在前人年内小屋的光伏发电总量产生的经济效益;表示第i月第j日t时刻小屋东面的辐射强度;表示第i月第j日t时刻小屋西面的辐射强度;表示第i月第j日t时刻小屋南面的辐射强度，表示第i月第j日t时刻小屋北面的辐射强度，表示第i月第j日t时刻小屋屋顶倾斜面的辐射强度，表示水平面的直接辐射量，表示法向直射辐射量，表示水平面的散射量，表示水平面的总辐射量，表示大气外层上太阳总辐射量，表示倾斜面上的太阳辐射总量表示光电板的倾角，表示地球表面的反射率（一般取P=0.2），表示倾斜面与水平面的直射量之比表示小屋东面表示小屋西面表示小屋南面表示小屋北面表示小屋顶面表示小屋东墙的可用面积bk2表示小屋西墙的可用面积bk3表示小屋南墙的可用面积表示小屋北墙的可用面积bk5Sxbk5Sxj表示小屋顶面的可用面积表示电池x的面积Syr表示电池七的面积表示电池x的端电压VyVyr表示电池七的端电压表示序号为t］的逆变器允许输入的电压范围UhP(UhP(k.)Xj1Py(k)表示序号为hl的逆变器允许输入的电压范围表示电池xj在匕面的转化功率表示电池七在匕面的的转化功率表示序号为t］的逆变器的配置容量（即额定功率）PhPhl额表示序号为ht的逆变器的配置容量（即额定功率）4模型建立与求解4.1问题一4.1.1斜面受到的辐射强度的计算由于光电板朝南（即朝向赤道方向），确定朝向赤道倾斜面上的辐射量通常采用Klein所提出的计算方法，认为天空散射辐射量是均匀分布，倾斜面上的太阳辐射总量Hf是由直接太阳辐射量H所，天空散射辐射量H出和地面反射辐射量H三部分组成。而Hay模型认为倾斜面上的天空散射辐射量H出是由太阳光盘的辐射量和其余天空穹顶均匀分布的散射辐射量两部分组成，可表达为

H广HdH｝—bIH0(1+cosaH广HdH｝—bIH0(1+cosa)(4.1)X(5))(X(5))(t)=HR+气(1+cosa)+2PH(1-cosa)(4.2)其中24

H=——H0兀sc(1+0.033cosI365360n/cos©cosPsin3+笠％sin©sin24

H=——H0兀sc(1+0.033cosI365360n/cos©cosPsin3+笠％sin©sinp'人0360(4.3)cos(©-a)cosPsin3+3-[^。sin(©-a)sinPb兀cos©cosPsin3+3-^^^sin©sinP(4.4)所以X((5),(t)的计算可转化为对%的计算。其中：©为当地的纬度3?,30是倾斜面和,一^…一...，一(284+n\水平面的日落时角，P=23.45sin2k-I365)为太阳的赤纬角，〃为一年中的第几天，3,=min妇,arccos(一tan(©-a)tanP)｝，H为太阳系数取1353W/m2，30为水平面的sc日落时角，可由①=arccos(-tan©tanP)确定。因为太阳在每刻的时角公式等于3,=15(t-12)，倾斜面的日落时刻的时角为3s=min妇,arccos(一tan(©-a)tanP)｝。时角3的计算公式为:3=15(t-12)(度)(4.5)太阳高度角中的计算公式为:sin中=sin©-sinP+cos©-cosP-cos3(4.6)太阳方位角(倾斜面方位角)y的计算公式为.-sin3-cosPsinY=(4.7)太阳光入射角0的计算公式为:co0=(sifnca-s巾osasijlycosPsin(4.8)+(cofcox+sfinain/ycosPc(+gocoasyinfsin3cossin当a=0,y=0时，可得水平面上的太阳入射角为cos0=sinfsinP+cosfcosPcos①则倾斜面与水平面上接受到的直射辐射分别为cos00HnHn为垂直于太阳光线的倾斜面的辐射强度，由上式可得:cos0Hb=cos0Hb=溢0THt-由式（5）可以计算出倾斜面上日出，日落的时刻为(4.9)Trs(i,j)(①)=121一I15)(①)=121+rI15)4.1.2电池将太阳能转化成电量的计算因为太阳光的辐射强度是随时间变化的连续函数，而附件4中只给出不同面所受太阳光辐射量的一些离散值，由此不妨假设已知在,时刻时小屋的某一面的辐射强度值为a，在'+1时刻的辐射强度值为b时用X（k）（t）表示太阳光某日在小屋的k（k为小屋的东西南北顶面五面中的一个面）的辐射强度，当时间tG（t0,t0+1）这段时间区域时认为太阳光在小屋的k面的辐射值满足线性关系式（4.10）X(k)(t)=(b-a)t+a(t0+1)-bt0则在第i月第j日小屋的k面一天的辐射总量e"、为（i,（4.10）e(k)=即+1x(k)(t%t+电-TX(k)(t%t+"Fx(k)(t)dt(i,j)奖j)经过简单的化简可得m=-Te(k)=即+1x(k)(t%t+电-TX(k)(t%t+"Fx(k)(t)dt(i,j)奖j)经过简单的化简可得m=-TC)-+1m(ij)j)]（4.11）zxb-ae(k)=(i,j)2Trs(ij)-+1(ij)「(j-+cTss+当（-(i,j)-2；Tss-(ij)(?j)+1-T(rs)+c^Tss--(ij),-)b-a-务)-Xm=-¥j)]+1-T.?(m+)1其中TrsI,\TssL(i,j)」，L(i,j)」则小屋一年受太阳辐射的总量可表示为E=Z以七(i=1,2.,.,jA-1-,2或，3031(4.13)i,j,k又根据附件给出的数据及建模的要求可以提出以下两个假设：①对于型号为A的单晶硅电池，若辐照强度在［80,200w/m2］时，A电池的转换效率为5%n-120X(k)(t)-805%n-120②对于薄膜电池c,c，若辐照强度在［200,1000w/m2］时，c,c电池的转换效率为1212「X］"一2

门—・1%，8其中［］为取整符号。光伏电池组件启动发电时其表面所应接受到的最低辐射量限值，单晶硅和多晶硅电池启动发电的表面总辐射量大于等于80W/m2、薄膜电池表面总辐射量大于等于30W/m2，可得不同类型的电池一天在k(i=1,2,3)面将太阳能转化成电量的表达式如下电池A:mWk^

A

m1(i,j)=£tj1nA•X(k1)(t)bdtAm(i,j)Am1tl0,X(k1)(t)g[0,80)w/m2(i,j)X(k)(t)—805%.120n

A

m(4.14)-■,X(k1)(t)e[80,200)w/m2m=1,2,--,6(i,j),电池A:mWk^

A

m1(i,j)=£tj1nA•X(k1)(t)bdtAm(i,j)Am1tl0,X(k1)(t)g[0,80)w/m2(i,j)X(k)(t)—805%.120n

A

m(4.14)-■,X(k1)(t)e[80,200)w/m2m=1,2,--,6(i,j),X(k1)(t)>200w/m2电池B:mWk2(i,j)=£'j1nA•X(k2)(t)。dtBmBm(,.,»B，m1t(4.15)An=<Bm0,X(k2)(t)g[0,80)w/m2

n,X(k2)(t)>80w/m2Bm(i,j)m'=1,2,•,7电池C:mAnCm"nCm，dtC,、Ctm(i,j)"mlt

l0,X(k3)(t)e[0,3)0(i,j)「XWt)]-2,、^^•1%,X(k3)(t)>2008(ij)n,当X(^3)(t)e[30,200)m〃=1,2(4.16)示产品型号为Am,Bm,C〃，的电池面积；所有光伏组件在0〜10年效率按100%,10〜25年按照90%折算，25年后按80%折算。对(4.14),(4.15),(4.16)式分别求和可得到电池的一年在该面上转化的总电量。电池A一年在k面上转化的总电量为WAk1=EWk1(i,j),电池B一年在k面上转化m1AAm2"i，jm的总电量Wk=EWk2(i,j)，电池C一年在k面上转化的总电量Wk=EWk3(i,j)。B,B,m3mmi,j4.1.3建立双目标函数模型根据所给附表及所知的数据和假设，可以按以下方法求解，具体思路如图4-1中所给出的流程图。TOC\o"1-5"\h\z假设在小屋k面有电池的串联阵列n个，并列阵列mk个，j£〈,2,...吃}，第j个iii串联阵列中的电池符号为X(k)共有a(k)个，lM,2,...,m}，第l个并联阵列中由T(k)jijikt1i个串联电路组成，对于re{(1,2...,T(")},第l个并联阵列中r个串联阵列的电池符号为y(k)，有b(k)个。构造目标函数为riri(4.17)总电量：Q=£a(k)W(k)+EZb(k)W(k),maxjix.iriyij=1jl=1r=1'

(4.18)费用：Z=gP(A)H+gP(B)H+gH+丈Fg，

min1Gm%2Gm3C(4.17)(4.18)f=1图4.1双目标函数建模过程TOC\o"1-5"\h\z其中，H^,^，W分别表示A类单晶硅电池在五个面上的个数和，B类多晶硅电池mm"m"在五个面上的个数和，C类薄膜电池在五个面上的个数和，Ff表示序号为f的逆变器的个数，g,g,g分别表示A类单晶硅电池的价格，B类多晶硅电池的价格，C类薄膜123电池价格，g,表示序号为f的逆变器的价格，P（A），P（B），P（C）产品型号为A，fGmGmGmmBm，匕,的电池的峰值功率（组建功率）。目标函数的其约束条件为式（4.19）——（4.25）。，分别为电池总数约束条件（4.19）H+H+H=££a(k)+£ZZb(k)，AmBmCm，jirij=1i=1i=1l=1r=1

面积约束条件/Sa(k).S+S区b(k)S<bj1Xr1yr7=1I=1（4.19）面积约束条件/Sa(k).S+S区b(k)S<bj1Xr1yr7=1I=1r=1Sa(k)-S+£Eb(k)Sj2xr2y.1j/11'j=1l=1r=1£a(k).S+SZb(k)Sj3Xjr3yr7=1l=1r=1Sa(k).S+S区)b(k)・Sj4%r4y,j=1l=1r=1Sa(k).S+S旧b(k)Sj5Xjr5yrIj=1l=1r=1<bk2<ck3<bk4<bk5(4.20)i小屋北面,k表示小屋顶面.S，S为电池x(k)和y(k)的面积。5x.yjiri根据附录中给出的要求我们可以得到串并联阵列与逆变器相连的约束条件，对于串联有a(k)-P(k)a(k)-P(k)<p，i=1,2,...5；j=1,2,...nji%itj额k,a(k)-VgU(k)；jiXj(、)tji对于并联有S'P(k)-b(k)<p，、ririh额r=1b(k)-VgU；"yrhlr,rg(1,2,...,T(k))，12lilg(1,2,...m)V*-V(k)

<1^2maxV(k),V(k)<10%(4.21)(4.22)(4.23)(4.24)(4.25)yr1Iyir2且t,hgTOC\o"1-5"\h\z(1,2,...,18)，x(k),y(k)g{A,A，…,A,B,...B,C,...,C}，即为所有电池的集合jiri12617111元素，P(k)，P(k)表示电池x(k)和y(k)在小屋k面上的转化功率且x.iyijiriiP(k)=n(x)X(k)(t)s，P(k)=n(y)X(k)(t)S且t,hgx1j(i,j)%yrir(i,j)yr门(x.)和门(y)示电池x(k)和y(k)的转换效率。利用软件求出产品型号为A，B，Cjrjirimmm"的电池在k面上电池的个数n[mnk.Bmnki在人年内小屋的太阳能光伏发电总量Cm，

①广〃**"心(1)+**的电池在k面上电池的个数n[mnk.Bmnki在人年内小屋的太阳能光伏发电总量Cm，①广〃**"心(1)+**吒W(1)+**AmAmm=1i=1,BmBm，m=1i=1nkWk「(1)[,0<X<10m〃=1i=1C"J①=(0.91叫**nkI一一"mm=1i=1WkX1)+**nkWk(1)+**nkW<(1)|,10<X<25Amm，=1i=1B".…CC|m"=1i=1%①广（0.813.5）普*nkiAmm=1i=1Wf(1)+**"(1)+**n#W『(1)m,mm„m"m"m'=1i=1„Cmm=1i=1=31.5【X£以皿（1）+E切以吭（1）+四*35-AmAmm=1i=1,BmBmm=1i=1nkW(1)〔，m"=1i=1Hl(4.26)C,，的电池使用的第一年在k面上转mf1其中W^k（.（1）,吟（1）,吟（1）表示型号为AC,，的电池使用的第一年在k面上转mf1化的总电能。从而，前人年内小屋的光伏发电总量产生的经济效益g=0.①入。（4.27）利用matlab软件，可求得总费用为300880.9,经24年时间能回收成本。各型号的电池板及逆变器的数量如表4.1所示，计算得到的逆变器的型号及数量表4.2所示，电池板、逆变器的选取及费用见表4.3。表4.1各面各型号电池板的数量表/个型号东面南面西面北面顶面总数A211A3126164377C81616C101515表4.2各面各型号逆变器的数量表/个型号东面南面西面北面顶面SN11SN41SN61SN71SN111SN131SN141

表4.3电池与逆变器的价格型号单价个数合计/元A24842.514842.5电池板A3298077229460C838.416614.4C1057.615864SN1290012900SN4690016900SN615000115000逆变器SN710200110200SN11450014500SN1310300110300SN1415300115300总计/元300880.9依据太阳能电池的型号及数量，按照逆变器连接要求，分析计算得出太阳能电池与逆变器的连接方式，在各面的铺设连接方式如下示意图：顶面铺设连接示意图东面铺设连接示意图:南面铺设连接示意图:

西面铺设连接示意图:三维铺设连接示意图:4.2问题二因为电池板的朝向与倾角均会影响到光伏电池的工作效率，选择架空方式安装光伏电池可以提高电池板的发电量，只考虑小屋屋顶的电池板的朝向与倾角。1、朝向赤道倾斜面的上的最佳倾角计算由（4.2）式，因为对于确定的地点，其太阳辐射量及地面反射率等均为常量，倾角为a的光伏板接收的辐射量是一个关于变量a的函数，所以将X户）、0）对a求导并取值。,j）为0可得4.2问题二因为电池板的朝向与倾角均会影响到光伏电池的工作效率，选择架空方式安装光伏电池可以提高电池板的发电量，只考虑小屋屋顶的电池板的朝向与倾角。1、朝向赤道倾斜面的上的最佳倾角计算由（4.2）式，因为对于确定的地点，其太阳辐射量及地面反射率等均为常量，倾角为a的光伏板接收的辐射量是一个关于变量a的函数，所以将X户）、0）对a求导并取值。,j）为0可得dRsina[H(H-H)-pHH-dot2(H「H：)%（4.28）又因为cosR=——

b(©-a)cosPsin①+①・-^-sin(©-a)sinP

180cos©cosPsin①+①--丸sin©sinP0180求导可得:dR=Rbsin(©-a)cosPsin①bda一①-1^cos(©-a)sinPcos©cosPsin①0兀+①-180sin©sinP（4.29）a=arctansin©sin①一①-Kcos©sinPss180u-cos©cosPsin①+①・-^sinI00180in©sinP+①.兀sin©sinP+cos©sin①s180s(4.30)其中咋［H以0-气）了HH°】.2（H0+H）气2、偏离赤道斜面上的太阳辐射量计算朝向赤道的倾斜面上能接收到更多的太阳辐射条件，然而在有些情况下，可能受地形、地物的限制，无法正对赤道安装，对于偏离赤道的倾斜面上太阳辐射的计算，普遍采用Klein推导，后经修正的计算方法，其中此时的Rb可以用下式团计算：

«ss-Trs)sinP(sinMosa-cos。sinacosy)180C,j)(i,j)J1c.兀•1.c)-2cos。cosPsin①+——sin。sinP"0180)-1+cosPtin/ss-sinTrsAcos^cosa+sin。sinacosy)CO,j、j)osTssJ1c.兀•1.c)-2cos。cosPsin①+——sin。sinP"0180)-1经计算，求的与水平面最佳倾斜角为37.6度。根据最佳倾斜角及第一问的计算结果可知，当把顶面的太阳能电板的倾斜角调整为最佳倾斜角，其它铺设方案不变时，效果最佳，图形如下；此时所需费用不变但是只需22年就能回本。4.3问题二在问题三中，需要重新设计一个小屋，并对所设计小屋的外表面优化铺设光伏电池，即要求所设计的小屋全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小。假设建筑平面体型长边为气，短边为气，最低净空高度为h，最高点高度为H，东西南北墙的面积分别用S,S,S,S表示，东西南北面上的窗户面积分别用&,&,&,&表12341234示，顶面的倾斜角仍用a表示，根据附件7给出的小屋建筑要求有S=S=hx+1(H-h)x=(H+h)x1222222H-h\S=xh+x311sinaS=xh+、；(H—h)2+tana)且有h>2.8,H<5.4,x-x<74,x<15,x>31212务伟-°-35&12p3<0.5(4.33)S3S=xh+、；(H—h)2+tana)要求所设计的小屋全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小，则可采用问题一中的双目标函数模型进行计算，其中(4.20)式中的b=S一6,c=S一；c=SYbS—goS一6(4.34)k111k222k333kq44^555根据要求，可设计以下图形：由Matlab计算可得到如下铺设方式:南面的铺设方法：5令A3^.•••••[A3^-4sN13]A3^••…

顶面的铺设方法:西面铺设连接示意图:东面铺设连接示意图:花费18203.5元，大约19年回本。5模型评价与改进5.1模型评价1模型优点：.变量之间的联系紧凑，误差范围减少。.建立的模型将问题分解为几个步骤来求，增大了能算出答案的概率。2模型缺点：文中涉及的量很多，在一定程度会使计算相对比较复杂，不易理解；模型未考虑到门和窗的不可拆开性(即门和窗的面积不能由分布在各处的小块面积组成)；5.2模型改进上述问题也可以看成一个矩形和圆形的装填问题，采用最大穴度算法来计算。首先根据问题一中的模型先判断出每个面的需要铺的电池种类和个数，假设计算后得到在ki面上需要铺设的电池为R,R,,R，每一块电池R的宽为①，长为l，让小屋的每个12njjj墙面的右下角和直角坐标平面的原点重合，墙面的短边为宽记为W,墙面的长边为宽记为L。现在的问题是要把这n块电池尽可能多的放到£墙上，也就是说，要使墙面的面积利用率尽可能大，并满足如下限制条件：任一放入容器的电池块的边必须和容器平行；任意两块放入墙面的电池块没有嵌入。几个概念：格局：在某一时刻，容器中已经互不嵌入地放了若干物体，还有若干物体在容器外面等待放入，这种状态，称为一个格局。若容器中还没有放任何物体。此时的格局称为初始格局；若全部物体都按放置规则放进了容器或容器中已放不下任何物体。此时的格局称为终止格局。占角动作：在某一格局下，若即将放进去的物体与容器中已有物体包中的某两个物体贴边.则称放此物体的这种作法为一个占角动作O两物体间的距离：在两电池A,A中各取一点中gA,中gA，设这两点的欧氏ijiijj距离为d，则两物体间的距离d=min(d)oij叩Ai，％GA,ij占角动作的穴度将电池R按占角动作放进矩形容器之后，R与形成这个角以外的其它物体的距离的最小值为dmin，则此占角动作的穴度Cd=l-±,q为一性价比常数。用引算法算法表小为Greedypackingprocedure(Con)此过程中的参数Con表示当前格局。在格局Con之下，枚举所有的占角动作，并计算每一个占角动作的穴度；While存在占角动作do按占角动作选择原则选择唯一一个占角动作；作这个占角动作，作完之后得到新的格局，仍记为Con;在格局Con之下，枚举所有的占角动作，并计算每一个占角动作的穴度；Endwhile返回格局Con;End参考文献日本太阳能学会，太阳能利用新技术[M],北京：科学出版社，2009。高虹，张爱黎，新型能源技术与应用[M],北京：国防工业出版社，2007。韩菲，潘玉良，苏忠贤，固定式太阳能光伏板最佳倾角设计方法研究[J],工程设计学报2009，16:348〜352.杨金焕，毛家俊，陈中华,不同方位倾斜面上太阳辐射及最佳倾角的计算[J]，上海交通大学学报，2002，36:1031~1036申政，吕建，杨洪兴,蒋英,太阳辐射角接受面最佳倾角的计算与分析[J]，天津城市建设学院学报2009，15:61〜64。陈端兵，黄文奇,求解矩形和圆形装填问题的最大穴度算法[J]，计算机工程与应用，2007，43:1〜3。附录Matlab程序：functionfankuijieguo(weizhi)num=xlsread('辐射强度.xls','逐时气象参数',weizhi);[n1total1]=geshu(num,0,30)[n2total2]=geshu(num,30,80)[n3total3]=geshu(num,80,200)[n4total4]=geshu(num,200,10000)function[ntotal]=geshu(num,shu1,shu2)n1=find(shu1<=num&num<shu2);n=length(n1);total=0;fori=1:1:ntotal=total+num(n1(i));endfunctionx=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);sin_a=sin(a);%求每天太阳赤纬角的正弦、余弦、正切值n=1:365;%表示天数b=23.4