实验五种群数量的状态转移微分方程.docx

课程 名称 数学建模 实验项目 名 称 实验五 种群数量的 状态转移 微分方 程 实验项目类型 验 证 演 示 综 合 设 计 其 /、 他 指导 教师 成绩 实验目的 归纳和学习求解常微分方程（组）的基本原理和方法； 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析； 熟悉MATLAB^件关于微分方程求解的各种命令； 通过范例学习建立微分方程方面的数学模型以及求解全过程； 基础实验 一、 实验内容 ?微分方程及方程组的解析求解法； ?微分方程及方程组的数值求解法——欧拉、欧拉改进算法； 3?直接使用MATLAB^令对微分方程（组）进行求解（包括解析解、数值解）； 4.利用图形对解的特征作定性分析； 5?建立微分方程方面的数学模型，并了解建立数学模型的全过程。 二、 实验过程（一般应包括实验原理或问题分析，算法设计、程序、计算、图表等， 实 验结果及分析） 1.求微分方程的解析解，并画出它们的图形， y = y + 2 x, y（0） = 1, 0 x1； 实验原理：y=dsolve（eqn1,eqn2, …,c1 , …,x ） 问题分析： y -Dy 图表： 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 x 实验结果： y =-2*x-2+exp(x)*C1 y1 =-2*x-2+3*exp(x) 实验分析： 由图表知结果正确。 2?用向前欧拉公式和改进的欧拉公式求方程 y = y - 2 x/y, y(0) = 1 (0 x 1,h = 0.1) 的数值解，要求编写程序，并比较两种方法的计算结果，说明了什么问题？ 实验原理： 1)向前欧拉公式：(y = f (x, y)) y (Xn+1)? y(Xn) + h f(Xn, y(Xn))(迭代式) yn+1 ? yn + h f (Xn, yn)( 近似式) 特点：f (X,y )取值于区间［Xn, Xn+1］的左端点 2) 向后欧拉公式 yn+1 ? y n + h f (Xn +1 , y n +1) 特点：①f (X,y )取值于区间［Xn, Xn+1］的右端点. ② 非线性方程，称‘隐式公式’。 问题分析：用不同的方法来求解。用向前欧拉公式时，用Xn代替，用向后欧拉公式，用x(n+1). 程序： x1(1)=0;y1(1)=1;y2(1)=1; for k=1:10 x1(k+1)=x1(k)+0.1; y1(k+1)=y1(k)+0.1*(y1(k)-2*x1(k)/y1(k)); y2(k+1)=y2(k)+0.1*(y2(k+1)-2*x1(k+1)/y2(k+1)); end x1,y1,y2, x=0:0.1:1; 结果x1 =0,0 0.3000 0.4000 0.50000.1000 0.20C 结果 x1 =0 ,0 0.3000 0.4000 0.5000 0.1000 0.20C 0.6000 0.7000 0.8000 0.9000 1.0000 y=sqrt(2*x+1) 1.8 plot(x,y,x1,y1,k:,x1,y2,r--) 1.7 图表： 1.6 y =1.0000 1.0954 1.18321.2649 1.34161.4142 y =1.0000 1.0954 1.1832 1.2649 1.3416 1.4142 1.2 y1 =1.0000 1.1000 1.1918 1.2774 1.3582 1.4351 10 01i5098.2 1.5803 d.?498^1%78 b.78480.8 0.9 y2 =1.0000 1.0811 1.1448 1.1928 1.2273 1.2501 1.2635 1.2694 1.2698 1.2662 1.2604 1.4832 1.5492 1.6125 1.6733 1.7321 分析：向前欧拉公式较接近真实结果，向后欧拉公式差距较大 应用实验（或综合实验） 一、实验内容 盐水的混合问题 一个圆柱形的容器，内装350升的均匀混合的盐水溶液。如果纯水以每秒 14升的速度 含盐量随时间而降低。 五、附录(程序等) y=zeros(1,1001);v=350;y(1)=7; dt=.1; for i=2:1001 v=350+14*dt-10.5*dt; y(i)=y(i-1)-y(i-1)/v*10.5*dt; end plot(1:1001,y) 总结与体会 通过该实验的学习，掌握微分方程（组）求解方法（解析法、欧拉法、梯度法、改进欧拉法等）， 对常微分方程的数值解法有一个初步了解，同时学会使用 MATLAB软件求解微分方程的基本 命令，学会建立微分方程方面的数学模型。 教师签名 年 月曰