1、2012高教社杯全國大學生數(shù)學建模競賽承諾書我們仔細閱讀了中國大學生數(shù)學建模競賽的競賽規(guī)則.我們完全明白，在競賽開始后參賽隊員不能以任何方式（包括電話、電子郵件、網(wǎng)上咨詢等）與隊外的任何人（包括指導教師）研究、討論與賽題有關的問題。我們知道，抄襲別人的成果是違反競賽規(guī)則的, 如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻的表述方式在正文引用處和參考文獻中明確列出。我們鄭重承諾，嚴格遵守競賽規(guī)則，以保證競賽的公正、公平性。如有違反競賽規(guī)則的行為，我們將受到嚴肅處理。我們授權(quán)全國大學生數(shù)學建模競賽組委會，可將我們的論文以任何形式進行公開展示（包括進行網(wǎng)上公示，在書

2、籍、期刊和其他媒體進行正式或非正式發(fā)表等）。我們參賽選擇的題號是（從A/B/C/D中選擇一項填寫）：B我們的參賽報名號為（如果賽區(qū)設置報名號的話）：S13006所屬學校（請?zhí)顚懲暾娜汉幽峡萍紝W院參賽隊員(打印并簽名) ：1.姜慧芳 2. 許朝陽 3.董文濤指導教師或指導教師組負責人 (打印并簽名)：數(shù)學建模指導小組 日期：2012年9月10日賽區(qū)評閱編號（由賽區(qū)組委會評閱前進行編號）：2012高教社杯全國大學生數(shù)學建模競賽編 號 專 用 頁賽區(qū)評閱編號（由賽區(qū)組委會評閱前進行編號）：賽區(qū)評閱記錄（可供賽區(qū)評閱時使用）：評閱人評分備注全國統(tǒng)一編號（由賽區(qū)組委會送交全國前編號）：全國評閱編

3、號（由全國組委會評閱前進行編號）：太陽能小屋的設計摘要太陽能是一種清潔的能源，它的應用正在世界范圍內(nèi)快速地增長。利用太陽光發(fā)電就是一種使用太陽能的方式，可是目前建一個太陽能發(fā)電系統(tǒng)的成本還是較高的。本題主要是在設計太陽能小屋的電池組件的分布時，在堅持節(jié)省成本的基礎上，為了更加充分有效地利用太陽能，建立優(yōu)化模型，以實現(xiàn)太陽能發(fā)電量較大。我們的想法如下：對于問題一，我們在合理假設的基礎上，建立了雙目標的整數(shù)線性規(guī)劃模型。具體是以年發(fā)電量最大和單位發(fā)電量費用最小為目標，再用面積和電池的連接加以約束。但是此模型計算繁瑣，故為了方便計算而簡化該模型，然后用lingo軟件對其進行求解，得出電池組件的數(shù)量分

4、布情況。通過人為地對兩目標進行加權(quán)，以發(fā)電量大和空間利用率大為原則，分組排列出較優(yōu)的電池組件鋪設方案。用逆變器的額定功率和允許輸入電壓的范圍對電池組件連接加以約束，求出較優(yōu)的電池組件連接方案。于是我們可以求出最大年發(fā)電量和最小投資，根據(jù)條件中給出的條件容易算出，小屋光伏電池35年壽命期內(nèi)的發(fā)電總量、經(jīng)濟效益及投資回收年限。對于問題二，通過比較東、南、西、北四個方位接受的年太陽輻射總量，我們發(fā)現(xiàn)南向面受到的年太陽輻射量遠遠多于其他三面，故我們?nèi)」夥M件的方位角為正南方。我們提出年最佳傾角這個概念，即是電池組件架空時與水平面的夾角，而太陽高度角是太陽相對于地平線的高度角，我們可以通過求解太陽高度角

5、來確定年最佳傾角。然后對問題一中的整數(shù)優(yōu)化模型進行改進，較之問題一，問題二中光伏電池的受輻射面積改變了，只把問題一的數(shù)學模型中光伏電池的受輻射面積改成架空后的受輻射面積，便能夠用類似的方法求接出架空安裝時的最大年發(fā)電量和最小投資，以及小屋光伏電池35年壽命期內(nèi)的發(fā)電總量、經(jīng)濟效益及投資回收年限。對于問題三，題中已給出要設計的太陽能小屋的建筑要求，問題二中也給出年最佳傾角和較好的方位角，以發(fā)電量最大為目標，使小屋的屋頂?shù)膬A角等于年最佳傾角，開門窗在年輻射強度較小的方位上。用和問題一同樣的方法對小屋外表面優(yōu)化鋪設光伏電池，給出鋪設方案及連接方案，易得小屋光伏電池35年壽命期內(nèi)的發(fā)電總量、經(jīng)濟效益及

