按自变量、因变量值绘图

这是老师教的方法。官方文档

这个方法的实现是:将输入的 n 个点中,相邻两个点都用直线连接。当点足够密时,可以认为曲线是平滑的了。需要散点图在最后加一个参数 .

已知表达式

1
2
3
xx = linspace(0, 2*pi, 100);
yy = sin(xx);
plot(xx,yy);

参数方程

1
2
3
4
t = linspace(0,2*pi,100);
x1 = r*cos(t);
y1 = r*sin(t);
plot(x1,y1)

plot 散点图

1
plot(1:5, 1:5, '.'); % 另可用 `o`, `*`

按符号函数绘图

官方文档

fplot 可以直接传函数进去(推荐用符号函数,用匿名函数会警告)。传一个函数即是普通函数,传两个就是参数方程。区间用一个二元向量 [1 2]表示。

1
2
3
4
fplot(f)
fplot(f,xinterval)
fplot(funx,funy)
fplot(funx,funy,tinterval)

ezplot(不知道是什么)

1
2
ezplot(sin(x)/x);
ezplot(t,t^2,[1,2]);

scatter 散点图

后来老师讲了带参数的 plot,真香。

1
2
3
for i = 1:10
scatter(i,rand()), hold on;
end

使用实心圆可以使用 scatter(i,rand(), 8, 'filled')

histogram 直方图

旧版使用 hist,新版使用 histogram

输入一个向量,绘出 x 的每个元素的直方图。

其中,默认把 x 的值分为 10 个区间。

三维线图

plot3,和 plot 用法类似,给定 n 个点的 (x, y, z) 坐标,plot3 绘制其连起来的线图。

三维面图

surf 三维面图要稍微麻烦一点,需要给定 X, Y, Z 三个矩阵(而不是向量),其中 X, Y 可以用 meshgrid 函数和 x, y 的 n 个坐标进行生成。

下面的代码将生成抛物线 $z=x^2+y^2$ 在 [-5,5] 的范围。

1
2
3
4
5
x = -5:0.1:5;
y = -5:0.1:5;
[X, Y] = meshgrid(x, y);
Z = X.^2 + Y.^2;
surf(X,Y,Z);

MATLAB 自带了生成绘制标准球(sphere)、椭球(ellipsoid)和标准圆柱侧面(cylinder)的 X,Y,Z 的函数。

如果想要生成其他大小、位置的球/圆柱,可以考虑坐标变换。下面的代码生成了半径为 R,球心在 [-R, -R, L] 的球。

1
2
3
4
5
[x1, y1, z1] = sphere;
x1 = x1 * R - R;
y1 = y1 * R - R;
z1 = z1 * R + L;
surf(x1, y1, z1);

也可以手算得到 X, Y, Z:

1
2
3
4
5
6
7
8
9
10
11
close all
a = 6000;
b = 5000;
[r,t]=meshgrid(linspace(0,a,50),linspace(0,2*pi,50));
x=r.*cos(t); y=r.*sin(t);
z = b*sqrt(1-(x.^2+y.^2)/(a*a));
z = real(z);
surf(x,y,z);
hold on;
surf(x,y,-z)
axis equal

如果给定函数和范围想要直接画出图形,可以改用 fsurf

mesh 散点图

动画程序框架

动画框架即是播放一张图片,然后暂停 0.05s。

显然,如果绘图过程耗时远小于 0.05s,则可以做出 20 帧的动画。

1
2
3
4
5
6
7
8
9
10
11
close all;
N = 100;
for i=1:N
hold off;
% ... 绘图语句
hold on;
% ... 绘图语句
% ... 绘图语句
% 完成绘图后暂停
pause(0.05);
end;

image 输出图片

1
2
background = imread('Snakes_And_Ladder.png');
image(background);

输出了以后还可以用 hold on; plot 在图上绘制几何图形,每个像素为一个单位长度。但需要注意,该图的原点在左上角,x、y 轴的正方向分别为右方、下方。

提示:还可以使用 axis equal off 来保持图片的长宽比,然后隐藏坐标轴。

输出视频

输出视频的原理是将绘图窗口(gcf 为当前绘图窗口)的图片插入视频,最后导出一个 30 帧的视频。

下面用到了类的写法。

1
2
3
4
5
6
7
mp4 = VideoWriter('highway.mp4','MPEG-4');
mp4.open();
for i=1:N
% 绘图
mp4.writeVideo(getframe(gcf));
end
mp4.close();

图的注释、美化

hold on

使用两个 plot (或其他画图语句)后,第二个会删掉前面的图并且重新画。

要想将两个图叠加在一起,需要在第一次绘图以后加一句 hold on。如:

1
2
3
4
plot(...);
hold on;
plot(...);
plot(...)

如下次绘图不想叠加了,就在这次绘图之后加一句 hold off 即可。

但请注意,hold onhold off 的原理是设置图/表的 NextPlot 属性为 addreplace,并非立即刷新。在下次绘图时,会检查 NextPlot 的属性,再决定是否覆盖掉,还是将图形进行叠加。

因此,下图的代码中的 hold off 其实无意义。

1
2
3
4
5
plot(...);
hold off;
hold on;
plot(...);
plot(...)

另外,很神奇的是,在 MATLAB R2020a 中,rectangle 的语法很诡异。(在 MATLAB R2017a 中没有这个问题)

1
2
3
4
5
6
plot(...);
hold off;
rectangle();
hold on;
plot(...);
plot(...)

多图并排

1
2
3
4
subplot(1,2,1);
scatter(rand(1,10000),rand(1,10000))
subplot(1,2,2);
scatter(rand(1,10000),rand(1,10000))

另外,从 R2019b 开始,您可以使用 tiledlayoutnexttile 函数显示平铺绘图。

1
2
3
4
5
6
7
8
9
10
11
12
x = linspace(0,10);
y1 = sin(x);
y2 = cos(x);
tiledlayout(2,1)

% Top plot
ax1 = nexttile;
plot(ax1,x,y1)

% Bottom plot
ax2 = nexttile;
plot(ax2,x,y2)

坐标轴格式调整

1
2
axis off     % 隐藏坐标轴
axis equal % 使坐标轴的 x, y, z 轴单位长度相同

图的注释

1
2
3
plot(x,y);
title('y=f(x)')
xlabel('x'), ylabel('y')

线图的美化

1
2
3
4
h = plot(x1,y1, 'r-');      % 选择曲线颜色、线型   r 红色,k 黑色,b 蓝色,o 
set(h, 'linewidth' ,2); % 调整线的属性
legend('X', 'Y');
axis([0 25 0 20]) % 调整坐标区的范围

散点图的美化

1
2
3
plot(d(:,1),y,'o','markersize',12)
set(gca,'fontsize',14)
set(gcf,'color','w')