数学建模大赛组题.pdf

《数学建模大赛组题.pdf》由会员分享，可在线阅读，更多相关《数学建模大赛组题.pdf（22页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、-20122012 高教社杯全国大学生数学建模竞赛高教社杯全国大学生数学建模竞赛承承诺诺书书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料） ，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的

2、论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等） 。我们参赛选择的题号是（从 A/B/C/D 中选择一项填写） ： B我们的参赛报名号为（如果赛区设置报名号的话） ：所属学校（请填写完整的全名） ：贵州师范大学参赛队员 (打印并签名) ：1.宋家贵 2.樊佐举 3.李红指导教师或指导教师组负责人 (打印并签名)：教练组日期： 2012年 09 月 10 日赛区评阅编号（由赛区组委会评阅前进行编号）：-20122012 高教社杯全国大学生数学建模竞赛高教社杯全国大学生数学建模竞赛编编 号号 专专 用用 页页评阅人评分备注赛区评阅编号（由赛区组委会评阅前

3、进行编号）：赛区评阅记录（可供赛区评阅时使用）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：-太阳能小屋的设计摘要随着太阳能在日常生活中的普及，研究太阳能光伏发电系统也因此具有重要的意义。为了解决问题一，首先对太阳能电池进行筛选，这是一个最优化问题，通过对小屋最大经济效率的计算，进而得到太阳能电池板的最优组合。为了使小屋的发电总量、经济效益最大， 从而决定了各电池组件之间的数量及连接方式， 选择逆变器的数量和容量。对于问题二，为了实现太阳能光伏电池板最佳倾角的设计及优化1，通过MATLAB建立倾斜放置的电池组件接收太阳辐射模型，计算得到光伏板上的辐射

4、能，进而建立倾斜面上总辐射的函数， 该函数是关于斜面倾斜角的一个函数， 通过求解该函数的极大值，从而得到光伏太阳能电池板的最佳倾角。对于问题三，为了使小屋的经济效益最高，采用最优化问题解决策略，制定了小屋的建造方式，因此使所限制的条件能够得以实现。关键词：最优化方法；光伏发电系统；太阳能辐射能；最佳倾角-一、 问题重述在设计太阳能小屋时，需在建筑物外表面（屋顶及外墙）铺设光伏电池，光伏电池组件所产生的直流电需要经过逆变器转换成 220V 交流电才能供家庭使用，并将剩余电量输入电网。不同种类的光伏电池每峰瓦的价格差别很大，且每峰瓦的实际发电效率或发电量还受诸多因素的影响，如太阳辐射强度、光线入射

5、角、环境、建筑物所处的地理纬度、地区的气候与气象条件、安装部位及方式（贴附或架空）等。因此，在太阳能小屋的设计中，研究光伏电池在小屋外表面的优化铺设是很重要的问题。附件 1-7 提供了相关信息。请参考附件提供的数据，对下列三个问题，分别给出小屋外表面光伏电池的铺设方案，使小屋的全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小， 并计算出小屋光伏电池 35 年寿命期内的发电总量、 经济效益 （当前民用电价按 0.5 元/kWh 计算）及投资的回收年限。在求解每个问题时，都要求配有图示，给出小屋各外表面电池组件铺设分组阵列图形及组件连接方式（串、并联）示意图，也要给出电池组件分组阵列容量

6、及选配逆变器规格列表。在同一表面采用两种或两种以上类型的光伏电池组件时，同一型号的电池板可串联，而不同型号的电池板不可串联。在不同表面上，即使是相同型号的电池也不能进行串、并联连接。应注意分组连接方式及逆变器的选配。问题 1：请根据山西省大同市的气象数据，仅考虑贴附安装方式，选定光伏电池组件，对小屋（见附件 2）的部分外表面进行铺设，并根据电池组件分组数量和容量，选配相应的逆变器的容量和数量。问题 2：电池板的朝向与倾角均会影响到光伏电池的工作效率，请选择架空方式安装光伏电池，重新考虑问题 1。问题 3：根据附件 7 给出的小屋建筑要求，请为大同市重新设计一个小屋，要求画出小屋的外形图，并对所