6、投資回收年限。在模型的推廣問題上，問題二中在確定年最佳傾角是，可以更加詳細地加上地勢環(huán)境、季節(jié)等影響因素，這個可作為深入研究的一個方向。關鍵詞： 最優(yōu)分布 整數(shù)線性規(guī)劃模型年最佳傾角 lingo算法一、問題的重述在太陽能小屋的設計中，研究光伏電池在小屋外表面的優(yōu)化鋪設是很重要的問題。光伏電池每峰瓦的實際發(fā)電效率受到很多因素的影響，如太陽輻射強度、環(huán)境、建筑物所處的地理緯度、安裝部位及方式（貼附或架空）等。參考附件提供的數(shù)據(jù)，我們主要解決下列三個問題，并分別給出小屋外表面光伏電池的鋪設方案，使小屋的全年太陽能光伏發(fā)電總量盡可能大，而單位發(fā)電量的費用盡可能小，并計算出小屋光伏電池35年壽命期內(nèi)的發(fā)

7、電總量、經(jīng)濟效益（當前民用電價按0.5元/kWh計算）及投資的回收年限。根據(jù)要求：在同一表面采用兩種或兩種以上類型的光伏電池組件時，同一型號的電池板可串聯(lián)，而不同型號的電池板不可串聯(lián)；在不同表面上，即使是相同型號的電池也不能進行串、并聯(lián)連接。我們要對以下三個問題進行考慮。問題1：根據(jù)山西省大同市的氣象數(shù)據(jù)，僅考慮貼附安裝方式，選定光伏電池組件，對小屋的部分外表面進行鋪設，并根據(jù)電池組件分組數(shù)量和容量，選配相應的逆變器的容量和數(shù)量。問題2：電池板的朝向與傾角均會影響到光伏電池的工作效率，我們選擇架空方式安裝光伏電池，重新考慮問題1。問題3：根據(jù)附件給出的小屋建筑要求，我們要為大同市重新設計一個小

8、屋，要求畫出小屋的外形圖，并對所設計小屋的外表面優(yōu)化鋪設光伏電池，給出鋪設及分組連接方式，選配逆變器，計算相應結(jié)果。二、合理的假設1 假設除門窗外小屋其他地方都可以鋪設光伏電池；2 假設額定功率是在輻射度為1000W/,光譜為AM1.5，溫度為25標準下測試結(jié)果；3 在實際發(fā)電量計算時可以不考慮AM值的影響，根據(jù)電池表面接收到的太陽總輻射量參數(shù)進行計算即可；4 假設小屋四面的門窗都是規(guī)則的；5 假設太陽高度角不受地球自轉(zhuǎn)的影響；6 假設光伏電池壽命約為35年。三、符號的說明m第m部分的編號；n第n組連接方式編號；i第i種電池的編號；j第j種逆變器的編號； 第n種連接方式中連接到逆變器j上的第i

9、種電池的數(shù)目； 第i種電池的轉(zhuǎn)換效率； 第j種逆變器的逆變效率（直流變交流時） 第i種電池的單位面積；第n組連接方式中第j種逆變器的數(shù)目;第i種電池的單價； 第j種逆變器的單價；小屋第部分所對應的單位面積的太陽輻射量；小屋第m部分的表面積；第j種逆變器允許輸入電壓的下限；第j種逆變器允許輸入電壓的上限；小屋第m部分上第i種電池數(shù)目；陽高度角；為時角；赤緯角，也稱為太陽赤緯，即太陽直射緯度；為當?shù)氐木暥?大同的緯度為);年最佳傾角；太陽時；Q光伏電池35年內(nèi)總發(fā)電量；z1光伏電池的年發(fā)電量。四、模型的建立與求解1問題一1.1求解電池組件數(shù)量分布建立整數(shù)線性規(guī)劃模型表1：各方位年輻射量（w/平方米

