基于分数阶微积分的Kelvin_Voigt流变模型_郭佳奇

基于分数阶微积分的Kelvin_Voigt流变模型_郭佳奇
基于分数阶微积分的Kelvin_Voigt流变模型_郭佳奇

分数阶微分方程-课件

分数阶微分方程 第三讲分数阶微分方程基本理论 一、分数阶微分方程的出现背景及研究现状 1、出现背景 分数阶微积分是关于任意阶微分和积分的理论,它与整数阶微积分是统一的,是整数阶微积分的推广。 整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受,很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题,其无论在理论分析还是数值求解方面都已有较完善的理论。但当人们进入到复杂系统和复杂现象的研究时,经典整数阶微积分方程对这些系统的描述将遇到以下问题: (1)需要构造非线性方程,并引入一些人为的经验参数和与实际不符的假设条件; (2)因材料或外界条件的微小改变就需要构造新的模型; (3)这些非线性模型无论是理论求解还是数值求解都非常繁琐。 基于以上原因,人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程,其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势,因而成为复杂力学与物理过程数学建模的重要工具之一。 2、研究现状 在近三个世纪里,对分数阶微积分理论的研究主要在数学的纯理论领域里进行,似乎它只对数学家们有用。然而在近几十年来,分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。分数阶微积分理论也受到越来越多的国内外学者的广泛关注,特别是从实际问题抽象出来的分数阶微分方程成为很多数学工作者的研究热点。随着分数阶微分方程在越来越多的科学领域里出现,无论对分数阶微分方程的理论分析还是数值计算的研究都显得尤为迫切。然而由于分数阶微分是拟微分算子,它的保记忆性(非局部性)对现实问题进行了优美刻画的同时,也给我们的分析和计算造成很大困难。 在理论研究方面,几乎所有结果全都假定了满足李氏条件,而且证明方法也和经典微积分方程一样,换句话说,这些工作基本上可以说只是经典微积分方程理论的一个延拓。对分数阶微分方程的定性分析很少有系统性的结果,大多只是给出了一些非常特殊的方程的求解,且常用的求解方法都是具有局限性的。 在数值求解方面,现有分数阶方程数值算法还很不成熟,主要表现为: (1)在数值计算中一些挑战性难题仍未得到彻底解决,如长时间历程的计算和大空间域的计算等; (2)成熟的数值算法比较少,现在研究较多的算法主要集中在有限差分方法与有限单元法; (3)未形成成熟的数值计算软件,严重滞后于应用的需要。

数学建模-微积分模型

