第三讲--绘图.ppt

其它修饰： a. 水线修饰（waterfall) waterfall(peaks(30)) b. 等高线修饰 二维 contour(Z,n)--------绘制n条等高线 C= contourc(Z,n)------计算n条等高线的坐标 Clable(c)------给等高线加标注 在二维平面上绘制peaks函数的10条等高线 contour(peaks,10); C=contourc(peaks,10);clabel(C) %三维 contour3(peaks,20) （六）动画效果 动画生成的步骤 创建帧矩阵 ————— moviein 对动画中的每一帧生成图形，并把它们放到帧矩阵中 ———— getframe 从帧矩阵中回放动画 函数Moviein()产生一个帧矩阵来存放动画中的帧； 函数getframe对当前的图象进行快照； 函数movie()按顺序回放各帧。 解 [x,y,z]=peaks(30); surf(x,y,z) axis([-3 3 -3 3 -10 10]) m=moviein(150); for i=1:150 view(-37.5+24*(i-1),30) m(:,i)=getframe; end movie(m) 彗星状轨迹图 n=0.5; t=n*pi*(0:0.0005:48); x=sin(t);y=cos(t); plot(x,y, r); hold on;comet(x,y,6);hold off 卫星返回地球的运动轨线示意。 R0=1; %以地球半径为一个单位 a=24*R0;b=15*R0;T0=8*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) %画全程轨线 [X,Y,Z]=sphere(30);X=R0*X;Y=R0*Y;Z=R0*Z; %获得单位球坐标 grid on,hold on,surf(X,Y,Z),shading interp %画地球 x1=-18*R0;x2=6*R0;y1=-12*R0;y2=12*R0;z1=-6*R0;z2=6*R0; axis([x1 x2 y1 y2 z1 z2]) %确定坐标范围 view([117 37]),comet3(x,y,z,0.02),hold off %设视角、画运动轨线 一,设在2×4m房间的顶部，距离地面高度为2m处安装两盏照明灯，两灯相距2m，距离墙面为1m. 现假设灯的光强为1个单位，空间任一点处的光强与它们之间的距离成反比。试分析地面上光强的分布情况. 距离地面1m的层面上光强分布又如何？距离地面1.5m呢？(surf,meshgrid) 分析：以地面为xoy平面，以中心点为坐标原点建立空间坐标系则两灯的位置为A(0,-1,2)、B(0,1,2)，则距离地面高为h的空间点p(x,y,h)处的光强根据光的迭加原理有： 程序： x=-1:.1:1; y=-2:.1:2; [X,Y]=meshgrid(x,y); h=1; q1=1./sqrt(X.^2+(Y-1).^2+(h-2)^2); q2=1./sqrt(X.^2+(Y+1).^2+(h-2)^2); q=q1+q2; surf(X,Y,q); subplot(4,2,6); %绘制魔方图 x=magic(6);area(x) subplot(4,2,7); %绘制饼图 x1=[1 2 3 4 5 6 7]; y1=[0 0 0 1 0 0 0]; pie(x1,y1) subplot(4,2,8); %绘制散点图 a=rand(200,1); b=rand(200,1); c=rand(200,1); scatter(a,b,100,c,p) 例：三维螺线 Matlab 空间曲线绘图举例 x=t, y=sin(t), z=cos(t), 0 t 20 1) 给出空间离散点的坐标 (x,y,z) 2) 将这些点按顺序连接即可 先画点，后连线 空间曲线作图举例 t=[0:0.5:20]; x=t; y=sin(t); z=cos(t); plot3(x,y,z,’.-’) t=[0:0.5:20]; x=t; y=sin(t); z