10、）方位東南西北年輻射量594213.541050165.66881238.3011261478.82由表1可知，北立面的年太陽輻射量較低，從投資收益角度上來看對北立面鋪設電池組件意義不大，所以我們舍棄對北立面進行電池組件鋪設的計劃。我們按附錄中給出的順序?qū)﹄姵亍⒛孀兤鬟M行編號，對給出的電池及逆變器的三組連接方式也進行編號，將小屋的外表面分割，并對每部分進行編號。首先按照題目中要求的年發(fā)電量最大和單位發(fā)電量費用最小為目標，即：建立約束條件：1、所有電池的面積之和不應小于太陽能小屋的有效面積，即：2、并聯(lián)的光伏組件端電壓相差不應超過10%，所以我們要求不同類型電池組件盡量不并聯(lián)，且兩并聯(lián)分支的所串

11、聯(lián)電池組件的數(shù)目應大致相等；則在三種連接方式里，應滿足連接好的電池組件的電壓應在所選擇的逆變器允許輸入電壓范圍內(nèi)，即： (j=1,2,318)13、在三種連接方式里，應滿足連接好的電池組件的功率和應不超過所選擇的逆變器的額定功率，即： （n=1,2,3; j=1,2,3,18）4、必須為整數(shù)；模型的簡化由于資源的有限，時間的緊迫，為了方便計算，對上述雙目標整數(shù)規(guī)劃模型進行簡化如下：目標：小屋的每部分外表面上的電池組件的年發(fā)電量（直流）最大；約束條件：1、小屋的該部分外表面上的電池組件的總面積不得大于該部分外表面的面積；2、為整數(shù)；模型的編程和求解用lingo軟件編程，求出小屋各個部分（除北立面

12、外）上第i種電池數(shù)目（i=1,2，24）。求解結(jié)果如下：表2：小屋各部分的第i中電池數(shù)目電池型號南向面東向面西向面前傾屋頂后傾屋頂A100000A211010A3201717185A412111A500010A600000B100000B200000B300000B400000B500000B600000B700000C100000C200000C300000C400000C500000C601017820C710100C800000C900000C1000000C1100000因為，該模型很簡略，很多地方考慮并不周全，故求解的結(jié)果不能嚴格執(zhí)行，只能作為參考，故問題二和問題三種不再考慮數(shù)量要求

13、。21.2求解電池組件的鋪設方案電池組件的鋪設模型考慮已建成模型的目標與約束要求，對小屋的每塊表面積進行鋪設電池組件，對于模型一中的兩個目標，我們考慮給它們加權(quán)，且我們決定發(fā)電量大這個目標權(quán)重較大；具體原則是光電轉(zhuǎn)換率大且面積較小的優(yōu)先。（圖1）通過對比分析圖1，我們發(fā)現(xiàn)A3較為理想，A4、B1次之；并且在鋪設電池組件過程中先使用A類、B類電池組件鋪設大片面積，之后再用C類電池組件填充小面積。電池組件的鋪設方案下面是小屋各個表面電池組件鋪設情況：前傾屋頂（南向）電池組件優(yōu)化鋪設圖如下：（圖2）后傾屋頂（北向）電池組件優(yōu)化鋪設圖如下：（圖3）南立面電池組件優(yōu)化鋪設圖如下：（圖4）東立面電池組件優(yōu)

14、化鋪設圖如下：（圖5）西立面電池組件優(yōu)化鋪設圖如下：（圖6）1.3求解電池組件的連接方案電池組件的連接模型下面要實現(xiàn)電池組件的串并聯(lián)和搭配逆變器的選擇；因為并聯(lián)電壓相差不能超過10%，所以我們盡量使同種類型電池組件間進行并聯(lián)，并使并聯(lián)的幾支中串聯(lián)的同種型號電池組件個數(shù)相同；在選擇逆變器上，我們要求逆變器利用率最大，即在逆變器自身客觀條件允許的情況下連接盡可能多的電池組件；另外在前面的那些前提下還要兼顧逆變器的價格，以使投資成本盡量小。3電池組件的連接方案綜上所述看，我們可以找到較好的一種組件連接方式如下：前傾屋頂（南向）電池組件優(yōu)化連接：1.型號為SN17的逆變器連接一個有三個分支的并聯(lián)，每個

15、分支均由17個型號為A3的電池組件串聯(lián)而成；2.型號為SN11的逆變器連接一個串聯(lián)，該串聯(lián)由8個型號為C6的電池組件串聯(lián)而成。后傾屋頂（北向）電池組件優(yōu)化連接：1.型號為SN17的逆變器連接一個有兩個分支的并聯(lián)，每個分支均由18個型號為C6的電池組件組成；2.型號為SN11的逆變器連接一個串聯(lián)，該串聯(lián)由6個型號為A3的電池組件串聯(lián)而成。南立面電池組件優(yōu)化連接：1.型號為SN17的逆變器連接一個串聯(lián)，該串聯(lián)由7個型號為A3的電池組件串聯(lián)而成；2.型號為SN11的逆變器連接一個串聯(lián)，該串聯(lián)由3個型號為C2的電池組件串聯(lián)而成；3.型號為SN17的逆變器連接一個串聯(lián)，該串聯(lián)由23個型號為C6的電池組件

