1、摘要本文對易拉罐的最優設計主要從節省用料的角度研究。建立模型時，在滿足體積以及其它設計要求的限制下，以用料最少為目標建立最優化模型。求解模型的主要方法包括：乘子法、條件極值法以及數學軟件(、)求解等。在對易拉罐形狀及尺寸設計進行研究時，需要選擇適當的工具，運用多次測量求平均值的方法確定出必要的數據。實際中易拉罐的設計考慮到多方面的因素，如美觀、實用、生產、運輸以及其自身各部分的抗壓能力等。本文主要在某些設計要求下，研究用料最省的形狀設計，將求解出的最優形狀和尺寸與相應實測數據作對比來衡量設計的合理性。當假設易拉罐的形狀為正圓柱體時，主要以圓柱體高度與半徑的比例關系確定易拉罐形狀是否符合用料最省

2、的最優設計。分別運用條件極值法、乘子法以及數學軟件求解的方法最終確定出高度與半徑的比值為4，與實際易拉罐高度與半徑的比值基本符合，能夠說明實際的易拉罐形狀設計符合用料最省的設計原則。當假設易拉罐由正圓柱體和圓臺兩部分組成時，分別假設易拉罐尺寸符合不同設計要求，運用逐步改進的方法，最終求得易拉罐各項尺寸與實測數據比較吻合，說明現實中易拉罐的設計滿足用料最省的要求。求解時主要運用軟件和乘子法并在軟件的輔助下分別求得各項尺寸。最后根據對易拉罐的觀察和想象，在保證易拉罐整體形狀變化不大的前提下，我們提出將易拉罐上端的圓臺改為球臺。同樣以用料最省作為其最優設計，求出新設計的易拉罐與現實中易拉罐相比大約能

3、節省0.49%的用料。在建立易拉罐用料最省模型的基礎上，對模型作出進一步推廣。聯系實際當中各種形狀比較規則的容器，對其最優設計的一般規律作出說明。 關鍵詞：乘子法 重積分 條件極值法 1 問題重述我們只要稍加留意就會發現銷量很大的飲料 (例如飲料量為355毫升的可口可樂、青島啤酒等) 的易拉罐(即易拉罐)的形狀和尺寸幾乎都是一樣的?？磥?，這并非偶然，這應該是某種意義下的最優設計。當然，對于單個的易拉罐來說，這種最優設計可以節省的錢可能是很有限的，但是如果是生產幾億，甚至幾十億個易拉罐的話，可以節約的錢就很可觀了。現在就請你們小組來研究易拉罐的形狀和尺寸的最優設計問題。具體說，請你們完成以下的任

4、務：1 取一個飲料量為355毫升的易拉罐，例如355毫升的可口可樂易拉罐，測量你們認為驗證模型所需要的數據，例如易拉罐各部分的直徑、高度，厚度等，并把數據列表加以說明；如果數據不是你們自己測量得到的，那么你們必須注明出處。2 設易拉罐是一個正圓柱體。什么是它的最優設計？其結果是否可以合理地說明你們所測量的易拉罐的形狀和尺寸，例如說，半徑和高之比，等等。3 設易拉罐是一回轉體，上面部分是一個正圓臺，下面部分是一個正圓柱體。什么是它的最優設計？其結果是否可以合理地說明你們所測量的易拉罐的形狀和尺寸。4 利用你們對所測量的易拉罐的洞察和想象力，做出你們自己的關于易拉罐形狀和尺寸的最優設計。5 用你們

5、做本題以及以前學習和實踐數學建模的親身體驗，寫一篇短文(不超過1000字，你們的論文中必須包括這篇短文)，闡述什么是數學建模、它的關鍵步驟，以及難點。2 問題分析在對易拉罐的形狀進行研究時，首先需要測量出實際中易拉罐形狀的各個尺寸，然后以易拉罐用料最少作為其最優設計，對題目給出的不同簡化形狀進行研究。對于問題一，取一個355毫升易拉罐（可口可樂），首先分析出驗證模型可能需要的數據，然后利用相應的工具對需要的數據進行多次測量取平均值確定易拉罐各項尺寸的大小。對于問題二，題目假設易拉罐的形狀為一正圓柱體，但并沒有對各部分的壁厚做出說明，在求解的過程中可分易拉罐各面厚度相同和不同的情況進行求解。最終

