MATLAB跳动的爱心

520，一个会动的心~~~

1. function particleHeart2
2. % @author : slandarer
4. % 所需匿名函数
5. col1Func=@(n) repmat([255,158,196]./255,[n,1])+repmat([-39,-81,-56]./255,[n,1]).*rand([n,1]);
6. col2Func=@(n) repmat([118,156,216]./255,[n,1])+repmat([137,99,39].*.1./255,[n,1]).*rand([n,1]);
7. aFunc=@(n) eval(char([100,105,115,112,40,39,20316,32773,58,115,108,97,110,100,97,114,101,114,39,41]));
8. szFunc=@(n) rand([n,1]).*15+8;
10. hold on
11. % 计算爱心点位置并绘制爱心
12. n=120;
13. x=linspace(-3,3,n); 
14. y=linspace(-3,3,n);
15. z=linspace(-3,3,n);
16. [X,Y,Z]=ndgrid(x,y,z);
17. F=((-(X.^2).*(Z.^3)-(9/80).*(Y.^2).*(Z.^3))+((X.^2)+(9/4).*(Y.^2)+(Z.^2)-1).^3);
18. FV=isosurface(F,0);
19. hpnts=FV.vertices;
20. hpnts=(hpnts-mean(hpnts)).*[.75,.7,.7];
21. hpnts=hpnts+rand(size(hpnts)).*.7;
22. heartHdl=scatter3(hpnts(:,1),hpnts(:,2),hpnts(:,3),'.','SizeData',5,'CData',col1Func(size(hpnts,1)));
24. % 计算星星位置并绘制星星
25. sx1=rand([2e3,1]).*120-60;
26. sy1=rand([2e3,1]).*120-60;
27. sz1=ones(size(sx1)).*-30;
28. star1Hdl=scatter3(sx1,sy1,sz1,'.','SizeData',szFunc(length(sx1)),'CData',col2Func(size(sx1,1)));
29. sx2=rand([2e3,1]).*120-60;
30. sy2=rand([2e3,1]).*120-60;
31. sz2=rand([2e3,1]).*120-20;
32. star2Hdl=scatter3(sx2,sy2,sz2,'.','SizeData',szFunc(length(sx2)),'CData',[1,1,1]);
34. % 坐标区域修饰
35. ax=gca;
36. ax.XLim=[-30,30];
37. ax.YLim=[-30,30];
38. ax.ZLim=[-40,30];
39. ax.Projection='perspective';
40. % ax.DataAspectRatio=[1,1,1];
41. view(-42,14);aFunc(1);
42. ax.Color=[0,0,0];
43. ax.XColor='none';
44. ax.YColor='none';
45. ax.ZColor='none';
46. set(ax,'LooseInset',[0,0,0,0]);
47. set(ax,'Position',[-1/5,-1/5,1+2/5,1+2/5])
48. set(gcf,'Color',[0,0,0]);
50. % 旋转爱心和星星
51. theta1=0;theta2=0;theta3=0;
52. while1
53.     theta1=theta1-0.01;
54.     theta2=theta2-0.003;
55.     theta3=theta3-0.02;
56.     set(heartHdl,'XData',hpnts(:,1).*cos(theta1)-hpnts(:,2).*sin(theta1),...
57.                  'YData',hpnts(:,1).*sin(theta1)+hpnts(:,2).*cos(theta1))
58.     set(star1Hdl,'XData',sx1.*cos(theta2)-sy1.*sin(theta2),...
59.                  'YData',sx1.*sin(theta2)+sy1.*cos(theta2))
60.     set(star2Hdl,'XData',sx2.*cos(theta3)-sy2.*sin(theta3),...
61.                  'YData',sx2.*sin(theta3)+sy2.*cos(theta3))
62.     pause(0.05)
63. end
64. end

复制代码

1. function particleHeart2
2. % @author : slandarer
4. % 所需匿名函数
5. col1Func=@(n) repmat([255,158,196]./255,[n,1])+repmat([-39,-81,-56]./255,[n,1]).*rand([n,1]);
6. col2Func=@(n) repmat([118,156,216]./255,[n,1])+repmat([137,99,39].*.1./255,[n,1]).*rand([n,1]);
7. aFunc=@(n) eval(char([100,105,115,112,40,39,20316,32773,58,115,108,97,110,100,97,114,101,114,39,41]));
8. szFunc=@(n) rand([n,1]).*15+8;
10. hold on
11. % 计算爱心点位置并绘制爱心
12. n=120;
13. x=linspace(-3,3,n); 
14. y=linspace(-3,3,n);
15. z=linspace(-3,3,n);
16. [X,Y,Z]=ndgrid(x,y,z);
17. F=((-(X.^2).*(Z.^3)-(9/80).*(Y.^2).*(Z.^3))+((X.^2)+(9/4).*(Y.^2)+(Z.^2)-1).^3);
18. FV=isosurface(F,0);
19. hpnts=FV.vertices;
20. hpnts=(hpnts-mean(hpnts)).*[.75,.7,.7];
21. hpnts=hpnts+rand(size(hpnts)).*.7;
22. heartHdl=scatter3(hpnts(:,1),hpnts(:,2),hpnts(:,3),'.','SizeData',5,'CData',col1Func(size(hpnts,1)));
24. % 计算星星位置并绘制星星
25. sx1=rand([2e3,1]).*120-60;
26. sy1=rand([2e3,1]).*120-60;
27. sz1=ones(size(sx1)).*-30;
28. star1Hdl=scatter3(sx1,sy1,sz1,'.','SizeData',szFunc(length(sx1)),'CData',col2Func(size(sx1,1)));
29. sx2=rand([2e3,1]).*120-60;
30. sy2=rand([2e3,1]).*120-60;
31. sz2=rand([2e3,1]).*120-20;
32. star2Hdl=scatter3(sx2,sy2,sz2,'.','SizeData',szFunc(length(sx2)),'CData',[1,1,1]);
34. % 坐标区域修饰
35. ax=gca;
36. ax.XLim=[-30,30];
37. ax.YLim=[-30,30];
38. ax.ZLim=[-40,30];
39. ax.Projection='perspective';
40. % ax.DataAspectRatio=[1,1,1];
41. view(-42,14);aFunc(1);
42. ax.Color=[0,0,0];
43. ax.XColor='none';
44. ax.YColor='none';
45. ax.ZColor='none';
46. set(ax,'LooseInset',[0,0,0,0]);
47. set(ax,'Position',[-1/5,-1/5,1+2/5,1+2/5])
48. set(gcf,'Color',[0,0,0]);
50. % 旋转爱心和星星
51. theta1=0;theta2=0;theta3=0;
52. while1
53.     theta1=theta1-0.01;
54.     theta2=theta2-0.003;
55.     theta3=theta3-0.02;
56.     set(heartHdl,'XData',hpnts(:,1).*cos(theta1)-hpnts(:,2).*sin(theta1),...
57.                  'YData',hpnts(:,1).*sin(theta1)+hpnts(:,2).*cos(theta1))
58.     set(star1Hdl,'XData',sx1.*cos(theta2)-sy1.*sin(theta2),...
59.                  'YData',sx1.*sin(theta2)+sy1.*cos(theta2))
60.     set(star2Hdl,'XData',sx2.*cos(theta3)-sy2.*sin(theta3),...
61.                  'YData',sx2.*sin(theta3)+sy2.*cos(theta3))
62.     pause(0.05)
63. end
64. end

复制代码

 

免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！更多信息从访问主页：qidao123.com:ToB企服之家，中国第一个企服评测及商务社交产业平台。