(1)求解:
输入:
dsolve('Dy=y+2*x','y(0)=1','x')
输出:
ans=-2*x-2+3*exp(x)
(3)作图:
输入:
>>x=0:
0.1:
1;
>>y2=-2*x-2+3*exp(x);
>>plot(x,y2)
输出:
图表1方程特解图形
分析:
注意dsolve的用法。
2.用向前欧拉公式和改进的欧拉公式求方程y’=y-2x/y,y(0)=1(0≤x≤1,h=0.1)的数值解,要求编写程序,并比较两种方法的计算结果,说明了什么问题?
(1)求解析解
输入:
dsolve('Dy=y-2*x/y','y(0)=1','x')
输出:
ans=(2*x+1)^(1/2)
(2)用向前欧拉公式和改进的欧拉公式求方程的数值解并与解析解作图比较
程序:
x1
(1)=0;
y1
(1)=1;
y2
(1)=1;
h=0.1;
fork=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*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:
0.1:
1;
y=(2*x+1).^(1/2);
plot(x,y,x,y1,'o',x,y2,'+')
结果:
x1=00.10000.20000.30000.40000.50000.60000.70000.80000.90001.0000
y1=1.00001.10001.19181.27741.35821.43511.50901.58031.64981.71781.7848
y2=1.00001.09591.18411.26621.34341.41641.48601.55251.61651.67821.7379
图表2向前欧拉公式和改进的欧拉公式所求方程数值解与解析解的比较
由图可得,改进后的欧拉公式求得的数值解更贴合解析解。
分析:
注意向前欧拉与改进后的欧拉公式的不同。
3.Rossler微分方程组:
当固定参数b=2,c=4时,试讨论随参数a由小到大变化(如a∈(0,0.65))而方程解的变化情况。
程序:
rossler.m:
functionxdot=rossler(t,x)
xdot=[0,-1,-1;1,0.1,0;x(3),0,-4]*x+[0,0,2]';
fangchengzu.m:
x0=[000.1];
[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))
gridon
结果:
a=0.1时:
a=0.25时:
a=0.5时:
a=0.6时:
上述图形表示了a由小到大变化时方程解的变化。
分析:
注意xdot的书写以及ode45的运用。
4.Apollo卫星的运动轨迹的绘制
程序:
apollo.m:
functionyp=apollo(t,x)
u=1/82.45;
u1=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)/r1^3-u*(x
(1)-u1)/r2^3;x(4);-2*x
(2)+x(3)-u1*x(3)/r1^3-u*x(3)/r2^3];
weixing.m:
x0=[1.2;0;0;-1.04935751];
[t,x]=ode45('apollo',[0,20],x0);
plot(x(:
1),x(:
3))
xlabel('x')
ylabel('y')
title('Apollo卫星运动轨迹')
结果:
图表3apollo卫星轨迹图
分析:
注意求数值解时,高阶微分方程必须等价的变为一阶微分方程组。
应用实验(或综合实验)
一、实验内容
盐水的混合问题
一个圆柱形的容器,内装350升的均匀混合的盐水溶液。
如果纯水以每秒14升的速度从容器顶部流入,同时,容器内的混合的盐水以每秒10.5升的速度从容器底部流出。
开始时,容器内盐的含量为7千克。
求经过时间t后容器内盐的含量。
二、问题分析
(1)已知:
水的密度为1kg/L,盐溶解度为36g。
可计算出7kg盐所需要的溶剂为194L水。
因此,由混合液体积即可知开始时刻的7kg盐是完全溶于水中的,并且没有饱和。
所以,整个过程为食盐水被再次稀释的过程,则不会出现有盐析出现象。
(2)由于容器的容积相对于单位时间内水的体积变化来说很大,所以可以忽略溶质盐在在不同浓度的水内扩散至均匀的时间。
根据在每个微小的时间段内,减少的盐加上容器内剩余的盐等于开始的盐量建立方程。
三、数学模型的建立与求解(一般应包括模型、求解步骤或思路,程序放在后面的附录中)
假设:
1)温度对盐在水中的溶解度变化影响不大。
2)任意时刻容器内混合的、流出的盐水都均匀。
3)水流入及盐水流出的速度均为匀速。
设注水时间为t,t时刻时容器内含盐量为P(t)、容器内混合盐水的体积为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=10.5
方程可化为:
dP/dt=-10.5*P(t)/(350+3.5*t),P(0)=7
用MATLAB求解该方程并作图。
四、实验结果及分析
求得方程的解析解为:
P(t)=7000000/(t+100)^3
曲线图像为:
图表4经过时间t后容器内盐的含量
五、附录(程序等)
y=dsolve('Dy=-14*y/(350+3.5*t)','y(0)=7','t')
ezplot('7000000/(t+100)^3',[0,100])
xlabel('t')
ylabel('P(t)')
总结与体会
通过该实验的学习,掌握微分方程(组)求解方法(解析法、欧拉法、梯度法、改进欧拉法等),对常微分方程的数值解法有一个初步了解,同时学会使用MATLAB软件求解微分方程的基本命令,学会建立微分方程方面的数学模型。
教师签名
年月日
小学二
(2)班班规
一、安全方面
1、每天课间不能追逐打闹。
2、中午和下午放学要结伴回家。
3、 公路上走路要沿右边走,过马路要注意交通安全。
4、 不能在上学路上玩耍、逗留。
二、学习方面
1、每天到校后,不允许在走廊玩耍打闹,要进教室读书。
2、每节课铃声一响,要快速坐好,安静地等老师来上课。
3、课堂上不做小动作,不与同桌说悄悄话,认真思考,积极回答问题。
4、养成学前预习、学后复习的好习惯。
每天按时完成作业,保证字迹工整,卷面整洁。
5、考试时做到认真审题,不交头接耳,不抄袭,独立完成答卷。
三、升旗排队和两操方面
1、升旗时,要快速出教室排好队,做到快、静、齐,安静整齐地排队走出课室门,班长负责监督。
2、上午第二节后,快速坐好,按要求做好眼保健操。
3、下午预备铃声一响,在座位上做眼保健操。
四、卫生方面
1、每组值日生早晨7:
35到校做值日。
2、要求各负其责,打扫要迅速彻底,打扫完毕劳动工具要摆放整齐。
3、卫生监督员(剑锋,锶妍,炜薪)要按时到岗,除负责自己的值日工作外,还要做好记录。
五、一日常规
1、每天学生到齐后,班长要检查红领巾。
2、劳动委员组织检查卫生。
3、每天负责领读的学生要督促学生学习。
4、上课前需唱一首歌,由文娱委员负责。
5、 做好两操。
6、 放学后,先做作业,然后帮助家长至少做一件家务事。
7、 如果有人违反班规,要到老师处说明原因。
班训:
坐如钟 站如松 快如风 静无声
班规:
课堂听讲坐如钟,精神集中认真听;
排队升旗站如松,做操到位展雄风;
做事迅速快如风,样样事情记得清;
自习课上静无声,踏实学习不放松;
个人努力进步快,团结向上集体荣;
我为领巾添光彩,标兵集体记我功。
加分标准
序号
考核项目
加分值
备注
1
单元考试满分
+2
2
单元考试85分以上
+1
3
课堂小测满分
+1
4
期中、期末考试满分
+3
5
在红领巾广播站投稿一次
+2
6
在校级活动中获奖
+5
7
作业十次全对得一颗星
+3
8
课堂上得到表扬
+1
9
班干部工作认真负责
+1
10
做好事、有利于班集体和学校的事
+2
11
进步比较明显
+2
12
连续一周该组值日卫生达标
本组值日生每人加2分
扣分标准
序号
考核项目
扣分值
备注
1
没交作业、不做晚作业
-1
2
忘带书本、学具
-1
3
迟到
-1
4
在课堂上被老师点名
-2
5
不穿校服,不戴红领巾
-1
6
吃零食、带钱、带玩具
-2
7
说脏话、打架
-3
请家长,写保证书
8
座位周围有垃圾
-2
9
课间操、眼保健操不认真做
-1
10
升旗时违反纪律
-2
11
来学校不进教室,在走廊聊天打闹
-1
12
体育课打闹说话、排队不整齐
-2
注:
每人基本分60分起,学期末核算总分,作为学期评先依据。