数学实验报告中学教育中学实验_中学教育-中学实验.pdf

《数学实验报告中学教育中学实验_中学教育-中学实验.pdf》由会员分享，可在线阅读，更多相关《数学实验报告中学教育中学实验_中学教育-中学实验.pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 DS1407 学 院 自动化 年级 2013 专业班 自动化02班 学 生 姓 名 侯刚 学 号 开 课 时 间 2014 至 2015 学年第 二 学期 总 成 绩 教师签名 数 学 与 统 计 学 院 制 开课学院、实验室：数统学院 DS1407 实验时间：2014年 4月 3日 课程 名称 数学实验 实验项目 名 称 种群数量的状态转移微分方程 实验项目类型 验证 演示 综合 设计 其他 指导 教师 龚劬 成 绩 实验目的 1 归纳和学习求解常微分方程(组)的基本原理和方法；2 掌握解析、数值解法，并学会用图形观

2、察解的形态和进行解的定性分析；3 熟悉 MATLAB 软件关于微分方程求解的各种命令；4 通过范例学习建立微分方程方面的数学模型以及求解全过程；基础实验 一、实验内容 1.微分方程及方程组的解析求解法；2.微分方程及方程组的数值求解法欧拉、欧拉改进算法；3.直接使用 MATLAB 命令对微分方程(组)进行求解(包括解析解、数值解)；4.利用图形对解的特征作定性分析；5.建立微分方程方面的数学模型，并了解建立数学模型的全过程。二、实验过程 1求微分方程的解析解,并画出它们的图形，y=y+2x,y(0)=1,0 x x=0:1;y2=-2*x-2+3*exp(x);plot(x,y2)输出：图表

