【原创】一颗变化的心

qibbxxt

单击图形外侧，图形会变化，之前自己发在人人网上，看的人太少，就写在这里，

1. 
   
2. main.m
   
3. clear;clc;close all
   
4. n=0.1;
   
5. figure('menubar','none','name','heart','NumberTitle','off')
   
6. [xx,yy]=heart_11(n);
   
7. hp=fill(xx,yy,'r');
   
8. title('Heart','Fontsize',22,'Fontname','华文行楷')
   
9. set(gca,'Xtick',[],'Ytick',[]);box on;
   
10. axis([-20 20 -20 20])
    
11. xlabel({'\copyright Qibinbin','qibbxxt@126.com'});
    
12. set(gcf,'ButtonDownFcn',['n=n+0.1;[xx,yy]=heart_11(n);',...
    
13.         'set(hp,''xdata'',xx,''ydata'',yy);'...
    
14.         'if n>4 set(hp,''facecolor'',''k'');end;'...
    
15.         'if n>5.3 close; end']);
    
16. heart_11.m
    
17. function [xx,yy]=heart_11(n)
    
18. if nargin==0
    
19.     n=1;
    
20. end
    
21. t=-pi:0.01:90*pi/180;
    
22. a=5*n;
    
23. b=3*n;
    
24. theta=45*pi/180;
    
25. x=a*cos(t);
    
26. y=b*sin(t);
    
27. xy=zeros(2,length(x));
    
28. for i=1:length(x)
    
29. xy(:,i)=expm([0,-theta;theta,0])*[x(i);y(i)];
    
30. end
    
31. xy(:,xy(1,:)<0)=[];
    
32. xy0=fliplr([-xy(1,:);xy(2,:)]);
    
33. xy=[xy,xy0];
    
34. xx=xy(1,:);
    
35. yy=xy(2,:);
    

复制代码