6、確定出高度與半徑的比值關系，并與實際測量數據進行分析比較，判斷易拉罐設計的合理性。對于問題三，結合第二問得出的結論并在易拉罐形狀為一正圓臺和圓柱組成的情況下，求解過程可仍以材料最省為最優設計，但同時要滿足上、下頂面的強度要求，還要滿足加工方面的要求，并建立一個廣泛的最優化模型。然后假設多種情況對模型進行逐步改進，最終確定一個滿足材料最省，同時又滿足其他方面（如強度、美觀、加工等）要求的易拉罐形狀和尺寸，與實際測量值進行比較，分析其設計的合理性。對于問題四，在易拉罐頂蓋厚和壁厚比例關系不變的前提下，以用料最省作為最優化設計，根據多面體中球體表面積與體積比值最小的基本原理，將易拉罐上部的圓臺設計為

7、球臺。將設計出的易拉罐用料與實際的相比較，并對其優點進行分析說明。對于問題五，結合前面學習和實踐數學建模的親身體驗，分別對數學建模的基本認識、建模過程中的關鍵步驟和難點進行討論分析。3 模型假設1) 所取易拉罐各面的厚度均勻2) 易拉罐的頂蓋和下底面都是規則的平面3) 易拉罐都是規則的多面體4) 不考慮溫度對測量儀器的影響5) 易拉罐用同一種材料制成4 符號說明符號含義單位易拉罐的表面積所用材料的總體積為罐體圓柱體部分圓的半徑為圓柱體的高度易拉罐圓柱部分的壁厚易拉罐的罐內體積為表示圓臺面的傾斜角度5 模型的建立與求解5.1 問題一5.1.1解題思路本題要求選取一個355毫升的易拉罐，并對其進行

8、數據測量，這些數據包括易拉罐的直徑、高度，厚度等。測量時，取一個可口可樂公司生產的易拉罐，然后利用相應的測量工具對其測量即可。5.1.2求解準備需要數據的確定 經過分析發現，模型中可能用到的數據種類有罐直徑、罐高、罐壁厚、頂蓋厚、圓臺高、頂蓋直徑、圓柱體高、罐底厚、罐內體積等。具體說明如下：圖1(1) 罐直徑：易拉罐圓柱體部分(罐體最胖部分)橫截面圓的直徑(2)罐高：從易拉罐的頂蓋到底面的高度(3)罐壁厚：圓柱體回轉面部分的瓶壁厚度(4)頂蓋厚：易拉罐的上頂蓋厚度(5)罐內體積：易拉罐內部的體積(6)罐底厚：易拉罐的下底面厚度(7)頂蓋直徑：圓錐臺的上表面圓的直徑即罐頂蓋的直徑(8)圓臺高：從

9、罐的頂蓋到圓柱體部分的高度(9)圓柱體高：圓柱體部分的高度各數據的測量方法(1)直接測量經過分析可得，罐桶直徑、罐高、圓臺高、頂蓋直徑、圓柱直徑這幾種數據類型屬于外部屬性，可以直接進行測量。測量時可選用以下兩種方法：用一條非常窄的薄紙條,環繞易拉罐相關部位一圈測得周長,然后再換算求得直徑、半徑、面積等。用游標卡尺（50分度）對相關部位進行直接測量，計算出直徑和高等。(2)間接測量由現實情況可知，易拉罐的罐壁、頂蓋和罐底有一面是在易拉罐的內部，不能直接進行測量，因此就需要對他們進行處理后再進行測量。對于厚度的測量都可用下面第一種方法，體積的測量可用第二種方法。首先用剪刀和鉗子對易拉罐進行刨切，由

10、于易拉罐厚度和頂蓋厚度較小，可利用螺旋測微器進行測量。取一個量筒()和空的易拉罐，首先將清水倒入易拉罐中直至與罐口相平；然后將易拉罐中的水倒入量筒中進行讀數，即得到了易拉罐的體積。5.1.3求解為確保數據的精確性，需要對所有數據進行多次測量求平均值，經多次測量求得所需的數據如下表： 表1易拉罐(可口可樂)各項尺寸列表數據種類實測數據平均值單位罐高686罐桶直徑2086.5861罐壁厚1253頂蓋厚16罐底厚0圓臺高821頂蓋直徑228圓柱體高286罐內體積364.9364.0364.8易拉罐體積的說明由現實情況可知，每瓶飲料都未完全裝滿，另外還有一部分空間余量，根據測量數據也可以得出，對于標注

