勾兑酒的酒精含量分析数学建模论文.docx

勾兑酒的酒精含量分析数学建模论文

勾兑酒的酒精含量分析Tom 要本文研究勾兑酒的酒精的相关含量变化的相关问题。在这个问题的研究中,我们采用了典型的微分模型。对恒速加入和抽取酒精在使得酒精浓度的调酒器内变化的情况进行分析,从而构建常微分方程来模拟求解。可以看出问题(1),(2)只需构造一个微分方程即可解出。并参考matlab对微分方程的求解,简化模型。随着勾兑过程不断进行,最终这会勾兑出一个相对固定浓度的酒,而不再变化,即酒精的相对比例不变。关键词:酒精含量变化常微分方程MATLA方程求解一、问题重述此问题为相关量分析问题,有一调酒器,且原始装有一定量的调制酒,现从甲管中向其注入一定浓度的勾兑酒。同时,从乙管中抽取固定量的勾兑酒。试分析调酒器内的酒精含量随时间的变化规律.调酒器内的酒精的变化问题是分析勾兑酒进入调酒器之后,在调酒器内均匀分布,而同时在乙管抽取固定量的勾兑酒之后,各个时间段内酒精含量在调酒器中变化的情况。讨论问题:1、在恒速注入和抽取勾兑酒的情况下,酒精含量在酒器有如何变化?2、半小时后,调酒器内的酒精含量有什么变化规律?二、问题分析由于此问题已简化,只考虑酒精对勾兑酒的变化影响。实际上,勾兑酒的酒精含量的变化是一个很复杂的问题,涉及到很多因素。在本文讨论的方法中,主要采用的是微分方程求解过程。已知恒速注入的酒精含量是一个常量,且随时间推移,调酒器中的总量在不断增加,即酒精含量在不断增加,而同时调酒器内的酒精浓度在逐渐降低。利用这些规律来问题进行模拟求解,建立基本的微分方程模型,模拟酒精含量在调酒器中的加入、抽取情况,且认为在酒精有调酒器中是均匀分布的。确定了以上这些基本问题之后,就通过高等数学中的常微分方程来模拟酒精在调酒器中的变化。三、模型假设1、假设忽略酒精的浓度,视酒精的密度与水的密度相同;2、假设在任一时刻酒精在水中都是均匀分布的;四、符号说明时刻调酒器中的酒精含量时刻调酒器中酒精注入量时刻调酒器中酒精抽取量调酒器中的酒精的增量五、模型建立与分析结合酒精在调酒器内均匀分布,建立酒精含量变化情况的模型,它假设:注入勾兑酒后在任一时刻调酒器中酒精仍均匀分布,那么每次增量就是,勾兑酒会以一定速率持续注入到调酒器内,则调酒器中的酒精含量的变化率会有变化。设时刻调酒器中的酒精总量为;系统处于一种动态平衡中,即成立有关系式:。(即此处认为:酒精含量 = 加入的总量 – 出抽取出的量)1、分析酒精含量变化含有酒精的勾兑酒以恒定的速率从甲中注入调酒器,则满足(为常数)。酒精的抽取量通常被认为是与调酒器中的酒精浓度乘当前调酒器中的总量,列出方程。随着时间的推移,最终的酒精浓度会逐渐下降,,若恒速注入则调酒器中的酒精含量可能最终恒速变化。调酒器中的酒精含量满足:;且时,调酒器内的酒精含量为30公斤,即此时为零时刻则:。组成微分方程组,利用Matlab的dsolve函数求解此微分方程,得: 利用Matlab编程模拟得出如下图形:2、半小时后的酒精含量情况分析 根据上一问求解得到的方程,利用matlab 进行分析,即可求解得到第30分钟的酒精含量为7.9471(L),即如下matlab 运行结果。x1 = 7.9471min_alcohol = 7.5378绘图并固定图像的目标区间有如下(第30分钟做标记):从模拟出的图形中可以看出,从第0分钟开始,调酒器的酒精含量会逐渐下降,并且最终会达到一个最小含量,即min_alcohol,为7.5378(L)。继续当前操作,则在50分钟后,最终酒精含量会以一个固定的比率匀速上升,此时这个含量显然是和勾兑酒的注入有关而和抽取的速率是无关的,因为此刻的酒精浓度保持不变。六、模型的综合评价模型的优点:上文的讨论都是将勾兑过程看作一个简单的微分模型,从而将问题简化而忽略勾兑具体过程,只考注入勾兑酒和抽取勾兑洒对调酒器内酒精含量的影响。这种方法有利于对主要因素的研究从而较易的观察总结规律。模型的推广:此模型可用于实时的,且相对单一的物质的浓度,及物质的量的检测。但本题而言,实际上影响酒的勾兑品质有很多方面,因此,要建立更接近实际情况的数学模型就必须正视影响酒精含量的诸多相关参数之间的关联关系,这就需要多数学模型进行检验。。七、参考文献[1]周义仓等,数学建模实验,西安交通大学出版社,第2版,28,2007.8;[2]朱旭等,MATLAB软件与基础数学实验,西安交通大学出版社,125-129,2008.1;[3]姜启源等,数学建模,高等教育出版社,第4版,p136,2010.8。八、附录酒精含量变化情况的程序:x = dsolve(Dx+(6*x)/(50+4*t)=0.2,x(0)=30,t);%求解所得到微分方程式t = 31;x1 = (2.*t)./25 + (3625.*2.^(1./2

文档评论(0)

1亿VIP精品文档

相关文档