數(shù)學(xué)建模——太陽(yáng)能小屋的設(shè)計(jì)模型.doc

太陽(yáng)能小屋的設(shè)計(jì)模型摘要本文通過(guò)附件所提供數(shù)據(jù)的分析，建立了太陽(yáng)能小屋的光伏電池的優(yōu)化鋪設(shè)的模型，得到了貼附和架空的兩種鋪設(shè)方案，并根據(jù)要求設(shè)計(jì)了一個(gè)新的更高效的太陽(yáng)能小屋，在此基礎(chǔ)上分別求出在不同的方案下，小屋鋪設(shè)的光伏電池在35年壽命期內(nèi)的發(fā)電總量，經(jīng)濟(jì)效益，以及投資回收年限。在對(duì)給定小屋僅考慮貼附鋪設(shè)時(shí)，根據(jù)可用的逆變器，結(jié)合同一種光伏組件不同的串聯(lián)組合方式，分別給出了光伏組件的可用的串聯(lián)組合。針對(duì)四個(gè)墻面和屋頂?shù)牟煌妮椛渌?，設(shè)計(jì)了一個(gè)判定光伏組件優(yōu)先順序的判定因子。將光伏組件的鋪設(shè)歸結(jié)為一個(gè)裝箱問(wèn)題，基于光伏組件的優(yōu)先因子設(shè)計(jì)了一個(gè)貪婪算法，對(duì)光伏組件分別在小屋的墻壁和屋頂進(jìn)行了鋪設(shè)。在給定的鋪設(shè)情況下，頂部可以在7年收回成本，且最后可以有13.7萬(wàn)的收益；東面可以在22年收回成本，且最后可以有0.8萬(wàn)的收益；西面可以在15年收回成本，且最后有2.5萬(wàn)的收益；南面可以在18年收回成本，且最后有1.5萬(wàn)的收益；對(duì)于整個(gè)房屋來(lái)說(shuō)，可以在11年之內(nèi)收回成本。在考慮對(duì)給定小屋進(jìn)行光伏電池的架空鋪設(shè)時(shí)，考慮墻壁鋪設(shè)的實(shí)際情況，采用了在墻壁上進(jìn)行貼附鋪設(shè)，而在屋頂進(jìn)行架空鋪設(shè)的方案。針對(duì)屋頂?shù)募芸?，利用太?yáng)角計(jì)算法向輻射最值的方法求解太陽(yáng)輻照的最佳接受角度為，在此基礎(chǔ)上，采用貪婪算法，給出了屋頂?shù)膬?yōu)化鋪設(shè)算法。房屋的四周沒(méi)有變化的情況下，屋頂改進(jìn)后可以有14.8萬(wàn)的收益，在6年收回成本。對(duì)于整體而言，可在10年收回成本。設(shè)計(jì)太陽(yáng)能小屋的時(shí)候，根據(jù)前面計(jì)算出的每個(gè)墻面和屋頂（架空）鋪設(shè)的優(yōu)化方案，首先計(jì)算了每一面墻和屋頂?shù)膯挝幻娣e的發(fā)電量，然后建立小屋設(shè)計(jì)的非線性規(guī)劃模型，通過(guò)Lingo9.0軟件的計(jì)算，給出了太陽(yáng)能小屋的優(yōu)化方案。最后，對(duì)所建模型進(jìn)行評(píng)價(jià)和改進(jìn)，并且就太陽(yáng)能小屋的設(shè)計(jì)和建造問(wèn)題給出了具體的建議。關(guān)鍵詞：光伏電池 優(yōu)先因子 貪婪算法 非線性優(yōu)化 最佳傾角一、 問(wèn)題重述在設(shè)計(jì)太陽(yáng)能小屋時(shí)，需在建筑物外表面（屋頂及外墻）鋪設(shè)光伏電池，光伏電池組件所產(chǎn)生的直流電需要經(jīng)過(guò)逆變器轉(zhuǎn)換成220V交流電才能供家庭使用，并將剩余電量輸入電網(wǎng)。不同種類(lèi)的光伏電池每峰瓦的價(jià)格差別很大，且每峰瓦的實(shí)際發(fā)電效率或發(fā)電量還受諸多因素的影響，如太陽(yáng)輻射強(qiáng)度、光線入射角、環(huán)境、建筑物所處的地理緯度、地區(qū)的氣候與氣象條件、安裝部位及方式（貼附或架空）等。因此，在太陽(yáng)能小屋的設(shè)計(jì)中，研究光伏電池在小屋外表面的優(yōu)化鋪設(shè)是很重要的問(wèn)題。針對(duì)該問(wèn)題，我們需要研究如下的問(wèn)題：（1）根據(jù)山西省大同市的氣象數(shù)據(jù)，僅考慮貼附安裝方式，選定光伏電池組件，對(duì)小屋的部分外表面進(jìn)行鋪設(shè)，并根據(jù)電池組件分組數(shù)量和容量，選配相應(yīng)的逆變器的容量和數(shù)量。（2）：電池板的朝向與傾角均會(huì)影響到光伏電池的工作效率，請(qǐng)選擇架空方式安裝光伏電池，重新考慮問(wèn)題1。（3）：根據(jù)附件7給出的小屋建筑要求，為大同市重新設(shè)計(jì)一個(gè)小屋，要求畫(huà)出小屋的外形圖，并對(duì)所設(shè)計(jì)小屋的外表面優(yōu)化鋪設(shè)光伏電池，給出鋪設(shè)及分組連接方式，選配逆變器，計(jì)算相應(yīng)結(jié)果。二、問(wèn)題分析對(duì)于在綠色小屋上鋪設(shè)太能電池板最優(yōu)解的問(wèn)題，我們可以設(shè)定為是多因素限制的一個(gè)多變量的規(guī)劃問(wèn)題。太陽(yáng)能電池由電池板和逆變器構(gòu)成，兩個(gè)相互影響，只允許相同型號(hào)的光伏組件進(jìn)行串聯(lián)，不同型號(hào)的光伏組件只能并聯(lián)。電池板本身大小固定，很難完全鋪設(shè)，這些都是限制因素，都要在模型中加以考慮。根據(jù)逆變器的參數(shù)調(diào)整設(shè)計(jì)電池組件分組陣列串并聯(lián)的方式以滿足相應(yīng)的輸出電壓和總功率，研究電池板型號(hào)對(duì)應(yīng)于逆變器的關(guān)系來(lái)確定電池板的組件；利用優(yōu)先因子的辦法來(lái)給組件排序，優(yōu)先因子的獲取依據(jù)光伏電池的大小，轉(zhuǎn)換效率以及價(jià)格等因素獲得，從而得到對(duì)于每個(gè)方向的組件的優(yōu)先順序。結(jié)合優(yōu)先因子，設(shè)計(jì)了一個(gè)貪婪算法來(lái)給每面墻進(jìn)行鋪設(shè)，以獲得比較好的鋪設(shè)方案。對(duì)于最佳太陽(yáng)輻射傾斜角的問(wèn)題，我們采用了循環(huán)的方式，結(jié)合大同市的數(shù)據(jù)，用逐個(gè)嘗試的方法來(lái)求出一個(gè)最優(yōu)角度。在問(wèn)題2的求解基礎(chǔ)上，給出了一個(gè)每面墻和屋頂?shù)膯挝幻娣e的產(chǎn)出電量，然后以此為基礎(chǔ)，給出了一個(gè)設(shè)計(jì)最優(yōu)太陽(yáng)能小屋的非線性優(yōu)化模型，并求出了結(jié)果。三、基本假設(shè)1.串聯(lián)的光伏電池水平放置。2.由于逆變器是一個(gè)將直流電轉(zhuǎn)化為交流電的電子設(shè)備，不宜放在陽(yáng)光直射的地方，也不易放在溫度較高，溫度變化較大的地方，所以我們假設(shè)逆變器的安裝都是在室內(nèi)，不會(huì)在外觀上影響電池板的排布3.為避免無(wú)序排列光伏電池，采用由逆變器控制光伏電池組件，以達(dá)到有序排列的目的。4. 光伏電池的鋪設(shè)建立在35年內(nèi)有經(jīng)濟(jì)收益的基礎(chǔ)上。5.光伏電池轉(zhuǎn)換率低于原有轉(zhuǎn)換率的5%時(shí)，光伏電池產(chǎn)生的電量忽略不計(jì)。6.太陽(yáng)能小屋的設(shè)計(jì)上，不考慮門(mén)的方向。四、變量說(shuō)明E36524-一年內(nèi)每天每小時(shí)的某一面的輻射量m35-35個(gè)組件對(duì)應(yīng)的轉(zhuǎn)化率s35-35個(gè)組件對(duì)應(yīng)的面積W -發(fā)電量-電池對(duì)應(yīng)的轉(zhuǎn)化率d-當(dāng)?shù)氐碾妰r(jià) 組件對(duì)應(yīng)逆變器的平均價(jià)格 組件對(duì)應(yīng)的成本太陽(yáng)時(shí)時(shí)角赤緯角太陽(yáng)高度角 -太陽(yáng)能小屋的長(zhǎng)寬高 房頂?shù)母叨?屋頂?shù)拈L(zhǎng)度 太陽(yáng)能小屋的門(mén)的長(zhǎng)寬，各個(gè)墻面的長(zhǎng)賀寬度五、模型建立和求解5.1 組件優(yōu)先因子模型首先將同類(lèi)型的光伏電池進(jìn)行適當(dāng)?shù)拇?lián)，基于電池的大小和墻面的尺寸，我們這里僅對(duì)行1、2、3個(gè)同類(lèi)型的光伏電池進(jìn)行串聯(lián)，然后對(duì)每一種串聯(lián)的組合求解優(yōu)先因子。逆變器存在最大允許電壓和最大電流，而每塊電池板本身也有輸出的額定電壓和額定電流，所以電流板在連接到逆變器上就會(huì)存在限制條件。通過(guò)計(jì)算最大的串聯(lián)電壓來(lái)設(shè)計(jì)組件。通過(guò)計(jì)算一個(gè)逆變器電流的限制來(lái)控制組件的并聯(lián)，實(shí)際上，對(duì)于電壓相同的逆變器，下一個(gè)逆變器的電流都比前一個(gè)大，更準(zhǔn)確的說(shuō)，對(duì)于電壓相同的逆變器，若干的逆變器可以下一個(gè)逆變器來(lái)取代，已達(dá)到降低成本的效果逆變器允許電壓一個(gè)串聯(lián)兩個(gè)串聯(lián)三個(gè)串聯(lián)SN12132c6c8c9c10SN22132c6c8c9c10SN34264a1a2a3a6b2b4b5c2c11SN44264a1a2a3a6b2b4b5c2c11SN54264a1a2a3a6b2b4b5c2c11SN64264a1a2a3a6b2b4b5c2c11SN799150c4c11a1a2a3a4a5a6b1b2b3b4b5b6b7c6SN899150c4c11a1a2a3a4a5a6b1b2b3b4b5b6b7c6SN999150c4c11a1a2a3a4a5a6b1b2b3b4b5b6b7c6SN1099150c4c11a1a2a3a4a5a6b1b2b3b4b5b6b7c6SN11180300c4c2c3c5SN12180300c4c2c3c5SN13180300c4c2c3c5SN14180300c4c2c3c5SN15180300c4c2c3c5SN16180300c4c2c3c5SN17250800c1c1SN18330800c4根據(jù)上表來(lái)設(shè)計(jì)組件出35個(gè)組件出來(lái)。如圖組件名稱(chēng)電池板長(zhǎng)寬對(duì)應(yīng)的逆變器組件1310355組件2615355組件3920355組件4818355組件51580808組件61956991組件71580808組件81956991組件91956991組件101640992組件111956992組件121321711組件131645712組件1414001100長(zhǎng)寬（）組件15329071216451424組件162800110014002200組件172600110013002200組件18474080815802424組件19586899119562973組件20495399216512976組件21495099116592973組件22586899119562973組件23495099116502973組件24586899119562973組件25444699214822976組件26492099216402976組件27586899219562976組件28586899219562976組件295004100016683000組件309303553101065組件31396371113212133組件324242111414143342組件334200110014003300組件343900110013003300組件354200110014003300右邊多出兩行是給出串聯(lián)時(shí)電池橫放和豎放的長(zhǎng)寬。優(yōu)先因子模型和求解方法方法如下：根據(jù)給出的35種組件，我們可以求出每個(gè)組件對(duì)應(yīng)的功率和面積， Y是關(guān)于性價(jià)比的目標(biāo)函數(shù)，代表某一個(gè)組件在一個(gè)方向上的平均面積純收益，代表上面程序算出的總發(fā)電量。由和構(gòu)建一個(gè)優(yōu)先因子模型，B是優(yōu)化因子，根據(jù)平均收益和單位發(fā)電量的比重來(lái)決定，F(xiàn)是優(yōu)化的目標(biāo)函數(shù)。（此處假設(shè)A=0.65，B=0.35）對(duì)于某一個(gè)組件來(lái)說(shuō)，五個(gè)方向都存在著這種優(yōu)化公式，我們通過(guò)循環(huán)程序來(lái)計(jì)算出某個(gè)組件在五個(gè)方向上的總功率，進(jìn)而求的單位面積上的發(fā)電功率（詳見(jiàn)附錄1）。然后求出每一個(gè)組件在每一個(gè)方向的目標(biāo)函數(shù)F，對(duì)于同一個(gè)方向的進(jìn)行排序，以確定在這個(gè)方向上優(yōu)先考慮哪一個(gè)組件，如圖（部分排序）組件名稱(chēng)電池板組件91956991組件275868992組件3439001100組件3342001100組件3242421114組件121321711頂部?jī)?yōu)先考慮以上組件；組件名稱(chēng)電池板組件3439001100組件1414001100組件3242421114組件3342001100組件121321711組件153290712組件91956991組件131645712組件3542001100東面優(yōu)先考慮以上組件；組件名稱(chēng)電池板組件91956991組件343900110013003300組件334200110014003300組件324242111414143342組件1414001100組件121321711組件172600110013002200組件354200110014003300西面優(yōu)先考慮以上組件組件名稱(chēng)電池板組件91956991組件343900110013003300組件334200110014003300組件324242111414143342組件121321711組件1414001100組件172600110013002200組件354200110014003300組件15329071216451424南面優(yōu)先考慮以上組件；對(duì)于北面和北屋頂，任何F目標(biāo)函數(shù)小于零，依次說(shuō)明在那里鋪設(shè)任何電板都是虧錢(qián)，不存在盈利，所以得出不用在北面和北屋頂鋪設(shè)電池板的結(jié)論。而且就北面鋪設(shè)電板所產(chǎn)生的電量也是極少，不足影響其他四個(gè)面的電量。5.2 裝箱模型可以將此問(wèn)題描述為, 給定n個(gè)組件, 分別具有n種類(lèi)型, 每個(gè)組件分別具有不同的特性(長(zhǎng)度和寬度); 有6個(gè)箱子（4面墻和屋頂）, 具有六種類(lèi)型, 每面墻和屋頂具有多個(gè)窗戶和門(mén)，在組件的鋪設(shè)上，需要規(guī)避這些門(mén)和窗戶。求將組件分鋪設(shè)在各個(gè)墻面上，使得發(fā)電量達(dá)到最大的最優(yōu)鋪設(shè)方案。這里稱(chēng)這種問(wèn)題為帶有多約束二維復(fù)雜裝箱模型，這里的多約束是指物品特性和類(lèi)型對(duì)于裝箱的約束, 以區(qū)別于傳統(tǒng)裝箱問(wèn)題空間或重量這一單一約束。由于普通的裝箱問(wèn)題為NP難問(wèn)題，這里的裝箱模型較之普通的裝箱問(wèn)題有更多的約束，它的求解難度只會(huì)比普通的裝箱問(wèn)題更難，這里采用貪婪算法，以求得該問(wèn)題的近似最優(yōu)解。5.2.1 貪婪算法 （1）首先計(jì)算每個(gè)組件的優(yōu)先因子 （2）根據(jù)優(yōu)先因子對(duì)組件進(jìn)行排序 （3）在需要鋪設(shè)的墻面上優(yōu)先鋪設(shè)優(yōu)先因子大的組件，直至該組件排不開(kāi)為止，然后選擇下一級(jí)優(yōu)先因子的組件進(jìn)行鋪設(shè)，直至整個(gè)墻面鋪滿為止。、 （3）根據(jù)墻面鋪設(shè)的組件的電壓的情況進(jìn)行并聯(lián)，以降低逆變器的價(jià)格。根據(jù)前面計(jì)算結(jié)果，基于貪婪算法，利用組件對(duì)每個(gè)面進(jìn)行鋪設(shè)，給出一些鋪設(shè)結(jié)果； 南面圖1 16個(gè)組件4和10個(gè)組件12 逆變器： SN1和SN3 南面圖2 30個(gè)組件12 和4個(gè)組件7 逆變器： SN3 南面圖3 19個(gè)組件4和7個(gè)組件13逆變器： SN1和SN3 東面圖1 2個(gè)組件9和15個(gè)組件 12逆變器： SN3 東面圖2 10個(gè)組件9逆變器： SN3 西面圖1 11個(gè)組件9逆變器： SN3 頂面圖1 25個(gè)組件9 逆變器： SN35.3 貼附的鋪設(shè)經(jīng)過(guò)計(jì)算比較，優(yōu)化解為：方向電池板構(gòu)成逆變器圖像東面10個(gè)組件92個(gè)sn4東面圖2西面11個(gè)組件9sn6西面圖1 南面19個(gè)組件4和7個(gè)組件13sn1和sn3南面圖3頂面25個(gè)組件9sn3，sn5,sn6頂面圖1在上圖給定的鋪設(shè)情況下，頂部可以有13.7萬(wàn)的收益，在7年收回成本；東面可以有0.8萬(wàn)的收益，在22年收回成本；西面可以有2.5萬(wàn)的收益，在15年收回成本；南面可以有1.5萬(wàn)的收益，在18年收回成本；對(duì)于整個(gè)房屋來(lái)說(shuō)，可以在11年之內(nèi)收回成本 5.4 架空的鋪設(shè)方式5.4.1最佳傾角問(wèn)題 參考文獻(xiàn)可以得知，對(duì)于理想狀態(tài)下，斜面太陽(yáng)能電池板在傾斜角等于所在地的緯度時(shí)會(huì)接收到最大輻射量，從而能產(chǎn)生最大的發(fā)電量。由資料可知，大同的緯度在40度左右。 實(shí)際情況要通過(guò)對(duì)大同市的輻射情況來(lái)進(jìn)行計(jì)算： 是相對(duì)于電池板的法向而言的太陽(yáng)角的分量，具體的三角函數(shù)：由以上公式，利用C語(yǔ)言設(shè)計(jì)循環(huán)算法來(lái)計(jì)算哪一個(gè)角度才能使架空的太陽(yáng)能電池板獲得最大的輻射量。（程序詳見(jiàn)附錄2）計(jì)算出的結(jié)果顯示，當(dāng)=30度時(shí)架空的電池板可以接受到最大的輻射量，屋頂一年之內(nèi)可發(fā)出479Kw的電量。（具體算法，詳見(jiàn)附錄3）房屋的四周沒(méi)有變化，屋頂改進(jìn)后可以有14.8玩的收益，在6年收回成本。對(duì)于整體而言，可在10年收回成本。5.5太陽(yáng)能小屋設(shè)計(jì)模型： 基于前面問(wèn)題的求解，我們?cè)O(shè)計(jì)一個(gè)單邊屋頂?shù)奶?yáng)能小屋，如下圖所示 為求得該房屋的具體數(shù)據(jù)，我們基于問(wèn)題2求解的結(jié)果，給出房屋的各個(gè)墻面的各個(gè)墻壁和屋頂?shù)膯挝幻娣e下的發(fā)電量；作為新設(shè)計(jì)房屋的每個(gè)墻面和屋頂?shù)膯挝幻娣e的發(fā)電量的參考依據(jù)，可以求出分別為，則得到一個(gè)非線性的規(guī)劃模型如下所示：5.41求解結(jié)果X1=24.66667； X2 = 3.000000； X3 = 5.400000； X4 = 0.000000 ； Y1 = 1.234568 ； Y2 = 1.234568； ； Y3 =； 1.234568； ；Y4 = 1.234568 ； Y5 = 1.234568 ； Y6 = 3.583568； Y7 = 1.234568； Y8 = 1.234568 。 由結(jié)果可知，所設(shè)計(jì)的太陽(yáng)能小屋是一個(gè)平房，屋頂?shù)墓夥姵貫榧芸珍佋O(shè)，房屋的門(mén)開(kāi)在北向。六、模型評(píng)價(jià)為求得光伏電池的合理的串并聯(lián)組合和鋪設(shè)方式，將該問(wèn)題轉(zhuǎn)化為一個(gè)二維的裝箱問(wèn)題，針對(duì)該裝箱問(wèn)題首先對(duì)相同類(lèi)型的光伏電池進(jìn)行串聯(lián)，基于房屋的大小和光伏電池的尺寸，給出了一個(gè)合適的串聯(lián)方式，組成了可以選用的組件。將組件鋪設(shè)的問(wèn)題轉(zhuǎn)化成一個(gè)二維的裝箱問(wèn)題，對(duì)于裝箱問(wèn)題的求解，由于該問(wèn)題的算法復(fù)雜性為NP-難，故而采用了貪婪算法。該貪婪算法的設(shè)計(jì)是基于一個(gè)組件的優(yōu)先因子，按照優(yōu)先因子的大小，優(yōu)先對(duì)值大的進(jìn)行排序。取得了很好的結(jié)果。對(duì)于架空的鋪設(shè)方法，首先根據(jù)問(wèn)題的實(shí)際情況，確立了僅在屋頂進(jìn)行架空鋪設(shè)的原則，然后建模計(jì)算出電池的在該地區(qū)的最佳傾角作為架空的依據(jù)。采用貪婪算法進(jìn)行了架空的排列。對(duì)于太陽(yáng)能小屋的設(shè)計(jì)，在問(wèn)題2求解的基礎(chǔ)上，給出一個(gè)每個(gè)墻面和屋頂?shù)膯挝幻娣e產(chǎn)出電量作為建模的依據(jù)，建立了一個(gè)非線性的優(yōu)化模型，設(shè)計(jì)出了一個(gè)太陽(yáng)能的小屋。七、建議總結(jié) 從模型中可以看出，如果考慮經(jīng)濟(jì)效益的話，屋頂鋪設(shè)光伏電池可以在較短的時(shí)間內(nèi)收回成本，而墻面的鋪設(shè)經(jīng)濟(jì)效益不高，回收年限時(shí)間長(zhǎng)。我們提出如下建議：1.對(duì)于太陽(yáng)能小屋的建設(shè)和光伏電池的鋪設(shè)，應(yīng)盡量采用屋頂鋪設(shè)的原則，并采用沿最佳的接受太陽(yáng)輻射的仰角架空鋪設(shè)的方法。2.應(yīng)加大科研的力度，提高光伏電池的轉(zhuǎn)換效率。3.發(fā)明可以切割的光伏電池板，以提高鋪設(shè)面的利用效率。4.發(fā)明透明的光伏電池板，以避免窗戶，天窗等對(duì)鋪設(shè)的干擾，在保證采光的前提下，提高鋪設(shè)面的利用效率。 八、參考文獻(xiàn)1 韓斐，潘玉良，蘇忠賢，固定式太陽(yáng)能光伏板最佳傾角設(shè)計(jì)方法研究，工程設(shè)計(jì)學(xué)報(bào)，第16卷；348頁(yè)-353頁(yè)，2009年10月2 蔣興波，呂自慶，劉成城，二維矩形條帶裝箱問(wèn)題的底部左齊擇優(yōu)匹配算法，軟件學(xué)報(bào)，第20卷；1528頁(yè)-1538頁(yè)，2009年6月3 姜啟元，謝金星，葉俊，數(shù)學(xué)模型，北京市：高等教育出版社，2003年8月第3版九、附件清單附錄1、計(jì)算發(fā)電量#include #include #include void main()int m,n,y35,i; float a36524,x35,z35,sum35=0.0,0.0;FILE *fp;fp=fopen(e:4.txt,r);if(fp=NULL)printf(can not open this file n);exit(0);for(m=0;m365;m+)for(n=0;n24;n+)fscanf(fp,%f,&amn);fclose(fp); fp=fopen(e:5.txt,r);if(fp=NULL)printf(can not open this file n);exit(0); for(m=0;m=0&m=35&m=70&m105)fscanf(fp,%f,&xm-70);fclose(fp);for(i=0;i35;i+)for(m=0;m365;m+)for(n=0;n24;n+)if(yi=1)if(amn=30) if(amn150|amn=80) sumi=sumi+amn*zi*xi; for (i=0;i35;i+)printf(%fn,sumi);fp=fopen(e:6.txt,a+);for(i=0;i35;i+)fprintf(fp,%.0fn,sumi);fclose(fp);附錄2、計(jì)算最佳傾斜角#include #include #include #define pi 3,1415void main()int m,n;float i,j,a36524,l=0.7,sum1=0.0,sum2=0.0,sum3=0.0,sum=0.0,w,q,x;FILE *fp;fp=fopen(g:1.txt,r);if(fp=NULL)printf(can not open this file n);exit(0);for(m=0;m365;m+)for(n=0;n24;n+)fscanf(fp,%f,&amn);fclose(fp);for(m=0;m365;m+)for(n=0;n24;n+)w=15*(n-11)*2*pi;w=w/360;q=2*pi*(284+m+1);q=q/365;q=23.45*2*pi*sin(q);q=q/360;x=sin(l)*sin(q)+cos(l)*cos(q)*cos(w); sum1=sum1+amn*cos(w)*x;sum2=sum2+amn*cos(w)*sqrt(1-x*x);for(i=-1;isum)sum=sum3;j=i;printf(%f %f %f %f,sum,j,sum1,sum2);附錄3、計(jì)算架空情況下的屋頂發(fā)電量#include #include #include #define pi 3.1415void main()int m,n;float a36524,b36524,s=0.1639,z=1.938,sum1=0.0,sum2=0.0,sum3=0.0,w,q,x,l=0.7;FILE *fp;fp=fopen(e:7.txt,r);if(fp=NULL)printf(can not open this file n);exit(0);for(m=0;m365;m+)for(n=0;n24;n+)fscanf(fp,%f,&amn);fclose(fp); fp=fopen(e:8.txt,r);if(fp=NULL)printf(can not open this file n);exit(0);for(m=0;m365;m+)for(n=0;n24;n+)fscanf(fp,%f,&bmn);fclose(fp); for(m=0;m365;m+)for(n=0;n24;n+) w=15*(n-11)*