用matlab绘制的漂亮图形

合集下载

陈泽灏 MATLAB巧绘七彩玫瑰线

陈泽灏 MATLAB巧绘七彩玫瑰线

-0.1
0 向日葵
0.1
0.2
0.3
0.4
0.5
荷花 2.5
2
1.5
1
0.5
0
-0.5
-1 -2
-1.5
-1
-0.5
0
0.5
1
1.5
2
狗尾巴草 2.5
2
1.5
1
0.5
0
-0.5
-1
0.6
-2
-1
0
1
2
0.3
0.4
0.2
0.2
0.1
0
0
-0.2ቤተ መጻሕፍቲ ባይዱ
-0.1
-0.4
-0.2
-0.6
-0.3
-0.8 -0.8
0.5
0.5
270 n=4
300
240
270 n=5
300
图2
对比图 1 和图 2,可以看出,当把 限制在 [0, ] 内时,玫瑰线就已经形成了 n 个花瓣,当 t 在 [ ,2 ] 内时,又新增了 n 个花瓣;只是若 n 为偶数,新增的 n 个花瓣不与原来的花瓣重合;若 n 为奇数,新增的 n 个花瓣与原来的花瓣重合罢 了。 进一步,令 n=3,我们在 [0, / 3],[ / 3,2 / 3],[2 / 3, ] 三个区间里分别用红、 绿、蓝三种颜色画玫瑰线,结果如图 3.
MATLAB 巧绘七彩玫瑰线
西工大附中高 2013 届 13 班 陈泽灏
摘要 研究了玫瑰线的性质及 matlab 的作图功能。 在认识玫瑰线性质的基础 上,以 matlab 为工具,用玫瑰线或其一部分作为基本线条,配以心形线以及直 线,画出了各种漂亮的图案,从而对数学的实用性有了更加深刻的认识。 关键词 玫瑰线;matlab;性质 先看看这些漂亮的图案吧!让我来告诉你,它们是用 matlab 画的,而且过 程非常简单!因为所有的这些图案虽然各不相同,但它们都是用一种线条画出来 的!它就是------ 玫瑰线。对,玫瑰线! 1.什么是玫瑰线 什么是玫瑰线呢?玫瑰线是一族平面曲线,又叫格拉金曲线,意大利数学 家格拉金是第一个研究玫瑰线的人(1728 年) 。玫瑰线的极坐标方程为

Matlab绘图系列之高级绘图

Matlab绘图系列之高级绘图

Matlab绘图系列之高级绘图一、目录1.彗星图二维彗星图三维彗星图2.帧动画3.程序动画4.色图变换5.V oronoi图和三角剖分V oronoi图三角剖分6.四面体7.彩带图彩带图三维流彩带图8.伪彩图9.切片图切片图切片轮廓线图10.轮廓图显示轮廓线显示围裙瀑布效果带光照模式的阴影图11.函数绘图轮廓线、网格图、曲面图、轮廓网格图轮廓曲面图、二维曲线、极坐标曲线图、自定义函数12.三维图形控制视点灯光效果色彩控制二、图形示例1.彗星图二维彗星图t=0:.01:2*pi;x=cos(2*t).*(cos(t).^2);y=sin(2*t).*(sin(t).^2);comet(x,y);title('二维彗星轨迹图')hold onplot(x,y)三维彗星图a=12;b=9;T0=2*pi;%T0是轨道的周期T=5*T0;dt=pi/100;t=[0:dt:T]';f=sqrt(a^2-b^2);%地球与另一焦点的距离th=12.5*pi/180;%未经轨道与x-y平面的倾角E=exp(-t/20);%轨道收缩率x=E.*(a*cos(t)-f);y=E.*(b*cos(th)*sin(t));z=E.*(b*sin(th)*sin(t));plot3(x,y,z,'g')%画全程轨线hold on,sphere(20);%画地球axis offtitle('卫星返回地球示例')x1=-18*T0;x2=6*T0;y1=-12*T0;y2=12*T0;z1=-6*T0;z2=6*T0;axis([x1 x2 y1 y2 z1 z2])% axis([-15 10 -15 10 -10 10])axis equalcomet3(x,y,z,0.02);%画运动轨线hold off2.帧动画Z=peaks;surf(Z)%绘制网格表面图axis tightset(gca,'nextplot','replacechildren');%设定axis覆盖重画模式title('帧动画播放示例')for j=1:20surf(sin(2*pi*j/20)*Z,Z)%重新绘制网格表面图,这里后面一个Z当成了颜色矩阵F(j)=getframe;%创建帧endmovie(F,20)%播放动画20次3.程序动画t=0:pi/50:10*pii=1;h=plot3(sin(t(i)),cos(t(i)),t(i),'*','erasemode','none');%设定擦除模式grid onaxis([-2 2 -2 2 -1 10*pi])title('程序动画示例')for i=2:length(t)set(h,'xdata',sin(t(i)),'ydata',cos(t(i)),'zdata',t(i));drawnowpause(0.01)end4.色图变换load spineimage(X)colormap coolspinmap(10)5.V oronoi图和三角剖分V oronoi图rand('state',5)x=rand(1,10);y=rand(1,10);subplot(131)voronoi(x,y);%绘制voronoi图形axis equalaxis([-0.2 1.6 -0.5 2.5])subplot(132)[vx,vy]=voronoi(x,y);plot(x,y,'r+',vx,vy,'b-');%应用返回值绘制axis equalaxis([-0.2 1.6 -0.5 2.5])subplot(133)rand('state',5);x=rand(10,2);[v,c]=voronoin(x);%返回值v参数维voronoi顶点矩阵,返回值c参数为voronoi元胞数组for i=1:length(c)if all(c{i}~=1)patch(v(c{i},1),v(c{i},2),i);%应用色图iendendaxis equalaxis([-0.2 1.6 -0.5 2.5])box on三角剖分[x,y]=meshgrid(1:15,1:15);tri=delaunay(x,y);z=peaks(15);trimesh(tri,x,y,z)6.四面体d=[-1 1];[x,y,z]=meshgrid(d,d,d);%定义一个立方体x=[x(:);0];y=[y(:);0];z=[z(:);0];%[x,y,z]分别为加上中心的立方体顶点X=[x(:) y(:) z(:)];Tes=delaunayn(X);%返回m×n的数组值tetramesh(Tes,X);%绘制四面体图camorbit(20,0);%旋转摄像目标位置7.彩带图彩带图[x,y]=meshgrid(-3:.5:3,-3:.1:3);z=peaks(x,y);ribbon(y,z)三维流彩带图load wind%打开保存的数据lims=[100.64 116.67 17.25 28.75 -0.02 6.86];%定义坐标轴范围[x,y,z,u,v,w]=subvolume(x,y,z,u,v,w,lims);%lims来定义数据子集[sx sy sz]=meshgrid(110,20:5:30,1:5);%定义网格点verts=stream3(x,y,z,u,v,w,sx,sy,sz,.5);%计算彩带顶点cav=curl(x,y,z,u,v,w);%计算卷曲角速度wind_speed=sqrt(u.^2+v.^2+w.^2);%计算流速h=streamribbon(verts,x,y,z,cav,wind_speed,2);%绘制流彩带图view(3)8.伪彩图n=6%定义轮数r=(0:n)'/n;%定义轮的半径theta=pi*(-n:n)/n;%定义轮的扇区角X=r*cos(theta);Y=r*sin(theta);%定义网格顶点C=r*cos(2*theta);%定义色图pcolor(X,Y,C)%绘制伪彩图axis equal tight 9.切片图切片图[x,y,z] = meshgrid(-2:.2:2,-2:.25:2,-2:.16:2);v = x.*exp(-x.^2-y.^2-z.^2);xslice = [-1.2,.8,2]; yslice = 2; zslice = [-2,0];slice(x,y,z,v,xslice,yslice,zslice)colormap hsv切片轮廓线图[x y z v]=flow;%打开水流数据h=contourslice(x,y,z,v,[1:9],[],[0],linspace(-8,2,10));%切片轮廓线view([-12 30])10.轮廓图显示轮廓线[x,y,z]=peaks;subplot(1,2,1)meshc(x,y,z);%同时画出网格图与轮廓线title('meshc 网格图与轮廓线')axis([-inf inf -inf inf -inf inf]);subplot(1,2,2)surfc(x,y,z);%同时画出曲面图与轮廓线title('surfc 曲面图与轮廓线')axis([-inf inf -inf inf -inf inf]);显示围裙[x y z]=peaks;meshz(x,y,z);瀑布效果[X,Y,Z]=peaks(30);waterfall(X,Y,Z)带光照模式的阴影图[x,y]=meshgrid(-3:1/8:3);z=peaks(x,y);surfl(x,y,z);shading interp%着色处理colormap(gray);%灰度处理axis([-3 3 -3 3 -8 8])11.函数绘图轮廓线、网格图、曲面图、轮廓网格图%图1绘制轮廓线、网格图、曲面图、轮廓网格图subplot(221)f=['3*(1-x)^2*exp(-(x^2)-(y+1)^2)-10*(x/5-x^3-y^5)*exp(-x^2-y^2)-1/3*exp(-(x+1)^2-y^2)'];%定义双变量x、y的函数式ezcontour(f,[-3,3],49)%x、y为[-3 3],网格为49×49 subplot(222)ezmesh('sqrt(x^2+y^2)');subplot(223)ezsurf('real(atan(x+i*y))')%经过滤波,如果相同数据surf绘图没有滤波subplot(224)ezmeshc('y/(1+x^2+y^2)',[-5,5,-2*pi,2*pi])%x、y的数值范围分别为[-5 5]、[-2*pi 2*pi]轮廓曲面图、二维曲线、极坐标曲线图、自定义函数%图2绘制轮廓曲面图、二维曲线、极坐标曲线图、自定义函数figure(2)subplot(221)ezsurfc('sin(u)*sin(v)')subplot(222)ezplot('x^2-y^4');subplot(223)ezpolar('1+cos(t)')subplot(224)fplot('myfun',[-20 20])function Y=myfun(x)Y(:,1)=200*sin(x(:))./x(:);Y(:,2)=x(:).^2;三维曲线图%绘制三维曲线图figure(3)ezplot3('sin(t)','cos(t)','t',[0,6*pi])12.三维图形控制视点View图形旋转subplot(121)surf(peaks);title('旋转前图形');subplot(122)h=surf(peaks);rotate(h,[1 0 1],180)title('旋转后图形');灯光效果%灯光效果(1)camlight(2)light(3)lightangle(4)lighting(5)materialsphere;camlight色彩控制%色彩控制(1)缺省设置colordef、whitebg(2)色图colormap(3)浓淡处理shadingload flujetimage(X)colormap(jet)subplot(131)sphere(16)axis squareshading flattitle('Flat Shading')subplot(132)sphere(16)axis squareshading facetedtitle('Faceted Shading')subplot(133)sphere(16)axis squareshading interptitle('Interpolated Shading')。

