高教社杯全國大學(xué)生數(shù)學(xué)建模競賽題目

1、電池剩余放電時間預(yù)測摘要 鉛酸電池的剩余放電時間是電池性能指標(biāo)中的一個重要參數(shù)，對電池的輸出性能和使用壽命有著至關(guān)重要的影響。但是在生活中不免遇到電池用完了卻沒有備用電池的情況，而且電池的使用時間在我們生活中與我們息息相關(guān)，有時候這會給我們的生活造成很大的困擾。本文就是對電池剩余電量和放電時間的矛盾進(jìn)行討論得出一個相對準(zhǔn)確的電池剩余放電時間預(yù)測的一個模型的建立。 針對問題一：電池的剩余放電量這一問題。首先，根據(jù)題目我們得知鉛酸電池在放電過程中電流時恒定的，并且鉛酸電池的額定保護(hù)電壓是9V，我們根據(jù)附件一給出的數(shù)據(jù)進(jìn)行分析，然后利用MATLAB軟件進(jìn)行圖像擬合，得到一個電壓與時間的放電曲線圖，并

2、且，根據(jù)MATLAB軟件擬合得出的圖像，經(jīng)過figure對圖像的精確處理，最終得出初等函數(shù)來表示各放電曲線。又根據(jù)平均相對誤差定義得出：平均相對誤差（MRE是預(yù)測誤差相對值的平均值，其計算公式為：）根據(jù)公式和附件一中給出的數(shù)據(jù)按照比例篩選出231組數(shù)據(jù)【附件二】，并對數(shù)據(jù)進(jìn)行處理，我們分別求出各放電曲線的平均相對誤差。并且根據(jù)測得電壓都為9.8V時，根據(jù)模型得出電池剩余放電時間分別是多少。針對問題二：根據(jù)數(shù)據(jù)擬合的曲線圖像得出三種情況進(jìn)行討論。用初等函數(shù)表示A20-A100的放電曲線，并根據(jù)放電曲線的平均相對誤差。然后分別求出它們剩余放電時間。然后我們利用求平均數(shù)的方法大概繪制出55A的曲線圖

3、。針對問題三：根據(jù)附件2中的數(shù)據(jù)，利用Excel表格數(shù)據(jù)整理得出圖形【附件三】，分析并計算得出模型1、模型2，通過優(yōu)化得出模型3。由于同一電池在不同衰減狀態(tài)下，在同一電流強(qiáng)度情況下從充滿電開始放電，時間隨電壓不斷變化，所以我們從電壓方面進(jìn)行具體分析。一.問題重述鉛酸電池作為電源被廣泛用于工業(yè)、軍事、日常生活中。在鉛酸電池以恒定電流強(qiáng)度放電過程中，電壓隨放電時間單調(diào)下降，直到額定的最低保護(hù)電壓（Um）。電池在當(dāng)前負(fù)荷下還能供電多長時間（即以當(dāng)前電流強(qiáng)度放電到Um的剩余放電時間）是使用中必須回答的問題。同時還要考慮到電池通過較長時間使用或放置，充滿電后的荷電狀態(tài)會發(fā)生衰減。問題1 附件1是同一生產(chǎn)

4、批次電池出廠時以不同電流強(qiáng)度放電測試的完整放電曲線的采樣數(shù)據(jù)。請根據(jù)附件1用初等函數(shù)表示各放電曲線，并分別給出各放電曲線的平均相對誤差（MRE，定義見附件1）。如果在新電池使用中，分別以30A、40A、50A、60A和70A電流強(qiáng)度放電，測得電壓都為9.8伏時，根據(jù)你獲得的模型，電池的剩余放電時間分別是多少？問題2 試建立以20A到100A之間任一恒定電流強(qiáng)度放電時的放電曲線的數(shù)學(xué)模型，并用MRE評估模型的精度。用表格和圖形給出電流強(qiáng)度為55A時的放電曲線。問題3 附件2是同一電池在不同衰減狀態(tài)下以同一電流強(qiáng)度從充滿電開始放電的記錄數(shù)據(jù)。試預(yù)測附件2中電池衰減狀態(tài)3的剩余放電時間。二問題分析針

