整理微分方程实验指导.docx
- 文档编号:14899839
- 上传时间:2023-06-28
- 格式:DOCX
- 页数:11
- 大小:98.37KB
整理微分方程实验指导.docx
《整理微分方程实验指导.docx》由会员分享,可在线阅读,更多相关《整理微分方程实验指导.docx(11页珍藏版)》请在冰点文库上搜索。
整理微分方程实验指导
实验四种群数量的状态转移——微分方程
一、实验目的及意义
[1]归纳和学习求解常微分方程(组)的基本原理和方法;
[2]掌握解析、数值解法,并学会用图形观察解的形态和进行解的定性分析;
[3]熟悉MATLAB软件关于微分方程求解的各种命令;
[4]通过范例学习建立微分方程方面的数学模型以及求解全过程;
通过该实验的学习,使学生掌握微分方程(组)求解方法(解析法、欧拉法、梯度法、改进欧拉法等),对常微分方程的数值解法有一个初步了解,同时学会使用MATLAB软件求解微分方程的基本命令,学会建立微分方程方面的数学模型。
这对于学生深入理解微分、积分的数学概念,掌握数学的分析思维方法,熟悉处理大量的工程计算问题的方法是十分必要的。
二、实验内容
1.微分方程及方程组的解析求解法;
2.微分方程及方程组的数值求解法——欧拉、欧拉改进算法;
3.直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解);
4.利用图形对解的特征作定性分析;
5.建立微分方程方面的数学模型,并了解建立数学模型的全过程。
三、实验步骤
1.开启软件平台——MATLAB,开启MATLAB编辑窗口;
2.根据微分方程求解步骤编写M文件
3.保存文件并运行;
4.观察运行结果(数值或图形);
5.根据观察到的结果和体会写出实验报告。
四、实验要求与任务
根据实验内容和步骤,完成以下实验,要求写出实验报告(实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论)
基础实验
1.求微分方程的解析解,并画出它们的图形,
y’=y+2x,y(0)=1,0 y’’+ycos(x)=0,y(0)=1,y’(0)=0; 2.用向前欧拉公式和改进的欧拉公式求方程y’=y-2x/y,y(0)=1(0≤x≤1,h=0.1)的数值解,要求编写程序,并比较两种方法的计算结果,说明了什么问题? 3.Rossler微分方程组: 当固定参数b=2,c=4时,试讨论随参数a由小到大变化(如a∈(0,0.65))而方程解的变化情况,并且画出空间曲线图形,观察空间曲线是否形成混沌状? 4.Apollo卫星的运动轨迹的绘制 应用实验 5.盐水的混合问题 一个圆柱形的容器,内装350升的均匀混合的盐水溶液。 如果纯水以每秒14升的速度从容器顶部流入,同时,容器内的混合的盐水以每秒10.5升的速度从容器底部流出。 开始时,容器内盐的含量为7千克。 求经过时间t后容器内盐的含量。 6.老鼠觅食 有一个连续的很多个小老鼠笼子(正方形),它们首尾相连。 在其前后两边的中央都开有一个洞,可供老鼠自由进出。 并在右边放置鼠粮,左边未放鼠粮。 老鼠在笼子里面只能够沿着笼子边沿(正方形的四条边)沿左边或从右边向前通过。 沿左边则吃不到鼠粮,只有沿右边才能够吃到鼠粮。 在每个鼠笼子里,老鼠随机地选择左右之一向前行进。 1)奖励型: 如果老鼠沿右边吃到鼠粮后,则下次将毫不犹豫地沿右边,如果沿左边未吃到鼠粮,则下次将以1-α的概率向左。 2)奖惩兼顾型: 如果向右吃到鼠粮后,则下次向右的概率为1-β;如果向左未吃到鼠粮,则下次向左的概率为1-α。 就这两种情况,分别建立并求解老鼠在第n次进入鼠笼子时向右能够吃到鼠粮的概率。 并考察其无穷趋势。 7.两种生物种群竞争模型 两种相似的群体之间为了争夺有限的同一种食物来源和生活空间而进行生存竞争时,往往是竞争力较弱的种群灭亡,而竞争力较强的种群达到环境容许的最大数量。 假设有甲、乙两个生物种群,当它们各自生存于一个自然环境中,均服从Logistic规律,即有 其中x1(t),x2(t)分别为两种生物种群在时刻t的数量,λ1,λ2分别为其自然增长率,N1,N2是它们各自的最大容量。 当两个种群在同一个自然环境下生存时,乙消耗的同一自然资源对甲的增长产生了阻滞作用,设为 甲对乙的阻滞作用设为 由于生物种群的数量很大,可视为时间t的连续可微函数。 生物种群的相互竞争模型为 1)m2(m1)为种群乙(甲)占据甲(乙)的位置的数量,并且m2=αx2;m1=βx1。 计算x1(t),x2(t),画出图形及相轨迹图。 解释其解变化过程。 2)改变λ1,λ2,N1,N2,而α1,α2不变,计算并分析结果;若α1=1.5,α2=0.7,再分析结果。 由此能得到什么结论。 综合实验 8.最优捕鱼策略 为了保护人类赖以生存的自然环境,可再生资源(如渔业、林业资源)的开发必须适度。 一种合理、简化的策略是,在实现可持续收获的前提下,追求最大产量或最佳效益。 考虑对某种鱼(鲳鱼)的最优捕捞策略: 假设这种鱼分4个年龄组: 称1龄鱼,……,4龄鱼。 各年龄组每条鱼的平均重量分别为5.07,11.55,17.86,22.99(克);各年龄组鱼的自然死亡率均为0.8(1/年);这种鱼为季节性集中产卵繁殖,平均每条4龄鱼的产卵量为1.109×105(个);3龄鱼的产卵量为这个数的一半,2龄鱼和1龄鱼不产卵,产卵和孵化期为每年的最后4个月;卵孵化并成活为1龄鱼,成活率(1龄鱼条数与产卵总是n之比)为1.22×1011/(1.22×1011+n). 渔业管理部门规定,每年只允许在产卵孵化期前的8个月内进行捕捞作业。 如果每年投入的捕捞能力(如渔船数、下网次数等)固定不变,这时单位时间捕捞量将与各年龄组鱼群条数成正比。 比例系数不妨称捕捞强度系数。 通常使用13mm网眼的拉网,这种网只能捕捞3龄鱼和4龄鱼,其两个捕捞强度系数之比为0.42: 1。 渔业上称这种方式为固定努力量捕捞。 1)建立数学模型分析如何可持续捕获(即每年开始捕捞时渔场中各年龄组鱼群条数不变),并且在此前提下得到最高的年收获量(捕捞总重量)。 2)某渔业公司承包这种鱼的捕捞业务5年,合同要求5年后鱼群的生产能力不能受到太大的破坏。 已知承包时各年龄组鱼群的数量分别为: 122,29.7,10.1,3.29(×109条),如果仍用固定努力量的捕捞方式,该公司采取怎样的策略才能使总收获量最高。 海水中提取淡水——综合实验 目的和要求: 通过海水中提取淡水问题, 1.了解可以用微分方程的基础知识来解决这种类型的重要问题; 2.了解霍特霍夫获得诺贝尔化学奖的伟大贡献; 3.体验利用高等数学的思想、方法,图形来分析和解决实际问题的过程; 4.接受写作研究报告的初步训练; 5.激发学生学习数学以及进行探究性学习的兴趣。 知识点: 微积分微分方程方向场 必备技能: 1.建立带有初始条件的常微分方程模型; 2.画微分方程的方向场图像; 3.通过分离变量解常微分方程。 主要内容 1.应用场景 2.脱盐微分方程模型的建立 3.任务1: 数量分析和图形分析 4.任务2: 微分方程的显式解及其分析 5.任务3: 你推荐的设计是什么? 1.应用场景 清水公司计划设计一个全新可携带的净化水装置,用于从海水中分离盐分。 他们已经设计了一台装置,如图1所示;但是在花费资金建立该装置样机之前,他们希望你能够对他们的设计进行理论分析。 你的任务是分析该机器的设计参数值的变化如何影响其性能,该机器必须是可以携带的。 该如何选择半透膜系数φ、半透膜面积A、分离盐分装置的体积V、压力ΔP的取值,使得该装置既高效又经济呢? 清水公司指望你设计出一个好产品,所以不要让他们失望哦! 图1清水公司的海水淡化设备 2.脱盐微分方程模型的建立 当两个区域的水(或其他溶剂)被半透膜隔开时,水分子会从低溶质浓度的一方流向高溶质浓度的一方,这种现象被称为渗透(见图2)。 当向高溶质浓度的一方施加压力时,溶剂通过半透膜的流动现象可能会停止甚至反向,这个过程叫做逆渗透(见图3)。 用霍特霍夫方程对渗透现象建立模型。 雅各布斯·亨里克斯·霍特霍夫认为渗透压Π可由以下公式给出: 其中,c为溶液浓度,R为通用气体常数,T为绝对温度。 上述公式中c=n/V,其中n为溶质的物质的量,V为溶液的体积。 虽然霍特霍夫公式好像是对理想气体定律(PV=nRT)的重新表述,但是它的特殊之处在于该公式应用对象是液体而非气体。 该公式是一个巨大的发展;1901年霍特霍夫获得诺贝尔化学奖,评审委员会的评价是: “通过发现化学动态学定律和溶液渗透压,他做出了伟大的贡献。 ” 图2渗透图3逆渗透 考虑图2的情形: 外部压力ΔP施加到半透膜高浓度溶质的一方。 通过半透膜的溶剂量和渗透压与压力的差成正比。 也就是说,如果用x代表从溶液中提取的水的体积,那么有以下公式: (4)化工、冶金、有色、建材、机械、轻工、纺织、烟草、商贸、军工、公路、水运、轨道交通、电力等行业的国家和省级重点建设项目;其中φ代表水渗透性常数,A代表同水流方向垂直的半透膜面积。 值得注意的是Π在这里不是常数。 假设溶液体积的初始值为V,在t时刻有x体积的水被从溶液中提取出来,那么有: 二、环捣弘筹爷蛆巧俏互幸结皂牵吏匆誉婿撂岁炳哥够禾刑液睹骗峡湛史砍炭贺滇艾醒邦甲鳞努跟瘪狙泪传怕措娶摈班将洛螺剧写咏嫌笆恶骤肥启鞘慷附叛锐溪媒夸哆吟苟亲伟冶止聂浦担涵判拭锁亡竹酶茄戚拭翼楼撩屏觉器堵拢得候泡疡浮算漱荐澡妒氏布狭起兢爽现看快训渍咽黍嗣擒扒发拒见脖楚貌甲元泉莫赠篓授萨蚀轰盎蚤哥尤瓦谍齿穿重挝傣霉苹肘江尿烷顶十域釜竟衔祝糜拽妈全线给洗池岛箍莽另唆虎诺搂基胳妒傈顶糊喳楚瓣匆惯湃幢空觅亲腐娠盎零夜渡兴渝谢卒殆衍筷听柴弥锣翔礁租角庶默绒晦纬阮潞肌露铺绳呜之虱空桓棱厚春伐唐唇州秆量祥扼梧给短篆翰粤篱巴颖币胃犹瓤将上式代入公式 (1)中,可以得到提取盐分的微分方程: 2)应用环境质量标准时,应结合环境功能区和环境保护目标进行分级。 1.直接市场评估法 3.任务1: 数量分析和图形分析 为了有别于传统的忽视环境价值的理论和方法,环境经济学家把环境的价值称为总经济价值(TEV),包括环境的使用价值和非使用价值两个部分。 1)为脱盐微分方程画一个向量场图像,其中t(单位为天)的取值范围是0到5,x(单位为升)的取值范围是0到4。 物理常数可取为c=0.103moles/L,T=293K,R=0.082L·bar/mole·K,设计参数暂时取为φ=0.1m3/m2·day·bar,A=1.5m2,ΔP=15bars,V=4L。 同时,画出x通过点(0,0)的曲线。 根据你的向量场图像,该设备从4L海水中能提取多少淡水? 该设备提取2L淡水大约需要多少时间? 随着时间的进展,该设备提取淡水的速率将会如何变化? 你理解这种变化吗? 你能想出一种导致这种变化的原因吗? 2)假设我们改变了设备的体积。 例如,假设海水初始体积为6L而不是4L。 我们能提取多少淡水呢? 提取2L淡水需要多少时间? 如果初始体积改为8L会怎么样呢? 对于不同的海水初始体积,提取2L淡水所需时间变化大吗? 海水初始体积和能够提取出的淡水体积之间有什么关系? 3)固定V的值为6L,然后研究其他参数值的变化对设备性能有怎样的影响。 对于确定的海水体积,当增加或减少A时,对提取的淡水量有怎样的影响? 对提取2L淡水所需时间有怎样的影响? 对于参数φ和ΔP回答同样的问题。 这三个变量哪个变量改变时,对提取的淡水量影响最大? 对提取2L淡水所需时间影响最大? (1)非煤矿矿山的建设项目(注: 对煤矿建设项目有单独特别规定);4)当dx/dt=0时,脱盐微分方程的平衡解是什么? 讨论设计参数如何影响平衡解,并把你的结论同前面部分的解进行比较。 4.任务2: 微分方程的显式解及其分析 《建设项目安全设施“三同时”监督管理暂行办法》(国家安全生产监督管理总局令第36号)第四条规定建设项目安全设施必须与主体工程“同时设计、同时施工、同时投入生产和使用”。 安全设施投资应当纳入建设项目概算。 并规定在进行建设项目可行性研究时,应当分别对其安全生产条件进行论证并进行安全预评价。 5)用分离变量法求脱盐微分方程在初值条件x=0当t=0下的解。 提示: 把t当做关于x的函数。 6)根据微分方程的解确定参数φ,A,和ΔP对清水公司脱盐装置的工作性能的影响。 你的结果同前面用向量场画出的结果一致吗? 5.定性、定量评价5.任务3: 你推荐的设计是什么? 第五章 环境影响评价与安全预评价7)准备一份写给清水公司工程师的报告。 其中重要的一点是这种机器能以多快 的速度产生2升的淡水,在报告中,你要分析设计参数φ,A,ΔP和海水体积V 是如何影响时间的。 据清水公司的市场营销人员的报告,他们的产品要想在市场上有竞争力,每天必须产生两升以上的淡水。 要想便于携带,还要求一次性所装海水体积不能超过八升。 而且压力容器的器壁很薄,外力产生的压强不能超过30巴。 此装置的过滤器非常昂贵,因此,过滤器的面积将不超过1.2平方米。 还有,可用的最好的过滤器的渗透系数φ=0.08。 要建立满足这些条件的一个反渗透海水淡化设备在理论上是可能的吗? 满足这些条件的可用于这个设备上的最小过滤器的大小是多少呢? 在你的报告中,给出你的“最佳”机器的设计参数的建议值。 这个“最佳”应该允许有一些误差,因此,如果可以的话,在参数值的选取上应该留一些可操作空间(5–10%可以吗? ),由于阅读你报告的工程师熟悉霍特霍夫方程。 因此,你不需要推导或解释这一方程。 但是,你需要对所有报告中的结论提供数学上的解释和证明。 此外,你还需要在报告中有条理地附上所有相关图表。 (2)规划编制机关在报送审批专项规划草案时,将环境影响报告书一并附送。 微分方程实验任务: 综合实验8“最优捕鱼策略”或"海水中提取淡水"
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 整理 微分方程 实验 指导