四讲Matlab绘图ppt课件精品文档

四讲Matlab绘图ppt课件精品文档

p l o t (. . ., s t r ) 使用字符串s t r指定的颜
色和线型进行绘图。表1 中列出了s t r可
以取的值。
2019/10/12
4
一、二维图形:表一


.

^
正三角
*
星号
v
倒三角
OO
s
正方形
++
d
菱形
X 乘号
p
五角星
<
左三角
h
六角星
>
右三角
none 无点
2019/10/12
5
2019/10/12
26
三、三维图形
m e s h ( X,Y,Z ) 将矩阵Z中的各个元素作 为矩形网格上的高度,对这些值绘图,
并且将相邻的点连接形成三维网格表面 图。颜色由高度,即Z中的元素指定。
s u r f ( X , Y , Z) 绘制出由坐标(Xi j, Yi j, Zi j)确定的表面图形。如果X和Y分别是长 度为m和n的向量,那么, Z必须为m×n 的矩阵,并且表面是由(Xj, Yi, Zi j)来定 义的。
一、二维图形:表一
线型
颜色
--. : none
实线 m 虚线 b 点划线 c 点线 w 无线 r
品红色 蓝色 灰色 白色 红色
颜色
k
g
绿色 y
黑色 黄色
2019/10/12
6
一、二维图形
1.2彗星图形
c o m e t ( x , y ) 绘制向量y对向量x的彗星 轨线。如果只给出一个向量,则用该向 量对其下标值绘图。
2019/10/12
10

用matlab绘制的漂亮图形

用matlab绘制的漂亮图形

用matlab绘制的漂亮图形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),[0 0 -1;0 1 0;1 0 0])); 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');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 onlambda=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.lambda=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',[200 31 700 696],'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份,然后一等分点为圆心,以等分点到该圆垂直直径%%的距离为半径画圆。

matlab画图

matlab画图
星号 小黑点
叉号 小正方形
菱形
标记符 ^ v > < p h
方式 向上三角形 向下三角形 向右三角形 向左三角形
五角星 六角星
7/ 38
fplot命令的使用
fplot(function,limits,LineSpec) 在limits指定的范围内绘制function函数的 图形. limits是一个矢量, 指定x轴上的范围[xmin, xmax], 或者x轴和y轴上的范围[xmin xmax ymin ymax]. function必须是M文件函数的名称或句柄, 或者含有变量x的字符串. LineSpec指定属性进行绘图.
meshz: 除了生成网格曲面外,
还在曲面下面加一长方形台柱。
meshz(Z);
colormap([0 0 1]) %设定色图参数
25/ 38
surf: 对网线图的网格块区域着色得到刻画面. surfc: 和meshc类似,在刻面图下方绘上等值线. surfl: 对刻画面中单元颜色进行平滑处理, 得到
27/38
饱和色
[0 0 0] — 黑色 [0 0 1] — 兰色 [0 1 0] — 绿色 [0 1 1] — 浅兰 [1 0 0] — 红色 [1 0 1] — 粉红 [1 1 0] — 黄色 [1 1 1] — 白色
调和色 [0.5 0.5 05] — 灰色 [0.5 0 0] — 暗红色 [1 0.62 0.4] — 铜色 [0.49 1 0.8] — 浅绿 [0.49 1 0.83] —宝石兰
例14. [X,Y]=meshgrid(-2:.2:2,-2:.2:2); Z=X.*exp(-X.^2-Y.^2); mesh(X,Y,Z)
24/38

第五章matlab绘图

第五章matlab绘图
例5-10 在图形窗口中,以子图形式同时绘制 多根曲线。
2024/3/28
21
5.2 其他二维图形
5.2.1 其他坐标系下的二维数据曲线图 1. 对数坐标图形 MATLAB提供了绘制对数和半对数坐标 曲线的函数,调用格式为: semilogx(x1,y1,选项1,x2,y2,选项2,…) semilogy(x1,y1,选项1,x2,y2,选项2,…) loglog(x1,y1,选项1,x2,y2,选项2,…)
2024/3/28
31
2. 绘制三维曲面的函数 surf函数和mesh函数的调用格式为:
mesh(x,y,z,c)
surf(x,y,z,c) 一般情况下,x,y,z是维数相同的矩阵。x,y 是网格坐标矩阵,z是网格点上的高度矩阵, c用于指定在不同高度下的颜色范围。
例5-17 绘制三维曲面图z=sin(x+sin(y))-x/10。
plotyy(x1,y1,x2,y2)
其中x1,y1对应一条曲线,x2,y2对应另一 条曲线。横坐标的标度相同,纵坐标有两 个,左纵坐标用于x1,y1数据对,右纵坐标 用于x2,y2数据对。
2024/3/28
8
例5-4 用不同标度在同一坐标内绘制曲线 y1=0.2e-0.5xcos(4πx) 和 y2=2e-0.5xcos(πx)。
程序如下: x=0:pi/100:2*pi; y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); plotyy(x,y1,x,y2);
2024/3/28
9
4. 图形保持
hold on/off命令控制是保持原有图形还是 刷新原有图形,不带参数的hold命令在两 种状态之间进行切换。

matlab画立体蝴蝶的代码

matlab画立体蝴蝶的代码

matlab画立体蝴蝶的代码Matlab是一款强大的工具,可以帮助我们完成许多数据分析、数学建模和图像处理等工作。

在Matlab中,我们可以用简单的代码绘制出许多精美的图像,例如立体蝴蝶。

下面就来介绍一下如何使用Matlab绘制立体蝴蝶的代码。

1.准备工作我们需要在Matlab中打开一个新的文件,然后输入下面这句代码创建一个3D画布:figure('units','normalized','outerposition',[0 0 1 1])这句代码将创建一个全屏的3D画布,供我们后续用来绘制立体蝴蝶。

