欢迎来到优知文库! | 帮助中心 分享价值,成长自我!
优知文库
全部分类
  • 幼儿/小学教育>
  • 中学教育>
  • 高等教育>
  • 研究生考试>
  • 外语学习>
  • 资格/认证考试>
  • 论文>
  • IT计算机>
  • 法律/法学>
  • 建筑/环境>
  • 通信/电子>
  • 医学/心理学>
  • ImageVerifierCode 换一换
    首页 优知文库 > 资源分类 > DOCX文档下载
    分享到微信 分享到微博 分享到QQ空间

    MATLAB作业5参考答案.docx

    • 资源ID:396350       资源大小:29.90KB        全文页数:4页
    • 资源格式: DOCX        下载积分:3金币
    快捷下载 游客一键下载
    账号登录下载
    微信登录下载
    三方登录下载: QQ登录
    二维码
    扫码关注公众号登录
    下载资源需要3金币
    邮箱/手机:
    温馨提示:
    快捷下载时,如果您不填写信息,系统将为您自动创建临时账号,适用于临时下载。
    如果您填写信息,用户名和密码都是您填写的【邮箱或者手机号】(系统自动生成),方便查询和重复下载。
    如填写123,账号就是123,密码也是123。
    支付方式: 支付宝    微信支付   
    验证码:   换一换

    加入VIP,免费下载
     
    账号:
    密码:
    验证码:   换一换
      忘记密码?
        
    友情提示
    2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,就可以正常下载了。
    3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
    4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰。
    5、试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。

    MATLAB作业5参考答案.docx

    MATLAB作业5参考答案1、试求出下面线性微分方程的通解。+13+64+152+176+80y(r)=e2zsin(2r+-)+cos(3r)drdrdrdrdt3假设上述微分方程满足已知条件XO)=1,WI)=3,y)=2,P(O)=1,(l)=2,试求出满足该条件的微分方程的解析解。【求解】先定义七为符号变量,求出等号右侧的函数,则可以由下面叮嘱求出方程的解析解,解的规模较大,经常能占数页。»symstexp(-2*t)*(sin(2*t+sym(pi)/3)+cos(3*t)ans=exp(-2*t)*(sin(2*t+13*pi)+cos(3*t)»y=dsolve(,D5y+13*D4y+64*D3y÷152*D2y+176*Dy+80*y=,exp(-2*t)*(sin(2*t+13*pi)+cos(3*t),'y(0)=,y(1)=3',y(pi)=2,.,Dy(0)=,Dy(1)=2,)略:事实上,细致阅读求出的解析解就会发觉,其中大部分表达式是关于系数的,所以假如能对系数进行近似则将大大减小解的困难度。»vpa(y)ans=.20576131687242798353909465020576e-2*exp(-2.*t)*cos*t)+.15538705805619602372728107411086e-1*exp(-2.*t)*sin(2.*t)+.76830587084294035590921611166287e-2*exp(-2.*t)*cos(2.*t)-106.24422608844727797303237726774*exp(-2.*t)*t2+98.159206062620455331994871615083*exp(-2.*t)*t+59.405044899367325888329709780356*e×p(-2.*t)*t3-30.741892776456442808809983330755*exp(-2.*t)+.20576131687242798353909465020576e-2*exp(-2.*t)*sin(3.*t)+31.732152104579289125415500223136*exp(-5.*t)2、试求解下面微分方程的通解以及满足MO)=I,%(乃)=2,y(0)=0条件下的解析解。x(r)+5x(r)+4x(r)+3y(t)=e'6tsin(4r)2y(r)+y(r)+4x()+6x(t)=e6tcos(4r)【求解】可以用下面的语句得出微分方程组的通解。»symstx,y=dsoIveCD2x+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)')解略。将已知初始条件代入,则可以得出下面的特解。»symstx,y=dsoIve('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),vpa(y)解略。3、试求出微分方程y(x)-(2-L)MX)+(1-L)y(x)=2"5的解析解通解,并求出满足XX边界条件XD=,y()=1的解析解。【求解】微分方程的通解可以由下面的函数干脆求出»symsXy=dsoIve('D2y-(2-1/x)*Dy+(1-1/x)*y=x2*exp(-5*x),x,)y=exp(x)*C2+e×p(x)*Iog(x)*C1÷11296*(6*exp(6*x)*Ei(1,6*x)+11÷30*x+36*x2)*exp(-5*x)若须要求取满足边界条件的特解,须要在求解时代入边界条件,这样就可以由下面的语句得出微分方程的特解。»symsXy=dsoIve('D2y-(2-1/x)*Dy+(1-1/x)*y=x2*exp(-5*x),.'y(1)=pi','y(pi)=,'x')y=-1/1296*exp(x)*(6*exp(1)*Ei(1,6)+77*exp(-5)-1296*sym(pi)/exp(1)+1/1296*exp(x)*Iog(x)*(6*Ei(1,6)*exp(6*sym(pi)+6)+77*exp(6*sym(pi)-1296*sym(pi)*exp(6*sym(pi)+5)+3*sqrt(-1)*pi*csgn(sym(pi)*exp(6*sym(pi)+6)-6*Ei(1,6*sym(pi)*exp(6*sym(pi)÷6)-3*sqrt(-1)*pi*exp(6*sym(pi)÷6)-3*sqrt(-1)*pi*csgn(6*sqrt(-1)*sym(pi)*exp(6*sym(pi)+6)-30*sym(pi)*exp(6)+3*sqrt(-1)*pi*csgn(sym(pi)*csgn(6*sqrt(-1)*sym(pi)*exp(6*sym(pi)+6)-36*sym(pi)2*exp(6)-11*exp(6)+1296*exp(5*sym(pi)+6)/1og(sym(pi)*exp(-6*sym(pi)-6)+1/1296*(6*e×p(6*x)*Ei(1,6*x)+11+30*x+36*x2)*e×p(-5*x)由于运用了Sym(Pi)这样精确的表示,其实还可以用数值解的方法对各个系数进行近似,这样可以得出如下结果。»vpa(y,10)ans=1.155578411*e×p(x)9717266142*exp(x)*Iog(x)+.7716049383e-3*(6.*exp(6.*x)*Ei(1,6.*x)+11.+30.*x+36.*x2)*exp(-5.*x)还可以对该结果进行图形显示,得出如图所示的解曲线,可见,该曲线通过两个给定点。»x1=0.5:0.01:4;y1=subs(y,x,x1);plot(x1,y1,1,pi,o',pi,1,o,)4、LOlka-Volterra扑食模型方程为区)=4所2?且初值为MO)="(ox,j(Z)=XOXO-3o试求解该微分方程,并绘制相应的曲线。【求解】用下面叮嘱可以立刻解出微分方程的模型,并绘制出解的时间响应曲线和相平面曲线,分别如图所示。»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);5、是给出求解下面微分方程的MATLAB叮嘱,y÷(yy+r2y2=e-ty,y(0)=2,y(0)=y(0)=0并绘制出y(t)曲线。试问该方程存在解析解吗?选择四阶定步长Runge-Kutta算法求解该方程时,步长选择多少可以得出较好的精度,MATLAB语言给出的现成函数在速度、精度上进行比较。【求解】因为该方程为含有非线性项的微分方程,所以一般来说该方程不能有解析解,只能求解其数值解。要求解该方程。,则须要引入状态变量x1,x2,x3t这样可以写出一阶微分方程组。就可以用下面的语句求解微分方程,并绘制出状态变量的时域响应曲线,可见,时变微分方程也可以简洁地求解。»f=inIine(,x(2);x(3);-t2*x(1)*x(2)-t2*x(2)*x(1)2+exp(-t*x(D)','t','x');t,x=ode45(f,0,10,2;0;0);6、试用解析解和数值解的方法求解下面的微分方程组x(t)=-2x(t)-3x(r)+e,x(0)=1,x(0)=2JQ)=2x(/)-3X0-4X0-4j(0-Sin)Mo)=3,y(0)=4【求解】该方程的解析解可以由下面的语句干脆得出»symstx,y=dsoIve('D2x=-2*x-3*Dx+exp(-5*t),D2y=2*-3*y-4*D-4*Dy-sin(t)',.'x(0)=,'Dx(0)=2','y(0)=3,Dy(0)=41)X=1/12*exp(-5*t)-103*exp(-2*t)+174*exp(-t)y=-715*exp(-3*t)-265/16*exp(-t)+11/48*exp(-5*t)+1003*exp(-2*t)+15*cos(t)-1/10*sin(t)+51/4*exp(-t)*t»f=inline(.'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)',.,t',x');t1,x1=ode45(f,0,10,1;2;3;4);ezplot(x,0,10),Iine(t1,x1(:,1)figure;ezplot(y,0,10),Iine(t1,×1(:,3)得出的x(t)和y*)曲线的数值解与解析解分别如图所示,可见,得出的数值解精度还是很高的,在曲线上无法区分出解析解和数值解。7、下面的方程在传统微分方程教程中经常被认为是刚性微分方程。运用常规微分方程解法和刚性微分方程解法分别求解这两个微分方程的数值解,并求出解析解,用状态变量曲线比较数值求解的精度。(1)y)=j2() = %=9y+24%+5cosr-sinr,y2=-24yl-5Iy2-9cosr+-sinf,yi=-0.1y1-49.9y2,yi(0)=1(2)y2=-502,必(°)=2¼=70y2-120y3,y3(0)=l解略,参见讲义。8、试求出习题3中给出的微分方程边值问题数值解,绘制出y()曲线,并和该习题得出的解析解比较精度。【求解】微分方程的解析解可以由下面的函数干脆求出»symsXy=dsoIve('D2y-(2-1/x)*Dy+(1-1x)*y=x2*exp(-5*x),.'y(1)=pi','y(pi)=,x')ezplot(y,0,10)2u2u.7÷T=Ox2y29、试用数值方法求解偏微分方程0力=1,un>=0,并绘制出函数曲面。=0fy>0IV=O,0x>0,y>0【求解】可以看出,该方程为椭圆型偏微分方程,且c=1,f=0,用Pdetool函数供应的界面求解该方程,选择并设置椭圆微分方程,再选择矩形区域,可以设置出横轴和纵轴的边界条件,求解微分方程则可以得出所示的解曲面。

    注意事项

    本文(MATLAB作业5参考答案.docx)为本站会员(王**)主动上传,优知文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知优知文库(点击联系客服),我们立即给予删除!

    温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    copyright@ 2008-2023 yzwku网站版权所有

    经营许可证编号:宁ICP备2022001189号-2

    本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。优知文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知优知文库网,我们立即给予删除!

    收起
    展开