一起来看一下matlab程序小彩蛋spy %显示二哈
xpbombs %扫雷小游戏
teapotdemo %水壶渲染图像
vibes %曲面波动动图
logo %显示MATLAB-logo
travel %旅行问题
why %tell me why!
truss %弯曲桁架
life %人生小游戏
fifteen %滑动拼图
xpquad %超二次曲面
wrldtrv %环球旅行
earthmap %世界地图
cplxdemo %复杂曲面
lorenz %洛伦兹吸引子
%彩蛋-彩蛋-彩蛋，图不适宜摆放，自行出图
% 画爱心
[X, Y] = meshgrid(0.01:0.01:1, 0.01:0.01:1);
Zfun =@(x,y)12.5*x.*log10(x).*y.*(y-1)+exp(-((25 …
*x – 25/exp(1)).^2+(25*y-25/2).^2).^3)./25;
Z = Zfun(X,Y);
figure;
surf(Y,Z,X,’FaceColor’,[1 0.75 0.65],’linestyle’,’none’);
hold on
surf(Y+0.98,Z,X,’FaceColor’,[1 0.75 0.65],’linestyle’,’none’);
axis equal;
view([116 30]);
camlight;
lighting phong; % 设置光照和光照模式
五星红旗% 五星红旗
function varargout=saxplaxliu(varargin)
x=[0 0 40 40 0];
y=[0 13 13 0 0];
k=40/13;
fill(x,y,’r’)
hold on
plot(x,y,’r’)
axis([0 40 0 13])
seta=pi/6;
R=1.5;R2=0.55;x0=7;y0=9;seta=0;
fliu(R,R2,x0,y0,seta,k)
R=0.6;R2=0.23;x0=15;y0=11.6;seta=-pi/20;
fliu(R,R2,x0,y0,seta,k)
R=0.6;R2=0.23;x0=20;y0=10;seta=pi/20;
fliu(R,R2,x0,y0,seta,k)
R=0.6;R2=0.23;x0=19;y0=7;seta=pi/25;
fliu(R,R2,x0,y0,seta,k)
R=0.6;R2=0.23;x0=16;y0=5;seta=-pi/10;
fliu(R,R2,x0,y0,seta,k)
axis off
function fliu(R,R2,x0,y0,seta,k)
x=zeros(1,10);
y=x;
for i=1:5
x(2*(i-1)+1)=R*cos(pi/2+2*pi/5*(i-1)+seta);
x(2*i)=R2*cos(pi/2+2*pi/5*(i-1)+pi/5+seta);
y(2*(i-1)+1)=R*sin(pi/2+2*pi/5*(i-1)+seta);
y(2*i)=R2*sin(pi/2+2*pi/5*(i-1)+pi/5+seta);
end
x(11)=x(1);
y(11)=y(1);
x=x*k+x0;
y=y+y0;
fill(x,y,[1 0.7 0])
plot(x,y,’y’)
心形线%%
subplot(2,3,1)
ezplot(‘x^2+(y-(x^2)^(1/3))^2=9′);
%%
subplot(2,3,2)
ezplot(‘-x^2*y^3+(x^2+y^2-1)^3=0′,[-1.5,1.5])
%%
subplot(2,3,3)
ezplot(’17*x.^2-16*abs(x).*y+17*y.^2=200′) ;
%%
subplot(2,3,4)
x=linspace(-2,2,1000);
y1=sqrt(2*sqrt(x.^2)-x.^2);
y2=-2.14*sqrt(sqrt(2)-sqrt(abs(x)));
plot(x,y1,’b’,x,y2,’b’);
axis([-2.5,2.5,-3,1.5]);
%%
subplot(2,3,5)
f=@(x,y,z)(x.^2+(9./4).*y.^2+z.^2-1).^3-x.^2.*z.^3-(9./80).*y.^2.*z.^3;
[x,y,z]=meshgrid(linspace(-3,3));
val=f(x,y,z);
[p,v]=isosurface(x,y,z,val,0);
patch(‘faces’,p,’vertices’,v,’facevertexcdata’,jet(size(v,1)),’facecolor’,’w’,’edgecolor’,’flat’);
view(3);grid on;axis equal;
%%
subplot(2,3,6)
f=@(x,y,z)(x.^2+ (9./4).*y.^2 + z.^2 – 1).^3 – x.^2.*z.^3 – (9./80).*y.^2.*z.^3;
[x,y,z]=meshgrid(linspace(-1.5,1.5));
val=f(x,y,z);
isosurface(x,y,z,val,0);
axis equal;view(3);colormap([1 0.2 0.2])
躁动的心gcp=figure;
M=moviein(100,gcp);
heartSpeed=15;
radiusRate=3;
for i = 1:100
axis([-20 20 -20 20]);
cla(gcp);
x=[];
y=[];
for t = 0:0.01:2*pi
r=(sin(t)*sqrt(abs(cos(t))))/(sin(t)+7/5)-2*sin(t)+2;
r = radiusRate * r;
x=[x,r*cos(t)*(0.5+0.5*HeartRate(i,heartSpeed))];
y=[y,r*sin(t)*(0.5+0.5*HeartRate(i,heartSpeed))];
end
hold on;
fill(x,y,’r’);
function ratio=HeartRate(i,heartSpeed)
i=mod(i,heartSpeed);
if i < heartSpeed / 2
ratio=sin( i/heartSpeed*pi);
elseif i >= heartSpeed / 2 && i < heartSpeed * 3/4
ratio=2*cos(i/heartSpeed*2/3*pi);
else
ratio = 0;
end
部分代码源自网络仅供娱快学习