11、為的可口可樂易拉罐，它的實際罐體容量為，所以在以下的各個問題中凡是對于易拉罐的體積在計算時都以為標準進行設計是與現實不符合的，實際所測易拉罐的容積為，所以在以下題目的求解過程中對于易拉罐的容積都以為標準進行計算。5.2問題二5.2.1模型分析對于本題來說，假設易拉罐是一個正圓柱體(見下圖2)，怎樣構思和設計才能達到用料最省的最優設計是關鍵所在。根據推理可知，易拉罐的形狀由罐直徑和罐高的比例關系決定的。罐直徑和罐高又是罐的主要尺寸，所以在對易拉罐進行設計時就主要以罐直徑和罐高的比例來作為檢驗設計是否合理的一個標準。前面我們已經對正在使用的易拉罐的相關數據進行了測量，因此我們可以將這些數據與自己認

12、同的最優設計得出的數據相比較并說明其合理性。由于現在是工業高度發達的現代社會，每個國家的原用料都是高度稀缺。因此省料省錢并能達到基本要求的制作方法是資本家追逐的最重要目標之一，下面我們以制作易拉罐用料最省作為它的最優設計目標來進行研究，并通過罐的直徑與罐高的比例關系來說明其合理性。5.2.2模型準備設計方案的確定以所用料最省作為目標對易拉罐進行最優設計時，當易拉罐各部位罐壁厚度不相同時對易拉罐尺寸的設計帶來較大影響。在各部位厚度不相同的情況下，相應的研究出了兩種研究方案。第一種方案：用手捏一下發現易拉罐非常的薄，在這種感性認識的支配下我們認為這是最理想的情況即易拉罐的壁厚均勻且無限薄(可以忽略

13、不計,見圖3)。在這種情況下主要考慮易拉罐的表面積來建立數學模型求解。第二種方案：用手往下按頂蓋能夠感覺到它的硬度要比其它部位的用料要硬,相比之下,硬度體現在同樣用料的厚度上；根據測量的數據可知，頂蓋厚度大約是其他部分的用料厚度的3倍(參考)。因此可以假設除易拉罐的頂蓋外,罐的厚度相同(見圖4)。在這種情況下制作易拉罐的用料就要通過各個部位體積來考慮，為求制作需要原用料最省，通過建立數學模型求解。各種情況的易拉罐側視圖如下：圖2圖3圖4各方案模型目標與約束條件的確定n 第一種方案(1)目標：假設易拉罐是一個正圓柱體，易拉罐各處厚度均勻且非常薄(可忽略其厚度)時。就不用具體考慮易拉罐用料的體積，

14、只以易拉罐的表面積最小為目標就可使用料最省。設易拉罐罐高為，罐體圓柱體部分圓的半徑為。即(2)主要約束：易拉罐的容積是一個固定的常量。在忽略罐壁厚的情況下我們可以認為易拉罐的體積與它的容積等價。設易拉罐的罐內體積為，即n 第二種方案(1)目標：本方案考慮到易拉罐的頂蓋同罐壁硬度不同，因此同種用料下易拉罐頂蓋和罐壁的厚度也各不相同。由于引入了易拉罐的厚度，需要研究易拉罐的罐殼體積而不是單方面的考慮易拉罐的表面積。所以說以易拉罐的用料體積最小為目標，可使制造易拉罐的用料最省。通過現實情況和對圖形4的觀察可知，易拉罐的用料體積主要分成三個部分：頂蓋所用的用料體積、罐底所用的用料體積、側面所用的用料體

15、積。符號的確定：假設除易拉罐的頂蓋外,罐的厚度相同，記作；頂蓋的厚度為()；易拉罐的半徑為，直徑為；罐高為；罐內體積為；所用用料的總體積為。下面是對易拉罐三部分用料體積的確定：易拉罐側面所用的用料體積為：易拉罐頂蓋所用的用料體積為：易拉罐底所用的用料體積為：綜上可得易拉罐用料的總體積為：因為，為簡化模型求解，所以的項可以忽略。所以注：由以上表達式可知罐的用料體積近似等于各部分罐面的外表面面積與相應厚度的乘積之和。因此，對下面關于罐的用料體積進行計算時，忽略厚度對半徑及高度的影響，直接對各部分管面的外表面面積與相應厚度的乘積進行求和。(2)主要約束：同第一種方案的約束條件一樣，易拉罐內部的體積為

16、一常量。在忽略罐壁厚的情況下我們可以認為易拉罐的體積與它的容積相等。5.2.3模型建立由于第一種方案是在感性認識下產生的方案，因此我們第一種方案建立的模型定義為感性模型；第二種方案是在以具體數據作為依托的情況下建立的方案，我們將這種模型定義為理性模型。兩種模型都是以需要制作易拉罐的用料最少為目標，都屬于最優化模型。根據兩種方案各自的特點分別建立模型如下：感性模型以易拉罐罐內體積和飲料容量相同為約束條件，以制作易拉罐需要的原用料最省為目標建立最優化模型。即符號說明：符號含義單位易拉罐的表面積易拉罐的罐內體積罐體圓柱部分的半徑易拉罐的罐高理性模型以易拉罐罐內體積一定為約束條件，以制作易拉罐需要的用

