1、降落傘的選購模型摘要本模型研究的是降落傘的選購方案問題，目的是在滿足空投要求的條件下，使費用最少。為了方便對降落傘進行受力分析，我們把降落傘和其負載的物資看做一個整體，忽略了傘和繩子的質(zhì)量，并假設(shè)降落傘只受到豎直方向上空氣阻力和重力的作用。通過對降落傘在空中的受力情況的分析建立起了高度與時間的方程，然后以高度與時間的方程作為擬合曲線與題中給出的時間與高度的數(shù)據(jù)進行擬合，得出阻力系數(shù)k的值。我們建立了速度與質(zhì)量的方程，并證明其為嚴格增函數(shù)（證明過程見建模與求解）。由于題中已限制降落傘的最大落地速度為20m/s，所以當速度為20m/s時，傘的承載量最大。建立高度與時間，速度與時間的方程組，代入最大

2、速度20m/s，高度500m,傘的半徑（題中已給出可能選購的每種傘的半徑），分別計算出每種傘的最大承載量。最后運用LINGO軟件進行線性規(guī)劃求解得：x1=0,x2=0,x3=6,x4=0,x5=0.即購買半徑為3m的降落傘6個時總費用最少為4932元。關(guān)鍵字：線性規(guī)劃、空氣阻力系數(shù)、擬合一、問題的重述為向災(zāi)區(qū)空投救災(zāi)物資共2000kg，需選購一些降落傘。已知空投高度為500m，要求降落傘落地時的速度不能超過20m/s。降落傘面為半徑r的半球面，用每根長 L, 共16根繩索連接的載重m的物體位于球心正下方球面處，每個降落傘的價格由三部分組成。傘面費用C1由傘的半徑r決定，見表1；繩索費用C2由繩

3、索總長度及單價4元/米決定；固定費用C3為200元。 表1 r（m）22.533.54費用（元）651703506601000降落傘在降落過程中受到重力作用外還受到的空氣阻力，可以認為與降落速度和傘的受力面積的乘積成正比。為了確定阻力系數(shù)，用半徑r=3m、載重m=300kg的降落傘從500m高度作降落試驗，測得各時刻的高度 ，見表2。表2 時刻t（s）036912151821242730高度h（m）500470425372317264215160108551試根據(jù)以上條件確定降落傘的選購方案，即共需多少個，每個傘的半徑多大（在表1中選擇），在滿足空投要求的條件下，使費用最低。二、模型的假設(shè)1、

4、 假設(shè)空投物資的瞬時傘已打開。2、 空投物資的總數(shù)2000kg可以任意分割。3、 空氣的阻力系數(shù)與除空氣外的其它因素無關(guān)。4、 降落傘和繩的質(zhì)量可以忽略不計。5、 假設(shè)降落傘只受到豎直方向上的空氣阻力作用。三、符號說明1、m 降落傘的負載重量2、g 重力加速度3、a 降落傘的加速度4、k 空氣的阻力系數(shù)5、S 降落傘的傘面面積6、 v 降落傘的速度7、 H 降落傘的位移8、 h 降落傘離地高度9、x1,x2,x3,x4,x5 分別為每種傘的個數(shù)四、問題的分析由題意可知每個傘的價格由三部分組成：三面費用C1、繩索費用C2、固定費用C3。傘面費用由傘的半徑r決定；繩索費用C2由繩索的長度及單價決定

5、，由圖一可知繩索的長度又由降落傘的半徑?jīng)Q定即；固定費用為定值200。因為題中已給出每種傘面的半徑，所以每種傘的價格為定值。要想確定選購方案，即共需半徑（在題中給出的半徑中選擇）為多大的傘的數(shù)量，在滿足空投物資要求的條件下使總費用最少。因此，我們需要確定每種傘的最大承載量。然后進行線性規(guī)劃，確定總費用和每種傘的個數(shù)。要確定最大載重量，我們需對降落傘進行受力分析（如圖二）。降落傘在降落過程中除受到豎直向下的重力作用外還受到豎直向上的空氣阻力的作用，而由題可知空氣阻力又與阻力系數(shù)、運動速度、傘的受力面積有關(guān)。運動速度和受力面積是已知的，所以要想確定每種傘的最大承載量，就必須先要確定空氣的阻力系數(shù)。圖

6、一 圖二 對圖二的分析可知降落傘的運動狀態(tài)是做加速度趨近于0的加速運動。因此，我們可以建立一個位移與時間的函數(shù)關(guān)系式，在根據(jù)題中所給的數(shù)據(jù)擬合出阻力系數(shù)k的值。然后再建立一個速度與時間的函數(shù)關(guān)系式，兩個關(guān)系式聯(lián)立求解出最大載重量（其中高度和速度由題目已經(jīng)給出）。最后用LINGO軟件進行線性規(guī)劃算出問題要的結(jié)果。五、建模與求解（1）首先確定阻力系數(shù)K為了方便對物資進行受力分析，我們把降落傘和物資看作一個整體如圖二。由假設(shè)5可知物體A只受到豎直向上的空氣阻力和豎直向下的重力作用。又由題可知空氣阻力與降落速度v和傘的受力面積S的乘積成正比。則物體A在豎直方向上受到的合外力為：由運動學方程：得由物體位

