目录1.不同坐标系下的图形对比 (4)2.球曲面的法线 (4)3.浪花—山峰 (5)4.色彩斑斓的圆筒 (7)5.分层的不明物 (8)6.马鞍面 (9)7.螺旋线 (11)8.光芒四射---矢量图+等势线 (13)9.山谷—山峰【线性】 (14)10.山谷—山峰【面】 (16)11.牛顿环（动态的） (16)12.衍射调制下的双孔干涉条纹 (18)13.双缝衍射图 (20)14.沙丘 (21)15.漂亮的尤物 (22)16.圆花饼 (23)17.火红的心—尤物 (24)18.神秘的罗盘 (25)19.山峰分析图 (26)20.小球 (27)21.波浪、涟漪---像不像装鸡蛋的 (28)22.（动画的）山峰 (28)23.瀑布图 (29)24.三宝 (30)25.涟漪四视图 (31)26.3D腰带 (32)27.彩皮球 (33)28.山崩地裂 (34)29.飘荡的柔纱 (35)30.波纹 (36)31.相交的椭圆 (37)32.飘落的线 (38)33.跳动的正弦线 (39)34.磁滞回线—尤物 (40)35.复制的美 (40)36.评议扫描图 (42)37.旋转扫描图 (43)38.混沌吸引子 (43)39.动画演示混沌吸引分子形成 (44)40.绘制Julia集图形 (44)41.擦除动画实例——卫星绕地球运动(注释很详细) (46)42.擦除动画实例——太阳｜地球｜月亮｜卫星，绕转演示动画(注释很详细) (48)43.流行划过天空 (50)44.旋转的山峰 (50)45.抛物运动 (51)1.不同坐标系下的图形对比theta=0:pi/20:4*pi;phi=theta.^2-theta;[t,p]=meshgrid(theta,phi);r=t.*p;subplot(1,2,1);mesh(t,p,r);ylabel('x');xlabel('y');zlabel('z'); [x,y,z]=sph2cart(t,p,r);subplot(1,2,2);mesh(x,y,z);ylabel('x');xlabel('y');zlabel('z');2.球曲面的法线[x,y,z]=sphere;Surfnorm(x,y,z)3.浪花—山峰x=rand(100,1)*16-8;y=rand(100,1)*16-8;r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;xlin=linspace(min(x),max(x),33);ylin=linspace(min(y),max(y),33); [X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y);mesh(X,Y,Z);axis tight;hold on;ylabel('x');xlabel('y');zlabel('z'); plot3(x,y,z,’r’,’MarkerSize’,15)x=rand(1000,1)*16-8;y=rand(1000,1)*16-8;r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;xlin=linspace(min(x),max(x),99);ylin=linspace(min(y),max(y),99); [X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y);mesh(X,Y,Z);axis tight;hold on;ylabel('x');xlabel('y');zlabel('z'); plot3(x,y,z,'r','MarkerSize',30);去掉网格后：4.色彩斑斓的圆筒fnplt(fncmb(rsmak(‘cylinder’,1,2),[00-1;010;100])); axis equal,axis off;shading interp5.分层的不明物c=8.01*10^(-6);T=9.8*0.940*10^(-7);for m=1:100a=m/10000000;for n=1:100b=n/10000000;omiga(m,n)=sqrt((10.6667*a^3*b^5+10.6667*a^5*b^3)/(2.3112*a^5*b^5*c)*T) /(2*pi);endendmesh((1:100)/10000000,(1:100)/10000000,omiga);ylabel('半短轴b');xlabel('半长轴a');zlabel('频率f');6.马鞍面x=-8:8;y=-8:8;[X,Y]=meshgrid(x,y);Z=(X.^2./4.^2-Y.^2./5.^2);meshz(X,Y,Z);ylabel('x');xlabel('y');zlabel('z');x=-8:8;y=-8:8;[X,Y]=meshgrid(x,y);Z=X.^2./4.^2-Y.^2./5.^2;surfc(X,Y,Z);ylabel('x');xlabel('y');zlabel('z');7.螺旋线t=0:pi/60:10*pi;x=sin(t);y=cos(t);plot3(x,y,t,'*-b');ylabel('x');xlabel('y');zlabel('z');t=0:pi/60:10*pi;x=sin(t);y=cos(t);fill3(x,y,t,'*-b');ylabel('x');xlabel('y');zlabel('z');动画版：%by dynamic t=0:pi/50:10*pi;x=30*sin(t);y=30*cos(t);z=t;plot3(x,y,z);hold on%axis equalcomet3(x,y,z,0.5)8.光芒四射---矢量图+等势线[x,y]=meshgrid([-2:0.2:2]); f=x.^2+y.^2;[fx,fy]=gradient(f,0.2,0.2); quiver(x,y,fx,fy);hold oncontour(x,y,f);9.山谷—山峰【线性】[x,y]=meshgrid([-4:0.1:4]); z=x.*exp(-x.^2-y.^2);plot3(x,y,z);10.山谷—山峰【面】[x,y]=meshgrid([-2:0.05:2]);z=x.*exp(-x.^2-y.^2);surf(x,y,z);11.牛顿环（动态的）clearR=1;N=400;lamda=200e-9;rr=0.0017; [x,y]=meshgrid(linspace(-rr,rr,N)); r=sqrt(x.^2+y.^2);phi=r.^2/R./lamda.*pi.*2+pi;I=4.*cos(phi./2).^2;H=imshow(I);t=0;k=1;while k;s=get(gcf,'currentkey');if strcmp(s,'space');clc;k=0;endt=t+0.01;pause(0.2);phi=phi+t;I=4.*cos(phi./2).^2;set(H,'CData',I);end12.衍射调制下的双孔干涉条纹clear alllambda=500e-9;D=1;d=1e-6;a=3.5*d;ymax=10*lambda*D/a;xmax=ymax;N=101;xs=linspace(-xmax,xmax,N);ys=linspace(-ymax,ymax,N);for i=1:N;for j=1:N;r1(i,j)=sqrt(D.^2+(ys(i)-a./2).^2+xs(i).^2);r2(i,j)=sqrt(D.^2+(ys(j)-a./2).^2+xs(j).^2);Dr=r2-r1;phi(i,j)=Dr(i,j)./lambda.*pi.*2;m(i,j)=pi.*d.*sin(phi(i,j))./lambda;I(i,j)=(4.*cos(phi(i,j)./2).^2).*real((2.*besselj(1,m(i,j))./m(i,j)).^2) ;endendNC=225;Ir=I*NC;subplot(1,2,1);image(xs,ys,Ir);axis equalcolormap(gray(NC));subplot(1,2,2);plot(I,m),hold onclear alllambda=600e-9;D=1;d=1e-6;a=2.5*d;ymax=10*lambda*D/a;xmax=ymax;N=101;xs=linspace(-xmax,xmax,N);ys=linspace(-ymax,ymax,N);for i=1:N;for j=1:N;r1(i,j)=sqrt(D.^2+(ys(i)-a./2).^2+xs(i).^2);r2(i,j)=sqrt(D.^2+(ys(j)-a./2).^2+xs(j).^2);Dr=r2-r1;phi(i,j)=Dr(i,j)./lambda.*pi.*2;m(i,j)=pi.*d.*sin(phi(i,j))./lambda;I(i,j)=(4.*cos(phi(i,j)./2).^2).*real((2.*besselj(1,m(i,j))./m(i,j)).^2) ;endendNC=225;Ir=I*NC;colormap(gray(NC));plot(m,I),hold on13.双缝衍射图clear alllambda=500e-9;D=1;d=1e-6;a=3.5*d;ymax=10*lambda*D/a;xmax=ymax;N=101;xs=linspace(-xmax,xmax,N);ys=linspace(-10000*ymax,10000*ymax,N);for i=1:N;for j=1:N;r1(i,j)=sqrt(D.^2+(ys(i)-a./2).^2+xs(i).^2); r2(i,j)=sqrt(D.^2+(ys(j)-a./2).^2+xs(j).^2); Dr=r2-r1;phi(i,j)=Dr(i,j)./lambda.*pi.*2;I(i,j)=real(4.*cos(phi(i,j)./2).^2);endendNC=225;Ir=I/3.*NC;image(xs,ys,Ir);colormap(gray(NC));hold on14.沙丘figure('Position',[20031700696],'color','w');hold on;axis equal; axis offtitle('沙丘','fontsize',16,'color','k');text(180,480,'——万祎龙','fontsize',12,'color','k');i=1;for py=1:2:360c=py*4*pi/360;b=pi*cos(c);a=0:pi/30:6*pi;x=100*a/pi+20;y=10*sin(a+cos(a)*pi-b)*cos(c)+py;plot(x,y,'color','k');Z(i,:)=sin(a+cos(a)*pi-b)*cos(c);i=i+1;end15.漂亮的尤物%存脚本文件heart.m%%一个半径为R的圆周等分成N份，然后一等分点为圆心，以等分点到该圆垂直直径%%的距离为半径画圆。