5、對問題一，我們對附表中的數(shù)據(jù)進(jìn)行整理，采用數(shù)據(jù)擬合的方法，根據(jù)數(shù)據(jù)擬合的曲線圖像得出三種情況進(jìn)行討論。用初等函數(shù)表示A20-A100的放電曲線，并根據(jù)放電曲線的平均相對誤差。然后分別求出它們剩余放電時間。針對問題二，根據(jù)附件1給出的數(shù)據(jù)，篩選數(shù)據(jù)后擬合出20A-100A之間任一恒定電流強(qiáng)度放電時的放電曲線的模型，三.模型假設(shè)1假設(shè)外界環(huán)境不發(fā)生變化2假設(shè)電池容量恒定五.符號定義及說明符號含義單位T電池剩余時間min(分鐘)U電壓V(伏)I電流A（安）Tn電池剩余時間預(yù)測min（分鐘）Um最低保護(hù)電壓V(伏)六.模型的建立與求解6.1 對問題一建立模型與求解多項式模型：首先對以20A恒定電流放電

6、數(shù)據(jù)進(jìn)行擬合。利用MATLAB建立電壓放電時間的線性模型：，圖 1該模型較大，并且觀察數(shù)據(jù)容易發(fā)現(xiàn)前10組數(shù)據(jù)殘差較大，因此將前10組數(shù)據(jù)舍去，并且采用多項式擬合，擬合模型如下：，圖 2該多項式模型相對線性模型較小，已經(jīng)達(dá)到0.99并且二次項系數(shù)已經(jīng)很小，因此采用二次多項式擬合。20A100A二次多項式模型表1：電流/A模型MRE20A2.2980.99190.035030.58%30A2.0590.98850.041111.1440A1.3110.98930.039211.39%50A0.97060.98830.038791.093%60A0.69090.98920.03672.16%70A

7、0.59130.98860.037430.846%80A0.43790.98960.0350790A0.34460.99010.03383100A0.24940.99130.03097表 1根據(jù)以上二次多項式模型，分別求出和的時間和，則 即為在新電池使用中，分別以30A、40A、50A、60A和70A電流強(qiáng)度放電，測得電壓都為9.8伏時，電池的剩余放電時間，如下表：電流20A30A40A50A60A70A80A90A100A剩余放電時間/min1387.1953.8655.2514.1418.4348.5303.9262.7238.6表 26.2對問題二建立模型與求解根據(jù)問題一建立的電壓與放電

8、時間的二次多項式模型，建立多項式系數(shù)隨電流變化的模型。二次多項式模型中多項式系數(shù)與電流的關(guān)系如下表：電流/A二次項系數(shù)a一次項系數(shù)b常數(shù)項c20-7.25E-08-0.0000758910.5630-1.48E-07-0.000167910.5840-3.35E-07-0.000163910.5350-5.72E-07-0.000177510.4860-8.97E-07-0.000200310.4570-1.4E-06-0.000157910.3880-1.92E-06-0.000181410.3490-2.68E-06-0.000170210.31100-3.35E-06-0.0002689

9、10.28表 3利用MATLAB分別建立系數(shù)與電流的關(guān)系如下： 圖 3 圖 4 圖 5根據(jù)以上模型求出電流為55A時電壓隨時間變化的二次模型的系數(shù)分別是： 則以55A恒定電流放電時電壓隨放電時間的模型為：55A恒定電流放電時放電曲線：圖 66.3對問題三建立模型與求解利用Excel建立不同衰減狀態(tài)下的放電時間與電壓的二次多項式模型，如圖7圖 7其中二次多項式系數(shù)呈現(xiàn)線性關(guān)系，如圖8可以得出不同衰減狀態(tài)下模型系數(shù)與衰減狀態(tài)的關(guān)系，從而預(yù)測出衰減狀態(tài)3下的二次多項式模型：圖 8根據(jù)衰減狀態(tài)3下的模型與已知衰減狀態(tài)3下的數(shù)據(jù)進(jìn)行比較；得到下圖：圖 9針對模型三：于同一電池在不同衰減狀態(tài)下，在同一電流