17、料最省為目標建立最優化模型。即符號說明：符號含義單位所用用料的總體積罐體圓柱部分的半徑易拉罐的罐高除易拉罐的頂蓋外,罐的厚度易拉罐頂蓋的厚度() 模型求解對于以上模型都可用兩種方法求解，感性模型可以通過條件極值法和利用數學軟件直接求解，理性模型可以通過乘子法、條件極值法和利用數學軟件直接求解。為了增加結果的準確性，下面分別利用這些方法對模型進行求解。感性模型求解(1)條件極值法模型中共有兩個變量和，體積的限制為一等式，即，通過等式變換可得將上面的表達式代入到目標函數中可得：此時目標函數中只含變量，對求導可得：由可得：對求二階導數可得：，由可得：即時，取極小值，且是唯一極值點。所以，時，取最小值

18、。由，可以確定：(2)借助數學軟件求解利用數學軟件求解得出：，根據半徑與高度的大小可確定：(3)結果檢驗。因此可得易拉罐的直徑同罐高之比為1：1關系，由此發現通過感性模型計算出來的數據同我們實際測量的數據罐直徑和罐高之比1：2相差甚遠。理性模型求解(1)乘子法將目標中的主要限制條件設為：要找目標函數在條件下的極值點，可以先作函數：于是可得到如下方程組：于是將問題化為求三元函數L的無條件極值的問題。根據2式解得：根據3式解得：將上面的兩式代入1式可得：將代入中得：即高度與半徑的關系式為：(2)求條件極值法由解出 ,代入得： 上式中只有為變量，對求導得：求極值點:令其導數為零得解得極值點為,因此計

19、算的二階導數因此使達到極小值，且是唯一極值點。所以，時，取最小值。另外，求極小的方法是應用算術幾何平均值不等式，即：當且僅當時等號成立。令,于是有當且僅當時等號成立,即所以，高度與的關系式為，此時用料體積最小。(3)結果檢驗通過兩種方法求得的結果相同，說明模型的正確性。經過以上兩種方法推導得到易拉罐的罐高,罐的半徑為，因此易拉罐的半徑與罐高之比。又因測量數據易拉罐的頂蓋厚大約是罐壁厚的3倍，即。代入可得罐的半徑與罐高之比為1：4，而實際測量的數據大概也是這個比例，根據半徑與高度的比值能夠說明易拉罐的形狀符合用料省的最優設計。問題三解題思路本題的最優設計仍然是用料最少，不過在求解的過程中還要考慮

20、到，上頂蓋和底蓋的強度要求，并且在最后還要兼顧到美觀的因素。在求解的時候，先建立一個最廣泛的模型，然后根據假設的不同，分別利用此模型進行求解實現模型的逐步改進，以達到最優的設計尺寸。模型建立對于易拉罐的簡化形狀，可將其分成兩部分考慮，上部分為一正圓臺，下部分為一正圓柱體，如右圖5所示。圖5圖6易拉罐容積的確定(1)正圓臺部分體積的確定正圓臺是對應圓錐的一部分，如右圖6所示。求解圖中圓臺部分的體積時，可先求出圖中小圓錐的體積，再求出大圓錐的體積，則圓臺的體積即為大圓錐體積減去小圓錐的體積。根據圓錐的體積公式：如圖中所示，其中小圓錐的體積為：大圓錐的體積為：則圓臺的體積為：(2)正圓柱部分體積的確

21、定綜上可得：易拉罐結構的總體積為：根據測量數據，易拉罐壁厚，；因此在確定易拉罐的容積時可近似看成體積，即：易拉罐用料體積的確定易拉罐用料主要包括四部分：上頂面、圓臺回轉面、圓柱回轉面、下底面。根據問題二可知，并且由于易拉罐壁厚，；為簡化計算，在求易拉罐用料體積時，可近似看成各個面的面積與其厚度乘積之和，忽略各個面由于相交產生的體積偏差。在求各個面的面積時，以外表面的測量值為準。設圓柱回轉面的厚度為一單位，上頂面、圓臺回轉面、下底面分別是圓柱回轉面厚度的、倍。上頂面的用料體積為：圓臺回轉面的用料體積為：下底面的用料體積為：圓柱回轉面的用料體積為：綜上可得，易拉罐用料體積為：于是，可根據易拉罐容積

