MATLAB入门与作图

MATLAB

一、变量与函数

1、变量

（1）变量名区分大小写

（2）变量名必须以字母开头，之后可以是任意字母、数字或下划线，不允许使用标点符号。

2、特殊变量

​ pi 圆周率

3、标点符号

​ 1、每条命令，若为逗号或无标点符号，则显示命令结果；若命令后为分号，则禁止显示结果。

​ 2、% 后面的文字为注释

​ 3、… 表示续行

4、数学函数

		sin(x) 正弦函数     cos(x) 余弦函数    tan(x) 正切函数

​ abs(x) 绝对值 sqrt(x)开平方 log(x)自然对数

​ exp(x)以e为底的指数 fix(x) 取整

5、数据显示格式

		format  short  小数点后4位，缺省显示    format  short e  5位科学记数法

​ format long 15位小数 format long e 15位科学记数法

​ format bank 小数点后两位

6、M文件

​ 第一行为特殊字符

​ 例 ： function 因变量名 = 函数名（自变量名）

​

function   f=fun(x)
f=100*(x(2)-x(1)^2)^2+(1-x(1))^2;
end

二、数组

1、创建数组

​ x = [a b c d e f] 创建指定元素的行向量

​ x= first：last 创建从first开始，加1计数，到last结束的行向量

​ x=first:increement:last 创建从first开始，加increment计数，last结束的行向量

​ x=linspace(first，last，n） 创建从first开始，到last结束，有n个元素的行向量

​ x=logspace(first，last，n） 创建从开始，到结束，有n个元素的对数分隔行向量

2、数组 元素的访问

​ （1）访问一个元素： x(i)表示访问数组x的第i个元素

​ （2）访问一块元素：x(a ：b ：c)表示访问数组x的从第a个元素开始，以步长为b到第c个元素（但不超过

​ c），b可以为负数，b缺损时为1.

​ （3）直接使用元素编址序号. x([a b c d]) 表示提取数组x的第a、b、c、d个元素构成一个新的数组
​ [x(a) x(b) x© x(d)]

3、数组的方向

​ 产生列向量的两种方法：

​ 直接产生 例： c=[1；2；3；4]

​ 转置产生 例 ：

 b=[1 2 3 4];  
 c=b’

​ 说明：以空格或逗号分隔的元素指定的是不同列的元素，而以分号分隔的元素指定了不同行的元素.

4、数组的运算

​ 1、标量-数组运算

​ 设：a=[a1,a2,…,an], c=标量

​ 则：a+c=[a1+c,a2+c,…,an+c]
​ a.c=[a1c,a2c,…,anc]
​ a./c= [a1/c,a2/c,…,an/c](右除）
​ a.\c= [c/a1,c/a2,…,c/an] (左除）
​ a.^c= [a1^c,a2c,…,an^c]
​ c.^a= [c^a1,ca2,…,c^an]

​ 2、数组-数组运算

​ 设：a=[a1,a2,…,an], b=[b1,b2,…,bn]

​ 则：a+b= [a1+b1,a2+b2,…,an+bn]
​ a.b= [a1b1,a2b2,…,anbn]
​ a./b= [a1/b1,a2/b2,…,an/bn]
​ a.\b=[b1/a1,b2/a2,…,bn/an]
​ a.^b=[a1b1,a2^b2,…,anbn]

​ 3、函数

​ length(x) 长度 mean(x) 平均值 median(x) 中值 prod(x) 乘积 sort(x) 从小到大排序

三、矩阵

1、矩阵的建立

​ 逗号或空格用于分隔某一行的元素，分号用于区分不同的行。

​ 特殊矩阵：

​ a=[ ] 产生一个空矩阵，当对一项操作无结果时，返回空矩阵，空矩阵的大小为零.

​ b=zeros(m，n) 产生一个m行、n列的零矩阵

​ c=ones(m，n) 产生一个m行、n列的元素 全为1的矩阵

​ d=eye(m，n) 产生一个m行、n列的单位矩阵

2、矩阵中的元素操作

​ （1）矩阵A的第r行：A（r，：）

​ （2）矩阵A的第r列：A（：，r）

3、矩阵的运算

​ （1）标量-矩阵运算 同标量-数组运算。

​ （2）矩阵-矩阵运算

​ 矩阵加法：A+B
​ 矩阵乘法：A*B
方阵的行列式：det（A）
​ 方阵的逆：inv（A）
​ 方阵的特征值与特征向量：[V，D]=eig[A]

四、控制流

1、for循环

      for  x=array
         commands;
      end

​ 例 对n=1,2,…,10,求Xn=
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Xl10hrUJ-1626608706958)(C:\Users\17117\AppData\Roaming\Typora\typora-user-images\image-20210711200745004.png)]的值。

for n=1:10
   x(n)=sin(n*pi/10);
end
x

2、while循环

​ while expression ​ commands; ​ end

​ 例 设银行年利率为11.25%。将10000元钱存入银行，问多长时间会连本带利翻一番？

clc;clear;
money=10000;
years=0;
while money<20000
   years=years+1
   money=money*(1+11.25/100)
end

3、If-Else-End结构

​ if expression ​ commands; ​ end

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-NjS7EVpG-1626608706959)(C:\Users\17117\AppData\Roaming\Typora\typora-user-images\image-20210711181703651.png)]

