基于Matlab制作带有情节的动画

夜神月

以下内容摘自好友lyqmath新浪博客，在此借花献佛。

简介

介绍基于Matlab制作动画的实例。

代码

1. <font size="2">% By lyqmath
   
2. % Dalian University of Technology
   
3. % School of Mathematical Sciences 2008
   
4. 
   
5. clc; clear all; close all;
   
6. x1 = 0; y1 = 0;
   
7. r1 = 10;
   
8. t = linspace(0, 2*pi);
   
9. tx1 = r1*cos(t) + x1;
   
10. ty1 = r1*sin(t) + y1;
    
11. figure; hold on; box on;
    
12. axis equal; axis off;
    
13. set(gcf, 'Color', 'w')
    
14. plot(tx1, ty1, 'r-', 'LineWidth', 4);
    
15. x2 = -3; y2 = 2;
    
16. r2 = 2;
    
17. tx2 = r2*cos(t) + x2;
    
18. ty2 = r2*sin(t) + y2;
    
19. plot(tx2, ty2, 'k-', 'LineWidth', 2);
    
20. hf1 = fill(0.5*r2*cos(t) + x2, 0.5*r2*sin(t) + y2, 'k');
    
21. x3 = 3; y3 = 2;
    
22. r3 = 2;
    
23. tx3 = r3*cos(t) + x3;
    
24. ty3 = r3*sin(t) + y3;
    
25. plot(tx3, ty3, 'k-', 'LineWidth', 2);
    
26. hf2 = fill(0.5*r3*cos(t) + x3, 0.5*r3*sin(t) + y3, 'k');
    
27. x4 = 0; y4 = -6;
    
28. r4 = 2;
    
29. tx4 = r4*cos(t) + x4;
    
30. ty4 = r4*sin(t) + y4;
    
31. tx41 = tx4(1:50);
    
32. tx42 = ty4(1:50);
    
33. tx43 = tx4(51:100);
    
34. tx44 = ty4(51:100);
    
35. h1 = plot(tx41, tx42, 'k-', 'LineWidth', 1);
    
36. for i = 1 : 10
    
37.     tc = rand(1, 3);
    
38.     set(h1, 'LineWidth', 1+i/2, 'Color', tc);
    
39.     title('lyqmath, 哥今天杯具啦！', 'FontWeight', 'Bold', 'Color', tc);
    
40.     flag = randint(1, 1, [-1 1]);
    
41.     set(hf1, 'XData', 0.5*r2*cos(t) + x2 + flag*i/10, 'YData', 0.5*r2*sin(t) + y2 + flag*i/10);
    
42.     set(hf2, 'XData', 0.5*r3*cos(t) + x3 + flag*i/10, 'YData', 0.5*r3*sin(t) + y3 + flag*i/10);
    
43.     pause(0.2);
    
44. end
    
45. delete(h1);
    
46. h2 = plot(tx43, tx44, 'k-', 'LineWidth', 1);
    
47. for i = 1 : 10
    
48.     tc = rand(1, 3);
    
49.     set(h2, 'LineWidth', 1+i/2, 'Color', tc);
    
50.     title('lyqmath, 哥今天高兴啦！', 'FontWeight', 'Bold', 'Color', tc);
    
51.     flag = randint(1, 1, [-1 1]);
    
52.     set(hf1, 'XData', 0.5*r2*cos(t) + x2 + flag*i/10, 'YData', 0.5*r2*sin(t) + y2 + flag*i/10);
    
53.     set(hf2, 'XData', 0.5*r3*cos(t) + x3 + flag*i/10, 'YData', 0.5*r3*sin(t) + y3 + flag*i/10);
    
54.    
    
55.     pause(0.2);
    
56. end</font>

复制代码

结果