2022年重庆大学数学实验报告_微分方程模型、求解及稳定性分析 .pdf

开课学院、实验室：实验时间：年月日课程名称数学实验实验项目名称微分方程实验项目类型验证演示综合设计其他指导教师成绩实验目的1 归纳和学习求解常微分方程(组)的基本原理和方法；2 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析；3 熟悉 MATLAB 软件关于微分方程求解的各种命令；4 通过范例学习建立微分方程方面的数学模型以及求解全过程；通过该实验的学习，使学生掌握微分方程(组)求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方程的数值解法有一个初步了解，同时学会使用MATLAB 软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。这对于学生深入理解微分、积分的数学概念，掌握数学的分析思维方法，熟悉处理大量的工程计算问题的方法是十分必要的。基础实验一、实验内容1 微分方程及方程组的解析求解法；2 微分方程及方程组的数值求解法欧拉、欧拉改进算法；3 直接使用MATLAB 命令对微分方程(组)进行求解(包括解析解、数值解)；4 利用图形对解的特征作定性分析；5 建立微分方程方面的数学模型，并了解建立数学模型的全过程。二、实验过程（一般应包括实验原理或问题分析，算法设计、程序、计算、图表等，实验结果及分析）基础实验1求微分方程的解析解,并画出它们的图形，(1)y=y+2x,y(0)=1,0 x1;程序：syms x y k k=dsolve(Dy=y+2*x,y(0)=1,x)ezplot(k,0,1)结果：ans=-2*x-2+3*exp(x)图像：名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 6 页 -(2)y+ycos(x)=0,y(0)=1,y(0)=0;程序：syms x y z z=dsolve(D2y+y*cos*(x)=0,y(0)=1,Dy(0)=0)ezplot(z)结果：z=cos(cos(1/2)*x(1/2)*t)图像：3.Rossler微分方程组当固定参数b=2,c=4 时，试讨论随参数a由小到大变化（如a(0,0.65)而方程解的变化情况，并且画出空间曲线图形，观察空间曲线是否形成混沌状？程序：建立 rossler.m文件：function r=rossler(t,x)global a;global b;global c;)(cxzbzayxyzyx名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 6 页 -r=-x(2)-x(3);x(1)+a*x(2);b+x(3)*(x(1)-c);建立 exp4-3.m 如下：global a;global b;global c;b=2;c=4;t0=0,200;for a=0:0.03:0.65 t,x=ode45(rossler,t0,0,0,0);subplot(1,2,1);plot(t,x(:,1),r,t,x(:,2),g,t,x(:,3),b);title(x(红色),y(绿色),z(蓝色)随 t 变化情况);xlabel(t);subplot(1,2,2);plot3(x(:,1),x(:,2),x(:,3)title(相图);xlabel(x);ylabel(y);zlabel(z);end 1.当a=0时，图像：050100150200-1-0.8-0.6-0.4-0.200.20.40.6x(红 色),y(绿 色),z(篮 色)随 t变 化 情 况t-0.500.5-1-0.5000.20.40.60.8x相 图yz所以当 a=0 时，(x,y,z)收敛于（0,0.5,0.5）2.当a=0.12 时，图像：050100150200-1.2-1-0.8-0.6-0.4-0.200.20.40.6x(红 色),y(绿 色),z(篮 色)随 t变 化情 况t-0.500.51-1.5-1-0.5000.20.40.60.8x相 图yz当a=0.27 时，(x,y,z)仍然收敛，但是收敛速度大大降低。3.当a=0.27 时，图像名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 6 页 -050100150200-5-4-3-2-101234x(红 色),y(绿 色),z(篮 色)随 t变 化情 况t-505-50500.511.522.53x相 图yz4.当a=0.39 时，图像050100150200-6-4-20246x(红 色),y(绿 色),z(篮 色)随 t变 化 情 况t-50510-10-5050123456x相 图yz5.当a=0.52 时，图像050100150200-10-5051015x(红 色),y(绿 色),z(篮 色)随 t变 化 情 况t-50510-10-505051015x相 图yz从这一系列变化的图像中可以看出，随着 a 的增大，(x,y,z)接近极限环的速度加快。4.Apollo卫星的运动轨迹的绘制名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 6 页 -程序：x0=1.2;0;0;-1.04935751;options=odeset(reltol,1e-8);t,y=ode45(appollo,0,20,x0,options);plot(y(:,1),y(:,3)title(Appollo卫星运动轨迹)xlabel(X)ylabel(Y)绘图如下：应用实验5盐水的混合问题一个圆柱形的容器，内装 350 升的均匀混合的盐水溶液。如果纯水以每秒14 升的速度从容器顶部流入，同时，容器内的混合的盐水以每秒10.5 升的速度从容器底部流出。开始时，容器内盐的含量为7 千克。求经过时间t后容器内盐的含量。做出如下假设：1.假设在不同浓度的水中的盐扩散速度都相同。2.假设任何时刻容器内的盐水都是均匀的。1133121331212222121()()2,2,1/82.45,1,(),()(0)1.2,(0)0,(0)0,(0)1.04935751xxxyxrryyyxyrrrxyrxyxxyy名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 6 页 -3.用 y（t）表示容器内t 时刻的盐的含量，用W（t）表示容器内t 时刻的水的总量，用O表示盐水流出的速度，用I 代表纯水流入的速度，时间变化t后容器内盐的含量为y（t+t）。考虑在t内流出的盐水的为Ot则其流出的盐为tOtWty)()(.通过以上假设可以得如下模型：tOtWtytytty)()()()(tOItWOdttytydyty)()()()()(0化简可得tOItWOtydtdyo)()()(程序：syms x y z k;k=dsolve(Dy=-(y*Y)/(T(t0)+(C-Y)*t),y(0)=7,t)结果：(7*exp(O*log(W(t0)/(I-O)/exp(O*log(W(t0)+t*(I-O)/(I-O)由题目可知W（t0）=350，O=10.5，I=14，从而 y(t)=7000000/(t+100)3。总结与体会通过本次的实验，使我对数学建模的过程有了一定的掌握，而使用matlab的使用方法和使用数学模型解决问题的思路有了进一步的了解和掌握，对以后数学实验的学习打下基础。Matlab的强大功能和简洁编程使平时一些很难得问题解决起来得心应手。教师签名年月日名师资料总结-精品资料欢迎下载-名师精心整理-第 6 页，共 6 页 -