16、串聯(lián)而成；4.型號為SN3的逆變器連接一個有三個分支的并聯(lián)，每個分支均由4個型號為C7的電池組件串聯(lián)而成。東立面電池組件優(yōu)化連接：1.型號為SN17的逆變器連接一個串聯(lián)，該串聯(lián)由15個型號為A3的電池組件串聯(lián)而成；2.型號為SN17的逆變器連接一個串聯(lián)，該串聯(lián)由11個型號為C6的電池組件串聯(lián)而成。西立面電池組件優(yōu)化連接：1.型號為SN17的逆變器連接一個串聯(lián)，該串聯(lián)由13個型號為A3的電池組件串聯(lián)而成；2.型號為SN17的逆變器連接一個串聯(lián)，該串聯(lián)由17個型號為C6的電池組件串聯(lián)而成。1.4計算35年壽命期內(nèi)的發(fā)電總量、經(jīng)濟效益及投資的回收年限數(shù)學公式模型在此已給定的優(yōu)化圖的基礎上，可計算出小

17、屋光伏電池35年壽命期內(nèi)的發(fā)電總量、經(jīng)濟效益（當前民用電價按0.5元/kWh計算）及投資的回收年限。根據(jù)以下公式：（1）35年壽命期內(nèi)的發(fā)電總量=年發(fā)電總量*（10+15*90%+10*80%）即Q=z1*（10+15*90%+10*80%）（2）年發(fā)電總量=小屋每塊表面積上每種電池的總面積*該面積對應的輻射強度*光電轉(zhuǎn)換率*交流電轉(zhuǎn)化率即（3）35年內(nèi)經(jīng)濟效益=35年壽命期內(nèi)的發(fā)電總量*交流電單價（0.15元/kwh）-投資費用Z=0.15*Q-z2（4）投資費用=電池組件購買費用+逆變器購買費用（5）令k年內(nèi)經(jīng)濟效益為0，即投資費用=年發(fā)電總量*交流電單價*k，若求出的k<=10，則

18、k即為所求投資的回收年限；若25>k>10,則再令k年內(nèi)經(jīng)濟效益為0，即投資費用=年發(fā)電總量*交流電單價*10+(k-10)*90%若求出k滿足25>k>10,則滿足題意，否則再令投資費用=年發(fā)電總量*交流電單價*10+15*90%+(k-10)*80%若求得的35>k>25,則滿足題意，否則就算投資失敗求解代入數(shù)據(jù)，經(jīng)計算得：35年壽命期內(nèi)的發(fā)電總量=(kw.h)35年內(nèi)經(jīng)濟效益=(元)k<1,即在一年內(nèi)即可收回投資成本。42問題二此題與第一題的差別在于要找一個最好的傾角并使電池組件朝向一個好的方位，我們用年太陽輻射強度來判斷這個好的傾角較好的方位。

19、由于小屋墻體上不易架空鋪設電池組件，后傾屋頂面積較小，為了簡化模型，我們只考慮前傾屋頂鋪設優(yōu)化問題，即四壁仍是貼附鋪設。2.1最佳方位角的選取分析附件中的資料，我們整理出東、南、西、北四個方位的年輻射強度如下：各方位年輻射強度（w/平方米）方位東南西北年輻射強度594213.541050165.66881238.3011261478.82通過對比分析，容易得到南向面的年輻射強度值較大，故我們選正南方為最佳方位角。2.2年最佳傾角的選取5傾斜角是太陽電池方陣平面與水平地面的夾角，并希望此夾角是方陣一年中發(fā)電量為最大時的最佳傾斜角度，稱之為年最佳傾角。 一年中的最佳傾斜角與當?shù)氐牡乩砭暥扔嘘P，當緯

