MATLAB.rar_heart simulation_心跳

clc clear all close all x = linspace(-2,2,100); [X,Y,Z] = meshgrid(x,x,x); I1 = (X.^2+9/4*Y.^2+Z.^2-1).^3-X.^2.*Z.^3-9/80*Y.^2.*Z.^3; p = patch(isosurface(X,Y,Z,I1,0)); set(p, 'FaceColor', 'red', 'EdgeColor', 'none'); view(3); axis equal ; axis off; light('Posi',[0 -2 3]); % 在(0,-2,3)点处建立一个光源 lighting phong set(gca,'nextplot','replacechildren'); % 记录电影 XX = get(p,'XData'); YY = get(p,'YData'); ZZ = get(p,'ZData'); for j = 1:20 bili = sin(pi*j/20); set(p,'XData',bili*XX,'YData',bili*YY,'ZData',bili*ZZ) F(j) = getframe; end % 放映10次 for i=1:10 movie(F,1) sound('heart.mp3',2000); end