2.绘制蝴蝶的翅膀接下来,我们需要定义蝴蝶的两个翅膀,具体代码如下:t = 0:pi/10:2*pi;x = sin(t).*(exp(cos(t))-2*cos(4*t)-sin(t/12).^5);y = cos(t).*(exp(cos(t))-2*cos(4*t)-sin(t/12).^5);z = sin(t/2).^2.*cos(5*t+pi/2);plot3(x,y,z,'r','LineWidth',2)这段代码的意思是先设置一个参数t,然后用该参数计算出x、y、z坐标分别对应蝴蝶翅膀上的点的位置。

最后使用plot3函数将这些点连接起来,形成一个立体的蝴蝶翅膀。

这时我们已经完成了蝴蝶的一只翅膀的绘制。

接下来我们需要将另一只翅膀按照一定的规律旋转、平移后绘制出来。

代码如下:hold onfor i=1:3xx = x*cos(pi*i/3) + y*sin(pi*i/3) - 5*(cos(pi*i/3)+1)/2; yy = -x*sin(pi*i/3) + y*cos(pi*i/3);zz = z+2*i;plot3(xx,yy,zz,'r','LineWidth',2)end这段代码首先使用for循环来遍历翅膀的三个位置。

MATLAB曲线绘制大全

MATLAB曲线绘制大全

一、二维数据曲线图1.1绘制单根二维曲线plot函数的基本调用格式为:plot(x,y)其中x和y为长度相同的向量,分别用于存储x坐标和y坐标数据。

例1-1在0x2p区间内,绘制曲线y=2e-0.5xcos(4x)程序如下:x=0:pi/100:2*pi;y=2*exp(-0.5*x).*cos(4*pi*x);plot(x,y)例1-2绘制曲线。

程序如下:t=0:0.1:2*pi;x=t.*sin(3*t);y=t.*sin(t).*sin(t);plot(x,y);plot函数最简单的调用格式是只包含一个输入参数:plot(x)在这种情况下,当x是实向量时,以该向量元素的下标为横坐标,元素值为纵坐标画出一条连续曲线,这实际上是绘制折线图。

1.2绘制多根二维曲线1.plot函数的输入参数是矩阵形式(1)当x是向量,y是有一维与x同维的矩阵时,则绘制出多根不同颜色的曲线。

曲线条数等于y矩阵的另一维数,x被作为这些曲线共同的横坐标。

(2)当x,y是同维矩阵时,则以x,y对应列元素为横、纵坐标分别绘制曲线,曲线条数等于矩阵的列数。

(3)对只包含一个输入参数的plot函数,当输入参数是实矩阵时,则按列绘制每列元素值相对其下标的曲线,曲线条数等于输入参数矩阵的列数。

当输入参数是复数矩阵时,则按列分别以元素实部和虚部为横、纵坐标绘制多条曲线。

2含多个输入参数的plot函数调用格式为:plot(x1,y1,x2,y2,,xn,yn)(1)当输入参数都为向量时,x1和yl,x2和y2,,xn和yn分别组成一组向量对,每一组向量对的长度可以不同。

每一向量对可以绘制出一条曲线,这样可以在同一坐标内绘制出多条曲线。

(2)当输入参数有矩阵形式时,配对的x,y按对应列元素为横、纵坐标分别绘制曲线,曲线条数等于矩阵的列数。

例1-3分析下列程序绘制的曲线。

x1=linspace(0,2*pi,100);x2=linspace(0,3*pi,100);x3=linspace(0,4*pi,100);y1=sin(x1);y2=1+sin(x2);y3=2+sin(x3);x=[x1;x2;x3]';y=[y1;y2;y3]';plot(x,y,x1,y1-1)3.具有两个纵坐标标度的图形在MATLAB中,如果需要绘制出具有不同纵坐标标度的两个图形,可以使用plotyy绘图函数。

十节MATLAB绘图课件

