*,/,25,安徽财经大学,Anhui University of Finance&Economics,1959,数模方法之五步法,*,科技论文中的方法,1.3、稳定性与稳健性,1.2、灵敏度分析,主讲:,朱家明,Mathematical model of five step method,数,学,1.1、五步方法,化,最,优,数,学,建,模,讲,座,循序渐进,追求卓绝,安徽财经大学,Anhui University of Finance&Economics,1959,建,模,电话:,18226682718,邮箱,:,2021/5/23,1.1、五步方法,1、,五步方法概要,2、,五步方法详解,1.2、灵敏性分析,1、问题的提出,2、最佳售猪时间,x,关于,价格下降速率,r,的灵敏性,3、最佳售猪时间,x,关于,生长率,g,的灵敏性,4、灵敏,性的相对改变量,1.3、稳定性与稳健性,1、关于稳键,性,2、,r,g不是常数时对模型结果的影响,1.4、小结,1.5、练习,题,数模方法之五步法,2021/5/23,1.1、五步方法,1、,五步方法概要,数学模型解决问题的一般过程分五步,称之为,五步方法。,定义:,五个步骤:,提出问题;,选择建模方法;,推导模型的数学表达式;,求解模型;,回答问题。,2021/5/23,1.1、五步方法,2、五步方法详解,例1.1、,一头猪重200磅,每天增重5磅,,饲养每天需花费45美分。猪的市场价格,为每磅65美分,但每天下降1%,求出售,猪的最佳时间,。(1磅=0.454kg),提出问题:,即如何用数学语言来表达问题。,列出问题涉及的变量,包括恰当的单位;,写出关于上述变量所做的假设,列出已知的或假设的这些变量之间的关系式(等式和不等式);,用明确的数学语言写出问题的目标的表达式。,变量、单位、等式、不等式、假设和目标表达式等构成,完整的问题,。,2021/5/23,1.1、五步方法,例1.1中,全部的变量包括:,猪的重量,w,(磅),,,从现在到出售猪期间经历的时间,t,(天),t,天饲养猪的花费,C,(美元),猪的市场价格,p,(美元/磅),售出生猪所获得的收益,R,(美元),我们最终获得的净收益,P,(美元),。,其他相关的参(非变)量:如猪的初始重量,(200磅),等。,写出关于上述变量所做的假设,考虑到参量在模型中的影响。,猪的重量从初始的,200磅,按每天5磅增加有,这里把变量的单位带进去,可以检查所列式子的意义.,该问题涉及到的其他假设包括:,2021/5/23,1.1、五步方法,售价,饲养成本,收益,利润,假设,t,0,目标:求利润或净收益,P,的最大值。,为了便于参考,下面对第一步所得的结果进行了如下的归纳(见下表),2021/5/23,1.1、五步方法,变量,:,t=,时间(天),w=,猪的重量(磅),p,=猪的价格(美元/磅),C,=饲养,t,天的花费(美元),R,=售出猪的收益(美元),P,=净收益(美元),假设,:,w=,200+5,t,p,=0.65-0.01,t,C,=0.45,t,R=p,w,P=R-C,t0,目标,:求的最大值,注意:第一部分三个阶段(变量、假设、目标)的确定不需要按特定的顺序。,图1-1 售猪问题的第一步的结果,2021/5/23,1.1、五步方法,选择建模方法:,即如何用数学方法来获得解。,许多问题都可表成一个已有有效方法的标准形式.,应用数学的多数研究,包含确定问题的一般类别,并提出解决该类问题的有效方法。,在应用数学领域中有许多的文献,并且不断取得许多新的进展。一般很少有学生对选择较好的建模方法有经验或熟悉参考文献。,注意:,下面除了极少例外,一般都给定所用的建模方法。如例1.1可定位为,单变量优化问题,或极大极小化问题,建模方法,为:设,y,=,f,(,x,),在,x,S,处是可微的,若,f,(,x,),在,x,处达到极大或极小,则,f,(,x,)=0,。详细可参阅微积分中导数应用部分的内容.,2021/5/23,1.1、五步方法,推导模型公式:,即要把第一步得到的问题应用于,第二步,写成所选建模方法需要的标准形式,以,于我们运用标准的算法过程求解。,如:,例1.1把问题中的变量名改换一下,在算法上就比较方便,。,P=R-C=pw-,0.45,t,=(0.65-0.01,t,)(200+5,t,),-,0.45,t,记,y,=,P,作为求最大值的目标变量,,x=t,作为自变量,我们的问题就化为在集合S=,x,:,x,0上求下面函数的最大值:,y,=,f,(,x,)=(0.65-0.01,x,)(200+5,x,),-,0.45,x,.,这是我们最熟悉不过的求一元函数极值问题。,2021/5/23,1.1、五步方法,利用第二步中确定的标准过程,求解,这个,模型,。,如本例中即对,y,=,f,(,x,)=(0.65-0.01,x,)(200+5,x,),-,0.45,x,在区间,x,0上,求最大值。,如图可知,y,=,f,(,x,),关于,x,是,二次的曲线图,易得,f,(,x,)=,-0.1,x,+0.8,则在点,x,=8,处,f,(,x,)=0.,由,f,在区间(,-,8,)上单升,而在区间(8,+,)上单减.,故点,x,=8,是整体最大值点.,且有,f,(8)=133.20,从而点(,x,y,)=(8,133.20,)是,f,在整个实轴上的整体,最大值点,也,是区间,x,0上,的最大值点。,图1-2 售猪问题的净收益,f,(,x,),关于时间,x,的曲线图,0,5,10,15,20,126,128,130,132,134,x,f,(,x,),y,=,-,0.05,x,2,+0.8,x+,130,2021/5/23,1.1、五步方法,回答问题:,回答第一步提问“何时售猪可以达到最大净收益.,由第四步我们得到的答案是在8天之后,可以获得净收益133.20美元。只要第一步假设成立,这一结果就是正确的。,相关的问题及其他不同的假设可以按照第一步中的做法调整得到。由于我们处理的是一个实际问题(,一个农民决定何时出售他饲养的生猪,),在第一步中会有一个风险因素存在,因此通常有必要研究一些不同的可能,这一过程称为,灵敏性分析,。我们将在下一节进行讨论。,本节主要介绍五步方法,下面将这一方法总结归纳成如下图表,以便以后参考.,2021/5/23,1.1、五步方法,第一步、提出问题.,列出问题涉及的变量,包括恰当的单位;,注意不要混淆了变量和常量;,列出你对变量所做的全部假设,包括等式和不等式;,检查单位从而保证你的假设有意义;,用准确的数学表达式给出问题的目标。,第二步、选择建模方法.,选择你问题的一个一般的求解方法;,一般地,这一步的成功需要经验、技巧的对相关文献有,一定的熟悉程度;,在本书中,我们通常会给定要用的建模方法。,第三步、推导模型的公式:,把第一步中得到的问题重新表达成第二步选定的建模,方法需要的形式;,图1-3 五步方法图,2021/5/23,1.1、五步方法,你可能需要将第一步中的一些变量名改成与第二步所用,的记号一致;,记下任何补充假设,这些假设是为了使在第一步中描述的问题与第二步中选定的数学结构相适应而做的。,第四步、求解模型.,将第二步中所选方法应用于第三步得到的表达式;,注意你的数学推导,检查是否有错误,答案是否有意义;,采用适当的技术,计算机代数系统、图形、数值计算的,软件等都能扩大你解决问题的范围,并减少计算错误.,第五步、回答问题.,用非技术性的语言将第四步中的结果重新表述;,避免数学符号和术语;,能理解最初提出问题的人就应该能理解你给出的解答.,图1-3 五步方法图(续),2021/5/23,1.2、灵敏性分析,1、问题的提出,灵敏性分析,是数学建模的一个重要方面,具体内容与所用的建模方法有关,关于它的讨论贯穿本书,下面仅对单变量优化问题进行灵敏性分析.,上用售猪说明五步法,图1-1列出了求解的所有假设,虽然数据和假设都有非常详细的说明,但还要再严格检查,由于,数据,是由,测量、观察,有时甚至完全是,猜测,得到的,故要考虑数据的不准确的可能性。,上概要介绍五步法,从假设开始,但难保证假设都正确.故要考虑结果对每一条假设的敏感程度即,灵敏性,.,可靠性高的数据:,生猪现在的重量、猪现在的价格、每天饲养的花费等易测量,确定性大;,可靠性低的数据:,猪的生长率,g,和价格的下降速率,r,.,2021/5/23,1.2、灵敏性分析,2、最佳售猪时间,x,关于价格下降速率,r,的灵敏性,粗分析,前面我们假定,r,=0.01,美元/天,现在假设,r,的实际值是不同的,对几个不同的,r,值,重复前面的求解过程,我们会对问题的解关于,r,的敏感程度有所了解.,即给定,r,对,y,=,f,(,x,)=(0.65-,r,x,)(200+5,x,),-,0.45,x,求导,令,f,(,x,)=0,,可得相应,x,值,下表1-1给出了选择几个不同的,r,值求出,x,的计算结果。,表1-1 售猪问题中最佳售猪时间,x,关于价格的下降速率,r,的灵敏性,r,(美元/天),x,(天),r,(美元/天),x,(天),0.008,0.009,0.01,0.011,0.012,15.0,11.1,8.0,5.5,3.3,2021/5/23,1.2、灵敏性分析,将上表1-1中的数据绘制在如下图1-4中。,图1-4 售猪问题中最佳售猪时间,x,关于价格的下降速率,r,的曲线,x,(天),r,(美元/天),2,4,6,8,10,12,14,16,0.008,0.009,0.010,0.011,0.012,我们可以看到售猪的最优时间,x,对参数,r,是很敏感的.,x,对价格下降速率,r,灵敏性的系统分析,将,r,作为未知的参数,仍按前面的步骤求解(见下页):,2021/5/23,1.2、灵敏性分析,出售价格:,p,=0.65-r,t,;,目标函数:,y,=,f,(,x,)=(0.65-,rx,)(200+5,x,),-,0.45,x,=130+2.8,x,-200,rx,-5,rx,2,;,求导,f,(,x,)=2.8-,200,r,-10,rx,;,使,f,(,x,)=0的点为,x=,(7-500,r,),/,25,r,.,若要,x,0,只要0,0.014,在0,+)上都有,f,(,x,)0,最佳售猪时间为,x=,0.图,1-5给出了,r=,0.015的情况.,图1-5 售猪问题的净收益,f,(,x,),在r=0.015关于时间,x,的曲线图,0,5,10,15,20,90,100,110,120,130,x,f,(,x,),y,=-0.075,x,2,-0.2,x+,130,2021/5/23,1.2、灵敏性分析,3、最佳售猪时间,x,关于生长率,g,的灵敏性,前面我们假定,g,=5,磅/天,一般地,我们有如下步骤,出售重量:,w,=200+g,t,;,目标函数:,y,=,f,(,x,)=(0.65-,0.01,x,)(200+g,x,),-,0.45,x,=130+0.65g,x,-2.45,x,-0.01,gx,2,;,求导,f,(,x,)=0.65g-,2.45-0.02,gx,;,使,f,(,x,)=0的点为,x=5,(13,g,-49),/,2,g,.,若要,x,0,最佳售猪时间可由,x=5,(13,g,-49),/,2,g,给出,图1-6 给出了最佳售猪时间和生长率,g,之间的关系.,图1-6 售猪问题中最佳售猪时间,关于生长率,g,的曲线图,3,4,5,6,7,-10,-5,0,5,10,g,x,x=5,(13,g,-49),/,2,g,15,2021/5/23,1.2、灵敏性分析,4、灵敏,性的相对改变量,意义:,相对改变量比绝对改变量更自然,、,更实用,例如,r,的10%下降导致了,x,的39%的增加,g,的10%下降导致了,x,的34%的下降.,x,对,r,的灵敏性:,对售猪问题中,由,x=,(7-500,r,),/