7、设计小屋的外表面优化铺设光伏电池，给出铺设及分组连接方式，选配逆变器，计算相应结果3。二、问题分析2.1 对问题一的分析根据附件 4 的气象数据初步判断出东、西、南、北、以及顶部的总的辐射强度。1）计算比较 A1A6，B1B7，C1C11 每一种电池在每一个平面上的发电量，从而算出每一个面上总的经济效益，以此为标准来确定铺设的太阳能电池组件的类型。根据附件 2 的每一面的几何数据确定每一面上可铺设电池的面积，由此确-2）-3）量。2.2 对问题二的分析由附件 4 给出的数据信息，根据Hay 提出的天空散射辐射各向异性模型，首先计算出倾斜面上总的太阳辐射强度（E） ：E定铺设电池的数量。根据附件

8、 1 电池的连接要求来选择逆变器的数量和容4 Eb EdEb为倾斜面上的直接太阳辐射总量，Ed为倾斜面上的散射辐射总量，此外，太阳的入射角在很大的程度上影响到倾斜面上对太阳辐射总量的吸收，因此斜面吸收的总的太阳辐射是倾斜角的函数，即：EN EN()对这个函数对进行求导其值为 0，即dEN() 0d解此方程算出的即为最佳倾斜角，进而求出太阳的入射角5。2.3 对问题三的分析对于问题三，由于其限制条件比较多，所以要同时实现其所给条件，我们运用了最优化模型对此进行解决，从而达到条件。其限定小屋使用空间高度为：建筑屋顶最高点距地面高度5.4m, 室内使用空间最低净空高度距地面高度为2.8m；建筑总投影

9、面积（包括挑檐、挑雨棚的投影面积）为74m2； 建筑平面体型长边应15m， 最短边应3m； 建筑采光要求至少应满足窗地比 （开窗面积与房间地板面积的比值，可不分朝向）0.2 的要求；建筑节能要求应满足窗墙比（开窗面积与所在朝向墙面积的比值）南墙0.50、东西墙0.35、北墙0.30。建筑设计朝向可以根据需要设计，允许偏离正南朝向。因此，根据上面条件，我们就能得出相应数据进行模型的建立2。三、模型假设-1、假设大气层对太阳的辐射总量没有削弱作用；2、假设附件 4 中的气象数据没有错误；3、太阳能电池板的铺设不考虑屋檐或不规则的区域；4、假设天空的散射辐射是均匀分布的；四、符号说明S：倾斜面与水平

10、面之间的夹角E： 倾斜面上太阳能电池板吸收的总的太阳能辐射Eb：倾斜面上的直接太阳辐射总量Ed：倾斜面上的散射辐射总量：太阳光线的入射角E(t)：倾斜辐射系数b： 太阳直辐射透明度系数d：太阳散射辐射透明度系数： 太阳时角： 太阳赤纬角Tss、Tst：倾斜面上的日出日落时刻ph/ p0：大气修正系数AMh ：表示一定地形高度下的大气量h： 大同市的海拔高度FZ ： 35 年的发电总量L ： 35 年的经济效益五、模型的建立及求解-5.1问题一的模型的建立与求解1）首先对大同市一年内的气象数据进行分析，得到东、西、南、北、以及顶部的总的辐射强度：东面总的辐射强594213.54度（W/m2）西面

11、881238.30南面北面1050165.66261478.82顶部4415600 2） 其次计算单位面积上每块太阳能电池板产生的经济效益 （JABC） ， 其计算公式为：转换率辐射强度0.510003）计算单位面积上太阳能电池板所需成本，其计算公式为：组件功率价格面积通过上述计算得出每一种电池板在每一面上所产生的经济效益，将所得的数据进行从高到低的排序，从而得到每一面上优先考虑铺设的电池顺序。紧接着计算出每一面上电池铺设的最佳组合数，求出在这个组合下这一面一年所产生的电能及经济效益，其计算公式分别为：组合数面积单位经济效益转换率组合数面积单位经济效益 转换率0.5-成本由一年产生的电能可以得