20、度較高時，相應的傾斜角也大；但是也要考慮到屋頂?shù)膬A斜角等方面的限制條件。我們?nèi)∫荒曛刑栞椛鋸姸茸畲蟮哪且惶旒澳翘炖锾栞椛鋸姸茸畲蟮哪且恍r來計算這個年最佳傾角；另外，我們知道太陽高度角是太陽相對于地平線的高度角，這是以太陽視盤面的幾何中心和理想地平線所夾的角度，且年最佳傾角和太陽高度角是互余的，所以可以通過求太陽高度角來求得這個年最佳傾角。（1）太陽高度角可以使用下面的算式，經(jīng)由計算得到很好的近似值：其中為太陽高度角，為時角，為當時的太陽赤緯，為當?shù)氐木暥?大同的緯度為)。（2）赤緯角也稱為太陽赤緯，即太陽直射緯度，其計算公式近似為其中為日期序號，例如，1月1日為，3月22日為。（3）時角

21、是以正午12點為0度開始算，每一小時為15度，上午為負下午為正，即10點和14點分別為-30度和30度。因此，時角的計算公式為其中為太陽時（單位：小時）。我們?nèi)∈率蝗帐藭r為一年中太陽輻射強度最大的時刻；代入公式，分析計算可得太陽高度角=。綜上，可求出年最佳傾角。；2.3前傾屋頂電池組件的鋪設可用與問題一同樣的方法對問題二中的模型進行闡述、求解。即據(jù)已建成模型的思路，對小屋前傾屋頂進行鋪設電池組件，對于模型中的兩個目標，我們考率給它們加權(quán)，且發(fā)電量大這個目標權(quán)重較大；具體原則是光電轉(zhuǎn)換率大且面積較小的優(yōu)先，并且在鋪設電池組件過程中先使用A類、B類電池組件鋪設大片面積，之后再用C類電池組件填

22、充小面積。下面是小屋前傾屋頂電池組件鋪設情況：前傾屋頂電池組件鋪設優(yōu)化圖（俯視圖）如下：（圖7）前傾屋頂電池組件鋪設優(yōu)化圖（側(cè)視圖）如下：（圖8）2.4前傾屋頂電池組件的連接下面實現(xiàn)電池組件的串并聯(lián)和搭配逆變器的選擇，可用同樣的方法對問題二中的模型進行闡述、求解。綜上所述看，我們可以找到較好的一種組件連接方式如下：前傾屋頂（南向）電池組件優(yōu)化連接：（1）型號為SN17的逆變器連接一個有三個分支的并聯(lián)，每個分支均由17個型號為A3的電池組件串聯(lián)而成；（2）型號為SN11的逆變器連接一個串聯(lián)，該串聯(lián)由8個型號為C6的電池組件串聯(lián)而成。2.5計算35年壽命期內(nèi)的發(fā)電總量、經(jīng)濟效益及投資的回收年限問題

23、一種已給出求小屋光伏電池35年壽命期內(nèi)的發(fā)電總量、經(jīng)濟效益（當前民用電價按0.5元/kWh計算）及投資的回收年限的各個計算公式，代入數(shù)據(jù)計算得：35年壽命期內(nèi)的發(fā)電總量= (kw.h)35年內(nèi)經(jīng)濟效益=(元)k<1,即在一年內(nèi)即可收回投資成本。3問題三3.1太陽能小屋的設計這是一個優(yōu)化問題，目標是獲得最大太陽輻射強度，為達到這一目標，我們設計的屋頂以問題二中求出的年最佳傾角作為屋頂?shù)膬A斜角，并朝向太陽輻射強度最大的方位正南方；在滿足太陽能小屋建筑要求的基礎上，盡量使南方少開或不開門窗，盡量使門窗分布在北面。6再綜合小屋各個建筑要求，我們作出如下設計：小屋的設計圖如下：（圖9）南立面?zhèn)纫晥D

24、如下：（圖10）北立面?zhèn)纫晥D如下：（圖11）東立面?zhèn)纫晥D如下： 西立面測試圖如下：（圖12）3.2電池組件鋪設方案用和問題一類似的方法，對電池組件進行鋪設，同樣為了簡化模型，減少計算量，合理地舍棄對北立面的鋪設計劃。下面是小屋各個表面電池組件鋪設情況：南立面電池組件鋪設優(yōu)化圖如下：（圖13）西立面電池組件鋪設優(yōu)化圖如下：東立面電池組件鋪設優(yōu)化圖如下： （圖14）屋頂電池組件鋪設優(yōu)化圖如下：（圖15）3.3電池組件連接方案用和問題一中同樣的方法，我們可以找到的新設計的小屋的一種較好的組件連接方式，如下：屋頂電池組件較優(yōu)連接方式：1.型號為SN17的逆變器連接一個有三個分支的并聯(lián)，每個分支均由17