第四章 微积分模型 今天人们不论从事什么活动都讲究高效益,即希望所采取的策略使某个或某些指标达到最优。商店订货要使订货、存贮等费用最小,体育比赛运动员要创造最好的成绩,工程设计要追求最佳方案。普遍存在的优化问题经常成为人们研究的对象,建立这类问题的模型,我们称为优化模型。 建立优化模型首先要确定所关心的优化指标的数量描述,然后构造包括这个指标及各种限制条件的模型,通过模型求解给出达到优化指标的所谓策略。本章仅考虑定常情况(即所给的策略不随时间改变)。 4.1 不允许缺货模型 某配送中心为所属的几个超市送配某种小电器,假设超市每天对这种小电器的需求量是稳定的,订货费与每个产品每天的存贮费都是常数。如果超市对这种小家电的需求是不可缺货的,试制定最优的存贮策略(即多长时间订一次货,一次订多少货)。 如果日需求量价值100元,一次订货费用为 5000元,每件电器每天的贮存费1元,请给出最 优结果。 模型假设: (1)每天的需求量为常数r ; (2)每次的订货费用为c 1,每天每件产品的存贮费为c 2 ; (3)T 天订一次货,每次订Q 件,且当存贮量 为0时,立即补充,补充是瞬时完成的; (4)为方便起见,将r ,Q 都视为连续量。 模型建立 将存贮量表示为时间的函数(),0q t t =时,进货Q 件这类小电器,储存量(0),()q Q q t =以需求r 的速率递减,直到q (T )=0。 易见 Q=rT (4.1) 一个周期的存贮费用 C 2= A c ds s q T 20 )(=? 一个周期的总费用 C =2 2 21rT c c + 每天平均费用

数学模型课后答案

数学模型课后答案

《数学模型》作业答案 第二章(1)(2012年12月21日) 1.学校共1000名学生,235人住在A宿舍,333人住在B宿舍,432人住在C宿舍.学生们要组织一个10人的委员会,试用下列办法分配各宿舍的委员数: (1). 按比例分配取整数的名额后,剩下的名额按惯例分给小数部分较大者; (2). §1中的Q值方法; (3).d’Hondt方法:将A、B、C各宿舍的人数用正整数n=1,2,3,……相除,其商数如下表:

将所得商数从大到小取前10个(10为席位数),在数字下标以横线,表中A 、B 、C 行有横线的数分别为2,3,5,这就是3个宿舍分配的席位.你能解释这种方法的道理吗? 如果委员会从10个人增至15人,用以上3种方法再分配名额,将3种方法两次分配的结果列表比较. 解:先考虑N=10的分配方案, , 432 ,333 ,235321 ===p p p ∑==3 1 . 1000i i p 方法一(按比例分配) , 35.23 1 11 == ∑=i i p N p q , 33.33 1 22 == ∑=i i p N p q 32 .43 1 33 == ∑=i i p N p q 分配结果为: 4 ,3 ,3321 ===n n n 方法二(Q 值方法) 9个席位的分配结果(可用按比例分 配)为: 4 ,3 ,2321===n n n 第10个席位:计算Q 值为

2. 试用微积分方法,建立录像带记数器读数n 与转过时间的数学模型. 解: 设录像带记数器读数为n 时,录像带转过时间为t.其模型的假设见课本. 考虑t 到t t ?+时间内录像带缠绕在右轮盘上的长度,可得,2)(kdn wkn r vdt π+=两边积分,得 ??+=n t dn wkn r k vdt 0 )(2π ) 2 2 2 n wk k(r n πvt +=∴ . 2 2 2n v k w n v rk t ππ+=∴ 《数学模型》作业解答 第三章1(2008年10月14日) 1. 在 3.1节存贮模型的总费用中增加购买货物本身的费用,重新确定最优订货周期和订货批量.证明在不允许缺货模型中结果与原来的一样,而在允许缺货模型中最优订货周期和订货批量都比原来结果减少.

清华大学微积分习题(有答案版)

第十二周习题课 一.关于积分的不等式 1. 离散变量的不等式 (1) Jensen 不等式:设 )(x f 为],[b a 上的下凸函数,则 1),,,2,1),1,0(],,[1 ==∈?∈?∑=n k k k k n k b a x λλΛ,有 2),(1 1≥≤??? ??∑∑==n x f x f k n k k k n k k λλ (2) 广义AG 不等式:记x x f ln )(=为),0(+∞上的上凸函数,由Jesen 不等式可得 1),,,2,1),1,0(,01 ==∈?>∑=n k k k k n k x λλΛ,有 ∑==≤∏n k k k k n k x x k 1 1 λλ 当),2,1(1 n k n k Λ==λ时,就是AG 不等式。 (3) Young 不等式:由(2)可得 设111,1,,0,=+>>q p q p y x ,q y p x y x q p +≤1 1 。 (4) Holder 不等式:设11 1, 1,),,,2,1(0,=+>=≥q p q p n k y x k k Λ,则有 q n k q k p n k p k n k k k y x y x 111 11?? ? ????? ??≤∑∑∑=== 在(3)中,令∑∑======n k q k n k p k p k p k y Y x X Y y y X x x 1 1,,,即可。 (5) Schwarz 不等式: 2 1122 1 121?? ? ????? ??≤∑∑∑===n k k n k k n k k k y x y x 。 (6) Minkowski 不等式:设1),,,2,1(0,>=≥p n k y x k k Λ,则有 ()p n k p k p n k p k p n k p k k y x y x 11111 1?? ? ??+??? ??≤??????+∑∑∑=== 证明: ()()() () () ∑∑∑∑=-=-=-=+++=+?+=+n k p k k k n k p k k k n k p k k k k n k p k k y x y y x x y x y x y x 1 1 1 1 1 1 1

分数阶微积分发展现状及展望教学文稿

分数阶微积分发展现状及展望 在数学领域中,大体分为五种研究方向:基础数学,应用数学,计 算数学,概率论与数理统计,统计学与控制论。这五个方向对数学在当 代的发展都有不可或缺的作用。从研究内容来讲,方程、算子、群论、 图论、代数、几何等等都是数学领域重要的研究对象。作为基础数学专 业分数阶微分方程方向的博士生,本文将从分数阶微分方程的发展的历 史及现状、本人对分数阶微分方程未来发展的看法来介绍分数阶微分的 基本知识。 (一)、发展历史及现状 牛顿和莱布尼兹发明的微积分是现代数学与古典数学的分水岭。分数阶微积分是关于任意阶微分和积分的理论,它与整数阶微积分是统一的,是整数阶微积分的推广。整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受,很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题,其无论在理论分析还是数值求解方面都已有了比较完善的理论。但当人们进入到复杂系统和复杂现象的研究时,经典整数阶微积分方程对这些系统的描述将遇到一些问题,如:需要构造非线性方程,并引入一些人为的经验参数和与实际不符的假设条件;因材料或外界条件的微小改变就需要构造新的模型等等。基于以上原因,人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程,其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势,因而成为复杂力学与物理过程数学建模的重要工具之一。 对大多数研究人员和工程师而言,分数阶微积分也许还是比较陌生的,但它实际上早在300多年前就被提出。1695年9月,洛必达 (L’Hospital)在给莱布尼兹的著名信件中就写到“对于简单的线性函数 f(x)=x,如果函数导数次数为分数而不是整数那会怎样”。这是公认的第一次提及分数阶微分。1832年,刘维尔(Liouville)成功的应用了自己提出的分数阶导数的定义,解决了势理论问题。之后刘维尔发表的一系列文 章使他成为分数阶微积分理论的实际级创始人。1974年,Oldham与Spanier出版了第一本关于分数阶微积分理论的专著。 在近三个世纪里,对分数阶微积分理论的研究主要在数学的纯理论领域里进行,但是从近几十年,分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。分数阶微积分理论也受到越来越多的国内

数学模型课后答案

《数学模型》作业答案 第二章(1)(2012年12月21日) 1. 学校共1000名学生,235人住在A 宿舍,333人住在B 宿舍,432人住在C 宿舍.学生们 要组织一个10人的委员会,试用下列办法分配各宿舍的委员数: (1). 按比例分配取整数的名额后,剩下的名额按惯例分给小数部分较大者; (2). §1中的Q 值方法; (3).d ’Hondt 方法:将A 、B 、C 各宿舍的人数用正整数n=1,2,3,……相除,其商数如下表: 将所得商数从大到小取前10个(10为席位数),在数字下标以横线,表中A 、B 、C 行有横线的数分别为2,3,5,这就是3个宿舍分配的席位.你能解释这种方法的道理吗? 如果委员会从10个人增至15人,用以上3种方法再分配名额,将3种方法两次分配的结果列表比较. 解:先考虑N=10的分配方案, ,432 ,333 ,235321===p p p ∑==3 1 .1000i i p 方法一(按比例分配) ,35.23 1 11== ∑=i i p N p q ,33.33 1 22== ∑=i i p N p q 32.43 1 33== ∑=i i p N p q 分配结果为: 4 ,3 ,3321===n n n 方法二(Q 值方法) 9个席位的分配结果(可用按比例分配)为: 4 ,3 ,2321===n n n

第10个席位:计算Q 值为 ,17.92043223521=?=Q ,75.92404333322=?=Q 2.9331544322 3=?=Q 3Q 最大,第10个席位应给C.分配结果为 5 ,3 ,2321===n n n 方法三(d ’Hondt 方法) 此方法的分配结果为:5 ,3 ,2321===n n n 此方法的道理是:记i p 和i n 为各宿舍的人数和席位(i=1,2,3代表A 、B 、C 宿舍). i i n p 是每席位代表的人数,取,,2,1 =i n 从而得到的i i n p 中选较大者,可使对所有的,i i i n p 尽量接近. 再考虑15=N 的分配方案,类似地可得名额分配结果.现将3种方法两次分配的结果列表如下: 2. 试用微积分方法,建立录像带记数器读数n 与转过时间的数学模型. 解: 设录像带记数器读数为n 时,录像带转过时间为t.其模型的假设见课本. 考虑t 到t t ?+时间内录像带缠绕在右轮盘上的长度,可得,2)(kdn wkn r vdt π+=两边积分,得 ?? +=n t dn wkn r k vdt 0 )(2π )22 2 n wk k(r n πvt +=∴ .2 22n v k w n v rk t ππ+=∴ 《数学模型》作业解答 第三章1(2008年10月14日)

微积分习题讲解与答案

习题8.1 1?指出下列微分方程的阶数,并指出哪些方程是线性微分方程: (3) x 2 y 4y (sin x)y = 0 ⑷^P p= sin 2 r d6 解(1)1阶非线性 (2) 1阶线性 (3) 3阶线性 (4) 1阶线性 2?验证下列函数是否是所给微分方程的解 /八 、亠 sinx (1) xy y = cosx, y = x (2) (4 - x 2)y ' xy = 2x,y = 2 ? C" - x 2 (C 为任意常数) (3) y 2y : y = 0, y 二 Ce x (C 为任意常数) (4) y" — (X , + 丸2 )y ' +餌丸2 y = 0, y = C 4e" + C 2e'2 x (C 1 ? 为任意常数) (5) (x -2y)y" =2x - y, x 2 - xy ? y 2 =C (C 为任意常数) (6) (xy -x)y xy 2 yy 1 -2y = 0, y = ln( xy) xcosx — sinx sin x 亠 解⑴是,左=x 2 cosx =右 x x (2) 是,左=(4 — X 2 )-^= + x(2 +C 訥—X 2) = 2x =右 訥-x 2 (3) 是,左=Ce x -2Ce x Ce x =0 =右 (4) 是,左= G :e i x C 2 2e 2 x )-(「-g re 4 x C 2 -e 2 x ) i 2(Se 4 x C 2e?0 =右 2x — y (5) 是,左=(x - 2y) 2x - y 二右 2 ⑴ x(y ) -2yy xy = 0 2 (2) x y - xy y = 0

第二章数学模型与定解问题

第二章数学模型与定解问题 2.1典型方程 三类基本的二阶偏微分方程是: (1)波动方程 0)(2 =++-zz yy xx tt u u u a u (2)热传导方程 0)(=++-zz yy xx t u u u k u (3)拉普拉斯方程 0=++zz yy xx u u u 许多数学物理问题都可归结为解偏微分方程的问题,特别是可归结为解上面所列举的三个偏微分方程的问题.我们将开始研究这些方程,首先仔细考察表示这些物理问题的数学模型. 2.2弦的振动 在数学物理中最重要的问题之一是拉紧的弦的振动问题.由于它较简单, 且经常出现在许多数学物理的分支中,所以在偏微分方程理论中把它作为一个典型的例子. 让我们考察一长为 l 的两端固定的拉紧的弦.我们的问题是要确定弦的运动方程,用它来描述在给定初始扰动后任一时刻t 的弦的位移u(x,t). 为了能.得出一个较简单的方程,我们作下面的一些假设: (1)弦是柔软与有弹性的,即它不能抵抗弯矩,因此在任何时刻弦的张力总是沿着弦的切线方向; (2)弦的每一段都不伸长,因此根据胡克(Hooke)定律,张力是常数; (3)弦的重量与其张力相比很小; (4)弦的偏移与其长度相比很小; (5)位移后的弦在任一点上的斜率与1相比很小; (6)弦只有横振动. 我们考察弦上一微小元素.设T 是如图2.1所示的两端点上的张力.作用在弦的这一微小元素上的垂直方向的力是: αβsin sin T T - 图(Figure )2.1

根据牛顿第二运动定律,合力等于质量乘以加速度.因此 tt su T T ?=-ραβsin sin (2.2.1) 其中ρ是弦的密度,s ?是这一小段位移后的弦的弧长.因为位移后的弦的斜率很小,所以有 x s ?≈? 因为角α和β都很小,所以 ααtan sin ≈, ββtan sin ≈ 于是等式(2.2.1)变成 tt u T x ?=-ραβtan tan (2.2.2) 但是,由微积分学我们知道,在时刻t 有 x x u )(tan ≈α 及 x x x u ?+≈)(tan β 于是等式(2.2.2)可以写成 tt x x x x x u t u u x ρ =-??+])()([1 令x ?趋于零取极限,得 xx tt u a u 2 = (2.2.3) 其中ρ T a = 2 。方程(2.2.3)称为一维波动方程. 如果在弦的每单位长度上有外力F 作用着,方程(2.2.3)具有下列形式: f u a u xx tt +=2 (2.2.4) Where ρ F f = ,而外力可以是压力、重力、阻力以及其他力等 2.3膜的振动 膜振动方程在数学物理的许多问题中出现.在我们导出膜振动方程前,像在弦振动的情形中一样,我们作下列一些简化的假设: (1) 膜是柔软与有弹性的,即它不能抵抗弯矩,因此在任何时刻它的张力 总是在膜的切平面内; (2) 膜的每一块元素都没有伸张变形, 因此根据胡克定律, 张力是常数;

经济类微积分课后习题答案解析

一.教材和大纲(3--6月) 教材往往容易被很多同学忽略,其实教材真的很重要,除非你的基础很好,比如我今年,就只看了一遍陈文灯复习指南,600题都没有做完。 就做了些模拟题。然后就上考场。但是同学们必须知道我08年是怎么复习的。我大概在6月份之前把4本教材和教材的课后习题全部都做了。其实 我想说,很多同学都说自己看了多少书,做了多少题为什么最后考得还是不好。我希望大家能够做到,不是你做了多少题,但是我们做题不能只 做不想,不懂脑。当然做题还是要一定的量,人家政治不也说“量变引起质变”吗。 我想说的是,大家如果有资源的话尽量用起来,有那种数学强人的话,尽量让他们给你们答疑。把你们不会的,全部问清楚,这一点真的很重要。 我男朋友是数学系的,我可以说即使计算不会都问他,因为说不定他们就能说出怎么样计算更简单,更不容易出错。 二。复习指南(7--10月) 其实我觉得复习指南的话,用谁的吧,我不细说,因为每个人的情况不一样,而且基础不同吧。但是还是给个建议吧,如果你的基础还可以的话 ,个人建议用陈文灯的,如果你觉得你的基础一般的话,那还是用李永乐的吧。 我的好多学弟学妹们常问我怎么用复习指南。我个人觉得复习指南吧,一般要看2遍吧。第一遍和第二遍,有一定的笔记差距。我看的时候一般是: 首先,我想说,同学们请你不要看一个题目是怎么做的,而是要你自己去做,因为咱们已经看过一遍教材了,所以我们看书时,把答案先盖住,然后 自己做,做完后看和答案有什么差距,然后调整下自己的思维,希望你在第二次或第三次的时候能会。 第一遍:如果这个题基本不怎么会的话,就用红色笔打上大大的问号,以便第二次的时候可以重点看看。如果是计算错误的话,还是用蓝色的笔标记吧。 也许很多同学都觉得我方法都对了,计算是小问题。那我告诉你,你错了。像我09年数学考134,就是因为忽略了计算。说实话,一般来说, 130和150的区别也许就是谁细心了,实力差距个人觉得不是很大,所以希望同学们不要忽略计算问题。 第二遍:其实做题还是和第一遍一样,盖住答案,多注意下第一遍画红色的部分。蓝色笔的部分,希望大家不要再计算错了。 三。600题和模拟题(11--12月) 希望大家买的600题是那种答案和题目分开比较远那种,不要前面是题,下面就是答案,这样的书不便于同学们去发现自己的弱点。 咱们怎么用这个600题呢,首先,咱们每天规定做30题吧,但是不是连续20天都做题。这里有个建议必须说一下,希望同学们,在做600题的时候, 不要再去翻复习指南了。如果你不会,说明这就是你的弱点了,你是不是该好好地补习下这部分呢。比如说,我先做的60题,发现我自己对间断点的类型 不是很清楚。咱们不会,没有关系,我用红笔在这页的上面写上,间断点的类型。说明这是你的弱点,然后你自己在第二天再看看,做点别的练习,然后 再继续600题。 其实是模拟题。我一般都是采取考试的形式来要求自己,我自己对自己的要求比较高,我

分数阶微分方程_课件

分数阶微分方程 一、 预备知识 1、 分数阶微积分经典定义回顾 作为分数阶微积分方程的基础,本书在第二章中对分数阶微积分的定义及性质做了系统的介绍,为了接下来讨论的需要,我们首先对其进行一个简要的回顾。 (1)分数阶微积分的主要思想 如上图所示,分数阶微积分的主要思想是推广经典的整数阶微积分,从而将微积分的概念延拓到整个实数轴,甚至是整个复平面。但由于延拓的方法多种多样,因而根据不同的需求人们给出了分数阶微积分的不同定义方式。然而这些定义方式不仅只能针对某些特定条件下的函数给出,而且只能满足人们的某些特定需求,迄今为止,人们仍然没能给出分数阶微积分的一个统一的定义, 这对分数阶微积分的研究与应用造成了一定的困难。 1、分数阶微分的定义 为了满足实际需要,下面我们试图从形式上对分数阶微积分给出一种统一的表达式。 分数阶微积分的主要思想是推广经典的累次微积分,所有推广方法的共同目标是以非整数参数p 取代经典微积分符号中的整数参数n ,实际上,任意的n 阶微分都可以看成是一列一阶微分的叠加: ()()n n n d f t d d d f t dt dt dt dt = (1) 由此,我们可以给出一种在很多实际应用中十分重要的分数阶微积分的推广方 式。首先,我们假设已有一种合适的推广方式来将一阶微分推广为α(01α≤≤) 阶微分,即d D dt α→是可实现的。那么类似地可得到(1)的推广式为: ()()n n D f t D D D f t αααα= (2) 这种推广方式最初是由..K S Miller 和.B Ross 提出来的,其中D α采用的是R L -分数阶微分定义,他们称之为序列分数阶微分。序列分数阶微分的其他形式可以通过将D α替换为G L -分数阶微分、Caputo 分数阶微分或其他任意形式

分数阶微积分发展现状及展望讲课讲稿

精品文档 分数阶微积分发展现状及展望 在数学领域中,大体分为五种研究方向:基础数学,应用数学,计算数学,概率论与数理统计,统计学与控制论。这五个方向对数学在当代的发展都有不可或缺的作用。从研究内容来讲,方程、算子、群论、图论、代数、几何等等都是数学领域重要的研究对象。作为基础数学专业分数阶微分方程方向的博士生,本文将从分数阶微分方程的发展的历史及现状、本人对分数阶微分方程未来发展的看法来介绍分数阶微分的基本知识。 (一)、发展历史及现状 牛顿和莱布尼兹发明的微积分是现代数学与古典数学的分水岭。分数阶微积分是关于任意阶微分和积分的理论,它与整数阶微积分是统一的,是整数阶微积分的推广。整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受,很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题,其无论在理论分析还是数值求解方面都已有了比较完善的理论。但当人们进入到复杂系统和复杂现象的研究时,经典整数阶微积分方程对这些系统的描述将遇到一些问题,如:需要构造非线性方程,并引入一些人为的经验参数和与实际不符的假设条件;因材料或外界条件的微小改变就需要构造新的模型等等。基于以上原因,人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程,其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势,因而成为复杂力学与物理过程数学建模的重要工具之一。 对大多数研究人员和工程师而言,分数阶微积分也许还是比较陌生的,但它实际上早在300多年前就被提出。1695年9月,洛必达(L’Hospital)在给莱布尼兹的著名信件中就写到“对于简单的线性函数f(x)=x,如果函数导数次数为分数而不是整数那会怎样”。这是公认的第一次提及分数阶微分。1832年,刘维尔(Liouville)成功的应用了自己提出的分数阶导数的定义,解决了势理论问题。之后刘维尔发表的一系列文章使他成为分数阶微积分 精品文档

微积分模型

第一篇 微积分模型 在微积分部分的应用实例中,通过对应用问题建模主要培养应用极限、连续、相对变化率、微元、无穷级数、最优化和微分与差分方程等思想解决实际应用问题的能力。 函数的性质包括分段性质、单调性、奇偶性等,由函数的基本性质可以产生对函数进行分类的方法。与函数基本特性相关的应用实例有:市话费是降了还是升了,外币兑换与股票交易中的涨跌停板,库存问题与库存曲线,“另类”的常量函数,蠓虫分类的初等数学模型,核军备竞赛问题等。 数列与函数的极限和函数连续性质是处理变量变化过程的工具,应用重要极限计算连续复利利率的计算,应用函数的连续性和介值定理解决特殊的应用问题。与极限和连续等内容相关的应用实例有:从科赫雪花谈起,复利、连续复利与贴现,出售相同产品的公司为什么喜欢扎堆,椅子为什么能放稳等。 导数、微分是函数的相对变化的极限过程,函数的特性和极值理论可以解决经济管理中的实际应用问题,导数、微分在经济管理中的应用反映为边际、弹性等。相关的应用实例有:影子为什么那么长,边际是什么?弹性是什么?商家应该怎样制定自己的价格策略?不同消费群体的需求弹性问题,机械与人工的调配问题,易拉罐的形状,这批酒什么时候出售最好,该不该接受供货商的优惠条件,作者与出版商的利益冲突等。 微元分析是微积分中一种重要的分析方法,特别是函数的连续求和归结为该函数的积分。与积分和微元分析内容相关的应用实例有:洛伦兹曲线与基尼系数,均匀货币流的总价值与投资回收期的计算,下雪时间的确定,第二宇宙速度是怎样计算出来的等。 离散变量的求和可以用无穷级数来表达,无穷级数的求和是一个极限过程。与无穷级数内容相关的应用实例有:最大货币供应量的计算,政府支出的乘数效应,运用现值计算进行投资项目的评估,谈谈龟兔赛跑悖论 等。 如果影响研究问题的主要因素有两个或者两个以上,则要用多元函数的微积分学来处理,涉及到多元函数偏导数、偏边际、偏弹性和交叉弹性、条件极值等内容。相关的应用实例有:空调销售量的预测,相互关联商品的需求分析,衣物怎样漂洗最干净,拉格朗日乘数与影子价格等。 变量的变化过程可以用微分方程或差分方程来描述,通过对微分方程或差分方程的建立与求解,可以研究变量的形态和变化规律。与微分方程和差分方程相关的应用实例有:人口模型,单种群动物模型,相对封闭环境中的传染病模型,江河污染物的降解系数,怎样计算固定资产的折旧,放射性元素衰变模型,市场上的商品价格是怎样波动的,再谈下雪时间的确定,溶液浓度模型,饲养物的最佳销售时机,信贷消费中每月还款金额的确定,资源的合理开发与利用,从诺贝尔奖谈起,蛛网模型,梵塔问题,平面内直线交点的个数,菲波那契数列的通项公式等。 1

微积分课后题答案高等教育出版社

习 题 六 (A ) 1.根据定积分的几何意义说明下列各式的正确性 (1)0d cos 20=? x x π (2) x x x x d )1(2d )1(22 22 2+=+? ? - (3) 0d 3 1 1 =?-x x (3)x x dx x d 4 21 1 1 ?? == 解:(1)该定积分的几何意义如右图所示阴影部分面积的代数和,由对称性可知正确. (2)该定积分的几何意义如右图所示阴影部分面积的代数和,且在) 2 , 2(-范围内对称,所以是正确的. (3)该定积分的几何意义如右图所示阴影部分面积的代数和,且关于原点对称,所以正确. (4)原式dx x ? -=1 1 2 等式左边的定积分的几何意义是右边图形阴影部分面积的代数和的2倍,且又因为阴影部分在1) , 1(-范围内关于轴对称,所以等式两边相等. 2.不计算积分,比较下列积分值的大小 (1)x x d 2 10?与x x d 3 10? (2)x x d 2 3 1?与x x d 3 3 1 ? (2)x x d ln 4 3 ? 与 x x d )(ln 2 4 3 ? (4)x x d sin 2 ? π 与 x x d 2 ? π 解:(1)由定积分的比较性可知在1) , 0(范围内32x x >,所以前者大于后者. (2)由定积分的比较性可知在3) , 1(范围内32x x <,所以前者小于后者. (3)由定积分的比较性可知在4) , 3(范围内2)(ln ln x x <,所以前者小于后者.1=a (4)由定积分的比较性可知在)2 , 0(π 范围x x

分数阶微积分的性质

分数阶微积分的性质 根据上述三种分数阶微积分的定义,可以得到分数阶微积分一些性质如下 [66] : (1) 记忆属性。当t 在时刻时,函数()f t 的分数阶微分值由初始时刻到t 时 刻的所有时刻的函数值取值。 (2) 当1a t D β算子的1β是整数时,整数阶微积分和分数阶微积分二者为等 同关系,1β为任意阶时,整数阶微积分被包含在分数阶微积分内。 (3) 分数阶微积分算子1a t D β是线性的,符合线性系统中的齐次特性和迭 加特性,即对任意常数,a b 均满足: 1110 00[()()]()()t t t D af t bg t a D f t b D g t βββ+=+ (4) 解析函数()f t 分数阶导数10()t D f t β对t 和a 都是可以解析的。 2.4 分数阶系统的模型描述 实际生活中,大多数的对象的内在特性都能通过整数阶微分方程的形式来表征,比如物理特性、化学特性等。但往往存在一些特别的对象其特性无法靠整数阶微分方程精确表征,但分数阶次的微分方程刚好能考虑到整数阶次微分方程所忽略的特性,所以,用分数阶微分方程描述的系统,其内在特性反应更真实、更全面。 一个典型的单输入单输出分数阶线性系统的微分方程可用如下形式来表示: 3123 1 2 123122()()()()=()()()() m n m n a D y t a D y t a D y t a D y t b D u t b D u t b D u t b D u t ααααββββ+++++++ + (2.10) 其中,(1,2, ,),(1,2, ,)i j a i m b j n ==分别表示输出和输入相应的系数, 12m ααα<<<,12n βββ<<<分别表示输出和输入分数阶的阶次,()() u t y t 、分别表示系统的输入和输出。 结合前面的式(2.6)和式(2.10)对系统进行拉普拉斯变换,得到系统的传递函数模型为: 121 2 1212()n m n m b s b s b s G s a s a s a s βββααα++=+++ (2.11) 若(1,2, ,)i i i m αα==,(1,2, ,)i i i n ββ==,该系统可称为“同源次”分 数阶系统,则上式进一步可表示为: 11 ()n j j j m i i i b s G s a s β α === ∑∑ (2.12)

数学建模微积分模型

第四章 微积分模型 今天人们不论从事什么活动都讲究高效益,即希望所采取的策略使某个或某些指标达到最优。商店订货要使订货、存贮等费用最小,体育比赛运动员要创造最好的成绩,工程设计要追求最佳方案。普遍存在的优化问题经常成为人们研究的对象,建立这类问题的模型,我们称为优化模型。 建立优化模型首先要确定所关心的优化指标的数量描述,然后构造包括这个指标及各种限制条件的模型,通过模型求解给出达到优化指标的所谓策略。本章仅考虑定常情况(即所给的策略不随时间改变)。 4.1 不允许缺货模型 某配送中心为所属的几个超市送配某种小电器,假设超市每天对这种小电器的需求量是稳定的,订货费与每个产品每天的存贮费都是常数。如果超市对这种小家电的需求是不可缺货的,试制定最优的存贮策略(即多长时间订一次货,一次订多少货)。 如果日需求量价值100元,一次订货费用为 5000元,每件电器每天的贮存费1元,请给出最 优结果。 模型假设: (1)每天的需求量为常数r ; (2)每次的订货费用为c 1,每天每件产品的存贮费为c 2 ; (3)T 天订一次货,每次订Q 件,且当存贮量 为0时,立即补充,补充是瞬时完成的; (4)为方便起见,将r ,Q 都视为连续量。 模型建立 将存贮量表示为时间的函数(),0q t t =时,进货Q 件这类小电器,储存量(0),()q Q q t =以需求r 的速率递减,直到q (T )=0。 易见 Q=rT (4.1) 一个周期的存贮费用 C 2= A c ds s q T 20 )(=? 一个周期的总费用 C =2 2 21rT c c + 每天平均费用 2 )(21rT c T c T c += (4.2) 模型求解 求T ,使)(T c 取最小值。 由 0=dT dc ,得 2 12 1 2,2c r c Q rc c T = = (4.3)

微积分习题讲解与答案

微积分习题讲解与答案Last updated on the afternoon of January 3, 2021

1?指出下列微分方程的阶数,并指出哪些方程是线性微分方程: (1)心)2-2少 + 和=0 (2) x2y-xy f + y = 0 (3)x2+ 4y n + (sinx)y = 0 (4) —+ P =sin26 de 解⑴1阶非线性 (2)1阶线性 ⑶3阶线性 (4)1阶线性 2.验证下列函数是否是所给微分方程的解 (1)xy f + y = cosx,j = ----------- x ⑵(l-x2)y f + xy = 2x,y = 2 + C^l-x2 (C 为任意常数) (3)y n-2y, + y=09y=Ce x (C 为任意常数) (4)y n-(^ + A2)j r += 0,J =+C2e^x (CiQ为任意常数) ⑸(x - 2y)y f = 2x-y9x2-xy + y2=C (C 为任意常数) (6) (xy-x)y H + xy f2 + = = ln(xj) xcosx-sinx sinx 解(1)是, = cosx 二右 左二X -- ;---- + X X (2)是,^=(l-x2) t X +x(2 + Cyll-x2 ) = 2x=右y/l-x2 ⑶是‘左=Ce x— 2Ce x +Ce x =0=右(4)是,左二

=右 2x — V ⑸是,左*-2刃口^2一尸右 ,+兀-^ +〉,亠_2亠 (xj-x) (xy-x) xy — x xy- x 二比二'尹+亠厶+(宀2丿)5—)“ (xy-x) (xy-x) (xy-x) 二右 3 ?求下列微分方程的解 (3) (l + y)dx —(1 一 y)dy = 0 (2) | = Jcosxdx,j r = 5111^ + ^ ⑶圧^訂张j%严峡皿 即一y + 21nll + y l=x + C ⑷估心侖必 解得 ln(l + j 2) = ln(l + x 2) + C^ 4?已知曲线y = f(x)经过原点,并且它在点(X 』)处的切线的斜率等于2,, 试求这条 曲线的方程。 解已知y f = 2x 2 ⑹是,左 ⑴加2; ⑵ 4-0SX ; (be

微积分课后题答案习题详解

第二章 习题2-1 1. 试利用本节定义5后面的注(3)证明:若lim n →∞ x n =a ,则对任何自然数k ,有lim n →∞ x n +k =a . 证:由lim n n x a →∞ =,知0ε?>,1N ?,当1n N >时,有 取1N N k =-,有0ε?>,N ?,设n N >时(此时1n k N +>)有 由数列极限的定义得 lim n k x x a +→∞ =. 2. 试利用不等式A B A B -≤-说明:若lim n →∞ x n =a ,则lim n →∞ ∣x n ∣=|a|.考察数列x n =(-1)n ,说明上 述结论反之不成立. 证: 而 n n x a x a -≤- 于是0ε?>,,使当时,有N n N ?> n n x a x a ε-≤-< 即 n x a ε-< 由数列极限的定义得 lim n n x a →∞ = 考察数列 (1)n n x =-,知lim n n x →∞不存在,而1n x =,lim 1n n x →∞ =, 所以前面所证结论反之不成立。 3. 利用夹逼定理证明: (1) lim n →∞2221 11(1)(2)n n n ??+++ ?+?? L =0; (2) lim n →∞2!n n =0. 证:(1)因为 222222111112(1)(2)n n n n n n n n n n ++≤+++≤≤=+L 而且 21lim 0n n →∞=,2lim 0n n →∞=, 所以由夹逼定理,得 222111lim 0(1)(2)n n n n →∞?? +++= ?+? ?L . (2)因为22222240!1231n n n n n < =<-g g g L g g ,而且4lim 0n n →∞=, 所以,由夹逼定理得

微积分课后习题答案

微积分课后习题答案 Document serial number【KKGB-LBS98YT-BS8CB-BSUT-BST108】

微积分第八章课后习题答案 习题8-1 1.(1)一阶;(2)二阶;(3)一阶;(4)三阶;(5)三阶;(6)一阶;(7)二阶;(8)一阶。 2.(1)、(2)、(3)、(4)、(5)都是微分方程的通解。 3.1 22 y x =+.4.将所给函数及所给函数的导数代人原方程解得: 21 ()(1)2 u x x dx x x C =+=++?. 习题8-2 1.(1)原式化为:ln dy x y y dx = 分离变量得: 11 ln dy dx y y x = 两边积分得:11 ln dy dx y y x =?? 计算得:()11ln ln d y dx y x =? ? 即:()1ln ln ln y x C =+ 整理:1ln y C x = 所以:原微分方程的通解为:Cx y e =; (2)原式化为:()()2211y x dy x y dx -=-- 分离变量得: ()()22 11y x dy dx y x -=-- 两边积分得:()() 22 11y x dy dx y x -=--?? 计算得: ()()() ()22 22 1111112211d y d x y x -=----?? 即:()()221ln 1ln 1y x C -=--+ 整理:22(1)(1)y x C --=

所以:原微分方程的通解为:22(1)(1)y x C --=; (3xydx =- 分离变量得: 1dy y = 两边积分得: 1dy y =? 计算得:()2 1ln 12y x = - 即:1ln y C = 整理:y = 所以:原微分方程的通解为:y = (4) 1y e Cx -=-; (5)sin 1y C x =-; (6)1010x y C -+=; (7)22ln 22arctan y y x x C -=-+; (8)当sin 02y ≠时,通解为ln |tan |2sin 42y y C =-;当sin 02 y =时,特解为2(0,1,2,)y k k π==±±; (9)222ln x y x C +-=; (10)22ln ln x y C +=。 2.(1)tan 2 x y e =;(2)(1)sec x e y +=;(3)2(1)22y x e y +-=;(4) 1 ln |1|1a x a y =--+;(5)24x y =;(6)323223235y y x x +--=;(7)sin y x =; (8)cos 0x y =。

分数阶微积分的定义

分数阶微积分的定义 分数阶微积分的研究对象是分数阶微分和分数阶积分,分数阶微积分定义是整合和统一分数阶微分和分数阶积分得到的。首先介绍常用的三种分数阶微分定义,具体为: (1)Grünwald -Letnikov 分数阶微分定义 若()f t 函数在区间[,]a t 存在1m +阶连续导数,当0α>时, m 至少取到[]α,则其次数为(1)m m αα≤<+的分数阶微分定义为: [()/] ()lim ()t a h a t i h i D f t h f t ih αα αω--→==-∑ (2.1) 其中,α表示阶次,h 为采样步长,a 表示初始时间,[]表示取整, = (-1)i i i ααω?? ??? 是多项式系数,(1)(2)(1) = ! i i i ααααα??---+ ??? ,我们可以用以下 递推公式直接求出该系数: 01+11,1,1,2,...,i i i n i α αααωωω-??==-= ??? (2.2) 进一步对式(2.1)求极限,可得到其详细定义: 0,0 ()lim ()()()1 ()()(1)(1)a t h nh t a i i m t m a i D f t h f t ih i f a t a t f d i i α α αααξξξαα-→=--+-=?? =- ??? -=+-Γ-++Γ-+∑? (2.3) 其中,()Γ?为欧拉gamma 函数,10 ()t z z e t dt ∞--Γ=?,当R α∈,上述定义也称为Grünwald -Letnikov 分数阶微积分定义。 若:()=0i f t ,,q p R ∈,则微分算子D 满足式(2.4): (2.4) (2)Riemann -Liouville 分数阶微分定义 对于1,m m m N α-<<∈,有 11() ()()()m t a t m m a d f D f t d m dt t α αττατ-+= Γ--? (2.5) 其中,当R α∈,上述定义也称为Riemann -Liouville 分数阶微积分定义。 通常情况下,为了方便使用Riemann -Liouville 分数阶微积分定义,要对其取拉普拉斯变换,假设()F s 表示()f s 的原函数,则式(2.5)经过拉普拉斯变换 +(())()q p q p a t a t a t D D f t D f t =

相关文档
最新文档