12、出 35 年所产生的电能及经济效益，其计算公式分别为：组合数面积单位经济效益 转换率(10+150.9+100.8)组合数面积单位经济效益 转换率(10150.9+100.8) 0.5-成本5.2 模型求解1）由附录（15）程序求解得到铺设部分结果如图所示：-j jc10c10c10c10c10c10c10c10c10c10c10c10c10c10c1c1c1c1c1c1c1c1k kc1c1c1c1空白空白c1c1大门大门c1c1c1c1c1c1东面电池板铺设情况东面电池板铺设情况（图 1）c10c10c10c10c10c10c10c10c10c10c10c10c10c10c10c10c7c

13、7c8c8j j 3 3个个C10C10c8c8c6c66 6窗窗大门大门c10c10个个C10C10c10c10c10c10c10c103 3个个C10C10c10c10c1c1c1c1c10c10窗窗c1c1c2c2南面电池板的铺设情况南面电池板的铺设情况（图 2）-2）太阳能电池的连接方式部分结果如下：共共103103对对C7C7C7C7C7C7C7C7C7C7C7C7SN2SN2( (逆变器逆变器) )C6C6C6C6（北面电池连接方式图）（北面电池连接方式图）（图（图 3 3）-共共4040对对C10C10A3A3A3A3C10C10C10C10SN4SN4（逆变器）（逆变器）共共1

14、313对对C10C10（顶部电池连接方式图）（顶部电池连接方式图）（图（图 4 4）3）发电量与经济效益结果如表一；4 )每一面上电池所产生的经济效益情况如表二 ；5.2模型二的建立及求解52.1 太阳赤纬角（）太阳赤纬角（） ：日地中心线与赤道的夹角。赤道以北为正、南为负，变化范围为15 23.523.5，可由 Cooper 近似公式得到223.45N 284sin(2) (N=1,365)3603655.2.2 太阳时角（）时角描述太阳在一天内的变化情况，在当地时间（T）正午时为 0，每走经度 1h-为15，上午为正，下午为负，即15ts1215度由此得到日出日落时刻Tss12(1)Tst

15、12(115)5.2.3 太阳直辐射和散射透明度系数（b）1在不考虑大气对太阳辐射的削弱作用时，计算直接辐射透明度系数（b）的经验公式：0.56 AMh0.0095 AMhb 0.56(ee)上式中，AMh 表示一定地形高度下的大气量，计算公式为PhAMh AM0p0AM0为海平面上的大气量，其值取决于太阳高度角，AM0 (1229(614sin) )614sin2 1/2ph / p0是大气修正系数，与当地的海拔高度 h 有关，ph2880.0065h5.256 ()p0288散射辐射透明度系数（d）与直接辐射透明度系数（b）之间存在线性关系d 0.2710.294b-5.2.4 倾斜面上太

16、阳能总辐射的计算-光伏太阳能电池板所吸收的太阳能总辐射（E）包括，倾斜面上的直接太阳辐射（Eb）加上倾斜面上的散射辐射（Ed） ，即然而E Eb EdEb E(t)bcoscos sin()sincos()coscos（式中的为大同市的地理纬度，其值为40.1）电池板上太阳直接辐射强度为cos2Ed E(t)d 2sin则一年内太阳能光伏电池组件吸收的总辐射为ENEdTN1Tst365Tss可知EN是一个关于太阳能光伏电池组件倾角的函数，即EN EN()则最佳倾角dEN() 0满足d求解此方程，即可得到最佳倾角，见附录六.5.3 模型三的建立与求解5.3.1 模型三的建立-设小屋最大高度为 H