10、強(qiáng)度情況下從充滿電開始放電，時間隨電壓不斷變化，所以我們從電壓方面進(jìn)行具體分析。7.1 對問題三建立模型與求解根據(jù)附件2中的數(shù)據(jù)，利用Excel表格圖表散點處理，生成圖像如圖1： 圖1 新電池狀態(tài)與放電時間的線性圖觀察圖1得出圖像逐漸遞增，新電池狀態(tài)不斷遞增，電壓逐漸減少為9V。根據(jù)此結(jié)論：新電池狀態(tài)不斷遞增時，電壓逐漸減少為定值9V。第一種情況：根據(jù)附件2中衰減狀態(tài)1的放電記錄得出放電曲線圖2： 圖2 衰減狀態(tài)1放電預(yù)測圖利用圖2所得的圖像進(jìn)行線性預(yù)測，并計算狀態(tài)1的剩余放電時間。從計算結(jié)果對比附件2中的數(shù)據(jù)存在誤差，跟實際不相符，為此不考慮采用圖2的預(yù)測曲線。第二種情況：根據(jù)附件2中衰減狀

11、態(tài)2的放電記錄得出放點曲線圖3： 圖3 衰減狀態(tài)2放電預(yù)測圖 利用圖3所得的圖像進(jìn)行分析，并驗證衰減狀態(tài)2放電所測量的數(shù)據(jù)，驗算結(jié)果與衰減狀態(tài)2的數(shù)據(jù)相差無幾，根據(jù)驗算結(jié)果得出圖3預(yù)測圖可以采用。第三種情況：根據(jù)附件2中衰減狀態(tài)3的部分放電記錄，利用圖3得出的結(jié)論推導(dǎo)得出： U=0-0.01229t2+7.1221t-69.753 模型1通過模型1的轉(zhuǎn)換與推導(dǎo)得出： U=4.4081t-23.492 模型2根據(jù)模型2得出圖4： 圖4 衰減狀態(tài)3的預(yù)測圖利用圖4所得的圖像進(jìn)行分析，并驗證衰減狀態(tài)3放電所測量的數(shù)據(jù)，驗算結(jié)果與衰減狀態(tài)3的數(shù)據(jù)與測量值接近，該預(yù)測圖不僅較好地擬合了電池衰減狀態(tài)3的變

12、化曲線，并較好的預(yù)測了電池衰減狀態(tài)3的剩余放電時間，證明圖4符合。通過Excel軟件得： U=-0.01075t2+5.7353t-47.814 模型3由模型3得： 剩余放電時間t=749(min)模型的優(yōu)點；1、 模型的建立便于理解，問題一的圖像更接近于實際值，且附件一大數(shù)據(jù)支持具有一定的可靠性和普遍性2、 模型的建立簡單明了，簡單化了問題，已于建立。3、 在建模過程中，采用MATLAB進(jìn)行擬合，進(jìn)行了大量數(shù)據(jù)的預(yù)測與篩選，更加優(yōu)化了整個模型，簡單高效。模型的缺點；1、模型的建立比較理想化，忽略實際存在的問題，算出的有結(jié)果有較大的誤差，數(shù)據(jù)由人工計算操作，會有一定的誤差。2、模型未考慮溫度等條件。模型推廣：本模型建立關(guān)于鉛酸電池的優(yōu)化模型，是解決很多實際問題的行之有效的方法，同時，我們還運(yùn)用了E