MATLAB作业.docx
- 文档编号:5330796
- 上传时间:2022-12-15
- 格式:DOCX
- 页数:9
- 大小:178.85KB
MATLAB作业.docx
《MATLAB作业.docx》由会员分享,可在线阅读,更多相关《MATLAB作业.docx(9页珍藏版)》请在冰豆网上搜索。
MATLAB作业
MATLAB作业5
1、试求出下面线性微分方程的通解。
假设上述微分方程满足已知条件,
试求出满足该条件的微分方程的解析解。
解:
>>symsty;
y=dsolve(['D5y+13*D4y+64*D3y+152*D2y+176*Dy+80*y=','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)
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)
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('D2y-(2-1/x)*Dy+(1-1/x)*y=x^2*exp(-5*x)','y
(1)=pi','y(pi)=1')
y=
(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语言给出的现成函数在速度、精度上进行比较。
该方程的解析解不存在
>>f=inline('[x
(2);x(3);-t^2*x
(1)*x
(2)-t^2*x
(2)*x
(1)^2+exp(-t*x
(1))]','t','x');
[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)','x(0)=1','Dx(0)=2','y(0)=3','Dy(0)=4')
x=
17/(4*exp(t))-10/(3*exp(2*t))+1/(12*exp(5*t))
y=
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];
>>[t,y]=ode45('lorenzeq',[0,t_final],y0);
>>plot(t,y)
>>opt=odeset;=1e-6;
>>[t,y]=ode15s('lorenzeq',[0,t_final],y0,opt);
>>plot(t,y)%刚性解法
解:
functionydot=lorenzeq(t,y)
ydot=[*y
(1)*y
(2);-50*y
(2);70*y
(2)-120*y(3)];
>>t_final=100;y0=[1;2;1];
[t,y]=ode45('lorenzeq',[0,t_final],y0);
>>plot(t,y)
>>opt=odeset;=1e-6;
>>[t,y]=ode45('lorenzeq',[0,t_final],y0,opt);
>>plot(t,y)%刚性解法
>>opt=odeset;=1e-6;
>>[t,y]=ode45('lorenzeq',[0,t_final],y0,opt);
>>plot(t,y)
8、试求出习题3中给出的微分方程边值问题数值解,绘制出
曲线,并和该习题得出的解析解比较精度。
9、试用数值方法求解偏微分方程
,并绘制出
函数曲面。
>>pdetool
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 作业