22、一定，用料體積最省的最優化設計建立以下模型：符號說明：表示圓臺上頂面的半徑；表示圓柱的半徑；表示圓臺面的傾斜角；表示圓柱體的壁厚；表示各面厚度與圓柱體壁厚比值；模型求解根據問題一中的數據，取易拉罐體積，在求解過程中，會發現對于各種壁厚，是不能通過求用料最省來求得的，因為對于上、下頂面是有一定的強度要求的，如果厚度達不到強度要求是容易發生危險的而強度達到了就勢必會造成用料費用的上升所以兩者是矛盾的，衡量其重要性，強度要求更重，所以在求解的時候對于各種壁厚要進行定量處理。根據問題一數據和問題二的求解過程，設上頂面、圓臺回轉面、下底面分別是圓柱回轉面厚度的，倍。該模型與問題二中的模型相比變量更多，求

23、解更加復雜。運用直接推導計算的方法是相當困難的。對于此種情況，考慮運用數學軟件求解。第一種方法可根據軟件直接求解，第二種方法可根據乘子法，將模型轉化為方程組的形式，運用軟件對方程組求解。下面分別用兩種方法對模型求解：n 方法一：軟件求解觀察本模型屬于最優化模型，可利用軟件直接求解（具體編程過程請參見附錄一），得到上頂蓋半徑圓柱體半徑圓臺傾斜角圓柱體高度所用材料0n 方法二： 乘子法和軟件求解首先構造函數分別對求偏導，并使之為零，與聯立得到如下方程組：把上述方程組，利用軟件中的函數求解得到（具體編程過程請參見附錄二）上頂蓋半徑圓柱體半徑圓臺傾斜角圓柱體高度乘子所用材料0圖7結果驗證通過觀察數據，

24、兩種方法求得的結果基本一樣，各項取其算數平均值，得到上頂蓋半徑，圓柱體半徑，圓臺傾斜角，圓柱體的高度，此時所用材料，所得形狀如圖7。根據以上兩種求解方法所求數據相同，但與易拉罐的實際形狀和尺寸相差很大。需要對模型作出進一步改進。5模型改進一上述模型中，在處理易拉罐底面厚度時是簡化為與圓柱體壁厚相同來解決的，但通過實際測量發現底面厚度與圓柱體壁厚是不相等的，并且底面厚度與頂面厚度相同都為，即此時，將修改后的值代入模型，利用上述的兩種方法，求得易拉罐形狀的各項尺寸如下表所示：上頂蓋半徑圓柱體半徑圓臺傾斜角圓柱體高度所用材料觀察數據會發現，所得形狀同圖7類似，只是尺寸略有不同，用所得數據與問題一中實

25、測數據進行比較，相差仍然比較大，所以還要對模型進行改進。模型改進二在模型改進一中，只對易拉罐的各種壁厚進行了分析研究，發現并不能滿足現實情況，所以要從其他分面考慮，通過對易拉罐的觀察發現，易拉罐頂蓋實際上并不是平面的, 略有上拱, 頂蓋實際上是半徑為3+0.4+0.2=的材料沖壓而成的, 從頂蓋到圓柱體部分的斜率為0.3,這些要求也許保證了和易拉罐的薄的部分的焊接(粘合)很牢固、耐壓. 所有這些都是物理、力學、工程或材料方面的要求，簡單通過求材料最省是得不到滿意的結果的，所以要對其進行改進，假設易拉罐的上頂蓋半徑是已知的，通過問題一數據分析得到，將此值代入本題模型，并利用兩種方法，求得易拉罐形

26、狀的各項尺寸如下表所示：上頂蓋半徑圓柱體半徑圓臺傾斜角圓柱體高度所用材料圖8利用上表數據，進行繪圖得到右圖8，觀察上表數據和所得圖形，再與實測數據作比較，發現所得數據與實測數據之間的差值已經非常小，因此可以說明易拉罐形狀和尺寸在滿足一定的設計要求下是符合用料最省的。5.4 問題四解題思路本題仍以用料最少為最優設計進行分析計算，首先提出自己設計的易拉罐形狀，然后對這種易拉罐用料及體積進行分析計算建立數學模型，并與問題三中簡化的易拉罐進行比較，最終用兩種方法對模型進行求解，得出自己易拉罐的尺寸標準。提出自己易拉罐的最優設計通過對所測量易拉罐的觀察分析，發現這種易拉罐還不是最省材料的，根據多面體中球