7、移H和時間的二次微分等于加速度建立方程得: 用MATLAB解微分方程得：(程序見附錄【1】)則題目已經(jīng)給t-h數(shù)據(jù)為：時刻t（s）036912151821242730高度h（m）500470425372317264215160108551對給定的數(shù)據(jù)以為擬合函數(shù)進行擬合，r=3m,m=300kg,g=9.8,得出k=2.9377 。（程序見附錄【2】）（2）求解最大承載量用速度對時間的微分等于加速度，且v0=0建立方程組得：用MATLAB解得(程序見附錄【3】)由前面的和函數(shù)建立方程組得：k=2.9377,g=9.8,r=2 2.5 3 3.5 4因為降落傘在下落過程中其質(zhì)量是不變的，所以我們

8、把關(guān)系式中t看做一個定值，則關(guān)于m的方程為從上式我們可以知道是關(guān)于m的單調(diào)遞增函數(shù)（證明見附件【7】），并且如果存在平衡狀態(tài)則必須滿足，那么 而又通過對 分析，只有在，這與實際矛盾，故降落傘是一直做加速度減小的加速運動，不存在平衡狀態(tài)。因此，求最大載重量取傘在下降到地面的瞬間達到最大速度,此時，由方程組調(diào)用MATLAB分別解得半徑為r的降落傘在滿足空投條件下的最大載重量如下表：（程序見附錄【5】）r（m）22.533.54最大承載（kg）150.6787235.4355339.0272461.4536602.7150取整(kg)150235339461602（3）線性規(guī)劃求解數(shù)量和費用由分析可

9、知每種傘的單價：由題可知為:r（m）22.533.54費用（元）651703506601000為： 為固定值即：由以上數(shù)據(jù)求得每種傘的單價見下表：r（m）22.533.54單價C446596.3821.51176.81562取整44659682211771562我們設(shè)每種傘分別取x1,x2,x3,x4,x5個,則其目標函數(shù)為：z=446x1+596x2+822x3+1177x4+1562x5對其進行優(yōu)化求解z的最小值，就是所需的最小費用。由分析可知其限制條件如下：s.t. 150x1+235x2+339x3+461x4+602x5>=2000;(x1，x2，x3，x4，x5);用LING

10、O求解得（程序見附件【6】）x1=0,x2=0,x3=6,x4=0,x5=0。最少總費用為4932元。六、模型的評價與改進優(yōu)點：1、 本模型的求解過程大量的運用了電腦軟件，使得計算更加精確。缺點：1、 本模型未考慮降落傘打開的時間，將其假設(shè)成在下降時傘就已經(jīng)打開。2、 由于在實際生活中降落傘還受到風向的影響，本模型假設(shè)的是理想的狀態(tài)下（無風）改進：由于本模型假設(shè)的是在物資拋落的瞬時傘已打開，而在實際情況中物資拋落后應(yīng)有一段自由落體運動。在模型的改進時應(yīng)考慮到這一點，以便讓模型更切合實際。七、參考文獻1、數(shù)學實驗蕭樹鐵主編 高等教育出版社 1999 7 1附錄【1】求解位移的程序H=dsolve

11、('m*D2H+k*S*DH=m*g','H(0)=0,DH(0)=0','t')解得：g/k2/S2*m2*exp(-k*S/m*t)+g/k/S*m*t-1/k2/S2*m2*g附錄【2】擬合k程序建立一個名為myfun的m文件function F=myfun(x,xdata)s=2*pi*32;m=300;g=9.8;F=500-m2*g/(x(1)2*s2)*exp(-x(1)*s*xdata/m)-m*g*xdata/(x(1)*s)+m2*g/(x(1)2*s2);在matlab command window中輸入下列命令：xdata

12、=0 3 6 9 12 15 18 21 24 27 30;ydata=500470425372317264215160108551 ;x0=1;x=lsqcurvefit(myfun,x0,xdata,ydata)附錄【3】求解速度程序v=dsolve('m*Dv+k*S*v-m*g=0','v(0)=0','t')解得：附錄【4】在v-t,m函數(shù)中對m求二階導(dǎo)數(shù)syms m t g S k f=g*m/(k*S)-g*m/(k*S)*exp(-k*S*t/m); diff(f,m2)求得：-g/m3*t2*k*s*exp(-k*s/m*t)附

13、錄【5】 求最大載重量在matlab中建立一個名為myfun的m文件，如下：function F=myfun(x)r=2.5;g=9.8;k=2.9377;s=2*pi*r2;F=x(1)2*g/(k2*s2)*exp(-k*s*x(2)/x(1)+x(1)*g*x(2)/(k*s)-x(1)2*g/(k2*s2)-500;g*x(1)/(k*s)-g*x(1)/(k*s)*exp(-k*s*x(2)/x(1)-20;在matlab中command window中輸入以下命令：x0 = 1; 1; % 初始點options=optimset('Display','ite

14、r'); % 顯示輸出信息x = fsolve(myfun,x0,options)在m文件中更改r的值，然后在命令窗口重復(fù)輸入以上命令就可分別求出不同半徑的降落傘的最大載重量。分別求解可得最大載重量如下表：r（m）22.533.54m150.6787235.4355339.0272461.4536602.7150附錄【6】優(yōu)化求解min=446*x1+596*x2+822*x3+1177*x4+1562*x5;150*x1+235*x2+339*x3+461*x4+602*x5>=2000;x1>=0;x2>=0;x3>=0;x4>=0;x5>=0;gin(x1);gin(x2);gin(x3);gin(x4);gin(x5);求解得：Global optimal solution found. Objective value: 4932.000 Extended solver steps: 0 Total solver iterations: 0 Variable Value Reduced Cost X1 0.000000 446.0000 X2 0.000000 596.00