function f=fun2(x)
if x>1
         f=x^2+1;
    else if x<=0
          f=x^3;
        else
          f=2*x;
         end
end

​ 先建立M文件fun2.m定义函数f（x），再在Matlab命令窗口输入fun2(2),fun2(0.5), fun2(-1)即可。

Matlab作图

一、二维图型

1、曲线图

Matlab作图是通过描点、连线来实现的，故在画一个曲线图形之前，必须先取得该图形上的一系列的点的坐标（即横坐标和纵坐标），然后将该点集的坐标传给Matlab函数画图.

曲线画图命令为：

plot(x,y,s)

x,y是数组,分别表示点集的横坐标、纵坐标

s对应的语言
·y  黄色   . 点   -  连线
·m  洋红   o 圈   :  短虚线
·c  蓝绿色     x  x-符号       -.  长短线     r  红色        +   加号      --  长虚线

例
plot(X,Y)   %S默认蓝色实线
plot(X1,Y1,S1,X2,Y2,S2,……,Xn,Yn,Sn)
            %将多条线画在一起

线型	符号	-	:	-.	–
	含义	实线	虚线	点划线	双划线
色彩	符号	b	g	r	c
	含义	蓝	绿	红	青
点形	符号	.	+	*	o
	含义	点	加号	星号	圈

设置曲线样式

1.线方式

-	：	-.	–
实线**（solid）**	点线（dotted）	虚点线**（dashdot）**	波折线**（dashed）**

2.点方式

.	+	*****	x	o	s	d	v	^	<	>	p	h
圆****点	加号	星号	x****形	小圆	正方形	菱形	下三角形	上三角	左三角形	右三角形	五角星	六角形

3.颜色

y	r	g	b	w	k	m	c
黄色	红色	绿色	蓝色	白色	黑色	洋红色	青色

图形标注与坐标控制

（1）图形标注

title (’图形名称’)	加图形标题
xlable （’x轴说明’)	加x轴说明
ylable (’y轴说明’)	加y 轴说明
text(x,y,’图形说明’)	在指定位置添加图形说明
legend(‘图例1’，‘图例2，…’)	加图例

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-WzMRYNwI-1626608706961)(C:\Users\17117\AppData\Roaming\Typora\typora-user-images\image-20210711185540931.png)]