27、體的表面積與其體積的比值最小的原理，提出將易拉罐中的圓臺設計成球臺，其形狀簡圖如下：易拉罐序號名稱備注頂蓋半徑為3cm，厚度為球臺對應的球半徑為,厚度為正圓柱體半徑為3cm，高度為罐壁厚度為罐底半徑為，厚度為模型的建立通過觀察設計形狀簡圖，可以看到這種易拉罐可以分為兩部分：一部分是上邊的球臺，一部分是下邊的圓柱體，在它們的體積和滿足要求的情況下，存在一種用料最省的尺寸標準。球臺部分的求解圖9球臺可近似看成是由圖9中的陰影部分繞其中軸線旋轉而成的立體。設球臺的邊緣弧線的半徑為，球臺下表面半徑為，球臺上表面半徑為，上邊緣端點與中軸線的夾角為，下邊緣端點與中軸線的夾角為，球臺邊緣弧形部分壁厚等于圓柱

28、體的壁厚為，球臺上頂蓋的厚度為。球臺的體積根據上圖可以看出球臺體積是球中半徑為 的圓面以上部分與半徑為的圓面以上部分的體積之差。首先確定出、的大小，的大小為圖中所示的半徑繞豎直軸線的旋轉面以上部分的多面體與其下半部分的錐形體的體積之差，運用對球面坐標積分的方法可求出的體積為：同理可求得體積為：所以可求得球臺體積為：圖10球臺上表面的面積為：球臺弧形部分表面積，根據二重積分用球面坐標對其表面積積分可求得球臺部分表面積為：圓柱體部分的求解設圓柱體的高為，底面半徑為，圓柱體的壁厚為，底面厚度為，如圖10 。圓柱體的體積為：側面表面積為：圓柱體的底面積為：易拉罐的容積由于，所以易拉罐的整體體積看成易拉

29、罐的容積，所以總的易拉罐容積為：易拉罐所用材料的體積由于薄片的體積等于面積乘以厚度，所以易拉罐所用材料的體積為：綜上所述，在易拉罐體積一定的條件下，以總用料最少為目標建立最優化模型如下：模型求解：模型簡化在本模型中有不少數據是不能通過材料最省來確定的，它們的尺寸與其他方面（如焊縫長度、工時少、運輸方便等）有關系，所以通過以上模型不能對其求解，要進行簡化，如下：（1）值的確定是根據易拉罐頂蓋和底蓋所需要的強度來確定的不能因為材料省而使值變小，通過問題一的測量數據可以得到；（2）頂蓋其實不是完全的平面形，而是向上拱起的，是由薄片擠壓而成，并且考慮到罐體的美觀、實用性、運輸方便等，其大小不能通過求材

30、料最省得到，根據問題一的測量數據得到；（3）對于本題假設易拉罐所要滿足的容量求解方法對于本模型在求解的時候有兩種方法：1 直接利用軟件求解2乘子法：構造函數，并對所有未知參數求偏導得到一組方程，然后利用軟件解非線性方程組得到，所有參數的值n 方法一軟件求解本模型屬于最優化模型，所以在求解時可利用軟件編程求解（具體編程過程請參見附錄五）。所得結果，圓柱體半徑，球臺的半徑圓柱體高，與弧形部分有關的角，此時所用材料為n 方法二乘子法和軟件求解構造函數分別對求偏導，并使之為零，得到如下方程組：把帶入以上方程組，利用中的函數編程的得到（具體編程過程請參見附錄六），圓柱體的半徑，球臺所在球體的半徑，圓柱體

31、的高，與弧形部分有關的兩個角度，所用材料為對兩種方法所得結果，求算術平均值得到易拉罐各項尺寸，圓柱體半徑為，圓柱體的高度為，圓柱體上頂蓋為，球臺上頂蓋半徑為，球臺弧形部分壁厚為，圓柱體底面厚度為結果檢驗通過觀察數據可以看到兩種解法求得的結果是比較接近的，但觀察兩種解法會發現，第一種解法明顯簡單與第二種解法，這也說明了數學軟件的在處理大型數學模型中的重要性，但第二種方法也是有道理的，并且在沒有計算機輔助下，它是求解多元最值問題的最基本、最常用的方法之一。1與原來的設計相比較此種設計易拉罐所用的總材料體積為，實際所用材料總體積為（問題三已求），所以這種設計能節省的材料，雖然每個節省的材料并不是很多

