第2l卷第1期 2018年1月 西安文理学院学报(自然科学版) Journal of Xi’an University(Natural Science Edition) Vo1.21 No.1 Jan.2018 文章编号:1008—5564(2018)01-0034-06 电源剩余放电时间预测 赵 燕 (山西省机电职业技术学院基础部,山西长治046001) 摘要:铅酸电池以恒定电流强度放电过程中,电压随放电时间降至保护电压.首先,根据测试数据 用MATLAB拟合出放电曲线模型为:u(t)=ae +ce ,计算出模型的平均相对误差;并根据模型计算 30 A一70 A不同电流强度下,新电池放电得到9.8V时的剩余放电时间.其次,根据九种电流强度下放电 曲线模型的系数,拟合出系数。、b、c、d随电流,变化的函数,得模型: (,,t)=a(,)e“”‘+c(,)e ” ,并绘 出电流强度为55 A时的放电曲线图形;最后,根据三种衰减状态下电压与时间的测试数据,绘出时间随 电压变化的散点图,根据新电池状态放电时间与衰减状态3放电时间的倍数关系,预测出衰减状态3的 剩余时间为194.6 min. 关键词:铅酸电池;放电曲线模型;放电时间 中图分类号:029 文献标志码:A Prediction of Residual Discharge Time of Power Supply ZHAO Yan (Department of Basic Courses,Shanxi Institute of Mechanical&Electrical Engineering,Changzhi 046001,China) Abstract:During the discharge process of a 1cad—acid battery with constant current intensity。 the voltage decreases with the discharge time to the protection voltage.Firstly,according to the test data,the discharge curve model is fitted by MATLAB as“(t)=ae +ce ,and the aver— age relative error of the model is calculated.According to the model,the residual discharge time of the new battery when the discharge of the new battery is 9.8 V is calculated under the different current intensity of 30 A一70 A.Secondly,according to the coeficient of tfhe dis— charge curve model under nine kinds of current intensity,the function of coeficifent a,b,C and d with current change is fitted,and the model A is obtained,and the diagram of discharge curve when the current intensity is 55A is drawn.Finally,according to the test data of voltage and time in three attenuating states,the scatter diagram of time varying with the voltage is drawn.According to the multiple relation between the discharge time of the new battery state and the 3 discharge time of the attenuation state,the remaining time under the attenuation state of 3 is predicted for 194.6 minutes. Key words:lead—acid battery;discharge curve model;discharge time 收稿日期:2017-05-26 作者简介:赵燕(1982一),女,山西长治人,山西机电职业技术学院基础部讲师,硕士,主要从事应用数学研究 第1期 赵燕:电源剩余放电时间预测 35 1 研究背景及问题重述 以2016年全国大学生数学建模竞赛c题为背景,主要解决三个问题: 问题一:同一批次电池出厂时以20 A到100 A不同电流强度下的放电测试数据赛题在图l中给出: 放电时『可(Iin) 2O^ 30A 40a 50^ 60A 7O^ 80^ 90A lOO矗 0 l1.1781 l1.0514 11.O650 11.0821 l1.1043 11.1243 l1.1536 11.1864 11.2l79 2 l0.8913 10.7179 10.6421 10.5650 10.5086 10.t257 10.3"/36 10.3Z5o 10.2750 哇 10.7419 1O.6i71 10.552J. 10.4793 10.,L221 1O.3450 1D.2979 1D.2529 10.2t50 6 10.6288 10.5886 l0.5321 10.,I,657 10.41 57 10.3507 10.3143 10.2843 10.2636 8 l0.5581 10.5807 10.5300 10.1 ̄714 l0.4z64 10.3707 10.3379 10.3107 10.2857 10 10.518l 10.5786 10.5300 10.4786 10.4379 lO。3850 10.3521 10.3236 10.295"t 12 10.4956 10.5793 1O.5321 10.4843 10.4493 10.39 ̄,3 10.3fi07 10.3300 10.3014 14 10.4850 10.5829 1O.5379 lO.4914 1O.45哇3 10.4007 l0.365"/ 10.3S43 10。3038 16 10.4819 10.58‘3 1 0.5414 10.4921 10.4586 10.4050 10.3右86 10.3364 10.3050 18 1O.485o lO.58"(1 lO.5429 10.J,979 10.4629 10.4071 lO.370"/ 10.3379 10.3O57 20 10.d,919 1O.5879 10.5457 10.5000 10.4664 10.4093 10.3721 10.33'/9 10.3050 22 10.5000 10.5900 10.5486 10.5029 10.4686 10.4100 l0.3"/21 l0.3379 10.3043 24 10.50"/5 10.5936 10.54B曷 t0.5057 l0.4693 l0.t107 t0.3721 10.3379 10.3029 26 10.513B 10.5964 10.5514 10.50 ̄d, 10.4707 10.4l14 10.3721 10.3364 1O.301哇 28 :tO.5200 10.5957 10.5536 J.O.5064 10.4700 10.411, ̄ 10。3"/14 lO.3357 10.3O00 30 10.5250 10.59"fl 10.5536 10.5071 10.6,707 10.4114 10.3714 :tO.3343 lO.2993 图1 赛题附件1部分数据的截图 请用初等函数表示九种不同电流强度下的放电曲 线;计算各放电曲线的平均相对误差MRE(从电池的最 低电压为9 V开始,按不超过0.005 V的最大间隔提取 231个电压样本点.这些电压值对应的模型已放电时间 与采样已放电时间的平均相对误差定义为MRE);新电 池使用中,分别以30 A、40 A、50 A、60 A和70 A电流强 度放电,测得电压都为9.8 V时,根据所获得的模型,计 算电池的剩余放电时间分别是多少? 问题二:建立20 A到100 A之间任意一个恒定电流 强度下,放电时的放电曲线模型,并用MRE评估模型精 度.用表格和图形给出电流强度为55 A时的放电曲线. 问题三:同一电池在新电池状态和3种不同衰减状 态下,以同一电流强度从充满电开始放电的记录数据给 出.试预测图2中电池衰减状态3的剩余放电时间. 2 问题分析 2.1 问题一分析 图2赛题附件2部分数据的截图 根据20 A到100 A不同电流强度下的放电测试数据画出电压随时问变化的散点图,用MATLAB工 具箱在散点图上拟合出各自的放电曲线模型;其次,根据MRE的定义,选2O A到100 A各电流强度下的 231个采样点,并用Excel计算九种电流的平均相对误差;根据拟合的放电模型,算出9.8 V时已放电时 间,再由测试数据得出总放电时间,最后相减求出9.8 V时电池的剩余放电时间估计值,可以和9.8 V时 电池的剩余放电时间实际值进行比较,进一步检验模型的精确度. 2.2 问题二分析 根据问题一的放电曲线模型,找规律并建立电压关于电流、已放电时间的二元函数模型u(,,t),由 MRE的定义,选出20 A到100 A不同电流强度下的231个采样点,并用Excel计算九种电流的MRE进 而确定模型精度;再将,=55 A代人二元函数u(,,£)模型中得55 A时的放电模型u(55,t),用MATLAB 画出其放电曲线;最后利用Excel列出放电时间所对应的电压值. 36 2.3 问题三分析 西安文理学院学报(自然科学版) 第21卷 首先,新电池、衰减状态1、2的散点图及衰减状态3的已知点散点图以一定规律分布,用Excel处理 放电时间的采样点,发现数据间的关系,根据关系补齐衰减状态3缺失的放电时间数据;最后,根据补齐 的数据计算出衰减状态3的剩余放电时间. 3 符号说明 4 模型建立与求解 4.1 问题一的模型及求解 4.1.1 不同电流强度下各放电曲线模型 利用MATLAB软件曲线拟合工具箱…,根据20 A一100 A的放电测试数据,以时间t为横轴,电压 u为纵轴绘出画出散点图,模型为: “(t)=口e +ce (1) 模型系数0、b、c、d如表l所示: 表1 电流强度下放电模型系数 图3为电流强度20 A的拟合效果图.从拟合图形看 出,该模型的拟合效果比较好,模型的大部分值与样本值 重合,其他电流强度的拟合图像与该图类似,不再列举. 4.1.2 计算平均相对误差(MRE) 从最低保护电压u 开始按不超过0.005 V的最大间 隔提取231个电压样本点. 平均相对误差: 1∑1/7l搓奎骂璺 垄 样本时间. 1算出各电流的平均相对误 …l…~。…一… 差如表2所示. 模型“(t)=0e +Ce 在各个电流强度下的平均相对 误差都没有超过1%,误差比较小,模型精度较高. 图3 电流强度为20 A的模型拟合图 (t) 第1期 赵燕:电源剩余放电时间预测 37 4.1.3 电压为9.8 V时电池的剩余放电时间估计 9.8 V时电池的剩余放电时间估计值=电池总放电时间一9.8 V时已放电时间 其中,电池总放电时间可从赛题附件1(图1)中查到,9.8 V时已放电时间可以根据 (t)的模型求出. 9.8 V时的剩余放电时间如表3所示. 表3 9.8 V时的剩余放电时间 说明:9.8V时电池的剩余放电时间实际值,可以根据赛题附件1的放电测试数据计算得到. 从表3可以看出,9.8 V时电池的剩余放电时间估计值与实际值很接近,进一步说明模型 (t)的精 确度较高. 4.2 问题二的模型建立与求解 4.2.1 电压随时间和电流变化的模型建立 观察问题一的九个放电曲线模型很类似,仅仅系数不同.在MATLAB中以电流强度为横轴、系数为 纵轴用多项式P1,+P2,+P3 +P4,+P5/4+P6,3+P7,2+P8,+P9进行拟合 ,得到系数a(I)、 b(,)、c(,)、d(,)(见表4).并代人模型(1)得电压随时间和电流变化的模型: u(1,t)=0(,)e ‘ ‘+c(i)e ‘ 表4 系数。(,)、6(,)、c(,)、d(,)多项式拟合系数表 (2) 4.2.2 电压随时间和电流变化的模型MRE评估 根据问题一筛选出的231个电压样本点,算出各自的平均相对误差,再求出MRE的平均值,如表5 所示. 表5 电压随时间和电流变化的模型各电流的平均相对误差 38 西安文理学院学报(自然科学版) 第21卷 模型 (,,t)=n(,)e“” +c(,)e 的平均相对误差的均值为0.43%,也非常小,说明该模型的精 度较高. 4.2.3 电流强度为55A时的放电曲线图形及电压值 将,=55 A代人u(,,t)即得出电流强度为55 A时的放电时间数学模型为: u(55,t)=(z(55)e +c(55)e 根据55 A的数学模型在Excel中计算并列出各个时间点的电压部分截图如图4所示. C E l:55^任童-|捌的I电压囊 量一j放电时间(min)电压(v)放电时间(m电压(v)放电时间 电压(v)放电n:,tfq·电压(v)放电时问-a电压(v)放电时问t电压(v) .量i 哇i 0 10.44948 2 10.44809 50 l0.4i457 52 10.41317 l00 l0.37936 102 10.37794 150 10.34356 152 10.34211 200 10.30668 202 10.30517 250 10.2678 ̄ 252 10.26627 一曼.~6『 4 10.4467 6 10.4453 54 10.41l76 56 10.41036 104 10.37652 106 10.3751 154 10.3 ̄066 156 10.3392 204 10 30366 206 10.30215 254 10.26466 256 10.26303 7 8 9 e l0 8 10.44391 10 10.44251 l2 10.44112 l4 10.43973 58 10.40896 6O 10.40756 62 10.40615 64 l0.40475 108 10。37368 110 l0.37226 i12 10.37083 l14 10.36941 158 10.33775 16O 10.33629 162 10.33483 164 10.33337 208 10.30063 210 10.29911 212 10.29759 214 10.29606 258 10.26141 2fi0 10.25978 262 10.2581哇 264 10.25649 l I 12; l6 l0.43833 18 10.43693 66 10.40334 68 10.4O194 1l6 10.36798 118 10.36656 166 10.33191 168 10,33044 216 10.29 ̄53 218 10.29299 266 10.25484I 268 l0,25319 13} 20 10.43554 70 10.40053 120 10.36513 l70 10.32897 220 10.29146 270 10.25153 i4 f 1§c l5 17 22 10.43414 24 10.43275 26 10.43135 28 10.42995 72 10.39912 74 10.39772 76 10.39631 78 l0.3949 122 10.3637 124 10.36227 126 10.36084 128 10.35941 172 10.32751 174 10.32603 176 10.32456 178 10.32308 222 10.28991 22哇10.28837 226 l0.28682 228 10.28526 272 10.24986 274 10.24818 276 10.2465 278 10.24 ̄81 l8 l9: 2O 21 30 10.42856 32 10.42716 34 10.42576 36 10.42436 80 10.39349 82 10.39208 84 10.39067 86 10.38926 l30 l0.35797 132 10.3565 ̄ 134 10.3551 136 10.35366 180 10.32161 182 10.32012 184 10.31864 186 10.31716 230 l0.28371 232 10.28214 234 10.28058 236 10.27901 280 10.24311 282 10.24141 284 l0.2397 286 10.23798 22 23 38 10.42297 40 10.42157 88 10.38785 90 10.38643 138 10.35223 140 10.35078 l88 10.31567 190 10.31418 238 10.27743 2哇0 10.27585 288 10.23625 290 10.23452 24 25 42 10.42017 44 10.41877 92 10.38502 94 10.3836 142 l0.34934 144 10.3479 192 10.31268 194 10.31119 242 l0.27427 244 10.27268 292 10.23277 294 10.23102 26 27 46 10.41737 48 10.哇1597 96 l0.¥8219 98 10.38077 146 10.34645 148 10.345O1 196 10.30969 198 10.30819 2哇6 10.27108 248 10.26948 296 10.22926 298 10.22749 图4 55 A放电时间所对应的电压值 电流强度为55 A时的放电曲线用Matlab绘图如图5所示. 图5 电流强度为55 A时的放电曲线图u(55。t) 4.3 问题三的模型建立与求解 仔细观察赛题附件2 r{1新电池状态和j种衰减状态的放电时间数据,通过Exeel计算发现绝大部分 数据满足以下规律:要是 一 . 6, .32, -.62,可以根据 新电池状态 1.62 哀减状态3,补齐衰减状态3的剩余放电时问的缺失数据见截图6. 以电压为横轴,以放电时间为纵轴画出新电池状态、衰减状态1、2及已知衰减状态3的散点图见 图7,发现图像中每条曲线之间的间距有一定的规律,补齐后并画出四种状态下的放电时间见图8. 第1期 赵燕:电源剩余放电时间预测 39 I 9.240 9.235 9.230 瓤 ~卫 A 殳王8—0一 他~ : B 600.4 9.630 9.625 9.620 9.615 9.610 9.605 D 660.1 663.5 665.8 E 殂 F : 卯 {G 9.370 9.365 9.360 H 741.2 743.2 7LI3.9 J 766.1 767.8 767.9 K 9.110 9.105 9.100 L 781.8 782.3 782.7 电压(v)衰减态3放电时间电压(v) 衰减态3电压(V)衰减态3电压( 衰减态3电压(v)衰减态3电压(v)衰减态3 9.1"60 9.755 9.750 9.745 9.740 9."/35 9.'730 9。"/25 9."/20 9.715 9.710 9."/05 9.700 9。695 9.690 9.685 9.680 9.675 9.670 9.665 9.660 9.655 9.650 9.500 9。495 9.哇90 706.6 708.3 711.0 603.3 606.3 607.6 610.2 612.7 617.2 619.4 666.5 669.6 671.6 673.0 6"/5.0 675.5 677.3 679.9 680.7 682.4 685.I 686.1 688.2 690.9 9.485 9.480 9.475 9.4"/0 9.465 9.460 9.455 9.450 9.445 9.440 9.435 9.43O 9.425 9.42O 711.2 713.3 713.6 715.1 718.0 719.2 719.6 720.3 722.1 724.8 725.5 727.8 728.8 729.8 9.355 9.350 9.345 9.340 9.335 9.330 9.325 9.320 9.315 9.310 9.305 9.300 9.295 9.290 745.2 745.2 746.9 748.0 749.0 749.9 750.0 751.5 752.7 753.5 755.0 "/55.7 756.5 756.9 9.225 9.220 9.215 9.210 9.205 9.200 9.195 9.190 9.185 9.180 9.175 9.170 9.165 9.160 768.8 768.8 769.4 770.2 770.6 771.6 771.8 '/72.2 774.1 774.3 774.3 775.9 776.0 776.2 9.095 9.090 9.085 9.080 9.075 9.070 9.065 9.060 9.055 9.050 9.045 9.040 9.035 9.030 783.5 783.7 784.1 784.6 785.1 785.6 786.4 786.5 786.5 "/86.5 787.1 787.2 "188.0 5'89.1 9.600 9.595 9.590 620.5 622.2 624.3 9.585 9.580 628.9 630.8 631.6 9.5"/5 9.570 9.565 9.560 635.2 63'7.8 638.9 642.7 64q.3 644.9 64"/.7 649.7 651.0 653.5 655.5 659.6 9.555 9.550 9.545 9.540 9.535 9.530 9.525 9.520 9.5l5 9.510 9.505 691.5 694.7 696.4 697.1 9.415 9.410 9.405 9.400 730.0 732.1 /'33。4 735.5 9.285 9.280 9.275 9.270 759.8 759.9 "/61。2 761.6 9.155 9.150 9.145 9.140 777.8 778.3 7"/9.0 779.4 9.025 9.020 9.O15 9.010 5'89.4 "/89.6 789.9 789.9 699.6 700.1 702.0 9.395 9.390 9.385 736.6 '/36.8 738.8 9.265 9.260 9.255 762.0 764.0 ''/64.0 9.135 9.130 9.125 779.8 779.9 780.1 9.005 9.000 790.4 790.8 9.645 9.640 9.635 703.3 704.5 9.380 9.3"/'5 739.3 741.0 9.250 9.245 765。1 766.0 9.120 9.115 780.1 781.2 图6 补齐衰减状态3的剩余放电时间的缺失数据 1 400 1 200 1 ooo 一1 400 l 200 l 000 电压/v —一电脏,V 新f乜池状态 800 600 —新电池状态衰减状态l 800 600 一衰减状态l ……400 …衰减状态2 衰减状态3 衰减状态2 衰减状态3 400 200 …200 O 0 一 g …罱宝宝= …品 一 g …罱昌 …蠹 rq I'M q r N r N N 图7 新电池状态及衰减状态l、2及衰减状态3 已知放电时间采样散点图 图8 补齐衰减状态3后。新电池状态及衰减 状态】、2、3的放电时间散点图 根据预测数据得出衰减状态3总的放电时间为790.8 rain,数据样本衰减状态3的已放电时间为 596.2 rain,故其剩余放电时问为194.6 rain. [参考文献] [1] 周凯,宋军全,邬学军.数学建模竞赛入门与提高[M].杭州:浙江大学出版社,2012:25—32 [2] 王新华.应用数学基础[M].北京:清华大学出版社,2010:28—33. [责任编辑王新奇]