clc,clear;
x=0:pi/50:2*pi;        %自变量散点坐标
y1=2*exp(-0.5*x);    %y1散点坐标
y2=cos(4*pi*x);        %y2散点坐标
plot(x,y1,'--ro',x,y2,'--b^');     %画图
title('函数图像');         %加图形标题
xlabel('自变量X’);      %加X轴说明
ylabel('因变量Y’);      %加Y轴说明
text(0.8,1.5,'指数函数'); 
%在指定位置添加图形说明
text(2.5,1.1,'余弦函数'); 
legend('y1:指数函数','y2:余弦函数’);    %加图例

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-27vy0Uac-1626608706963)(C:\Users\17117\AppData\Roaming\Typora\typora-user-images\image-20210711185632526.png)]

也可以用鼠标来确定字符串的位置，方法是输入命令：gtext

2.符号函数(显函数、隐函数和参数方程)画图

(1) ezplot

ezplot(‘f(x)’,[a,b]) 表示在a<x<b绘制显函数f=f(x)的函数图

ezplot(’f(x,y)’,[xmin,xmax,ymin,ymax])，表示在区间xmin<x<xmax和 ymin<y<ymax绘制隐函数f(x,y)=0的函数图

ezplot(’x(t)’,’y(t)’,[tmin,tmax]) ，表示在区间tmin<t<tmax绘制参数方程 x=x(t),y=y(t)的函数图

例3 在[0,pi]上画y=sin(x)的图形

输入命令 ezplot(‘sin(x)’,[0,pi])

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-IQCTEfgd-1626608706964)(C:\Users\17117\AppData\Roaming\Typora\typora-user-images\image-20210711194643049.png)]

输入命令
ezplot(’cos(t)^{3’,’sin(t)}3’，[0,2*pi])

例6 ezplot 绘图应用举例

subplot(2,2,1);
ezplot('x^2+y^2-9’);
axis equal;
subplot(2,2,2);
ezplot('x^3+y^3-5*x*y+1/5’);
subplot(2,2,3);
ezplot('cos(tan(pi*x))',[ 0,1]);
subplot(2,2,4);
ezplot('8*cos(t)','4*sqrt(2)*sin(t)',[0,2*pi]); 

(2) fplot

fplot(‘fun’,lims) 表示绘制字符串fun指定的函数在lims=[xmin,xmax]的图形。

注意：
[1] fun必须是M文件的函数名或是独立变量为x的字符串.
[2] fplot函数不能画参数方程和隐函数图形，但在一个图上可以画多个图形。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-rD6aB8qa-1626608706965)(C:\Users\17117\AppData\Roaming\Typora\typora-user-images\image-20210711194921077.png)]

先自定义函数文件myfun1.m
   function y=myfun1(x)
   y=exp(2*x)+sin(3*x.^2);
    end

在命令行窗口输入命令：
fplot(‘myfun1’,[-1,2])；

例 在[-2,2]范围内绘制函数tanh的图形

  fplot(‘tanh’,[-2,2])；

3、对数坐标图

在很多工程问题中,通过对数据进行对数转换可以更清晰地看出数据的某些特征,在对数坐标系中描绘数据点的曲线,可以直接地表现对数转换.对数转换有双对数坐标转换和单轴对数坐标转换两种.用loglog函数可以实现双对数坐标转换,用semilogx和semilogy函数可以实现单轴对数坐标转换.

loglog(Y) 表示 x、y坐标都是对数坐标系
semilogx(Y) 表示 x坐标轴是对数坐标系
semilogy(…) 表示y坐标轴是对数坐标系
plotyy 有两个y坐标轴，一个在左边，一个在右边

例 用方形标记创建一个简单的loglog

输入命令:
      x=logspace(-1,2);
      loglog(x,exp(x),’-s’)
      grid on     %标注格栅

二、三维图型

1、空间曲线

1、 一条曲线

plot3(x,y,z,s)

y：n维向量，分别表示曲线上点集的横坐标、纵坐标、函数值

s：指定颜色、线型等

例 在区间[0，10*pi]画出参数曲线x=sin(t),y=cos(t), z=t.

