[转帖]心形线动态生成演示

suffer

1. function [xcar,ycar]=drawcardioid(varargin)
   
2. % A demonstration to show the form of the cardioid
   
3. 
   
4. error(nargchk(0,3,nargin));
   
5. nin = nargin;
   
6. if nin==0
   
7. cx1=0; cy1=0; % The coordinate of the center point.
   
8. r=1; % The radius of the circle.
   
9. elseif nin==1
   
10. cx1=0; cy1=0;
    
11. r=varargin{1};
    
12. elseif nin==2
    
13. cx1=0;
    
14. cy1= vararin{1};
    
15. r = varargin{2};
    
16. else
    
17. cx1=vararin{1};
    
18. cy1=vararin{2};
    
19. r=vararin{3};
    
20. end
    
21. 
    
22. if nargout==0
    
23. flag=1;
    
24. else
    
25. flag=0;
    
26. end
    
27. 
    
28. theta = linspace(0,2*pi,120);
    
29. x0 = r*cos(theta);
    
30. y0 = r*sin(theta);
    
31. x1 = x0+cx1;
    
32. y1 = y0+cy1;
    
33. cx2 = 2*x0+cx1;
    
34. cy2 = 2*y0+cy1;
    
35. x3=x0.*cos(theta)-y0.*sin(theta);
    
36. y3=x0.*sin(theta)+y0.*cos(theta);
    
37. for k=1:120
    
38. x2(k,:) = cx2(k)+x0;
    
39. y2(k,:) = cy2(k)+y0;
    
40. xx(k)=cx2(k)+x3(k);
    
41. yy(k)=cy2(k)+y3(k);
    
42. end
    
43. xcar = xx;
    
44. ycar = yy;
    
45. if flag
    
46. plot(x1,y1);
    
47. axis([cx1-3.2*r cx1+3.2*r cy1-3.2*r cy1+3.2*r]);
    
48. axis manual
    
49. hold on
    
50. daspect([1 1 1]);
    
51. set(gcf,'doublebuffer','on');
    
52. for k = 1:120
    
53. plot(x2(k,:),y2(k,:),'g',xx(k),yy(k),'r*');
    
54. pause(0.05);
    
55. %plot(x2,y2,'color','w');
    
56. end
    
57. end

复制代码

来自山城棒棒儿

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