MATLAB作业Word格式文档下载.docx
《MATLAB作业Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《MATLAB作业Word格式文档下载.docx（9页珍藏版）》请在冰豆网上搜索。
exp（-2*t）*（sin（2*t+pi/3）+cos（3*t））'
],'
y（0）=1'
y
（1）=3'
y（pi）=2'
Dy（0）=1'
Dy
（1）=2'
）;
vpa（y,20）
ans=
.2057654e-2*exp（-2.*t）*cos（3.*t）+.9602373e-1*exp（-2.*t）*sin（2.*t）+.0587e-2*exp（-2.*t）*cos（2.*t）+*exp（-2.*t）*t+*exp（-2.*t）*t^*exp（-2.*t）*t^*exp（-2.*t）+.2057654e-2*exp（-2.*t）*sin（3.*t）+*exp（-5.*t）
2、试求解下面微分方程的通解以及满足
条件下的解析解。
[x,y]=dsolve（'
D2x+5*Dx+4*x+3*y=exp（-6*t）*sin（4*t）'
2*Dy+y+4*Dx+6*x=exp（-6*t）*cos（4*t）'
x（0）=1'
x（pi）=2'
y（0）=0'
）；
vpa（x,10）
ans=
*exp（t）-（exp*t）^（1/4）*exp（t）^（13/4））+*exp*t）^（1/4））/exp（t）^（13/4）+*cos*t））/exp*t）-*sin*t））/exp*t）
vpa（y,10）
-*exp（t）+（exp*t）^（1/4）*exp（t）^（13/4））+*exp*t）^（1/4））/exp（t）^（13/4）-*cos*t））/exp*t）+*sin*t））/exp*t）
3、试求出微分方程
的解析解通解，并求出满足边界条件
的解析解。
symsxy;
y=dsolve（'
D2y-（2-1/x）*Dy+（1-1/x）*y=x^2*exp（-5*x）'
）
y=
C3*exp（t）+C2*exp（（t*（x-1））/x）+x^3/（exp（5*x）*（x-1））
y=dsolve（'
y
（1）=pi'
y（pi）=1'
（exp（t）*（exp（（x-1）/x）-x*exp（（x-1）/x）-pi*exp（（pi*（x-1））/x）-（x^3*exp（（pi*（x-1））/x））/exp（5*x）+（x^3*exp（（x-1）/x））/exp（5*x）+pi*x*exp（（pi*（x-1））/x）））/（exp（pi）*exp（（x-1）/x）-exp
（1）*exp（（pi*（x-1））/x）-x*exp（pi）*exp（（x-1）/x）+x*exp
（1）*exp（（pi*（x-1））/x））-（exp（（t*（x-1））/x）*（exp
（1）-x*exp
（1）-pi*exp（pi）-（x^3*exp（pi））/exp（5*x）+（x^3*exp
（1））/exp（5*x）+pi*x*exp（pi）））/（exp（pi）*exp（（x-1）/x）-exp
（1）*exp（（pi*（x-1））/x）-x*exp（pi）*exp（（x-1）/x）+x*exp
（1）*exp（（pi*（x-1））/x））+x^3/（exp（5*x）*（x-1））
4、Lotka-Volterra扑食模型方程为
，且初值为
，试求解该微分方程，并绘制相应的曲线。
symsxyt;
f=inline（'
[4*x
（1）-2*x
（1）*x
（2）;
x
（1）*x
（2）-3*x
（2）]'
t'
x'
[t,x]=ode45（f,[0,10],[2;
3]）;
plot（t,x）
5、是给出求解下面微分方程的MATLAB命令，
并绘制出
曲线。
试问该方程存在解析解吗选择四阶定步长Runge-Kutta算法求解该方程时，步长选择多少可以得出较好的精度，MATLAB语言给出的现成函数在速度、精度上进行比较。
该方程的解析解不存在
[x
（2）;
x（3）;
-t^2*x
（1）*x
（2）-t^2*x
（2）*x
（1）^2+exp（-t*x
（1））]'
[t,x]=ode45（f,[0,10],[2;
0;
0]）;
plot（t,x）
6、试用解析解和数值解的方法求解下面的微分方程组
解析解：
symstxy
[x,y]=dsolve（'
D2x=-2*x-3*Dx+exp（-5*t）'
D2y=2*x-3*y-4*Dx-4*Dy-sin（t）'
Dx（0）=2'
y（0）=3'
Dy（0）=4'
x=
17/（4*exp（t））-10/（3*exp（2*t））+1/（12*exp（5*t））
100/（3*exp（2*t））-265/（16*exp（t））-71/（5*exp（3*t））+11/（48*exp（5*t））+cos（t）/5-sin（t）/10+（51*t）/（4*exp（t））
数值解：
functiondx=apolloeq（t,x）
dx=[x
（2）;
-2*x
（1）-3*x
（2）+exp（-5*t）;
x（4）;
2*x
（1）-3*x（3）-4*x
（2）-4*x（4）-sin（t）];
x0=[1;
2;
3;
4];
[t,y]=ode45（'
apolloeq'
[0,20],x0）;
plot（y（:
1）,y（:
3））
7、下面的方程在传统微分方程教程中经常被认为是刚性微分方程。
使用常规微分方程解法和刚性微分方程解法分别求解这两个微分方程的数值解，并求出解析解，用状态变量曲线比较数值求解的精度。
functionydot=lorenzeq（t,y）
ydot=[9*y
（1）+24*y
（2）+5*cos（t）-1/3*sin（t）;
-24*y
（1）-51*y
（2）-9*cos（t）+1/3*sin（t）];
t_final=100;
y0=[1/3;
2/3];
lorenzeq'
[0,t_final],y0）;
plot（t,y）
opt=odeset;
=1e-6;
[t,y]=ode15s（'
[0,t_final],y0,opt）;
plot（t,y）%刚性解法
ydot=[*y
（1）*y
（2）;
-50*y
（2）;
70*y
（2）-120*y（3）];
y0=[1;
1];
8、试求出习题3中给出的微分方程边值问题数值解，绘制出
曲线，并和该习题得出的解析解比较精度。
9、试用数值方法求解偏微分方程
，并绘制出
函数曲面。
pdetool
升级会员 