[转帖]（来自罗卜驿站）气泡在液体中运动

huright

1. clear;
   
2. density=2.7e3;
   
3. g=9.8;
   
4. Pa=1.0e5;
   
5. tension=1.0;
   
6. h=20;
   
7. viscosity=100.0e-3;
   
8. 
   
9. axis([0 10 0 20])%如何同一单位尺寸？
   
10. set(gcf,'DoubleBuffer','on');
    
11. lbefore=0;
    
12. deltt=1;
    
13. r=1e-4;
    
14. circle=inline('exp(i*linspace(0,pi*2,200))*r+x+y*i','r','x','y');
    
15. while(lbefore<h)
    
16. v=density*g*(r^2)/(3*viscosity);
    
17. llast=lbefore+1000*deltt*v;
    
18. R=((Pa+1000*2*tension/r+density*g*(h-lbefore)/1000)/...
    
19.   (Pa+1000*2*tension/r+density*g*(h-llast)/1000))^(1/3)*r;
    
20. lbefore=llast;
    
21. r=R;
    
22. z=circle(1000*r,5,lbefore);
    
23. plot(z);axis([0 10 0 20]);
    
24. pause(0.3);
    
25. end

复制代码

[ 本帖最后由 suffer 于 2006-10-9 20:45 编辑 ]

hustbill

<P>不哦错，还搞不懂如何实现的<BR></P>