24小时热门版块排行榜    

查看: 430  |  回复: 2

pl03818948

新虫 (初入文坛)

[求助] 求助把这段动画保存为GIF,很简单的。急求! 已有1人参与

dx=0.1;x=(-8+dx:dx:8)';
nx=length(x);
k=dx^3;nsteps=2/k;
u=-5*(sech(x.^2)).^3
set(gcf,'doublebuffer','on');
for ii=1:nsteps
    k1=k*kdvequ(u,dx);
    k2=k*kdvequ(u+k1/2,dx);
    k3=k*kdvequ(u+k2/2,dx);
    k4=k*kdvequ(u+k3,dx);
    u=u+k1/6+k2/3+k3/3+k4/6;
    if mod(ii,10)==0
        plot(x,-u,'LineWidth',2);
        axis([-8,8,-2,12])
        
      drawnow
    end
end




function dudt=kdvequ(u,dx)
u=[u(end-1:end);u;u(1:2)];
dudt=-6*(u(3:end-2)).*(u(4:end-1)-u(2:end-3))/2/dx+(u(5:end)-2*u(4:end-1)+2*u(2:end-3)-u(1:end-4))/2/dx^3;

end
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老明古寺

铜虫 (职业作家)

2楼2016-05-27 18:18:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlg154003231

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
pl03818948: 金币+15, ★★★★★最佳答案 2016-07-30 12:08:35
GIF 图如下:
----------------------

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : mygif.gif
  • 2016-05-27 20:01:33, 221.6 K
3楼2016-05-27 20:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pl03818948 的主题更新
信息提示
请填处理意见