17、，最小高度为 h，长为 x，宽为 y，设小屋窗的面积为 S，则：2.8H5.4;2.8h5.4;3 x 15;3y15;0.2s/x*y 1;0=200);Ac=a(find(80=a=80);Bs=0.1621 0.1639 0.1598 0.1480 0.1598 0.1520 0.1499;Be=zeros(1,7);for j=1:7 Be(j)=(sum(Bs(j)*Bb)*0.5/1000;endBeCb=a(find(a=200);Cc=a(find(30=a=30);Ce(3)=(sum(0.0635*Cd)*0.5/1000;-Ce(4)=(sum(0.0584*Cd)*0.

18、5/1000;Ce(5)=(sum(0.0649*Cd)*0.5/1000;Ce(6)=(sum(0.0363*Cd)*0.5/1000;Ce(7)=(sum(0.0363*Cd)*0.5/1000;Ce(8)=(sum(0.0366*Cd)*0.5/1000;Ce(9)=(sum(0.0366*Cd)*0.5/1000;Ce(10)=(sum(0.0413*Cd)*0.5/1000;Ce(11)=(sum(0.0427*Cd)*0.5/1000;JABC=Ae Be Ce(1:11);CA=1580 1956 1580 1651 1650 1956;KA=808 991 808 992 99

19、1 991;ZA=215 325 200 270 245 295;MA=CA.*KA;JA=(ZA./MA)*14.9*106;CB=1650 1956 1482 1640 1956 1956 1668;KB=991 991 992 992 992 992 1000;ZB=265 320 210 240 280 295 250;MB=CB.*KB;JB=(ZB./MB)*12.5*106;CC=1300 1321 1414 1400 1400 310 615 615 920 818 1645;KC=1100 711 1114 1100 1100 355 180 355 355 355 712;

20、ZC=100 58 100 90 100 4 4 8 12 12 50;MC=CC.*KC;JC=(ZC./MC)*4.8*106CGB=JA JB JC; XD=JABC-CGB;Q=sort(XD);F1=MA MB MC;F=0 0 0 0 0 0 0 0 0 0 0 0 0 10 1 0 0 0 0 0 0 0 7 0;FZ=sum(F.*F1.*JABC)*(10+15*0.9+10*0.8)/0.5L=FZ*0.5-sum(F.*CGB)附件三a1=;%这是一个87601的矩阵，由于矩阵很大故不列出Ab=a1(find(a1=200);Ac=a1(find(80=a1=80);Bs

21、=0.1621 0.1639 0.1598 0.1480 0.1598 0.1520 0.1499;Be=zeros(1,7);for j=1:7 Be(j)=(sum(Bs(j)*Bb)*0.5/1000;endBeCb=a1(find(a=200);Cc=a1(find(30=a1=30);Ce(3)=(sum(0.0635*Cd)*0.5/1000;Ce(4)=(sum(0.0584*Cd)*0.5/1000;Ce(5)=(sum(0.0649*Cd)*0.5/1000;Ce(6)=(sum(0.0363*Cd)*0.5/1000;Ce(7)=(sum(0.0363*Cd)*0.5/10

22、00;Ce(8)=(sum(0.0366*Cd)*0.5/1000;Ce(9)=(sum(0.0366*Cd)*0.5/1000;Ce(10)=(sum(0.0413*Cd)*0.5/1000;Ce(11)=(sum(0.0427*Cd)*0.5/1000;JABC=Ae Be Ce(1:11);CA=1580 1956 1580 1651 1650 1956;KA=808 991 808 992 991 991;ZA=215 325 200 270 245 295;MA=CA.*KA;JA=(ZA./MA)*14.9*106;CB=1650 1956 1482 1640 1956 1956

23、1668;KB=991 991 992 992 992 992 1000;ZB=265 320 210 240 280 295 250;MB=CB.*KB;JB=(ZB./MB)*12.5*106;CC=1300 1321 1414 1400 1400 310 615 615 920 818 1645;KC=1100 711 1114 1100 1100 355 180 355 355 355 712;ZC=100 58 100 90 100 4 4 8 12 12 50;MC=CC.*KC;JC=(ZC./MC)*4.8*106;CGB=JA JB JC; XD=JABC-CGB;Q=sor

24、t(XD);F1=MA MB MC;F=0 0 0 0 0 0 0 0 0 0 0 0 0 2 1 0 0 0 1 1 2 0 26 0;FZ=sum(F.*F1.*JABC)*(10+15*0.9+10*0.8)/0.5L=FZ*0.5-sum(F.*CGB)附件四a2=;%这是一个87601的矩阵，由于矩阵很大故不列出Ab=a2(find(a2=200);Ac=a(find(80=a2=80);Bs=0.1621 0.1639 0.1598 0.1480 0.1598 0.1520 0.1499;Be=zeros(1,7);for j=1:7 Be(j)=(sum(Bs(j)*Bb)*0.

25、5/1000;endBe-Cb=a2(find(a2=200);Cc=a2(find(30=a2=30);Ce(3)=(sum(0.0635*Cd)*0.5/1000;Ce(4)=(sum(0.0584*Cd)*0.5/1000;Ce(5)=(sum(0.0649*Cd)*0.5/1000;Ce(6)=(sum(0.0363*Cd)*0.5/1000;Ce(7)=(sum(0.0363*Cd)*0.5/1000;Ce(8)=(sum(0.0366*Cd)*0.5/1000;Ce(9)=(sum(0.0366*Cd)*0.5/1000;Ce(10)=(sum(0.0413*Cd)*0.5/100

26、0;Ce(11)=(sum(0.0427*Cd)*0.5/1000;JABC=Ae Be Ce(1:11);CA=1580 1956 1580 1651 1650 1956;KA=808 991 808 992 991 991;ZA=215 325 200 270 245 295;MA=CA.*KA;JA=(ZA./MA)*14.9*106;CB=1650 1956 1482 1640 1956 1956 1668;KB=991 991 992 992 992 992 1000;ZB=265 320 210 240 280 295 250;MB=CB.*KB;JB=(ZB./MB)*12.5*

27、106;CC=1300 1321 1414 1400 1400 310 615 615 920 818 1645;KC=1100 711 1114 1100 1100 355 180 355 355 355 712;ZC=100 58 100 90 100 4 4 8 12 12 50;MC=CC.*KC;JC=(ZC./MC)*4.8*106;CGB=JA JB JC; XD=JABC-CGB;Q=sort(XD);F1=MA MB MC;F=0 0 0 0 0 0 0 0 0 0 0 0 0 12 0 0 0 0 1 0 0 0 11 0;FZ=sum(F.*F1.*JABC)*(10+1

28、5*0.9+10*0.8)/0.5L=FZ*0.5-sum(F.*CGB)附件五 a3=;%这是一个87601的矩阵，由于矩阵很大故不列出Ab=a3(find(a3=200);Ac=a(find(80=a3=80);Bs=0.1621 0.1639 0.1598 0.1480 0.1598 0.1520 0.1499;-Be=zeros(1,7);for j=1:7 Be(j)=(sum(Bs(j)*Bb)*0.5/1000;endBeCb=a3(find(a3=200);Cc=a3(find(30=a3=30);Ce(3)=(sum(0.0635*Cd)*0.5/1000;Ce(4)=(su

29、m(0.0584*Cd)*0.5/1000;Ce(5)=(sum(0.0649*Cd)*0.5/1000;Ce(6)=(sum(0.0363*Cd)*0.5/1000;Ce(7)=(sum(0.0363*Cd)*0.5/1000;Ce(8)=(sum(0.0366*Cd)*0.5/1000;Ce(9)=(sum(0.0366*Cd)*0.5/1000;Ce(10)=(sum(0.0413*Cd)*0.5/1000;Ce(11)=(sum(0.0427*Cd)*0.5/1000;JABC=Ae Be Ce(1:11);CA=1580 1956 1580 1651 1650 1956;KA=808

30、 991 808 992 991 991;ZA=215 325 200 270 245 295;MA=CA.*KA;JA=(ZA./MA)*14.9*106;CB=1650 1956 1482 1640 1956 1956 1668;KB=991 991 992 992 992 992 1000;ZB=265 320 210 240 280 295 250;MB=CB.*KB;JB=(ZB./MB)*12.5*106;CC=1300 1321 1414 1400 1400 310 615 615 920 818 1645;KC=1100 711 1114 1100 1100 355 180 3

31、55 355 355 712;ZC=100 58 100 90 100 4 4 8 12 12 50;MC=CC.*KC;JC=(ZC./MC)*4.8*106;CGB=JA JB JC; XD=JABC-CGB;Q=sort(XD);F1=MA MB MC;F=0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 206 0 0 0 0;FZ=sum(F.*F1.*JABC)*(10+15*0.9+10*0.8)/0.5L=FZ*0.5-sum(F.*CGB)附件六a3=;%这是一个87601的矩阵，由于矩阵很大故不列出S=atan(3/16); K_1=cos(ata

32、n(3/16)/2)2;K_2=cos(atan(12/7)/2)2; I_1=I_h+I_d*K_1 ;I_2=I_h+I_d*K_2 ; I=I_1+I_2Ab=I(find(I=200);Ac=I(find(80=I=80);Bs=0.1621 0.1639 0.1598 0.1480 0.1598 0.1520 0.1499;Be=zeros(1,7);for j=1:7 Be(j)=(sum(Bs(j)*Bb)*0.5/1000;endBeCb=I(find(I=200);Cc=I(find(30=I=30);Ce(3)=(sum(0.0635*Cd)*0.5/1000;Ce(4)=

33、(sum(0.0584*Cd)*0.5/1000;Ce(5)=(sum(0.0649*Cd)*0.5/1000;Ce(6)=(sum(0.0363*Cd)*0.5/1000;Ce(7)=(sum(0.0363*Cd)*0.5/1000;Ce(8)=(sum(0.0366*Cd)*0.5/1000;Ce(9)=(sum(0.0366*Cd)*0.5/1000;Ce(10)=(sum(0.0413*Cd)*0.5/1000;Ce(11)=(sum(0.0427*Cd)*0.5/1000;JABC=Ae Be Ce(1:11);CA=1580 1956 1580 1651 1650 1956;KA=

34、808 991 808 992 991 991;ZA=215 325 200 270 245 295;MA=CA.*KA;JA=(ZA./MA)*14.9*106;CB=1650 1956 1482 1640 1956 1956 1668;KB=991 991 992 992 992 992 1000;ZB=265 320 210 240 280 295 250;MB=CB.*KB;JB=(ZB./MB)*12.5*106;CC=1300 1321 1414 1400 1400 310 615 615 920 818 1645;KC=1100 711 1114 1100 1100 355 18

35、0 355 355 355 712;ZC=100 58 100 90 100 4 4 8 12 12 50;MC=CC.*KC;JC=(ZC./MC)*4.8*106CGB=JA JB JC; XD=JABC-CGB;Q=sort(XD);F1=MA MB MC;F=0 0 40 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 26 0;FZ=sum(F.*F1.*JABC)*(10+15*0.9+10*0.8)/0.5L=FZ*0.5-sum(F.*CGB)附件七j=1;-t=0;bata=0.1853;for bata =0.1853:0.1:pi/2for t

36、=9:17 derta=23.45*sin(2*pi*(284+j)/365*pi/180); womiga=15*(t-12)*pi/180; I=1367*(1+0.033*cos(2*pi*j/365); aerfa=pi-(asin(sin(40.1*pi/180)*sin(23.45. *sin(2*pi*(284+j)/365*pi/180)+. cos(40.1*pi/180)*cos(23.45*sin(2*pi*(284+j). /365*pi/180)*. cos(15*(t-12)*pi/180); taob=0.56*(exp(-0.56*(1229+(614*sin(

37、aerfa). 2)(1/2)-614*sin(aerfa)*. (288-0.0065*1000)/288)5.256)+exp(-0.095*. (1229+(614*sin(aerfa)2)(1/2)-614*sin(aerfa)*. (288-0.0065*1000)/288)5.256); taod=0.271-0.294*taob;endend I_d=I*taod*cos(bata)2/2*sin(aerfa) I_b=I*taob*sin(40.1*pi/180-bata*pi/180)*. sin(23.45*sin(2*pi*(284+j)/365*pi/180)+. co

38、s(40.1*pi/180-bata*pi/180)*. cos(23.45*sin(2*pi*(284+j)/365*pi/180)*. cos(15*(t-12)*pi/180); F=I_b+I_d; T_1=12*(1-acos(-tan(40.1*pi/180-bata*pi/180). *tan(23.45*sin(2*pi*(284+1)/365*pi/180)/15); T_2=12*(1+acos(-tan(40.1*pi/180-bata*pi/180). *tan(23.45*sin(2*pi*(284+1)/365*pi/180)/15); diff(int(F,bata)=0-