十节MATLAB绘图课件
怎样画具有不同半径旳同心圆?
MATLAB绘图
函数相对另一函数旳图形绘制 theta=linspace(0,2*pi,101); plot(sin(theta),sin(2*theta+pi/4)) % 里萨如图形 axis off %隐去坐标轴
MATLAB绘图
曲线族旳绘制
例, 抛物线族 y =a^2-x^2, -5<x<5, a=1,2…5.
Fplot 绘制函数在指定函数区间旳图形 fplot(fun,lims,tol) %lims=[xmin,xmax]
tol 默以为0.002. fplot('sin(1 / x)', [0.01 0.1],1e-3)
MATLAB绘图
fplot(fun,lims,tol) % fun 必须是m文件或者独立变 量为x旳字符串,假如fun为[f1(x),f2(x),f3(x)],则 当输入[x1,x2]后,返回旳是矩阵 [f1(x1),f2(x1),f3(x1)
MATLAB绘图
图形注释和可视化提升 1. 使用坐标轴,图形标题,曲线标注,填充
区域和添加文本 2. 变化轴、曲线和文本旳属性 3. 使用希腊字母、数学符号、上下标等
MATLAB绘图
例 x=0:0.05:6; plot(x,cos(x),'k',x,1./cosh(x),'k',[4.73 4.73],[-1
1],'k'); xlabel('x'); ylabel('Value of Function') title('Visualization of two intersectting curves') text(4.8,-.5,'x=4.73','fontname','times','fontsize',14) text(2.1,.3,'1/cosh(x)','fontsize',16) text(1.2,-.4,'cos(x)','fontsize',16,'fontname','times') set(gca,'fontsize',14,'LineWidth',2) PropertyName={'LineWidth','LineWidth','LineWidth'}; PropertyValue={2.5,2.5,2.5;7,7,7;1,1,1}; set(h,PropertyName,PropertyValue) [legendhandle

利用Matlab绘制有趣图像的示例代码

利用Matlab绘制有趣图像的示例代码

利⽤Matlab绘制有趣图像的⽰例代码⽬录1.随机樱花树2.苹果绘制3.南⽠绘制4.⼀堆三⾓形绘制5.⽉饼绘制6.⼤钻⽯绘制7.有趣曲线18.有趣曲线29.有趣曲线3——蝴蝶曲线10.有趣曲线4——⼼形曲线11.有趣曲线512.会害羞的含羞草13.随机雪景1.随机樱花树function sakura% @author:slandarer% 随机形状⼆叉树樱花树绘制hold on,axis equalaxis(0.5+[-10,50,0,50])set(gca,'xtick',[],'ytick',[],'xcolor','w','ycolor','w')set(gca,'color',[0.5020 0.5020 0.5020])length_trunk=6;width_trunk=4;k1=0.9;k2=0.8;number_branch=15;alp=pi/10;length_branch=k1*length_trunk;width_branch=k2*width_trunk;trunk=[12,0;12,length_trunk];plot(trunk(:,1),trunk(:,2),'color',[0 0 0],'Linewidth',width_trunk)begins=[trunk(2,:),pi/2,1];grow=begins;plotdata=[0 0 0 0 0 0 0 0];plotdata(1,:)=[];for i=1:number_branchcontrol=randi(25,[length(grow(:,1)),1])>=10;ag=grow(:,3);l=length(ag);parta=[length_branch.*k1.^grow(:,4).*cos(ag+ones(l,1)*alp),length_branch.*k1.^grow(:,4).*sin(ag+ones(l,1)*alp),ones(l,1)*alp,ones(l,1)]; partb=[length_branch.*k1.^grow(:,4).*cos(ag-ones(l,1)*alp),length_branch.*k1.^grow(:,4).*sin(ag-ones(l,1)*alp),-ones(l,1)*alp,ones(l,1)]; parta2=[0.8.*length_branch.*k1.^grow(:,4).*cos(ag),0.8.*length_branch.*k1.^grow(:,4).*sin(ag),zeros(l,1),ones(l,1)];partb2=[0.8.*length_branch.*k1.^grow(:,4).*cos(ag),0.8.*length_branch.*k1.^grow(:,4).*sin(ag),zeros(l,1),ones(l,1)];parta=control.*parta+(1-control).*parta2;partb=control.*partb+(1-control).*partb2;parta=parta+grow;partb=partb+grow;congress=[parta;partb];grow=[grow;grow];judge=[grow,congress];judge=unique(judge,'rows');grow=judge(:,5:end);plotdata=[plotdata;judge];endfor i=1:number_branchtemp_w=width_branch*0.8^i;temp_branch=plotdata(plotdata(:,4)==i,:);plx=[temp_branch(:,1),temp_branch(:,5)];ply=[temp_branch(:,2),temp_branch(:,6)];plx=plx';ply=ply';plot(plx,ply,'color',[0 0 0]+i*[0.3020 0.3020 0.3020]./number_branch,'Linewidth',temp_w)endbloom_pos=plotdata(plotdata(:,8)==number_branch+1,[5,6]);scatter(bloom_pos(:,1),bloom_pos(:,2),10,'CData',[0.8549 0.6824 0.6824]) bloom_pos=plotdata(plotdata(:,8)==number_branch,[5,6]);scatter(bloom_pos(:,1),bloom_pos(:,2),8,'CData',[0.7451 0.5961 0.5961].*0.97) end2.苹果绘制function appleG[u,v]=meshgrid(linspace(-pi,pi,100));XMesh=(1+cos(v)).*cos(u)+0.085.*cos(5.*u)+(0.994.*v./pi).^100;YMesh=(1+cos(v)).*sin(u);ZMesh=(6.*sin(v)+2.*cos(v)-2.9.*log(1-v./pi))./4.2;ax=gca;ax.DataAspectRatio=[1 1 1];hold on;grid onsurf(XMesh,YMesh,ZMesh,'FaceColor',[140,192,30]./255,'EdgeColor','none')view(45,25)h=light;h.Color=[1 1 1];h.Style='infinite';end3.南⽠绘制% 构造⽹格[t,p]=meshgrid(linspace(0,2*pi,200),linspace(0,pi-.05,200));% ⼆元分段函数foutline=@(t,p)(p<=.14).*90.*cos(7.*p)+(p>.14).*(cos(20.*t)+70).*(p.*(pi-p+.1)).^.2; % 球⾯坐标转化为X,Y,Zh=cos(p).*foutline(t,p);R=sin(p).*foutline(t,p);R=R;X=cos(t).*R;Y=sin(t).*R;% 颜⾊矩阵构造CMap=ones([size(t),3]);tMap=ones(size(t));c1=[253,158,3]./255;c2=[76,103,86]./255;for i=1:3tMap(:,:)=c1(i);tMap(p<=.14)=c2(i);CMap(:,:,i)=tMap;endfigure()surf(X,Y,h,'CData',CMap,'EdgeColor',[0,0,0],'EdgeAlpha',.2)4.⼀堆三⾓形绘制function tritritri% @author:slandarerhold on,axis equalaxis([-10,10,-10,10])set(gca,'xtick',[],'ytick',[],'xcolor','w','ycolor','w')set(gca,'color','w')for N=1:5drawTri([0,0],[cos(pi/2+2*N*pi/5),sin(pi/2+2*N*pi/5)].*5,...[cos(pi/2+pi/5+2*N*pi/5),sin(pi/2+pi/5+2*N*pi/5)].*5,10)drawTri([0,0],[cos(pi/2+2*N*pi/5),sin(pi/2+2*N*pi/5)].*5,...[cos(pi/2-pi/5+2*N*pi/5),sin(pi/2-pi/5+2*N*pi/5)].*5,10)endL5=5/sin(pi/2-pi/10)*sin(pi/5);L5s=(5+L5)/sin(3*pi/10)*sin(3*pi/5);for N=1:5drawPent([cos(pi/2+2*N*pi/5),sin(pi/2+2*N*pi/5)].*5,...[cos(pi/2+2*N*pi/5),sin(pi/2+2*N*pi/5)].*(5+L5),...[cos(pi/2+pi/10+2*N*pi/5),sin(pi/2+pi/10+2*N*pi/5)].*L5s,...[cos(pi/2+pi/5+2*N*pi/5),sin(pi/2+pi/5+2*N*pi/5)].*(5+L5),...[cos(pi/2+pi/5+2*N*pi/5),sin(pi/2+pi/5+2*N*pi/5)].*5,15)drawPent([cos(pi/2+2*N*pi/5),sin(pi/2+2*N*pi/5)].*5,...[cos(pi/2+2*N*pi/5),sin(pi/2+2*N*pi/5)].*(5+L5),...[cos(pi/2-pi/10+2*N*pi/5),sin(pi/2-pi/10+2*N*pi/5)].*L5s,...[cos(pi/2-pi/5+2*N*pi/5),sin(pi/2-pi/5+2*N*pi/5)].*(5+L5),...[cos(pi/2-pi/5+2*N*pi/5),sin(pi/2-pi/5+2*N*pi/5)].*5,15)end% 功能函数================================================================= function drawTri(A,B,C,times)LW=1.2;COLOR=[0 0 0];k=1/6;k2=0.8;plot([A(1);B(1)],[A(2);B(2)],'color',COLOR,'Linewidth',LW)plot([B(1);C(1)],[B(2);C(2)],'color',COLOR,'Linewidth',LW)plot([C(1);A(1)],[C(2);A(2)],'color',COLOR,'Linewidth',LW)for i=1:timesb=k*(C-B)+B;plot([A(1);b(1)],[A(2);b(2)],'color',COLOR,'Linewidth',LW);B=b;pause(0.01)c=k*(A-C)+C;plot([B(1);c(1)],[B(2);c(2)],'color',COLOR,'Linewidth',LW);C=c;pause(0.01)a=k*(B-A)+A;plot([C(1);a(1)],[C(2);a(2)],'color',COLOR,'Linewidth',LW);A=a;pause(0.01)LW=LW*k2;endendfunction drawPent(A,B,C,D,E,times)LW=1.2;COLOR=[0 0 0];k=1/4;k2=0.9;plot([A(1);B(1)],[A(2);B(2)],'color',COLOR,'Linewidth',LW)plot([B(1);C(1)],[B(2);C(2)],'color',COLOR,'Linewidth',LW)plot([C(1);D(1)],[C(2);D(2)],'color',COLOR,'Linewidth',LW)plot([D(1);E(1)],[D(2);E(2)],'color',COLOR,'Linewidth',LW)plot([E(1);A(1)],[E(2);A(2)],'color',COLOR,'Linewidth',LW)for i=1:timesb=k*(C-B)+B;plot([A(1);b(1)],[A(2);b(2)],'color',COLOR,'Linewidth',LW);B=b;pause(0.01)c=k*(D-C)+C;plot([B(1);c(1)],[B(2);c(2)],'color',COLOR,'Linewidth',LW);C=c;pause(0.01)d=k*(E-D)+D;plot([C(1);d(1)],[C(2);d(2)],'color',COLOR,'Linewidth',LW);D=d;pause(0.01)e=k*(A-E)+E;plot([D(1);e(1)],[D(2);e(2)],'color',COLOR,'Linewidth',LW);E=e;pause(0.01)a=k*(B-A)+A;plot([E(1);a(1)],[E(2);a(2)],'color',COLOR,'Linewidth',LW);A=a;pause(0.01)LW=LW*k2;endendend5.⽉饼绘制function moonCake% @author:slandarerax=gca;hold(ax,'on');axis equalax.XLim=[-15,15];ax.YLim=[-15,15];CSet=[0.92 0.51 0.11;1 0.7 0.09;0.87 0.41 0.05];for i=[1:7,9,8]if i==1tt=linspace(0,-pi/16,100);elseif i==9tt=linspace(-pi+pi/16,-pi,100);elsett=linspace(-pi/16-(i-2)*pi/8,-pi/16-(i-1)*pi/8,100);endxSet=cos(tt).*(10+abs(cos(tt.*8)));xMin=find(xSet==min(xSet));tt(xMin)xMax=find(xSet==max(xSet));t1=min([xMin(1),xMax(1)]);t2=max([xMin(1),xMax(1)]);xSet=cos(tt(t1:t2)).*(10+abs(cos(tt(t1:t2).*8)));ySet=sin(tt(t1:t2)).*(10+abs(cos(tt(t1:t2).*8)))-3;fill([xSet(1),xSet,xSet(end)],[ySet(1)+3,ySet,ySet(end)+3],CSet(mod(i,2)+1,:),'EdgeColor','none') endt=linspace(0,2*pi,640);fill(cos(t).*(10+abs(cos(t.*8))),sin(t).*(10+abs(cos(t.*8))),CSet(1,:),'EdgeColor','none')plot(cos(t).*(9+abs(cos(t.*8))),sin(t).*(9+abs(cos(t.*8)))-0.3,'Color',CSet(3,:),'LineWidth',6)plot(cos(t).*8.7,sin(t).*8.7-0.3,'Color',CSet(3,:),'LineWidth',4)plot(cos(t).*(9+abs(cos(t.*8))),sin(t).*(9+abs(cos(t.*8))),'Color',CSet(2,:),'LineWidth',6)plot(cos(t).*8.7,sin(t).*8.7,'Color',CSet(2,:),'LineWidth',4)plot([0 0],[-7 7]-0.3,'Color',CSet(3,:),'LineWidth',4)plot([-7 7],[0 0]-0.3,'Color',CSet(3,:),'LineWidth',4)plot([0 0],[-7 7],'Color',CSet(2,:),'LineWidth',4)plot([-7 7],[0 0],'Color',CSet(2,:),'LineWidth',4)t4=linspace(0,pi/2,100);xSet4=[cos(t4).*6+1,1,6,cos(t4(1:end-12)).*5+1,2,4.8,cos(t4(17:end-40)).*3.9+1];ySet4=[sin(t4).*6+1,1,1,sin(t4(1:end-12)).*5+1,2,2,sin(t4(17:end-40)).*3.9+1];plot(xSet4,ySet4-0.3,'Color',CSet(3,:),'LineWidth',4)plot(-xSet4,ySet4-0.3,'Color',CSet(3,:),'LineWidth',4)plot(xSet4,-ySet4-0.3,'Color',CSet(3,:),'LineWidth',4)plot(-xSet4,-ySet4-0.3,'Color',CSet(3,:),'LineWidth',4)plot(xSet4,ySet4,'Color',CSet(2,:),'LineWidth',4)plot(-xSet4,ySet4,'Color',CSet(2,:),'LineWidth',4)plot(xSet4,-ySet4,'Color',CSet(2,:),'LineWidth',4)plot(-xSet4,-ySet4,'Color',CSet(2,:),'LineWidth',4)end6.⼤钻⽯绘制function diamondx2% @author:slandarerhold onax=gca;axis equalgrid onax.GridLineStyle='--';ax.LineWidth=1.2;ax.XLim=[-2,2];ax.YLim=[-2,2];ax.ZLim=[-1.8,1.5];ax.XColor=[1,1,1].*0.4;ax.YColor=[1,1,1].*0.4;ax.ZColor=[1,1,1].*0.4;ax.DataAspectRatio=[1,1,1];ax.DataAspectRatioMode='manual';ax.CameraPosition=[-67.6287 -204.5276 82.7879];light('Position',[5,6,5])light('Position',[5,6,5])faceColor=[0 71 177];facePara={faceColor./255,'FaceAlpha',0.4,'EdgeColor',faceColor./255.*0.8,'EdgeAlpha',0.2,...'FaceLighting','gouraud','EdgeLighting','gouraud',...'SpecularStrength',0.9,'DiffuseStrength',0.8,'AmbientStrength',0.7}; % ⾯参数drawDiamond(facePara)faceColor=[237 64 153];facePara={faceColor./255,'FaceAlpha',0.4,'EdgeColor',faceColor./255.*0.8,'EdgeAlpha',0.2,... 'FaceLighting','gouraud','EdgeLighting','gouraud',...'SpecularStrength',0.9,'DiffuseStrength',0.8,'AmbientStrength',0.7}; % ⾯参数drawDiamond(facePara)function drawDiamond(facePara)thetaSet8=linspace(0,2*pi-pi/4,8)';thetaSet16=linspace(0,2*pi-pi/8,16)';% 绘制最上⽅⼋边形pntSet1=[cos(thetaSet8),sin(thetaSet8),ones(8,1)];fill3(pntSet1(:,1),pntSet1(:,2),pntSet1(:,3),facePara{:})% ⼀圈8个三⾓形pntSet2=[cos(thetaSet8+pi/8).*1.3,sin(thetaSet8+pi/8).*1.3,ones(8,1).*0.7];for i=1:8p1=i;p2=mod(i,8)+1;p3=i;fill3([pntSet1([p1,p2],1);pntSet2(p3,1)],...[pntSet1([p1,p2],2);pntSet2(p3,2)],...[pntSet1([p1,p2],3);pntSet2(p3,3)],facePara{:})end% ⼀圈8个四边形pntSet3=[cos(thetaSet8).*1.7,sin(thetaSet8).*1.7,ones(8,1).*0.3];for i=1:8p1=mod(i,8)+1;p2=i;p3=p1;p4=p1;fill3([pntSet1(p1,1);pntSet2(p2,1);pntSet3(p3,1);pntSet2(p4,1)],...[pntSet1(p1,2);pntSet2(p2,2);pntSet3(p3,2);pntSet2(p4,2)],...[pntSet1(p1,3);pntSet2(p2,3);pntSet3(p3,3);pntSet2(p4,3)],facePara{:})end% ⼀圈16个三⾓形pntSet4=[cos(thetaSet16).*1.7,sin(thetaSet16).*1.7,ones(16,1).*0.3];for i=1:8p1=i;p2=2*i-1;p3=2*i;fill3([pntSet2(p1,1);pntSet4([p2,p3],1)],...[pntSet2(p1,2);pntSet4([p2,p3],2)],...[pntSet2(p1,3);pntSet4([p2,p3],3)],facePara{:})p2=2*i;p3=mod(2*i,16)+1;fill3([pntSet2(p1,1);pntSet4([p2,p3],1)],...[pntSet2(p1,2);pntSet4([p2,p3],2)],...[pntSet2(p1,3);pntSet4([p2,p3],3)],facePara{:})end% ⼀圈16个四边⾓形pntSet5=[cos(thetaSet16).*1.7,sin(thetaSet16).*1.7,ones(16,1).*0.15];for i=1:16p1=i;p2=mod(i,16)+1;p3=mod(i,16)+1;p4=i;fill3([pntSet4(p1,1);pntSet4(p2,1);pntSet5(p3,1);pntSet5(p4,1)],...[pntSet4(p1,2);pntSet4(p2,2);pntSet5(p3,2);pntSet5(p4,2)],...[pntSet4(p1,3);pntSet4(p2,3);pntSet5(p3,3);pntSet5(p4,3)],facePara{:})endpntSet6=[cos(thetaSet8+pi/8).*0.4,sin(thetaSet8+pi/8).*0.4,ones(8,1).*(-1.15)];% ⼀圈16个长三⾓形for i=1:8p1=i;p2=2*i-1;p3=2*i;fill3([pntSet6(p1,1);pntSet5([p2,p3],1)],...[pntSet6(p1,2);pntSet5([p2,p3],2)],...[pntSet6(p1,3);pntSet5([p2,p3],3)],facePara{:})p2=2*i;p3=mod(2*i,16)+1;fill3([pntSet6(p1,1);pntSet5([p2,p3],1)],...[pntSet6(p1,2);pntSet5([p2,p3],2)],...[pntSet6(p1,3);pntSet5([p2,p3],3)],facePara{:})end% ⼀圈8个长四边形pntSet7=[cos(thetaSet8).*1.7,sin(thetaSet8).*1.7,ones(8,1).*0.15];for i=1:8p1=mod(i,8)+1;p2=i;p4=mod(i,8)+1;fill3([pntSet7(p1,1);pntSet6(p2,1);0;pntSet6(p4,1)],...[pntSet7(p1,2);pntSet6(p2,2);0;pntSet6(p4,2)],...[pntSet7(p1,3);pntSet6(p2,3);-1.5;pntSet6(p4,3)],facePara{:})endendend7.有趣曲线1% 鱼形曲线hold onaxis equalgrid onX=0:1:1023;Y=0:1:1023;[gridX,gridY]=meshgrid(X,Y);FishPatternFcn=@(x,y)mod(abs(x.*sin(sqrt(x))+y.*sin(sqrt(y))).*pi./1024,1);contour(gridX,gridY,FishPatternFcn(gridX,gridY),[0.7,0.7])8.有趣曲线2% ⽣成⼀个40x40 DFT矩阵% 并在复平⾯将点连接起来plot(fft(eye(40)));axis equal9.有趣曲线3——蝴蝶曲线x=linspace(0,2*pi,1000);hdl=polarplot(x,sin(x)+sin(3.5.*x).^3,'LineWidth',1.5);ax=gca;ax.RTick=[];10.有趣曲线4——⼼形曲线function loveFunc% @author:slandarerLoveFunchdl=@(x,a)(x.^2).^(1/3)+0.9.*((3.3-x.^2).^(1/2)).*sin(a.*pi.*x);hold ongrid onaxis([-3 3,-2 4])x=-1.8:0.01:1.8;text(0,3.3,'$f(x)=x^{\frac{2}{3}}+0.9(3.3-x^2)^{\frac{1}{2}}\sin(\alpha\pi x)$',...'FontSize',13,'HorizontalAlignment','center','Interpreter','latex');txt2=text(-0.35,2.9,'','FontSize',13,'HorizontalAlignment','left','Interpreter','latex','tag','alphadata'); hdl=plot(x,LoveFunchdl(x,1),'color',[0.2118 0.4667 0.9961],'LineWidth',1.5);for a=1:0.01:20AlphaString=['$\alpha=',num2str(a),'$'];Color=([1.0000 0.4902 0.6627]-[0.2118 0.4667 0.9961]).*(a/20)+[0.2118 0.4667 0.9961];set(txt2,'string',AlphaString)hdl.YData=LoveFunchdl(x,a);hdl.Color=Color;drawnowendend11.有趣曲线5[x,y]=meshgrid(-10:.01:10,-10:.01:10);cm=cos(cos(min(sin(x)+y,x+sin(y))))-cos(sin(max(sin(y)+x,y+sin(x))))>0;imshow(cm)12.会害羞的含羞草function mimose% @author:slandarer% ========================================================================= % <构造含羞草叶⼦轮廓曲线>% ⼤体是画⼀个⼀个压扁倾斜的半圆弧拼在⼀起% 然后乘以⼀个类似弧形曲线让每⽚⼩叶⼦不⼀样长% 计算⼀个压扁倾斜的半圆弧曲线t1=0:pi/100:pi+2*pi/100;y=5*abs(sin(t1).^(1/2));y(t1>pi)=-y(t1>pi);t1=t1.*cos(pi/9)-y.*sin(pi/9);y=t1.*sin(pi/9)+y.*cos(pi/9);% 将多个压扁倾斜的半圆弧曲线拼接T=[];Y=[];for i=1:20T=[T,(i-1)*(pi+2*pi/100)+t1];Y=[Y,y];end% 乘以⼀个类似弧形曲线让每⽚⼩叶⼦不⼀样长L=length(T);t2=linspace(pi/8,pi-2*pi/5,L);Y=Y.*sin(t2);% 旋转并调整⽐例X=T.*cos(pi)-Y.*sin(pi);Y=T.*sin(pi)+Y.*cos(pi);X=X-min(X);X=X.*0.8;Y=Y.*1.5;% 创建figure及axes=========================================================fig=figure('units','pixels',...'position',[500 100 500 500],...'Numbertitle','off',...'Color',[1 1 1],...'resize','off',...'visible','on',...'name','mimose',...'menubar','none');ax=axes('Units','pixels',...'parent',fig,...'Color',[0.2667 0.4039 0.2745],...'Position',[0 0 500 500],...'XLim',[0 100],...'YLim',[0 100],...'XColor',[1 1 1],...'YColor',[1 1 1]);hold(ax,'on')%plot(ax,X+20,Y+50)% 调⽤后⾯定义的函数绘图===================================================% 绘制树枝drawBranch([20,0],[50,10])drawBranch([20,0],[35,12])drawBranch([0,0],[15,45])drawBranch([4,0],[0,20])% 绘制树叶leaf.l1=drawLeaf([50,10]+1.*[cos(pi/1.7),sin(pi/1.7)],pi/1.7,X,Y,0.8,1);leaf.l2=drawLeaf([50,10]+1.*[cos(-pi/8),sin(-pi/8)],-pi/8,X,Y,0.8,1);leaf.l3=drawLeaf([50,10]+1.*[cos(pi/12),sin(pi/12)],pi/12,X,Y,1,1);leaf.l4=drawLeaf([50,10]+1.*[cos(pi/3),sin(pi/3)],pi/3,X,Y,1,1);leaf.l5=drawLeaf([0,45]+1.*[cos(pi/3),sin(pi/3)],pi/3,X,Y,1,1);leaf.l6=drawLeaf([0,30]+1.*[cos(pi/5),sin(pi/5)],pi/5,X,Y,1,1);% 绘制花朵flower.fl1=drawFlower([35,12]);flower.fl2=drawFlower([20,5]);% ⿏标移动回调============================================================= fps=25;mitimer=timer('ExecutionMode', 'fixedRate', 'Period',1/fps, 'TimerFcn', @miMove);start(mitimer)% 花朵树叶回缩函数function miMove(~,~)for ii=1:6if leaf.(['l',num2str(ii)]).ratio<1 %如果收缩率⼩于1则每隔⼀段时间+0.05,以下同理leaf.(['l',num2str(ii)]).ratio=leaf.(['l',num2str(ii)]).ratio+0.05;endresetH(leaf.(['l',num2str(ii)]),X,Y); % 依据收缩率重新绘制叶⼦endfor ii=1:2if flower.(['fl',num2str(ii)]).ratio<1flower.(['fl',num2str(ii)]).ratio=flower.(['fl',num2str(ii)]).ratio+0.05;endresetL(flower.(['fl',num2str(ii)]));endendset(gcf,'WindowButtonMotionFcn',@whilemovefcn)% ⿏标移动调⽤的函数function whilemovefcn(~,~)% 检测⿏标位置xy=get(gca,'CurrentPoint');xw=xy(1,1);yw=xy(1,2);pos=[xw,yw];% ⿏标附近的叶⼦收缩率-0.2for ii=1:6LF=leaf.(['l',num2str(ii)]);if closeLeaf(LF.pos1,LF.pos2,pos)% 如果距离叶⼦距离⾜够近if leaf.(['l',num2str(ii)]).ratio>0.2leaf.(['l',num2str(ii)]).ratio=leaf.(['l',num2str(ii)]).ratio-0.2;endendend% ⿏标附近的花收缩率-0.2for ii=1:2FL=flower.(['fl',num2str(ii)]);if closeFlower(FL.pos,pos)% 如果距离花距离⾜够近if flower.(['fl',num2str(ii)]).ratio>0.6flower.(['fl',num2str(ii)]).ratio=flower.(['fl',num2str(ii)]).ratio-0.2;endendendend% 距离检测函数,叶⼦形状近乎于椭圆,当距离俩焦点和⼩于⼀定值认为⾜够近% 若⾜够近返回TRUE 否则返回FALSEfunction bool=closeLeaf(pos1,pos2,pos3)pL=pos1+(pos2-pos1).*0.1;pR=pos1+(pos2-pos1).*0.9;lTotal=norm(pL-pos3)+norm(pR-pos3);bool=lTotal<=norm(pos1-pos2);end% 距离检测函数,花形状近乎于圆,当距离圆⼼⼩于⼀定值认为⾜够近% 若⾜够近返回TRUE 否则返回FALSEfunction bool=closeFlower(pos1,pos2)bool=norm(pos1-pos2)<=5;end% 依据收缩率重新计算并绘制叶⼦function resetH(lf,X,Y)X=X.*lf.Len;Y=Y.*lf.h.*lf.ratio;xl1=X.*cos(lf.alf)-Y.*sin(lf.alf)+lf.pos1(1);yl1=X.*sin(lf.alf)+Y.*cos(lf.alf)+lf.pos1(2);xl2=X.*cos(lf.alf)+Y.*sin(lf.alf)+lf.pos1(1);yl2=X.*sin(lf.alf)-Y.*cos(lf.alf)+lf.pos1(2);XL=[xl1,xl2(end:-1:1)];YL=[yl1,yl2(end:-1:1)];set(lf.LHdl,'XData',XL,'YData',YL);end% 依据收缩率重新计算并绘制花朵function resetL(fl)xf=fl.x.*fl.ratio+fl.pos(1);yf=fl.y.*fl.ratio+fl.pos(2);xxf(2:2:2*length(xf))=xf;yyf(2:2:2*length(yf))=yf;xxf(1:2:2*length(xf))=fl.pos(1);yyf(1:2:2*length(yf))=fl.pos(2);%set(fl.plHdl,'XData',xxf,'YData',yyf);set(fl.scHdl,'XData',xf,'YData',yf);set(fl.plHdl,'XData',xxf,'YData',yyf);end% 绘制花朵(圆⼼)function fl=drawFlower(pos)% 随机⽣成⼀些点作为花蕊并绘图theta=rand([1,120]).*2.*pi;r=rand([1,120]).*3+5;xf=r.*cos(theta)+pos(1);yf=r.*sin(theta)+pos(2);xxf(2:2:2*length(xf))=xf;yyf(2:2:2*length(yf))=yf;xxf(1:2:2*length(xf))=pos(1);yyf(1:2:2*length(yf))=pos(2);plHdl=plot(ax,xxf,yyf,'Color',[0.7608 0.4863 0.7216],'LineWidth',1);scHdl=scatter(ax,xf,yf,3,'filled','CData',[0.9569 0.9804 0.9333]);% 存储基础信息到实类fl.x=xf-pos(1);fl.y=yf-pos(2);fl.pos=pos;fl.scHdl=scHdl;fl.plHdl=plHdl;fl.ratio=1;% 计算可伸缩部分花蕊坐标,绘图,存储信息theta2=rand([1,80]).*2.*pi;r2=rand([1,80]).*1+3;xf2=r2.*cos(theta2)+pos(1);yf2=r2.*sin(theta2)+pos(2);xxf2=[ones([1,length(xf2)]).*pos(1);xf2];yyf2=[ones([1,length(yf2)]).*pos(2);yf2];plot(ax,xxf2,yyf2,'Color',[0.7804 0.6078 0.5020],'LineWidth',1)end% 绘制树枝(端点1,端点2)function drawBranch(pos1,pos2)% 依据起始及末尾两点计算⼀个⼀头宽⼀头窄四边形坐标数据并绘图dir=(pos2-pos1)./norm(pos2-pos1);len=norm(pos2-pos1);xb=[0 1 1 0].*len;yb=[len.*0.02,len.*0.012,-len.*0.012,-len.*0.02];xxb=xb.*dir(1)-yb.*dir(2)+pos1(1);yyb=xb.*dir(2)+yb.*dir(1)+pos1(2);fill(ax,xxb,yyb,[0.8157 0.6431 0.6078],'EdgeColor',[0.6157 0.5529 0.4510],'LineWidth',1.5)end% 绘制叶⼦(起始位置,叶⼦⽅向,轮廓x数据,轮廓y数据,轮廓x⽅向放缩⽐例,轮廓y⽅向放缩⽐例) function lf=drawLeaf(pos,alf,X,Y,l,h)X=X.*l; % 放缩叶⼦长度Y=Y.*h; % 放缩叶⼦宽窄% 叶⼦数据旋转并把两个叶⼦数据并在⼀起xl1=X.*cos(alf)-Y.*sin(alf)+pos(1);yl1=X.*sin(alf)+Y.*cos(alf)+pos(2);xl2=X.*cos(alf)+Y.*sin(alf)+pos(1);yl2=X.*sin(alf)-Y.*cos(alf)+pos(2);XL=[xl1,xl2(end:-1:1)];YL=[yl1,yl2(end:-1:1)];% 绘图lf.LHdl=fill(XL,YL,[0.4980 0.7216 0.3176]);% 基础数据存到实类⾥lf.pos1=pos;lf.pos2=[pos(1)+l*51*cos(alf),pos(2)+l*51*sin(alf)];lf.alf=alf;lf.Len=l;lf.h=h;lf.ratio=1;% 绘制⼀下叶⼦中线plot(ax,[pos(1),pos(1)+l*51*cos(alf)],[pos(2),pos(2)+l*51*sin(alf)],'Color',[0.6627 0.6980 0.4431],'LineWidth',3); endend13.随机雪景function snowCover% @author:slandarerMainFig=figure('units','pixels','position',[300 80 500 500],...'Numbertitle','off','menubar','none','resize','off',...'name','snow covered landscape');axes('parent',MainFig,'position',[0 0 1 1],...'XLim', [0 500],...'YLim', [0 500],...'NextPlot','add',...'layer','bottom',...'Visible','on',...'Color',[0 59 129]./255,...'XTick',[], ...'YTick',[]);hold on%画雪地==================================================================== layerBEPos=[150;140;100;90];layerColor=[146 210 245;101 183 231;0 101 181;0 59 129];excursion=15;interval=10;xSep=25;layerPos=zeros(size(layerBEPos,1),length(xSep:xSep:500));layerPos(:,1)=layerBEPos(:,1);for i=1:size(layerBEPos,1)for j=xSep:xSep:500listPos=round(j./xSep+1);tempRandi=randi([-excursion,excursion]);yPos=tempRandi+layerPos(i,listPos-1);if i>1&&yPos>=layerPos(i-1,listPos)-5yPos=layerPos(i-1,listPos)-interval;endyPos(yPos<0)=0;layerPos(i,listPos)=yPos;endendfor i=1:size(layerBEPos,1)XData=0:xSep:500;YData=layerPos(i,:);Yq=interp1(XData,YData,0:1:500,'spline');Xq=[0,0:1:500,500];Yq=[0 Yq 0];fill(Xq,Yq,layerColor(i,:)./255,'EdgeColor','none')end%画星星==================================================================== XRandiS=randi([50 450],[1,15]);YRandiS=randi([300 460],[1,15]);scatter(XRandiS,YRandiS,10,'o','filled','CData',[252 241 0]./255)XRandiB=randi([50 450],[1,15]);YRandiB=randi([300 460],[1,15]);scatter(XRandiB,YRandiB,18,'o','filled','CData',[252 241 0]./255)%画雪花====================================================================snowXpos=[randi([0,30],[1,6]),randi([470,500],[1,6]),randi([0,500],[1,12])];snowYpos=[randi([220,470],[1,6]),randi([220,470],[1,6]),randi([470,500],[1,12])];snowSize=randi([25,50],[1,24]);snowAngle=2*pi*rand([1,24]);snowWidth=1+0.6*rand([1,24]);for i=1:24drawSnow(snowXpos(i),snowYpos(i),snowSize(i),snowAngle(i),snowWidth(i))end%画⽉亮====================================================================moonPos=[320 300]+[randi([0 50]),randi([0 80])];moonSize=randi([28 32]);blankPos=randi([floor(0.3*moonSize) ceil(0.5*moonSize)],[1,2]);drawMoon(moonPos(1),moonPos(2),moonSize,moonPos(1)-blankPos(1),moonPos(2)+blankPos(2))%相关函数==================================================================function drawMoon(x,y,R,bx,by)t=0:pi/50:2*pi;X=x+cos(t).*R;Y=y+sin(t).*R;BX=bx+cos(t).*R;BY=by+sin(t).*R;fill(X,Y,[255 251 219]./255,'EdgeColor','none')fill(BX,BY,[0 59 129]./255,'EdgeColor','none')endfunction drawSnow(x,y,len,angle,width)for theta=0:pi/3:2*pi-pi/3xTail=cos(theta+angle)*len+x;yTail=sin(theta+angle)*len+y;plot([x,xTail],[y,yTail],'color',[0 134 207]./255,'lineWidth',width)for branchRatio=[0.4 0.55 0.7]BX=x+cos(theta+angle)*len*branchRatio;BY=y+sin(theta+angle)*len*branchRatio;LX=BX+cos(theta+angle+pi/3).*len.*0.4;LY=BY+sin(theta+angle+pi/3).*len.*0.4;RX=BX+cos(theta+angle-pi/3).*len.*0.4;RY=BY+sin(theta+angle-pi/3).*len.*0.4;plot([BX,LX],[BY,LY],'color',[0 134 207]./255,'lineWidth',width.*0.8)plot([BX,RX],[BY,RY],'color',[0 134 207]./255,'lineWidth',width.*0.8)endendendend以上就是利⽤Matlab绘制有趣图像的⽰例代码的详细内容,更多关于Matlab绘制图像的资料请关注其它相关⽂章!。

MATLAB第三节绘图

MATLAB第三节绘图
Note: (1) fun为函数名,以字符串形式出现。lims为变量取 值范围,tol为相对允许误差,其系统默认值为2e-3。 (2) fun必须是M文件的函数名或是独立变量为x的字 符串. (3) fplot函数不能画参数方程和隐函数图形,但在一 个图上可以画多个图形.
【例】 在[-2,2]范围内绘制函数tanh的图形.
最基本的三维图形函数为plot3,它是将二维函数plot 的有关功能扩展到三维空间,用来绘制三维图形。
Plot3 ( x,y,z,c)
Note:函数功能:以向量x,y,z为坐标,绘制三维曲 线,其中x,y,z表示三维坐标向量,表示线形或颜色。
4.4 三维图形
【例】 绘制三维螺旋曲线:
t=0:pi/50:10*pi; y1=sin(t),y2=cos(t); plot3(y1,y2,t); title('helix'); text(0,0,0,‘origin’); %在(0,0,0)点添加文本框 xlabel('sin(t)'),ylabel('cos(t)'),zlabel('t');
4.3 其它图形函数
3、填充图形
fill(x,y,’c’) — 绘制并填充二维多边图形
Note:x和y为二维多边形顶点坐标向量。字符 ’c’ 规 定填充颜色,其取值前已叙述。
如:绘制一正方形并以黄色填充: x=[0 1 2 4 0]; %正方形顶点坐标向量 y=[0 0 1 1 0]; fill(x,y,‘b’) %绘制并以蓝色填充正方形图
Note:表示在区间tmin<t<tmax绘制参数方程 x=x(t), y=y(t)的函数图.
【例】 在[0, 2 ]π上画 y sin3 t 的x星形co图s3 t.

matlab三维绘图ppt课件

matlab三维绘图ppt课件
MATLAB 绘图
xx
1
三维曲线绘图 三维曲面绘图
2
xx
二、MATLAB三维曲线绘图
plot3——三维曲线绘制指令 plot3的调用格式:
plot3(X,Y,Z) plot3(X,Y,Z,’String’) plot3(X1,Y1,Z1,’ String1’,X2,Y2,Z2,’ String2’,…) plot3与plot的 用法相同
➢ shading faceted
以平面为单位进行着色,在flat用色基础上,在贴片的 四周勾出黑色网线。
23
subplot(131),surf(peaks(40)); subplot(132),surf(peaks(40));shading flat subplot(133),surf(peaks(40));shading interp
subplot(122),h = surf(ones(10,10));rotate(h,[0 0 1],i,[1 0 0]) end
➢ rotate3d
动态旋转命令,可以让用户使用鼠标来旋转视角
19
背景颜色控制命令colordef
➢ colordef white 将图形的背景颜色设置为白色 ➢ colordef black ➢ colordef none 将图形背景和图形窗口的颜色设置
28
对象和句柄
MATLAB把构成图形的各个基本要素称为图形对象, 包括图形窗口、坐标轴、线条、曲面、文本和其它对 象。
每一个对象都有一个数字来标识,叫做句柄。当每次 创建一个对象时,MATLAB就为它建立一个唯一的句 柄。
29
作业ห้องสมุดไป่ตู้P79
5,7,8,9
30
22

Matlab绘图技巧与实例

Matlab绘图技巧与实例

Matlab绘图技巧与实例绘图在科学和工程领域中起着重要的作用,而Matlab作为一种功能强大的数学软件,具有丰富的绘图功能。

本文将介绍一些Matlab的绘图技巧,并通过一些实例来展示其用法和优势。

一、基本的绘图命令Matlab提供了一系列用于绘图的基本函数,最常用的是plot和scatter。

plot函数用于绘制曲线图,而scatter函数则用于绘制散点图。

这两个函数都可以接受多组数据,并且具有丰富的参数设置,可以对图形进行自定义。

例如,我们可以设置线条的颜色、线型和线宽,还可以添加标签和图例等。

二、特殊图形的绘制除了常见的曲线图和散点图外,Matlab还可以绘制一些特殊的图形,如柱状图、饼图和雷达图等。

这些图形可以用于展示不同类型的数据,从而更直观地呈现结果。

例如,柱状图可以用于比较不同组的数据,饼图则可以用于显示百分比等。

在绘制这些特殊图形时,Matlab提供了相应的函数,如bar、pie和polar等,使用这些函数可以轻松实现各种图形的绘制。

三、绘制3D图形Matlab还支持绘制3D图形,通过将数据在三维坐标系中表示,可以更全面地展示数据的分布和关系。

Matlab提供了许多用于绘制3D图形的函数,如plot3、scatter3和surf等。

使用这些函数可以绘制出各种复杂的3D图形,如曲面图、散点云和体积渲染等。

在绘制3D图形时,我们可以设置视角、光照和颜色等参数,从而使图形更加生动逼真。

四、图形的美化与字体设置除了绘图功能外,Matlab还提供了一些功能用于美化图形和设置字体。

通过设置标签和标题的字体、大小和颜色等,可以让图形更加清晰和美观。

此外,Matlab 还支持设置坐标轴的刻度、标签和范围,以及图形的背景颜色和边框样式等。

这些设置可以提高图形的可读性和视觉效果,从而更好地传达数据和结果。

五、图形的输出与保存Matlab不仅可以在软件中生成图形,还可以将图形输出为不同的格式,如图片文件和矢量图等。

四讲Matlab绘图ppt课件

四讲Matlab绘图ppt课件
第四讲 Matlab绘图
2019/12/15
2019/12/15
1
第四讲 Matlab绘图
• 3.1 二维图形 • 3.2 数据分析图 • 3.3 三维图形
2019/12/15
2
第四讲 Matlab的图形
M AT L A B拥有大量简单、灵活、易用的 二维和三维图形命令,许多图形命令都 在M AT L A B所带的演示程序d e m o中 给出。
2019/12/15
12
一、二维图形
• 1.m a n u a l 固定坐标轴刻度。如果当前图形 窗口为h o l d打开状态,则后面的图形将采用 同样的刻度
• 2. auto 把坐标轴刻度重新设置为缺省状态值。 • 3.e q u a l 设置x轴和y轴为同样的刻度增量。 • 4.t i g h t 以数据的大小为坐标轴的范围。
2019/12/15
17
一、二维图形
• t e x t ( x , y , t x t , ’ s c ’ ) 在图形窗口的 (x, y)处输出字符串t x t,给定左下角的 坐标为(0.0, 0.0),右上角的坐标则为(1.0, 1.0)。
• g t e x t ( t x t ) 通过使用鼠标或方向键, 移动图形窗口中的十字光标,让用户将 字符串t x t放置在图形窗口中。当十字光 标走到所期望的位置时,用户按下任意
2019/12/15
11
一、二维图形
a x i s(…)用行向量中给出的值,设置坐 标轴的最大和最小值。对于二维图形, 该向量中含有元素: [xmin, xmax, ymin, ymax]。对于三维图形,是[xmin, xmax, ymin, ymax,zmin, zmax]。

如何用matlab画出高质量的图

如何用matlab画出高质量的图

在我们写论文的时候难免会需要画一些图,但我们总是感觉matlab 画的图差强人意,在介绍如何用matlab 画出高质量的图之前,我们先看4张正弦图:XY图一图三图二图四图一是通过Copy Figure 得到的,图二是通过Export Setup 得到的,图三是通过截图得到的,图四也是通过Export Setup 得到的。

如果你觉得图一、图二、图三比图四要好,那么你可以不用往下面看了,免得浪费时间。

如果你觉得图四是最好的,那请继续看下去。

首先说一下利用Export Setup 得到图片,一般保存为两种格式比较好,一种是*.tif 格式,一种是*.eps 格式。

这两种格式的区别在于,*.tif 格式用于word 比较好,如果你是用word 写论文,请将图片保存为此种格式;*.eps 格式用于Latex 编写的文档比较好(我没有试过,网上别人是这么说的),如果你是用Latex 写论文,请将图片保存为此种格式。

不过我们可以看图二,图二就是*.eps 格式的图片,在word 中显示确实不好。

下面来讲如何通过Export Setup 获得*.tif 和*.eps 格式的图片,获取此两种格式的图片,前期处理都是一样的,只是到最后保存时有所区别。

第一步:我们在M 文件中写好程序,按下F5就会出现Figure 的对话框。

第二步:点击,然后点击Export Setup,会出现Export Setup对话框在里面有Size、Rendering、Fonts、Lines四个选项,Size可以调节输出图片大小,Rendering是最关键的,点击Rendering我们会看到(其他两个自己上网查一下如何使用)将的值直接设置为1000,虽然它本身的选择值只有300、600和auto三个选项,但可以自己直接设置,用键盘直接输出1000。

这个设置值是设置像素,像素高画出的图会感觉很光滑。

设置完后可以点击下图中的Save按钮,这样可以在下次输出图片时,保持原有设置的值。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

目录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份,然后一等分点为圆心,以等分点到该圆垂直直径%%的距离为半径画圆。

相关文档
最新文档