25、個型號為A3的電池組件串聯(lián)而成；2.型號為SN17的逆變器連接一個串聯(lián)，該串聯(lián)由15個型號為C6的電池組件串聯(lián)而成。南立面電池組件優(yōu)化連接：1.型號為SN17的逆變器連接一個有兩個分支的并聯(lián)，每個分支均由9個型號為A3的電池組件串聯(lián)而成；2.型號為SN17的逆變器連接一個串聯(lián)，該串聯(lián)由11個型號為C2的電池組件串聯(lián)而成。3.型號為SN11的逆變器連接一個串聯(lián)，該串聯(lián)由10個型號為C6的電池組件串聯(lián)而成。東立面電池組件優(yōu)化連接：1. 型號為SN11的逆變器連接一個串聯(lián)，該串聯(lián)由5個型號為A3的電池組件串聯(lián)而成。西立面電池組件優(yōu)化連接：1. 型號為SN11的逆變器連接一個串聯(lián)，該串聯(lián)由5個型號為A

26、3的電池組件串聯(lián)而成。3.5計算光伏電池35年壽命期內(nèi)的發(fā)電總量、經(jīng)濟效益及投資的回收年限問題一種已給出求小屋光伏電池35年壽命期內(nèi)的發(fā)電總量、經(jīng)濟效益（當前民用電價按0.5元/kWh計算）及投資的回收年限的各個計算公式，代入數(shù)據(jù)計算得：35年壽命期內(nèi)的發(fā)電總量=(kw.h)35年內(nèi)經(jīng)濟效益=(元)k<3,即在3年內(nèi)即可收回投資成本。五、模型的評價問題一中建立的雙目標規(guī)劃模型，結(jié)構(gòu)、算法比較嚴謹，可以得到很好地計算結(jié)果；不足之處是計算工作過于繁瑣復雜；光伏電池鋪設和連接方法簡單實用且樣式靈活。問題二中認真分析了并確定了最佳方位和年最佳傾角，繁簡得當，充分體現(xiàn)了建模思想靈活性。問題三種結(jié)合

27、給定條件，成功設計出簡單實用的太陽能小屋，方法簡單，得出的結(jié)果卻很好。在確定年最佳傾角上，可以更加詳細地加上地勢環(huán)境、季節(jié)等影響因素，這個可作為深入研究的一個方向。六、參考文獻1 胡運權(quán)，運籌學基礎及應用（第五版），北京：高等教育出版社，2008。2王冬琳，數(shù)學建模及實驗，北京：國防工業(yè)出版社，2004。3劉鋒,數(shù)學建模,南京：南京大學出版社，2005。4 范玉妹，徐爾,數(shù)學規(guī)劃及其應用(第二版), 北京：冶金工業(yè)出版社，2003。5 H. P. Williams， 數(shù)學規(guī)劃模型建立與計算機應用,北京：國防工業(yè)出版社, 1991。6姜啟源，謝金星，葉俊，數(shù)學模型（第三版），北京：高等教育出版社

28、，2005。七、附錄1.使用Lingo編寫的程序：model:max=1050165.66*x11*1.27664*0.1684+1050165.66*x12*1.938396*0.1664+1050165.66*x13*1.27664*0.187+1050165.66*x14*1.637792*0.165+1050165.66*x15*1.63515*0.1498+1050165.66*x16*1.938396*0.1511+1050165.66*x17*1.63515*0.1621+1050165.66*x18*1.938396*0.1639+1050165.66*x19*1.470144

29、*0.1598+1050165.66*x20*1.62688*0.148+1050165.66*x21*1.940352*0.1598+1050165.66*x22*1.940352*0.152+1050165.66*x23*1.668*0.1499+1050165.66*x24*1.43*0.0699+1050165.66*x25*0.939231*0.0617+1050165.66*x26*1.575196*0.0635+1050165.66*x27*1.54*0.0584+1050165.66*x28*1.54*0.0649+1050165.66*x29*0.11005*0.0363+1

30、050165.66*x30*0.1107*0.0363+1050165.66*x31*0.218325*0.0366+1050165.66*x32*0.3266*0.0366+1050165.66*x33*0.29039*0.0413+1050165.66*x34*1.17124*0.0427;x11*1.27664+x12*1.938396+x13*1.27664+x14*1.637792+x15*1.63515+x16*1.938396+x17*1.63515+x18*1.938396+x19*1.470144+x20*1.62688+x21*1.940352+x22*1.940352+x23*1.668+x24*1.43+x25*0.939231+x26*1.575196+x27*1.54+x28*1.54+x29*0.11005+x30*0.1107+x