t=0:pi/50:10*pi;
      plot3(sin(t),cos(t),t);
      title('螺旋线');   

2、多条曲线

plot3(x,y,z)

其中x，y，z是都是m*n矩阵，其对应的每一列表示一条曲线.

画多条曲线观察函数Z=(X+Y).^2.

x=-3:0.1:3;
y=-4:0.1:5;
[X,Y]=meshgrid(x,y);
Z=(X+Y).^2;
plot3(X,Y,Z);
title('z=(x+y)^2');

##### 2、空间曲面

（这里meshgrid(x,y)的作用是产生一个以向量x为行、向量y为列的矩阵）

2、空间曲面

(1) surf(x,y,z)画出数据点（x，y，z）表示的曲面

y：数据矩阵。分别表示数据点的横坐标、纵坐标、函数值

例 画函数Z=(X+Y).^2的图形.

clc,clear;
x=-3:0.1:3;
y=-4:0.1:5;
[X,Y]=meshgrid(x,y);
Z=(X+Y).^2;
surf(X,Y,Z);
title('z=(x+y)^2');    
shading  flat    %将图形变得平滑

（2） Mesh(x,y,z) 画网格曲面

y：数据矩阵。分别表示数据点的横坐标、纵坐标、函数值

例 画出曲面Z=(X+Y).^2在不同视角的网格图.

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

(3)meshz(X,Y,Z) 在网格周围画一个curtain图(如,参考平面)

例 绘函数的网格图

输入命令:
clc;
clear;
[X,Y]=meshgrid(-3:0.125:3);
Z=(X+Y).^3;
meshz(X,Y,Z);
xlabel('x');
ylabel('y');
zlabel('z');
title('z=(x+y)^3');

三、图形处理

1、在图形上加格栅、图例和标注

（1）GRID ON: 加格栅在当前图上
GRID OFF: 删除格栅

（2）hh = xlabel(‘string’):
在当前图形的x轴上加图例string

​ hh = ylabel(‘string’):
​ 在当前图形的y轴上加图例string

​ hh = zlabel('string’):
​ 在当前图形的z轴上加图例string

​ hh = title(‘string’):
​ 在当前图形的顶端上加图例string

例 在区间[0,2*pi]画sin(x)的图形，并加注图例“自变量
X”、“函数Y”、“示意图”, 并加格栅.

clc;
clear;
x=linspace(0,4*pi,101);
y=sin(x);
plot(x,y,'-r');
xlabel('自变量X');
ylabel('函数Y');
title('函数图像');
legend('y=sin(x)’);
grid on;

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-GBEpimQ2-1626608706965)(C:\Users\17117\AppData\Roaming\Typora\typora-user-images\image-20210711200957316.png)]

（3） hh = gtext(‘string’)

命令gtext(‘string’)用鼠标放置标注在现有的图上.
运行命令gtext(‘string’)时，屏幕上出现当前图形，在
图形上出现一个交叉的十字，该十字随鼠标的移动移动，
当按下鼠标左键时，该标注string放在当前十交叉的位
置.

例 在区间[0,2*pi]画sin(x)，并分别标注“sin(x)”
”cos(x)”.

x=linspace(0,2*pi,30);
 y=sin(x);
 z=cos(x);
 plot(x,y,x,z)
 gtext(‘sin(x)’);
 gtext(’cos(x)’)

2、定制坐标

Axis([xmin xmax ymin ymax zmin zmax]) 定制图形坐标

分别表示x、y、z的最大值最小值

Axis auto将坐标轴返回到自动缺省值

例 在区间[0.005,0.01]显示sin(1/x)的图形。

解  clc;clear;
x=linspace(0.005,0.01,2000);
y=sin(1./x); 
plot(x,y,'--ro');
title('y=sin(1/x)');xlabel('x');
ylabel('y');legend('sin(1/x)')
axis([0.005 0.01 -1 1]);

3、图形保持