3、1 方程特解图形 分析：注意 dsolve 的用法。2用向前欧拉公式和改进的欧拉公式求方程y=y-2x/y,y(0)=1(0 x1,h=的数值解，要求编写程序，并比较两种方法的计算结果，说明了什么问题？（1）求解析解 输入：dsolve(Dy=y-2*x/y,y(0)=1,x)输出：ans=(2*x+1)(1/2)（2）用向前欧拉公式和改进的欧拉公式求方程的数值解并与解析解作图比较 程序：x1(1)=0;y1(1)=1;时间至学年第二学期总成绩教师签名数学与统计学院制开课学院实验室数统学院实验时间年月日课程名称指导教师实验目的数学实验龚劬实验项目名称成绩种群数量的状态转移微分方程实验项目类型验

4、证演示综合设计其他归纳和学件关于微分方程求解的各种命令通过范例学习建立微分方程方面的数学模型以及求解全过程基础实验一实验内容微分方程及方程组的解析求解法微分方程及方程组的数值求解法欧拉欧拉改进算法直接使用命令对微分方程组进行求解实验过程求微分方程的解析解并画出它们的图形求解输入输出作图输入输出图表方程特解图形分析注意的用法用向前欧拉公式和改进的欧拉公式求方程的数值解要求编写程序并比较两种方法的计算结果说明了什么问题求解析解输入y2(1)=1;h=;for k=1:10 x1(k+1)=x1(k)+h;y1(k+1)=y1(k)+h*(y1(k)-2*x1(k)/y1(k);k1=y2(k)-2

5、*x1(k)/y2(k);k2=y2(k)+h*k1-2*x1(k+1)/(y2(k)+h*k1);y2(k+1)=y2(k)+h*(k1+k2)/2;end x1,y1,y2 x=0:1;y=(2*x+1).(1/2);plot(x,y,x,y1,o,x,y2,+)结果：x1=0 y1=y2=图表 2 向前欧拉公式和改进的欧拉公式所求方程数值解与解析解的比较 由图可得，改进后的欧拉公式求得的数值解更贴合解析解。分析：注意向前欧拉与改进后的欧拉公式的不同。3Rossler微分方程组：当固定参数 b=2,c=4 时，试讨论随参数a由小到大变化（如a(0,)而方程解的变化情况。程序：:functi

6、on xdot=rossler(t,x)xdot=0,-1,-1;1,0;x(3),0,-4*x+0,0,2;:x0=0 0;t,x=ode45(rossler,0,10,x0);plot(t,x(:,1),-,t,x(:,2),.,t,x(:,3),+)pause plot3(x(:,1),x(:,2),x(:,3)grid on 结果：a=时：a=时：)(cxzbzayxyzyx时间至学年第二学期总成绩教师签名数学与统计学院制开课学院实验室数统学院实验时间年月日课程名称指导教师实验目的数学实验龚劬实验项目名称成绩种群数量的状态转移微分方程实验项目类型验证演示综合设计其他归纳和学件关于微分方

7、程求解的各种命令通过范例学习建立微分方程方面的数学模型以及求解全过程基础实验一实验内容微分方程及方程组的解析求解法微分方程及方程组的数值求解法欧拉欧拉改进算法直接使用命令对微分方程组进行求解实验过程求微分方程的解析解并画出它们的图形求解输入输出作图输入输出图表方程特解图形分析注意的用法用向前欧拉公式和改进的欧拉公式求方程的数值解要求编写程序并比较两种方法的计算结果说明了什么问题求解析解输入 a=时：a=时：上述图形表示了 a 由小到大变化时方程解的变化。分析：注意 xdot 的书写以及 ode45 的运用。卫星的运动轨迹的绘制 程序：:function yp=apollo(t,x)u=1/;u

8、1=1-u;r1=sqrt(x(1)+u)2+x(3)2);r2=sqrt(x(1)-u1)2+x(3)2);yp=x(2);2*x(4)+x(1)-u1*(x(1)+u)/r13-u*(x(1)-u1)/r23;x(4);-2*x(2)+x(3)-u1*x(3)/r13-u*x(3)/r23;.m:x0=;0;0;t,x=ode45(apollo,0,20,x0);plot(x(:,1),x(:,3)xlabel(x)ylabel(y)title(Apollo 卫星运动轨迹)结果：1133121331212222121()()2,2,1/82.45,1,(),()(0)1.2,(0)0,(0

9、)0,(0)1.04935751xxxyxrryyyxyrrrxyrxyxxyy&时间至学年第二学期总成绩教师签名数学与统计学院制开课学院实验室数统学院实验时间年月日课程名称指导教师实验目的数学实验龚劬实验项目名称成绩种群数量的状态转移微分方程实验项目类型验证演示综合设计其他归纳和学件关于微分方程求解的各种命令通过范例学习建立微分方程方面的数学模型以及求解全过程基础实验一实验内容微分方程及方程组的解析求解法微分方程及方程组的数值求解法欧拉欧拉改进算法直接使用命令对微分方程组进行求解实验过程求微分方程的解析解并画出它们的图形求解输入输出作图输入输出图表方程特解图形分析注意的用法用向前欧拉公式和改

10、进的欧拉公式求方程的数值解要求编写程序并比较两种方法的计算结果说明了什么问题求解析解输入 图表 3 apollo卫星轨迹图 分析：注意求数值解时，高阶微分方程必须等价的变为一阶微分方程组。应用实验（或综合实验）一、实验内容 盐水的混合问题 一个圆柱形的容器，内装 350 升的均匀混合的盐水溶液。如果纯水以每秒 14 升的速度从容器顶部流入，同时，容器内的混合的盐水以每秒 10.5 升的速度从容器底部流出。开始时，容器内盐的含量为 7 千克。求经过时间t后容器内盐的含量。二、问题分析（1）已知：水的密度为 1kg/L,盐溶解度为 36g。可计算出 7kg 盐所需要的溶剂为 194L水。因此，由混

11、合液体积即可知开始时刻的 7kg 盐是完全溶于水中的,并且没有饱和。所以，整个过程为食盐水被再次稀释的过程，则不会出现有盐析出现象。（2）由于容器的容积相对于单位时间内水的体积变化来说很大，所以可以忽略溶质盐在在不同浓度的水内扩散至均匀的时间。根据在每个微小的时间段内，减少的盐加上容器内剩余的盐等于开始的盐量建立方程。三、数学模型的建立与求解(一般应包括模型、求解步骤或思路，程序放在后面的附录中)假设：1）温度对盐在水中的溶解度变化影响不大。2）任意时刻容器内混合的、流出的盐水都均匀。3）水流入及盐水流出的速度均为匀速。设注水时间为 t，t 时刻时容器内含盐量为 P(t)、容器内混合盐水的体积

12、为 V(t)，纯水流入容器的速度为 v1,混合液流出的速度为 v2。可列出方程组：P(t+t)=P(t)-P(t)*v2*t/V(t)V(t)=V(t0)+(v1-v2)*t V(t0)=350,P(0)=7,v1=14,v2=方程可化为：dP/dt=*P(t)/(350+*t),P(0)=7 用 MATLAB 求解该方程并作图。四、实验结果及分析 求得方程的解析解为：P(t)=7000000/(t+100)3 曲线图像为：图表 4 经过时间t后容器内盐的含量 时间至学年第二学期总成绩教师签名数学与统计学院制开课学院实验室数统学院实验时间年月日课程名称指导教师实验目的数学实验龚劬实验项目名称成

13、绩种群数量的状态转移微分方程实验项目类型验证演示综合设计其他归纳和学件关于微分方程求解的各种命令通过范例学习建立微分方程方面的数学模型以及求解全过程基础实验一实验内容微分方程及方程组的解析求解法微分方程及方程组的数值求解法欧拉欧拉改进算法直接使用命令对微分方程组进行求解实验过程求微分方程的解析解并画出它们的图形求解输入输出作图输入输出图表方程特解图形分析注意的用法用向前欧拉公式和改进的欧拉公式求方程的数值解要求编写程序并比较两种方法的计算结果说明了什么问题求解析解输入五、附录（程序等）y=dsolve(Dy=-14*y/(350+*t),y(0)=7,t)ezplot(7000000/(t+1

14、00)3,0,100)xlabel(t)ylabel(P(t)总结与体会 通过该实验的学习，掌握微分方程(组)求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方程的数值解法有一个初步了解，同时学会使用MATLAB 软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。教师签名 年 月 日 时间至学年第二学期总成绩教师签名数学与统计学院制开课学院实验室数统学院实验时间年月日课程名称指导教师实验目的数学实验龚劬实验项目名称成绩种群数量的状态转移微分方程实验项目类型验证演示综合设计其他归纳和学件关于微分方程求解的各种命令通过范例学习建立微分方程方面的数学模型以及求解全过程基础实验一实验内容微分方程及方程组的解析求解法微分方程及方程组的数值求解法欧拉欧拉改进算法直接使用命令对微分方程组进行求解实验过程求微分方程的解析解并画出它们的图形求解输入输出作图输入输出图表方程特解图形分析注意的用法用向前欧拉公式和改进的欧拉公式求方程的数值解要求编写程序并比较两种方法的计算结果说明了什么问题求解析解输入