1、数学实验三实验三 种群数量的状态转移微分方程一、实验目的及意义1 归纳和学习求解常微分方程(组)的基本原理和方法;2 掌握解析、数值解法,并学会用图形观察解的形态和进行解的定性分析;3 熟悉MATLAB软件关于微分方程求解的各种命令;4 通过范例学习建立微分方程方面的数学模型以及求解全过程; 通过该实验的学习,使学生掌握微分方程(组)求解方法(解析法、欧拉法、梯度法、改进欧拉法等),对常微分方程的数值解法有一个初步了解,同时学会使用MATLAB软件求解微分方程的基本命令,学会建立微分方程方面的数学模型。这对于学生深入理解微分、积分的数学概念,掌握数学的分析思维方法,熟悉处理大量的工程计算问题的
2、方法是十分必要的。二、实验内容1 微分方程及方程组的解析求解法;2 微分方程及方程组的数值求解法欧拉、欧拉改进算法;3 直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解);4 利用图形对解的特征作定性分析;5 建立微分方程方面的数学模型,并了解建立数学模型的全过程。三、实验步骤1开启软件平台MATLAB,开启MATLAB编辑窗口; 2根据微分方程求解步骤编写M文件3保存文件并运行;4观察运行结果(数值或图形);5根据观察到的结果和体会写出实验报告。四、实验要求与任务根据实验内容和步骤,完成以下实验,要求写出实验报告(实验目的问题数学模型算法与编程计算结果分析、检验和结论)基
3、础实验1求微分方程的解析解, 并画出它们的图形, y= y + 2x, y(0) = 1, 0x0.15时,可以看到图像逐渐发散,并且随着a值的变大,出现周期性的振荡。从右边的相图可以看见,当a增大时,单个极限环的稳定性逐渐失去,进入第二周期极限环,当a值继续变大时,第二周期的极限环的稳定性也失去,当a一直变大时,任意周期的极限环都将失去其稳定性,这就是产生混沌的必要条件。4.Apollo卫星的运动轨迹的绘制实验过程:首先创建一个M文件如下:function f=fun(x,y)u=1/82.45;u1=1-u;r1=sqrt(y1)+u)2+y(3)2),r2=sqrt(y(1)-u1)2+
4、y(3)2);f=y(2);2*y(4)+y(1)-u1*(y(1)+u)/(r1)3-u*(y(1)-u1)/(r2)3; y(4);-2*y(2)+y(3)-u1*y(3)/(r1)3-u*y(3)/(r2)3;然后,绘图程序如下:x,y=ode45(fun,0,20,1.2,0,0,-1.04935751);plot(y(:,1),y(:,3),grid ontitle(Appollo卫星运动轨迹),xlabel(X),ylabel(Y)实验结果如下:应用实验盐水的混合问题一个圆柱形的容器,内装350升的均匀混合的盐水溶液。如果纯水以每秒14升的速度从容器顶部流入,同时,容器内的混合的盐
5、水以每秒10.5升的速度从容器底部流出。开始时,容器内盐的含量为7千克。求经过时间t后容器内盐的含量。过程:分析:由于容器的容积相对于单位时间内水的体积变化来说很大,所以可以忽略溶质盐在在不同浓度的水内扩散至均匀的时间。根据在每个微小的时间段内,减少的盐加上容器内剩余的盐等于开始的盐量建立方程。实验过程:假设:容器内盐的含量在每一时刻都是均匀的,记y(t)为每时刻容器内盐的含量,T(t)为每时刻容器内水的体积,盐水流出的速率为v1,纯水流入的速率为v2则,在t内有:y(t+t)=y(t)-y(t)/T(t)*v1t又T(t)=T(t0)+(v2-v1)t,原方程可化为:y(t)+dy=y(t)
6、-y(t)v1dt/(T(t0)+(v2v1)t),化简得:dy/dt=yt*v1/(T(t0)+(v2v1)t)又由题目已知条件可得:T(t0)=350,v1=10.5,v2=14,y(0)=7,解微分方程得:y(t)= 7000000/(t+100)3实验结果及分析:当t=0时,y(0)=7,当t趋向于无穷大时,y(t)也趋向于0,但总有y(t)0,与实际情况相符的较好.总结体会:1.通过本次实验,我进一步熟悉了MatLab的工作环境,学会了用MatLab解方程的方法,掌握了MatLab数组的操作,强化了用MatLab作图的能力。2.使我对数学建模的过程有了一定的掌握,而使用matlab的使用方法和使用数学模型解决问题的思路有了进一步的了解和掌握,对以后数学实验的学习打下基础。Matlab的强大功能和简洁编程使平时一些很难得问题解决起来得心应手。3.学会怎样解微分方程。还要学会怎样建立合理的数学模型,以便使我们的实验结果更加正确4.;掌握解析、数值解法,并学会用图形观察解的形态和进行解的定性分析;熟悉MATLAB软件关于微分方程求解的各种命令.5.了解和熟悉了微分方程及方程组的解析求解法. 6.直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解)以及利用图形对解的特征作定性分析。
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1