(1) hold on 保持当前图形, 以便继续画图到当前图上
hold of 释放当前图形窗口

例 将y=sin(x)、y=cos(x)分别用点和线画出在同一屏幕上。

解  x=linspace(0,2*pi,30);
z=cos(x);
y=sin(x);
plot(x,z,'--rp')
hold on;
plot(x,y,'--bo');

(2) figure(n) 新建窗口，激活图形使其可见，并把它置于其它图形之上

例 区间[0,2*pi]新建两个窗口分别画出y=sin(x)；
z=cos(x)。

解  x=linspace(0,2*pi,100);
y=sin(x);
z=cos(x);
plot(x,y);
title('sin(x)');
xlabel('x'),ylabel('y');
figure(2);
plot(x,z);
title('cos(x)');
xlabel('x'),ylabel('y');

4、分割窗口

h=subplot(mrows,ncols,thisplot)

划分整个作图区域为mrows*ncols块（逐行对块访问）并激活第thisplot块，其后的作图语句将图形画在该块上。

subplot(mrows,ncols,thisplot)

激活已划分为mrows*ncols块的屏幕中的第thisplot块，其后的作图语句将图形画在该块上。

subplot(1,1,1)

命令Subplot(1,1,1)返回非分割状态。

例 将屏幕分割为四块，并分别画出y=sin(x)，z=cos(x)，a=sin(x)*cos(x),b=sin(x)/cos(x)。

clc,clear;
x=linspace(0,2*pi,200);
y1=sin(x); 
y2=cos(x);   
y3=sin(x).*cos(x);
y4=sin(x)./(cos(x)+eps);
subplot(2,2,1);plot(x,y1),title('sin(x)');
subplot(2,2,2);plot(x,y2),title('cos(x)');
subplot(2,2,3);plot(x,y3),title('sin(x)cos(x)');
subplot(2,2,4);plot(x,y4),title('sin(x)/cos(x)');

5、缩放图形

zoom on 为当前图形打开缩放模式。 单击鼠标左键，则在当前图形窗口中，以鼠标点中的点为中心的图形放大2倍；单击鼠标右键，则缩小2倍

zoom off 关闭缩放模式

例 缩放y=sin(x)的图形

x=linspace(0,2*pi,30);
 y=sin(x);
 Plot(x,y);
 zoom on; 

6、改变视角view

（1）view(a,b)
命令view(a,b)改变视角到(a,b),a是方位角,b为仰角。缺省视角为（-37.5，30）。

（2）view（[x，y，z]）

view用空间矢量表示的，三个量只关心它们的比例，与数值的大小无关，x轴view（[1，0，0]），y轴view（[0，1，0]），z轴view（[0，0，1]）。

例 画出曲面Z=(X+Y).^2在不同视角的网格图.

clc,clear;x=-3:0.1:3;  y=1:0.1:5;
[X,Y]=meshgrid(x,y);Z=(X+Y).^2;
subplot(2,2,1);mesh(X,Y,Z);
subplot(2,2,2);mesh(X,Y,Z);view(50,-34);
subplot(2,2,3);mesh(X,Y,Z);view(-60,70);
subplot(2,2,4);mesh(X,Y,Z);view([0,1,1]);

7、动画

Moviein(),getframe,movie()
函数Moviein()产生一个帧矩阵来存放动画中的帧；函数getframe对当前的图象进行快照；函数movie()按顺序回放各帧。

peaks（峰值函数）是MATLAB自带的函数，peaks 函数是典型的多元函数，其本质是二元高斯分布的概率密度函数

例 将曲面peaks做成动画

[x,y,z]=peaks(30);
surf(x,y,z);
axis([-3 3 -3 3 -10 10]);
m=moviein(20);
for i=1:20
   view(-37.5+18*(i-1),30);
   m(:,i)=getframe;
end
movie(m);

四、特殊二三维图形

1、特殊的二维图形函数

（1）极坐标图：polar (theta,rho,s)

