| 查看: 434 | 回复: 2 | ||
[求助]
求助把这段动画保存为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 |
回复此楼» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有6人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有4人回复
-
今年春晚有几个节目很不错，点赞！
已经有6人回复
-
球磨粉体时遇到了大的问题，请指教！
已经有15人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
-
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复
老明古寺
铜虫 (职业作家)
- 应助: 0 (幼儿园)
- 金币: 7492.6
- 散金: 124
- 红花: 7
- 帖子: 3383
- 在线: 179.7小时
- 虫号: 2988654
- 注册: 2014-02-23
- 专业: 意识障碍
2楼2016-05-27 18:18:51
hlg154003231
金虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 1881
- 红花: 1
- 帖子: 280
- 在线: 276.7小时
- 虫号: 1763257
- 注册: 2012-04-18
- 性别: GG
- 专业: 植物遗传学
3楼2016-05-27 20:01:35