MATLAB函数图形绘制

MATlAB命令	功能
fplot(fun,lims)	绘制符号函数fun在区间lims=[xmin,xmax]间的图像
plot(x,y,s)	绘制由向量x和向量y给定的离散数据连接起来的图像, s用来定义函数曲线的颜色和线型
ezplot(fun,lims)	绘制符号函数图像的简易方法,变量的变化范围lim可以省略, 表示-2pi<x<2pi, 如Fun为二元函数f(x,y),则绘制隐函f(x,y)=0的图像
plot3(X,Y,Z)	绘制三维空间的线点
mesh(Z) mesh(X,Y,Z,C)	绘制着色的三维网纹曲面, 颜色由C决定
meshgrid(x,y)	由向量x和y生成网格点(x,y),与mesh()配合使用
ezmesh(f)	3-D网格图的简单绘制方法, f是一个符号函数
cylinder(R,N)	绘制基于用向量R表示的曲线绕x轴旋转的旋转曲面, 与surf配合使用
subs(S,NEW)	以NEW替代表达式S中的默认变量
subs(S,OLD,NEW)	以NEW替代符号表达式S中的变量OLD OLD 是一个符号变量，字符串量，或字符串表达式； NEW 是符号变量，数值变量或数值表达式

函数plot中参数s的含义
.	point
o	circle
v	x-mark
+	plus
*	star
s	square
d	diamond
v	triangle(down)
^	triangle(up)
<	triangle(left)
>	triangle(right)
p	pentagram
h	hexagram
-	solid
:	dotted
-.	dashdot
--	dashed
(none)	no line
b	blue
g	green
r	red
c	cyan
m	magenta
y	yellow
k	black

在区间[−𝜋, 𝜋]中绘制函数𝑓(𝑥 )= sin2𝑥 + cos$x^2$ 的图形。

x=-pi:0.1:pi;
y=sin(2x)+cos(x.x); 由于x是向量，所以 x.*x
plot(x,y)

在同一个坐标系中绘制函数𝑓(𝑥)= 𝑒−𝑥2和其导函数𝑓′(𝑥)的图形

系统自动求函数的导数
syms x
y=exp(-x*x)
dy=diff(y,x)
x1=-3:0.1:3
y1=subs(y,x1)
dy1=subs(dy,x1)
plot(x1,y1,'r',x1,dy1,'b')

绘制函数𝑓(𝑥)=$x^3$ + 2$x^2$+ $e^x$的图形

先将函数定义为内连函数再使用ezplot将其调用
f=inline('x^3+x2+exp(x)')
ezplot('x^3+x2+exp(x)')

绘制空间曲线

$$
\left{
\begin{array}{lr}
x= & t^3 \
y= & cost & t\in[0,6] \
z= & sin2t
\end{array}
\right.
$$

clear
t=0:0.1:6;
x=t.^3;
y=cos(t);
z=sin(2*t);
plot3(x,y,z)

用红、绿、蓝三种颜色在同一个空间坐标系中, 绘制如下3条曲线：

$$
\left{
\begin{array}{lr}
x= & 4cost \
y= & 4sint & 0 \leq t \leq 20\pi \
z= & 2t
\end{array}
\right.
$$
$$
\left{
\begin{array}{lr}
x= & 6cost \
y= & 4sint & 0 \leq t \leq 20\pi\
z= & 2t
\end{array}
\right.
$$
$$
\left{
\begin{array}{lr}
x= & 4cost \
y= & 6sint & 0 \leq t \leq 20\pi\
z= & 2t
\end{array}
\right.
$$
clear
t=0:0.1:20pi;
x1=4cos(t);
y1=6sin(t);
z1=2t;
x2=6cos(t);
y2=6sin(t);
z2=2t;
x3=8cos(t);
y3=6sin(t);
z3=2t;
plot3(x1,y1,z1,'r', x2,y2,z2,'g', x3,y3,z3,'b')

绘制曲面𝑧 =$\sqrt{x^2 + y^2 }$

clear
s=-10:0.1:10;
t=-10:0.1:10;
[x,y]=meshgrid(s,t);
z=sqrt(x.^2+y.2);
mesh(x,y,z);

绘制由下面二张曲面围成的立体

（1）$z=x^ 2 + y^2$
（2）$z= 3 + \sqrt{9 − (x^2 + y^2)}$
clear
s=-2.3:0.1:2.3;
t=-2.3:0.1:2.3;
s1=-2.3:0.1:2.3;
t1=-2.3:0.1:2.3;
[x,y]=meshgrid(s,t);
[x1,y1]=meshgrid(s1,t1);
z=x.^2+y.2;
mesh(x,y,z);
hold on
z1=abs(3+sqrt(9-(x1.^2+y1.2))); 加abs防止报错无法显示
mesh(x1,y1,z1);
axis equal

旋转单叶双曲面的绘制

s=-pi/5.5:pi/66:pi/5.5
[X,Y,Z]=cylinder(1./(sqrt(cos(s).^2 -sin(s).^2)),60)
surf(X,Y,Z)