几个MATALB的小程序 （收集的）

xieyubing734

MATLAB花心形
clc;
const=0;
x=-5:0.05:5;y=-5:0.05:5;z=-5:0.05:5;
[x,y,z]=meshgrid(x,y,z);
f=(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9/80)*y.^2.*z.^3-const;
p=patch(isosurface(x,y,z,f,0));
set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
daspect([1 1 1])
view(3)
camlight; lighting phong



地球模拟
ap=0;
load topo         
figure('colormap',topomap1,'Color',[.8 .8 .8]);        %        <6>
[x,y,z] = sphere(50);
azzy.DataAspectRatio = [1 1 1];azzy.PlotBoxAspectRatioMode = 'auto';
fa = axes('Visible','off', azzy);        
szzy1.AmbientStrength = 0.1;szzy1.DiffuseStrength = 1;
szzy1.SpecularColorReflectance = .5; szzy1.SpecularExponent = 20;
szzy1.SpecularStrength = 1;
surface(x,y,z,szzy1,'FaceLighting','phong','FaceColor','texture',...
       'EdgeColor','none','Cdata',topo,'Parent',fa);%        <13>
if ap==1,set(fa,'CameraViewAngle',0.1*get(fa,'CameraViewAngle'));end   
light('position',[-1 0 1],'color',[0.5 1 0.5]);
light('position',[-1.5 0.5 -0.5],'color',[.6 .2 .2]);
light('Position',[1.5 1.5 -1]);
light('Position',[0 -1.5 0],'color',[0.6 0.6 1]);
view([-17 26]) %变化这个view值 可以改变观看视角



%%% 螺旋面
clear all
u=linspace(0,2*pi);
v=linspace(0,6*pi);
[U,V]=meshgrid(u,v);
x=(4+sin(20*U)).*cos(V);
y=(4+sin(20*U)).*sin(V);
z=2*V;
mesh(x,y,z)
axis equal
%%%===================%%%
%%% 螺旋卷筒
clear all
R=2;
b=2;
u=linspace(0,2);
v=linspace(0,6*pi);
[U,V]=meshgrid(u,v);
x=U.*cos(V);
y=U.*sin(V);
z=U.*V;
mesh(x,y,z)
axis equal
axis off
%%%====================%%%
%%% 圆环面
clc
clear
u=linspace(0,2*pi,100);
v=linspace(0,2*pi,100);
[U,V]=meshgrid(u,v);
X=(4+sin(U)).*cos(V);
Y=(4+sin(U)).*sin(V);
Z=cos(U);
mesh(X,Y,Z)
axis equal
axis off

%%% 螺旋管，俗称“弹簧”
clc
clear
u=linspace(0,2*pi,100);
v=linspace(0,6*pi,100);
[U,V]=meshgrid(u,v);
X=(4+sin(U)).*cos(V);
Y=(4+sin(U)).*sin(V);
Z=cos(U)+V;
mesh(X,Y,Z)
axis equal
axis off

%%% 圆柱上的缎带
clc
clear
u=linspace(0,2);
v=linspace(0,6*pi);
[U,V]=meshgrid(u,v);
X=2*cos(V);
Y=2*sin(V);
Z=U+V;
mesh(X,Y,Z)
axis equal
axis off

%%% 纽带
clc
clear
u=linspace(0,2);
v=linspace(0,6*pi);
[U,V]=meshgrid(u,v);
X=U.*cos(V);
Y=U.*sin(V);
Z=5*V;
mesh(X,Y,Z)
axis equal
axis off

675977947

呵呵 谢谢楼主奉献 我试着玩玩

npu250

去玩玩，谢谢楼主！

sun150

不错的程序，谢谢楼主

dqlxn2008

不错的程序，谢谢楼主:)

wyc412721

不错的程序，谢谢楼主分享

zf5566

感谢分享呵呵

冰海深处

有意思，好好玩~~~

bainhome

有几个是振动论坛版主发的帖子，本来很有意思，以前原作者在simwe也常出现，不过让这么转来转去，连个出处都没有，再加上屁股后面跟一大堆没营养的白开水回帖，就透着十足的没品+无聊，招人烦。
版内好几个这样的主题，有等于没有，感觉还不如关掉。