32、，但當大量生產的時候，能夠節省的經濟價值也是相當可觀的。2 驗證該設計符合審美標準 通過本題求解發現這種設計，得到的圓柱體的高度為，直徑為，直徑與高度的比值為0.66，這也是基本符合“黃金分割”的。5.5問題五(短文)關于數學建模問題 結合三晝夜的競賽和兩個多月辛苦培訓所獲得的經驗和平時學習過程中得到的知識，然后認真深入的探討和研究題目中給我們的提問，經過反復揣摩得到了以下真諦。所謂數學模型就是用數學符號對一類實際問題或實際系統發生的現象的做近似的描述，而數學建模則是獲得該模型、求解該模型并得到結論以及驗證結論是否正確的全過程,數學建模不僅是了解系統基本規律的強有力工具,而且從應用的觀點來看更

33、重要的是預測和控制所建模系統的行為的強有力的工具.許多重要的物理現象,常常是從某個實際問題的簡化數學模型的求解中發現,并給予明確的數學表述。在對實際現象中利用數學知識建立能夠更好地了解該現象,并且可以應用數學方法來解決的數學模型是比較麻煩的，因為實際現象通常都是極為復雜的,所以不經過系統化理想化的化簡是很難進行研究的。為此，數學建模的系統大體上可分為以下7個關鍵步驟：(1)對選定的實際問題進行觀察、分析(2)對實際問題進行必要的抽象或簡化，作出合理的假設(3)確定出建立的模型中需要的變量和參數(4)根據某種數學方法、各學科中的定律、已推證的算法、相關的經驗(都是正確的)等建立變量和參數間確定的

34、數學關系(數學模型)。實際當中這個關系是非常難以發現(5)解析或近似地求解該數學問題，可以以復雜的數學理論和方法為基礎然后利用計算機編程或者啟發式算法對結果進行求解。對結果進行計算，方法也是很難確定的。(6)得出結果后，用計算出的結果對實際問題中出現的現象進行對比或用某種方法 (例如：已存數據、實驗數據或現場測試數據等)來驗證結果是否正確。(7)如果第 6 步的結果是可行的，那么就可以試用于現實; 如果是否定的，那就要回頭對第(1)(7)進行仔細分析和檢驗，如果還是無果就重復上述建模過程從新開始。如果要對數學建模下定義，上述7個步驟的多次重復的過程就是它最好的詮釋。認真研讀和通過實踐上面的7個

35、步驟后可以發現，其中的三個要素是數學建模過程中的主干。其它部分都是對他們的一個銜接和延伸。然而這三個要素也是數學建模過程中的重點難點，它們分別是：(1) 如何從實際情況著手做出合理的假設,并通過各種手段(例如：數學方法、物力知識、計算機知識、生物知識等)設計出符合題意并且能夠進行求解的合理數學模型。(2) 如何求解模型中出現的數學問題,它可能是非常困難的問題，因為求解用到的知識太多，算法也太廣。(3) 如何驗證模型是正確、可行的。這個步驟也是比較難操作的，因為得出的結果都有隨機性。檢驗的方法不可能都十分相似，所以怎樣才能選取一個合適的方法進行檢驗也較困難。對于本題來說，它的建模步驟在整個論文中

36、已經體現出來了。而它建模的難點分別從以上三個要素進行說明。1. 本題對易拉罐的形狀進行了多種假設，而每種假設都可能因為形狀中某種屬性(例：易拉罐的頂蓋和罐壁的厚度)的變化而對數學模型的建立產生影響。在對模型的建立時要考慮多方面因素，比較麻煩。2. 當根據各種假設建立出數學模型后，同時涌現出了多種數學方法(例：拉格朗日乘子法、求條件極值法、數學軟件等)；而從中選取一個或者多個較準確的方法進行求解也是相當困難的。3. 通過某些算法將結果求出后，是利用其它數學方法算出來的結果進行檢驗還是利用歷史數據或者實測數據進行檢驗?對檢驗的方法也是相當困難。6 模型推廣對于本題中所建立的模型，對于那些形狀與易拉

37、罐相似的容器同樣適用。當容器的各個面的材料不同時，可以根據材料價格的不同，以材料費用最少為目標對原模型進行改進。當容器不需要頂蓋時，同樣可在上面的模型的基礎上進行改進。下面分別針對不同的情況進行說明：1 第一種情況將容器視為正圓柱體，當各個面所用材料不相同時，需要考慮每種材料的價格，此時需要以容量一定，容器的材料成本最低為目標建立最優化模型。由實際情況可知，生產一個正圓柱形容器所需要的總材料費用由易拉罐的頂蓋、圓柱回轉面和下底面的材料費用組成。將各個子費用分別求出相加從而簡化了模型中目標的難度。(1) 目標確定參數確定：設圓柱體回轉面單位面積材料價格為一個單位，上頂材料價格為，下底材料價格為；

