欢迎来到淘文阁 - 分享文档赚钱的网站! | 帮助中心 好文档才是您的得力助手!
淘文阁 - 分享文档赚钱的网站
全部分类
  • 研究报告>
  • 管理文献>
  • 标准材料>
  • 技术资料>
  • 教育专区>
  • 应用文书>
  • 生活休闲>
  • 考试试题>
  • pptx模板>
  • 工商注册>
  • 期刊短文>
  • 图片设计>
  • ImageVerifierCode 换一换

    数值分析作业8.docx

    • 资源ID:69680305       资源大小:84.61KB        全文页数:7页
    • 资源格式: DOCX        下载积分:15金币
    快捷下载 游客一键下载
    会员登录下载
    微信登录下载
    三方登录下载: 微信开放平台登录   QQ登录  
    二维码
    微信扫一扫登录
    下载资源需要15金币
    邮箱/手机:
    温馨提示:
    快捷下载时,用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)。
    如填写123,账号就是123,密码也是123。
    支付方式: 支付宝    微信支付   
    验证码:   换一换

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

    数值分析作业8.docx

    9-1 使用四阶Adams预测校正算法和经典RK算法。 RK 算法function yrk4=rk4(fxy,a,b,y0,n)format long;h=(b-a)/n;x=a:h:b;yrk4=zeros(n,1);yrk4(1)=y0;for i=1:n-1 k1=fxy(x(i),yrk4(i); k2=fxy(x(i)+h/2,yrk4(i)+h*k1/2); k3=fxy(x(i)+h/2,yrk4(i)+h*k2/2); k4=fxy(x(i)+h,yrk4(i)+h*k3); yrk4(i+1)=yrk4(i)+h*(k1+2*k2+2*k3+k4)/6;endfunction z=dm(x,y)z=-y+2*cos(x);>> yrk4=rk4(dm,0,pi,1,10)yrk4 = 1.000000000000000 1.260037091060021 1.396736267543184 1.396715573348376 1.259976385352878 0.999903214083919 0.641953487423769 0.221165565757292 -0.221271082024303 -0.642047809096457精确值x=0:pi/10:pi;y=cos(x)+sin(x)>> deey = Columns 1 through 4 1.000000000000000 1.260073510670101 1.396802246667421 1.396802246667421 Columns 5 through 8 1.260073510670101 1.000000000000000 0.642039521920206 0.221231742082474 Columns 9 through 11 -0.221231742082474 -0.642039521920206 -1.0000000000000009-2. function re(a,b,c,d,max)f1=(t,x,y,z)-y-z;f2=(t,x,y,z)x+a*y;f3=(t,x,y,z)(x-c)*z+b;x(1)=d;y(1)=d;z(1)=d;h=0.1;i=1;for t=0:h:max-h k1(1)=f1(t,x(i),y(i),z(i); k1(2)=f2(t,x(i),y(i),z(i); k1(3)=f3(t,x(i),y(i),z(i); k2(1)=f1(t+h,x(i)+h*k1(1),y(i)+h*k1(1),z(i)+h*k1(1); k2(2)=f2(t+h,x(i)+h*k1(2),y(i)+h*k1(2),z(i)+h*k1(2); k2(3)=f3(t+h,x(i)+h*k1(3),y(i)+h*k1(3),z(i)+h*k1(3); x(i+1)=x(i)+h/2*(k1(1)+k2(1); y(i+1)=y(i)+h/2*(k1(2)+k2(2); z(i+1)=z(i)+h/2*(k1(3)+k2(3); i=i+1;endfigure (1)plot3(0:h:max,x,y,'-');figure (2)plot(x,y);endre (0.2,0.8,4.5,1,1000)9-3 function de(h,max)f1=(t,x,y,z,m,n,p)m;f2=(t,x,y,z,m,n,p)n;f3=(t,x,y,z,m,n,p)p;f4=(t,x,y,z,m,n,p)2*x*z*m+3*x2*y*t2;f5=(t,x,y,z,m,n,p)exp(y)*n+4*x*z*t2;f6=(t,x,y,z,m,n,p)2*t*p+2*t*exp(x*y);x(1)=3;y(1)=3;z(1)=3;m(1)=2;n(1)=2;p(1)=2;i=1;for t=1:h:max-h k1(1)=f1(t,x(i),y(i),z(i),m(i),n(i),p(i); k1(2)=f2(t,x(i),y(i),z(i),m(i),n(i),p(i); k1(3)=f3(t,x(i),y(i),z(i),m(i),n(i),p(i); k1(4)=f4(t,x(i),y(i),z(i),m(i),n(i),p(i); k1(5)=f5(t,x(i),y(i),z(i),m(i),n(i),p(i); k1(6)=f6(t,x(i),y(i),z(i),m(i),n(i),p(i); k2(1)=f1(t+h,x(i)+h*k1(1),y(i)+h*k1(1),z(i)+h*k1(1),m(i)+h*k1(1),n(i)+h*k1(1),p(i)+h*k1(1); k2(2)=f2(t+h,x(i)+h*k1(2),y(i)+h*k1(2),z(i)+h*k1(2),m(i)+h*k1(2),n(i)+h*k1(2),p(i)+h*k1(2); k2(3)=f3(t+h,x(i)+h*k1(3),y(i)+h*k1(3),z(i)+h*k1(3),m(i)+h*k1(3),n(i)+h*k1(3),p(i)+h*k1(3); k2(4)=f4(t+h,x(i)+h*k1(4),y(i)+h*k1(4),z(i)+h*k1(4),m(i)+h*k1(4),n(i)+h*k1(4),p(i)+h*k1(4); k2(5)=f5(t+h,x(i)+h*k1(5),y(i)+h*k1(5),z(i)+h*k1(5),m(i)+h*k1(5),n(i)+h*k1(5),p(i)+h*k1(5); k2(6)=f6(t+h,x(i)+h*k1(6),y(i)+h*k1(6),z(i)+h*k1(6),m(i)+h*k1(6),n(i)+h*k1(6),p(i)+h*k1(6); x(i+1)=x(i)+h/2*(k1(1)+k2(1); y(i+1)=y(i)+h/2*(k1(2)+k2(2); z(i+1)=z(i)+h/2*(k1(3)+k2(3); m(i+1)=m(i)+h/2*(k1(4)+k2(4); n(i+1)=n(i)+h/2*(k1(5)+k2(5); p(i+1)=p(i)+h/2*(k1(6)+k2(6); i=i+1;endfigure (1)plot3(x,y,z,'-');ii=1;for t=1:h:max fprintf('t=%.6f,x=%.6f,y=%.6f,z=%.6fn',t,x(ii),y(ii),z(ii) ii=ii+1;endendde(1e-5,1.03325)t=1.033080,x=3.280523,y=3.151594,z=45.355426t=1.033090,x=3.280953,y=3.151698,z=45.523123t=1.033100,x=3.281384,y=3.151802,z=45.695588t=1.033110,x=3.281816,y=3.151906,z=45.873292t=1.033120,x=3.282249,y=3.152011,z=46.056803t=1.033130,x=3.282684,y=3.152115,z=46.246814t=1.033140,x=3.283120,y=3.152220,z=46.444196t=1.033150,x=3.283557,y=3.152324,z=46.650070t=1.033160,x=3.283996,y=3.152429,z=46.865926t=1.033170,x=3.284437,y=3.152534,z=47.093846t=1.033180,x=3.284878,y=3.152639,z=47.336910t=1.033190,x=3.285321,y=3.152744,z=47.600105t=1.033200,x=3.285766,y=3.152849,z=47.892665t=1.033210,x=3.286212,y=3.152954,z=48.236208t=1.033220,x=3.286659,y=3.153059,z=48.713819t=1.033230,x=3.287108,y=3.153165,z=51.049529t=1.033240,x=3.287558,y=3.153270,z=77662871239455918000000.000000t=1.033250,x=3.288010,y=3.153376,z=Inf

    注意事项

    本文(数值分析作业8.docx)为本站会员(asd****56)主动上传,淘文阁 - 分享文档赚钱的网站仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知淘文阁 - 分享文档赚钱的网站(点击联系客服),我们立即给予删除!

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




    关于淘文阁 - 版权申诉 - 用户使用规则 - 积分规则 - 联系我们

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

    工信部备案号:黑ICP备15003705号 © 2020-2023 www.taowenge.com 淘文阁 

    收起
    展开