Matlab基础常用绘图例子
Matlab基础常用绘图例子
学习练习时的一些作业结果分享。
aa.m文件
v=[’-1’;’+1’];
t=[0:0.05:2pi]; %初值为0，末值为2π，步长为0.05的一维矩阵
plot(t,sin(t)); %绘制sin(t)，t在[0,2π]的函数图像
xlabel(‘time(0–2\pi)’); %在横轴标记time(0–2\pi)
ylabel(‘value’); %在纵轴标记value
text(pi/2,0.9,[’\downarrowwsin(\pi1/2)=’,v(2,:)]);
text(pi3/2,-0.9,[’\downarrowwsin(\pi3/2)=’,v(1,:)]);
text(0,-0.6,[‘Date:’,date]); %在（0，-0.6）处添加当前日期
运行结果
clf %清除当前图形
t=0:0.1:4pi;
y=exp(-0.1t).sin(t);
figure(1) %创建1个图形窗口
subplot(2,2,1) , %将当前图形窗口分成22个绘图区，将下图绘于第一个区
stem(t,y) %绘制以t为横坐标，y为纵坐标的杆图
title(‘stem(t,y)’), %杆图标题标记stem(t,y)
pause %暂停
subplot(2,2,2), %在第一行第二列位置绘图
stairs(t,y) %绘制阶梯图
title(‘stair(t,y)’),
pause
subplot(2,2,3), %在第二行第一列绘图
bar(t,y)
title(‘bar(t,y)’),
pause
subplot(2,2,4),
fill([0,t,4pi],[0,y,0],‘r’) %按[0,t,4pi]，[0,y,0]向量下标渐增次序依次用红色直线连接t,y的数据点
title(‘fill(t,y,’‘r’’)’)
t=0:0.1:2pi;
x=sin(t);
y=cos(t);
z=sin(2t);
plot(t,x,‘k*’,t,y,‘rs’,t,z,‘bo’)
xlabel(‘t(0–2π)’)
legend(‘x=sin(t)’,‘y=cos(t)’,‘z=sin(2t)’)
t=0:10pi;
x=2*(cos(t)+t.sin(t));
y=2(sin(t)-t.cos(t));
z=1.5t;
plot3(x,y,z,‘b’);
xlabel(‘x’);
ylabel(‘y’);
zlabel(‘z’);
title(‘三维螺旋线’);
grid on;
x=-1.5pi:1.5pi;
y=-1.5pi:1.5pi;
z=cos(x).sin(y)./y;
figure(1);
subplot(3,1,1);
plot3(x,y,z);
title(‘三维曲线图’);
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
grid on;
[X,Y]=meshgrid(x);
Z=cos(X).sin(Y)./Y;
subplot(3,1,2);
mesh(X,Y,Z);
title(‘三维网格图’);
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
grid on;
subplot(3,1,3);
surf(X,Y,Z);
shading interp;
title(‘三维光滑曲面图’);
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
grid on;
figure(2);
Z1=Z;
i=find(Z<-0.1);
Z1(i)=NaN;
surf(X,Y,Z1);
title(‘镂空图’);
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
colormap(gray);
R=6;r=2;
u=0:2pi;
v=0:2pi;
X=(R+rcos(u)).cos(v);
Y=(R+rcos(u)).sin(v);
Z=rsin(u);
ezmesh(’(6+2cos(u)).cos(v)’,’(6+2cos(u)).sin(v)’,'2sin(u)’,[0,2pi,0,2pi]);
axis equal;
grid on;
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
hold on;
ezmesh('2cos(u).cos(v)’,'2cos(u).sin(v)’,'2sin(u)’,[0,2pi,0,2*pi]);
axis equal;
title(‘环面及球面’);
山区地貌：在某山区测得一些地点的高程如下表：(平面区域1200<=x<=4000,1200<=y<=3600)，试先利用数据作出该山区的地貌图和等高线图，再利用几种插值方法作出该山区的地貌图和等高线图进行比较。
x=1200:400:4000;
y=1200:400:3600;
z=[ 1130 1250 1280 1230 1040 900 500 700;
1320 1450 1420 1400 1300 700 900 850;
1390 1500 1500 1400 900 1100 1060 950;
1500 1200 1100 1350 1450 1200 1150 1010;
1500 1200 1100 1550 1600 1550 1380 1070;
1500 1550 1600 1550 1600 1600 1600 1550;
1480 1500 1550 1510 1430 1300 1200 980]
[X,Y]=meshgrid(x,y);
subplot(1,2,1)
meshz(X,Y,z);
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
title(‘地貌图’)
subplot(1,2,2)
meshc(X,Y,z);
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
title(‘等高线图’)
xi=linspace(1200,4000,20);
yi=linspace(1200,3600,20);
[x1,y1]=meshgrid(xi,yi);
z1=interp2(x,y,z,x1,y1,‘nearest’);
z2=interp2(x,y,z,x1,y1,‘linear’);
z3=interp2(x,y,z,x1,y1,‘cubic’);
figure(1)
meshc(x1,y1,z1);
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
title(‘带等高线地貌图最邻近插值’)
figure(2)
meshc(x1,y1,z2)
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
title(‘带等高线地貌图双线性插值’)
figure(3)
meshc(x1,y1,z3)
xlabel(‘X’);ylabel(‘Y’);zlabel(‘Z’);
title(‘带等高线地貌图三次多项式插值’)
插值的结果比原本的数据绘图结果精细些，比较接近地貌。三种插值，其中，3次多项式的插值效果最好，地表、等高线均比较平滑；线性插值效果次之，但也比原始图美观；最近点插值效果最差。