38、設圓柱體底面半徑為，高度為。則上頂面的材料費用為：圓柱回轉面的材料費用：下底面的材料費用為：生產一個正圓柱形容器所需要的總材料費用為：圓柱的體積為：(2) 模型建立根據生產材料費用最少可建立如下模型：(3)模型求解在求解模型時可運用求解易拉罐材料最省時的模型的方法，如變量代換求極值、軟件直接求解、乘子法。具體求解方法不再詳述?？汕蟮门c的關系式如下：第二種情況將容器視為無頂蓋的圓柱體，可考慮制造該容器的材料相同，各面的厚度不同，以容量一定，所用材料最少為目標建立最優化模型。由題意可知，正圓柱形容器所需要的總用料體積分為圓柱回轉面的用料體積和下底面的用料體積。因此(1) 目標確定參數確定：設圓柱體

39、回轉面的厚度為一個單位，下底厚為。圓柱回轉面的用料體積為：下底面的用料體積為：則生產一個正圓柱形容器所需要的總用料體積為圓柱的體積為：(2)模型建立于是可建立模型如下：(3)模型求解對本模型求解時，也是運用上面所說的方法進行求解，具體求解過程不再贅述。求得與的關系式如下： 總結綜合本論文中的所有模型可發現一普遍規律，在容器形狀及容積固定的前提下，只要以用料費用最少為目標，根據容器用料費用關系與容器形狀各尺寸的關系式及容積與容器各尺寸的關系式，即可建立相應模型如下：模型說明：表示容器形狀的所有尺寸變量；表示容器用料費用；表示容器的容積。對于該種模型可運用軟件直接求解的方法，也可運用乘子法將模型轉

40、化為方程組的形式在求解。對于尺寸變量為2的情況可直接運用變量代換求極值的方法比較簡單，但對于變量較多時，運用軟件求解是一種很快捷的方法。當要求容器中的某幾個尺寸為常量時，可在模型中直接限制；當要求容器的某幾個尺寸的比例關系一定時，可運用列方程的方法對變量之間的關系進行求解。綜上可以看出對于此類問題具有很強的規律性，由簡單的易拉罐設計問題可以推廣到其它各種容器的優化設計。參考文獻1 姜啟源著.數學模型.北京：高等教育出版社，2002 2 葉其孝.數學建模教育與國際數學建模競賽.中國工業與應用數學學會工科數學雜志社，1994年3 飛思科技產品研發中心.MATLAB7基礎與提高.北京：電子工業出版社

41、，2005年4 同濟大學應用數學系.高等數學(第五版·下冊).高等教育出版社，2002年5 葉其孝.，2005年附錄附錄一 問題三模型求解軟件程序365=1/3*3.1415926*TAN(E)*(R13-R23)+3.1415926*R12*H;26*R12;！E：代表圓臺的傾斜角R1>R2;！R1：代表圓柱體的半徑E<=1.57; ！R2：代表上頂蓋的半徑附錄二問題三模型求解軟件程序1 在文件編輯區建立待求方程組文件并保存function y=fun(x)y=2*pi*x(1)/cos(x(3)+2*pi*x(5)+2*pi*x(1)+x(4)*(pi*tan(x(3

42、)*x(1)2+2*pi*x(1)*x(5),6*pi*x(2)-2*pi*x(2)/cos(x(3)-x(4)*pi*tan(x(3)*x(2)2,pi*(x(1)2-x(2)2)*sin(x(3)/(cos(x(3)2+x(4)*pi/3*(x(1)3-x(2)3)/(cos(x(3)2, 2*pi*x(1)+x(4)*pi*x(1)2, 1/3*pi*tan(x(3)*(x(1)3-x(2)3)+pi*x(1)2*x(5)-3552 在的命令窗口輸入以下命令求解clearx0=1，1，1，1，1；fsolve(fun,x0,optimset)(注：輸出結果中分別代表圓柱體半徑、上頂蓋半徑、圓臺傾斜角、拉格朗日乘子、圓柱體的高度)附錄三 問題三模型改進一軟件程序MIN=3*3.1415926*R22+3.1415926*(R12-R22)/COS(E)+2*3.1415926*R1*H+3*3.1415926*R12;365=1/3*3.1415926*TAN(E)*(R13-R23)+3.1415926*R12*H;R1&