用角度theta（弧度表示）和极半径rho作极坐标图，用s指定线型。

例 r=sin2θ×cos2θ的极坐标图形

theta=linspace(0,2*pi,500);                        
        rho=sin(2*theta).*cos(2*theta);
        polar(theta,rho,'g');
        title('Polar plot of sin(2*theta).*cos(2*theta)');

（2）气泡散点图: scatter（X,Y,S,C）

在向量X和Y的指定位置显示彩色圈．X和Y必须大小相同．

例 绘制seamount散点图

输入命令：
     load seamount
     scatter(x,y,5,z)

（3）平面等值线图： contour (x,y,z,n)

绘制n个等值线的二维等值线图

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-nxapELFQ-1626608706966)(C:\Users\17117\AppData\Roaming\Typora\typora-user-images\image-20210713162358035.png)]

输入命令：
[X,Y]=meshgrid(-5:0.1:5,-5:0.1:5);
Z=X.^2+Y.^2;
[C,h]=contour(X,Y,Z);
clabel(C,h)
colormap cool

（4）统计分析图

MATLAB中，二维统计分析图形很多，常见的有条形图、阶梯图、杆图和填充图等

bar（x,y,’选项’）	条形图
stairs （x,y,’选项’）	阶梯图
stem（x,y,’选项’）	杆图
fill （x1,y1,’选项1’,x2,y2,’选项2’,…）	填充图

2、特殊的三维图形函数

（1）空间等值线图： contour 3(x,y,z,n)

其中n表示等值线数。

例 山峰的三维和二维等值线图。

[x,y,z]=peaks;
        subplot(1,2,1);       
        contour3(x,y,z,16,'s');   
        grid;xlabel('x-axis'); ylabel('y-axis');
        zlabel('z-axis');
        title('contour3 of peaks'); 
        subplot(1,2,2);
        contour(x,y,z,16,'s');
        grid; xlabel('x-axis'); ylabel('y-axis');
        title('contour of peaks');

（2）三维散点图 scatter3（X,Y,Z,S,C）

在向量X,Y和Z指定的位置上显示彩色圆圈.
向量X,Y和Z的大小必须相同.

例 绘制三维散点图

输入命令:
[x,y,z]=sphere(16); % 生成 16 x 16个面的球面
X=[x(:)*.5 x(:)*.75 x(:)];
Y=[y(:)*.5 y(:)*.75 y(:)];
Z=[z(:)*.5 z(:)*.75 z(:)];
S=repmat([1 .75 .5]*10,prod(size(x)),1);
C=repmat([1 2 3],prod(size(x)),1);
scatter3(X(:),Y(:),Z(:),S(:),C(:),'filled'),view(-60,60)

其中n表示等值线数。

例 山峰的三维和二维等值线图。

[x,y,z]=peaks;
        subplot(1,2,1);       
        contour3(x,y,z,16,'s');   
        grid;xlabel('x-axis'); ylabel('y-axis');
        zlabel('z-axis');
        title('contour3 of peaks'); 
        subplot(1,2,2);
        contour(x,y,z,16,'s');
        grid; xlabel('x-axis'); ylabel('y-axis');
        title('contour of peaks');

（2）三维散点图 scatter3（X,Y,Z,S,C）

在向量X,Y和Z指定的位置上显示彩色圆圈.
向量X,Y和Z的大小必须相同.

例 绘制三维散点图

输入命令:
[x,y,z]=sphere(16); % 生成 16 x 16个面的球面
X=[x(:)*.5 x(:)*.75 x(:)];
Y=[y(:)*.5 y(:)*.75 y(:)];
Z=[z(:)*.5 z(:)*.75 z(:)];
S=repmat([1 .75 .5]*10,prod(size(x)),1);
C=repmat([1 2 3],prod(size(x)),1);
scatter3(X(:),Y(:),Z(:),S(:),C(:),'filled'),view(-60,60)