（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软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。
教师签名
年月日
-7-
升级会员 