MATLAB中绘图命令介绍
- 格式:doc
- 大小:386.50 KB
- 文档页数:11
matlab里面plot用法MATLAB是一款功能强大的编程软件,它提供了很多实用的功能,如绘制图形、数据分析等。
在MATLAB中,plot命令是最常用的绘制图形命令之一,它可以帮助我们绘制各种类型的图形,如折线图、散点图、柱状图等。
本文将详细介绍MATLAB中plot命令的用法。
一、基本用法1. 绘制折线图要使用plot命令绘制折线图,首先需要确定x和y两个变量。
x变量表示横坐标,y变量表示纵坐标。
例如,我们要绘制一组数据的折线图,可以使用以下命令:```matlabx = [1, 2, 3, 4, 5];y = [2, 4, 6, 8, 10];plot(x, y);```这将绘制一条经过点(1, 2)、(2, 4)、(3, 6)、(4, 8)和(5, 10)的折线。
2. 绘制散点图要使用plot命令绘制散点图,同样需要确定x和y两个变量。
例如,我们要绘制一组数据的散点图,可以使用以下命令:```matlabx = [1, 2, 3, 4, 5];y = [2, 4, 6, 8, 10];plot(x, y,'o');```这将绘制一个包含点(1, 2)、(2, 4)、(3, 6)、(4, 8)和(5, 10)的散点图。
'o'表示以点'o'的形式绘制数据点。
3. 绘制柱状图要使用plot命令绘制柱状图,需要确定x和y两个变量,以及一个表示柱子宽度的参数。
例如,我们要绘制一组数据的柱状图,可以使用以下命令:```matlabx = [1, 2, 3, 4, 5];y = [2, 4, 6, 8, 10];bar(x, y);```这将绘制一个包含柱子(宽度分别为2、4、6、8和10)的柱状图,横坐标为[1, 2, 3, 4, 5]。
二、其他用法1. 添加标题和标签要为图形添加标题和标签,可以使用'title'和'xlabel'、'ylabel'命令。
matlab画图常用命令clc 清理命令窗口历史内容clear 清除所有内存存储的变量值clf 清除图形whos 显示各变量信息sqrt 开方edit 开编辑窗口linspace(a,b,N) 定义等差数列,a初值,b末值,N步数(即数据个数)logspace(a,b,N) 定义等比数列,初值10^a,末值10^b,N步数(即数据个数)A.*B 矩阵点乘,对应项相乘A./B 矩阵点除A.^B 矩阵点方(指数相同也要用点方)A=[a:n:b] 定义以a为开始,步长为n的等差数列,最后一个数不超过b(n省略代表步长为1)A' 矩阵转置A=[B,C;D] 矩阵拼凑e *10^exp e^format long 后续数据显示小数点后15位format short 后续数据显示小数点后4位format bank 后续数据显示小数点后2位(不适用于复数)format long/short e 后续数据科学技术法显示,并且小数点后15位/4位format long/short eng 后续数据类似科学技术法显示,但指数保持为3的整数倍,并且有效位数(15位+1/4位+1)format + 矩阵中各元素只显示正负,零为空格format rat 以分数形式显示有理数format long/short g Matlab自定最优显示load/save +文件名载入/储存工作区数据rem(a,b) a/b的余数size(A) A矩阵的大小[行数列数]ylim([0,1])help 打开帮助界面help+帮助界面中对应标题查看对应函数的使用nthroot(x,n) x的n次实数根sign(x) x大于零输出1;x等于零输出0;x小于零输出-1log10(x) lg(x)log(x) ln(x)【注:logb(a)=ln(a)/ln(b)】fix(x) 取整round(x) 对x四舍五入floor(x) 对x向负取整ceil(x) 对x向正取整factor(x) 对x因式分解gcd(a,b) 求a,b最大公约数lcm(a,b) 求a,b最小公倍数rats(x) 用分数表示xfactorial(x) x!nchoosek(n,k) 组合数n选kprimes(x) 找出小于x的素数isprime(x) x是素数,返回1sin(),cos(),tan()... 自变量为弧度asin(),acos(),atan()... 结果为弧度max(x),min(x) x适量中的最大、最小值[a,b]=max(A) A为一行时,a为最大值,b为最大值单行位置A为m行n列时,a为m行向量,对应各列最大值,b为m行向量,对应各列最大值在该列位置多个最大值时,位置默认第一个max(A,B) A、B同大,结果为A,B中对应位置最大值的汇总矩阵mean()/median()/mode() 求平均值/中位数/众数(众数选最小值)cumsum/cumprod(A) 求A的累加/累乘结果,生成与A同大小矩阵,(列运算)单矩阵运算sum/prod(A) A矩阵列求和/求积或求行向量和/积sort(A)/sort(A,'descend') 将A升序/降序排列(行向量自身升序/降序,矩阵列升序/降序)sortrows(A,n) 按第n列排列各行,n正升序,n负降序,n省略第一列升序size(x)/[a,b]=size(x) 返回[行数,列数]/给a、b赋值length(A) 矩阵A的最大长度(行数和列数的最大值)std(A) 求A的标准差(行向量自身求解,矩阵列求解)var(A) 求A的方差(行向量自身求解,矩阵列求解)rand/randn(m,n) 生成(0,1)m×n随机数矩阵/生成均值为0,标准差为1的高斯随机数矩阵(正态分布)【通过randn(m,n)*std+mean可得到均值为mean,标准差为std的正态分布随机数矩阵】A+B*i(complex(A,B)) 生成复数或复数矩阵real(A)/imag(A) 求实部/虚部isreal(A) 实数返回1coni(A) 求共轭【或用A'也可,但会发生行列互换】x为复数时abs(x)/angle(x) 求复数的模,与水平方向的夹角realmax/realmin 返回MATLAB能够使用的最大/最小浮点数intmax/intmin 返回MATLAB能够使用的最大/最小整数pi/i/j 圆周率/虚数/虚数clock 当前时间(一般使用fix(clock)增加可读性)date 返回日期,以字符串形式eps 返回MATLAB最小间隔矩阵A(n,:)/(:,m)【A(n,end)/(end,m)】A矩阵的第n行【最后一列】/第m列【最后一行】[A,B]=meshgrid(a,b),A.*B a,b为行向量,运行结果得a*b的m*n 维矩阵【meshgrid(x)等价于meshgrid(x,x)】zeros(m)/(m,n) m*m/m*n全零矩阵ones(m)/(m,n) m*m/m*n全一矩阵diag(A) 取对角元素为列向量diag(x) 若x为行向量或列向量,结果为对角阵其他元素为零diag(A,n/-n) 对角线右上/左下第n斜线上的元素fliplr(A)/flipud(A) A矩阵列/行进行对称翻转magic(m) 创建m*m维魔方矩阵作图xlabel/ylabel('') 添加x/y轴坐标title('') 添加表头grid 使图像出现网格figure(x) 创建或打开figure x窗口,之后作图均在该窗口进行hold on 保持图像窗口中之前的图像,进而在此作图不会清除之前图像(hold off取消)plot(x1,y1,x2,y2) 同时做两个图像plot(x) x为行向量,则以点数1至n为横轴,x为纵轴作图,按顺序依次连线plot(A) A为m*n矩阵,则图像为那条曲线,每条曲线横轴均为1至m,纵轴为相应列对应值plot(x,A) 以x为横轴,A的每一列为纵轴作图(x与A同维)plot(A,B) A与B需同维,对应列分别作为横轴和纵轴作图plot('标识符') 线型:-实线:点-.点画线--虚线点型:.点o圆圈xx形状+加号*星号s方形d菱形v下三角^上三角<左三角>右三角p五角星h六角星颜色:b蓝色g绿色r红色c青色m洋红色y黄色k黑色w白色【注】,多重输出可多重设定axis([a,b,c,d]) 限制图像x轴在[a,b],y轴在[c,d]legend('string1','string2',etc) 按照作图顺序添加图注text(x,y,'string') 在(x,y)处添加文本‘string’gtext('string') 添加文本‘string’,位置由鼠标点击确定【注】(适用于string形式)输入希腊字母需要'\'+希腊字母读法;^ 可出现上标,_ 可出现下标若想输出_或^,可用\+相应符号subplot(m,n,k) 将图形窗口划分成m行n列,所有的绘图操作都在一行一行数的第k个子图中进行【注】clf针对消除一个figure窗口内的内容,而plot等一系列操作针对一个子图中,且hold on/off被限于特定一个子图中,不影响其他子图polar(x,y) 绘制极图semilogx/semilogy(x,y) x轴对数,y轴线性/x轴线性,y轴对数作图loglog(x,y) 双对数坐标作图bar(x)/barh(x) x为矢量时,按x绘制垂直/水平条形图x为矩阵时,按各行分组绘制垂直/水平条形图bar3(x)/bar3h(x) 同上,绘制三维条形图pie(x)/pie3(x) 绘制(三维)饼状图。
一、二维数据曲线图1、MATLAB 最常用的画二维图形的命令是plot, plor 函数的基本调用格式为:plot(x.y)其 中x 和y 为长度相同的向豈,分别用于存储x 坐标和y 坐标数据。
例 1:在[0,2 7T ]画 Sill(.v) 0生成的图形如下图1所示:图1说明:(1) plot 函数的输入参数是矩阵形式时A 、 当x 是向量,y 是有一维与x 同维的矩阵时,则绘制出多根不同颜色的曲线。
曲线 条数等于y 矩阵的另一维数,x 被作为这些曲线共同的横坐标。
B 、 当x,y 是同维矩阵时.则以x,y 对应列元素为横、纵坐标分别绘制曲线,曲线条数 等于矩阵的列数。
C 、对只包含一个输入参数的plot 函数,当输入参数是实矩阵时,则按列绘制每列元素 值相对其卜.标的曲线,曲线条数等于输入参数矩阵的列数:当输入参数是复数矩阵时,则按 列分别以元素实部和虚部为横、纵坐标绘制多条曲线。
(2) 含多个输入参数的plot 函数 调用格式为:plot(xl,yl.x2,y2,"--.xn.yn)A, 当输入参数都为向量时,xl 和yl, x2和y2, xn 和yn 分别组成一组向量对,每一 组向量对的长度可以不同。
每一向量对可以绘制出一条曲线,这样可以在同一坐标内绘制岀 多条曲线。
B.当输入参数有矩阵形式时,配对的x_y 按对应列兀素为横、纵坐标分别绘制曲线,曲线 条数等于矩阵的列数。
例2:如卜所示的程序:x 1 =liuspace(0,2 *pi,l 00);x2=luispace(0.3 *pi,l 00);x3=linspace(0.4*pi,100);yl=sin(xl); y2=l+sin(x2);y3=2+sin(x3);x=[xl;x2;x3]';0.80.60.40.2-0.2-0.4-0.6-0.8y=[yl;y2;y3「plot(x,y,xl,yl-l) 其图形如图2所示:图2(3)plot函数最简单的调用格式是只包含一个输入参数:plot(x),在这种情况卜,当x是实向量时,以该向量元素的下标为横坐标,元素值为纵坐标画出一•条连续曲线,这实际上是绘制折线图。
∙图形绘制命令∙∙title ——给图形加标题∙xlable ——给x轴加标注∙ylable ——给y轴加标注∙text ——在图形指定的任意位置加标注∙gtext ——利用鼠标将标注加到图形任意位置∙grid on ——打开坐标网格线∙grid off——关闭坐标网格线∙legend ——添加图例∙axis ——控制坐标轴刻度∙hold on; %后续图形叠加显示plot函数绘制二维曲线,常用格式有:plot(x):缺省自变量的绘图格式,x可为向量或矩阵。
plot(x, y):基本格式,x和y可为向量或矩阵。
plot(x1, y1, x2, y2,…):多条曲线绘图格式,在同一坐标系中绘制多个图形。
plot(x, y, ‘s’):开关格式,开关量字符串s设定了图形曲线的颜色、线型及标示符号。
cleart=0:0.1:10;y1=sin(t);y2=cos(t);plot(t,y1,'r',t,y2,'b--');x=[1.7*pi;1.6*pi];y=[-0.3; 0.7];s=['sin(t)';'cos(t)'];text(x, y, s); %指定位置加标注title('正弦和余弦曲线'); %标题legend('正弦','余弦') %添加图例注解xlabel('时间') %x坐标名ylabel('正弦&余弦') %y坐标名grid on %添加网格axis square %将图形设置为正方形采用图形窗口分割方法进行比较显示cleart=0:pi/10:2*pi;y1=sin(t);y2=cos(t);y3=cos(t+pi/2);y4=cos(t+pi);%将图形窗口分割成两行两列,要画的图形为第1行第1列subplot(2,2,1);plot(t,y1);%将图形窗口分割成两行两列,要画的图形为第1行第2列subplot(2,2,2);plot(t,y2);%将图形窗口分割成两行两列,要画的图形为第2行第1列subplot(2,2,3);plot(t,y3);%将图形窗口分割成两行两列,要画的图形为第2行第2列subplot(2,2,4);plot(t,y4);。
Matlab绘图代码以及代码说明文档1.绘制椭圆曲线1)clear:指令,用于清空工作空间2)clc用于清空命令窗口.3)color=’gbkymcrgb’;表示一串字符,可以理解为一个字符的数组(或向量).4)a=4。
5:-0。
5:0.5;上述的a为一个向量(或数组),其取值从4.5开始,每间隔—0。
5取一个数,直到0.5为止。
即a=4.5,4,3.5…0.5;5)for——end是一个循环体,以end结束。
for i=1:1:length(a)表示循环的次数,i从1开始,每次加1,直到length(a)为止,length(a)表示数组(向量)a的长度。
6)a(i)表示a的第i个元素,a(1)=4.5;7)x。
*x表示向量x对应的元素相乘,由于x是向量,因此称号前面有一点。
8)((a(i))^2)表示a的第i个元素的平方。
9)sqrt(a),是一个函数,对a进行开方.10)color(i)是字符数组color的第i个元素,数组前面有定义。
11)hold on是图像保持,就是绘制下一个椭圆时,上一个已经绘制的仍在图形界面上,不会消失。
2.多图形绘制1)sin(t)表示对t求它的正弦,是一个正弦函数.2)subplot(2,2,1),是一个函数,第一个参数2表示将绘图的窗口分割成两行显示;第二个参数2表示将绘图的界面分为两列显示,因此为两行两列4块显示。
第三个参数1表示在第1块(从左向右,从上向下)绘制图像。
3)plot(t1,y1,'.r’)其中,"。
”表示图像的形状为“点”,r为颜色。
4)axis([0,3。
1427,—1,1])表示固定坐标轴,只显示x轴的0到3.1427的区间;y轴的—1到1的区间.5)plot(t1,y1,'b'),由于b前面没有一点,因此是坐标点之间的连线。
6)title(’子图(3)’)表示该图的标题为“子图(3)”,是一个函数,参量为字符串。
Matlab绘图强大的绘图功能是Matlab的特点之一,Matlab提供了一系列的绘图函数,用户不需要过多的考虑绘图的细节,只需要给出一些基本参数就能得到所需图形,这类函数称为高层绘图函数。
此外,Matlab还提供了直接对图形句柄进行操作的低层绘图操作。
这类操作将图形的每个图形元素(如坐标轴、曲线、文字等)看做一个独立的对象,系统给每个对象分配一个句柄,可以通过句柄对该图形元素进行操作,而不影响其他部分。
本章介绍绘制二维和三维图形的高层绘图函数以及其他图形控制函数的使用方法,在此基础上,再介绍可以操作和控制各种图形对象的低层绘图操作。
一.二维绘图二维图形是将平面坐标上的数据点连接起来的平面图形。
可以采用不同的坐标系,如直角坐标、对数坐标、极坐标等。
二维图形的绘制是其他绘图操作的基础。
一.绘制二维曲线的基本函数在Matlab中,最基本而且应用最为广泛的绘图函数为plot,利用它可以在二维平面上绘制出不同的曲线。
1. plot函数的基本用法plot函数用于绘制二维平面上的线性坐标曲线图,要提供一组x坐标和对应的y 坐标,可以绘制分别以x和y为横、纵坐标的二维曲线。
plot函数的应用格式plot(x,y) 其中x,y为长度相同的向量,存储x坐标和y坐标。
例51 在[0 , 2pi]区间,绘制曲线程序如下:在命令窗口中输入以下命令>> x=0:pi/100:2*pi;>> y=2*exp(-0.5*x).*sin(2*pi*x);>> plot(x,y)程序执行后,打开一个图形窗口,在其中绘制出如下曲线注意:指数函数和正弦函数之间要用点乘运算,因为二者是向量。
例52 绘制曲线这是以参数形式给出的曲线方程,只要给定参数向量,再分别求出x,y向量即可输出曲线:>> t=-pi:pi/100:pi;>> x=t.*cos(3*t);>> y=t.*sin(t).*sin(t);>> plot(x,y)程序执行后,打开一个图形窗口,在其中绘制出如下曲线以上提到plot函数的自变量x,y为长度相同的向量,这是最常见、最基本的用法。
MATLAB中的绘图函数介绍概述:MATLAB是一种非常强大的数值计算和科学绘图软件,在各个领域中都得到广泛的应用。
在MATLAB中,绘图函数是其中一个非常重要的功能,它可以帮助我们将数据可视化,并进行分析和解释。
在本文中,我们将详细介绍一些常用的MATLAB绘图函数及其功能。
一、plot函数:plot函数是MATLAB中最基本的绘图函数之一,它可以绘制线性图。
通过将一系列的点连接起来,我们可以绘制出数据的变化趋势。
下面是plot函数的一个简单示例:```matlabx = 0:0.1:10;y = sin(x);plot(x, y);```在这个例子中,我们首先定义了x的取值范围为0到10,间隔为0.1。
然后通过使用sin函数计算出对应的y值。
最后,调用plot函数将x和y的数值传入,即可得到一条关于sin函数的图形。
除了基本的线性图,plot函数还可以绘制不同颜色和线型的曲线,并添加标题、标签等。
它是进行简单数据可视化的利器。
二、scatter函数:相比于plot函数,scatter函数可以绘制散点图,用于展示多个不同数据点之间的分布关系。
通过scatter函数,我们可以方便地比较不同变量之间的相关性。
以下是scatter函数的一个示例:```matlabx = randn(100,1);y = 0.5*x + randn(100,1);scatter(x, y);```在这个例子中,我们首先生成了两组随机数x和y。
然后使用scatter函数将它们绘制成散点图。
通过观察散点图的分布,我们可以判断出x和y之间是否存在线性相关性。
scatter函数还支持设置散点的颜色、大小和透明度等参数,以满足不同的需求。
它是进行多变量分析的重要工具之一。
三、bar函数:bar函数可以用于绘制柱状图,常用于展示各个类别的数据之间的差异。
通过柱状图,我们可以清晰地比较不同类别之间的数值大小。
以下是bar函数的一个示例:```matlabx = categorical({'A', 'B', 'C', 'D'});y = [10, 15, 8, 12];bar(x, y);```在这个例子中,我们首先定义了四个类别,分别是'A'、'B'、'C'和'D'。
绘图指令1 二维曲线图1.1 绘制折线图plot调整坐标范围:axis axis([0,300,0,2])1.2 绘制自定义函数1.3 绘制符号函数1.4 绘制自定义函数1.5 图形修饰设置颜色 y m c r g b w k 设置线型 - : -. --设置标记 . o x + *在指定坐标处,书写文字:1.6 更多类型的二维图e=0.2*rand(size(X));1.7 数值函数的二维图可绘制系统函数,也可绘制自定义函数的图形。
2 三维曲线图2.1 三维曲线plot32.2 三维面填色fill33 曲面图形3.1 网格点坐标的表示3.2 三维网格mesh、meshc、meshz消影开关:hidden on / hidden off利用peaks(50)作为模拟数据矩阵;3.3 着色表面图surf、surfc3.4 二元函数的伪彩色图pcolor3.5 等高线contour3.6 矢量场图quiverZ=X.*exp(-X.^2-Y.^2);3.7 视角控制view视点控制方式及效果:方位角、仰角控制方式及效果:缺省为(-37.5,30)。
3.8 多视区控制subplot3.9 制作、播放动画[x,y,z]=peaks(30); surf(x,y,z)% 制作动画m=moviein(5); % 5帧画面的动画变量for i=1:5view([1 1 i]) % 不断调整视点m(:,i)=getframe; % 将当前画面作为帧保存到序列中end% 播放动画movie(m,2,1); % 每隔1秒播放1帧,循环播放2次3.10 示例:地形图3.11 示例:汶川地形图load dem.txtmesh(dem)4 数据特征的图形展示以pdf结尾的函数,计算概率密度或概率;以cdf结尾的函数,计算累计分布函数。
如:示例:显示t分布与标准正态分布之间的关系。
示例:观察数据的分布特征。
绘图的指令1、plot函数MATLAB中最常用的绘图函数是plot( )。
其基本调用格式为:plot(x, y ),其中x, y均为向量,该函数表示以x向量作为X轴,以y向量作为Y轴。
示例3:绘制从0 ~ 4范围的正弦函数曲线。
t = 0 : pi/20 : 4*pi;y = sin(t);plot(t, y);2、stem函数对于离散序列,MATLAB用stem( )命令实现其绘制。
示例4:绘制从0 ~ 4范围的正弦函数序列。
t = 0 : pi/10 : 4*pi;y = sin(t);stem(t, y);3、subplot函数如果要在一个绘图窗口中显示多个图形,可用subplot函数实现。
其基本调用格式为:subplot(m, n, k) 或subplot(m n k),其中m, n, k取值为1 ~ 9。
该函数表示将绘图窗口划分为m×n个子窗口(子图),并在第k个子窗口中绘图。
示例5:在一个绘图窗口中绘制下列函数的图形。
t = 0 : 0.1 : 2;y1 = 2*exp(-3*t); % (a) 指数函数figure;subplot(221);plot(t,y1);xlabel('(a)');y2 = 2*t.*exp(-3*t); % (b) t乘指数函数subplot(222);plot(t,y2);xlabel('(b)');t1 = -4 : 0.1 : 4;y3 = 1/(2^0.5)*exp(-0.5*t1.^2); % (c) 正态分布函数subplot(223);plot(t1,y3);xlabel('(c)');t2 = -5 : 0.1 : 5;y4 = sinc(t2); % (d) 取样函数subplot(224);plot(t2,y4);xlabel('(d)');ylabel('sinc(t)');axis([-5 5 -0.25 1.1])grid on二维图形的修饰1、坐标轴名称标识函数xlabel、ylabel、title绘制 的nyquist 图和bode 图102.110)(2++=s s s G解MATLAB编程如下:»num=[10];»den=[1,1.2,10];»w=logspace(-2,2,1000); » nyquist(num,den,w)»grid»bode(num,den,w)»grid极小化»fx='10*exp(-x)*cos(x)';»fplot(fx,[2,5])»xmin=fmin(fx,2,5)xmin =2.35619746669214求零点»xzero=fzero(fx,5)xzero =4.71238898038469»xzero=fzero(fx,2)xzero =1.57079632679490曲线拟合»x=0:0.1:1;»y=[-.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2];»p=polyfit(x,y,2);»xi=0:0.01:1;»yi=polyval(p,xi);»plot(x,y,xi,yi) 三维图形»x=-4:0.1:4;»y=x;»[x,y]=meshgrid(x,y);»z = 3*(1-x).^2.*exp(-(x.^2) - (y+1).^2) ...- 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2) ...- 1/3*exp(-(x+1).^2 - y.^2);»mesh(x,y,z)二维图形»x=0:0.01:pi*3;»y=sin(x);»plot(x,y)»title('First Figure Example')»xlabel('Time(s)')»ylabel('Value(v)')»grid on»gtext(‘sinx’)»legend(‘sinx’)管理命令和函数help 在线帮助文件doc 装入超文本说明what M、MA T、MEX文件的目录列表 type 列出M文件lookfor 通过help条目搜索关键字which 定位函数和文件Demo 运行演示程序Path 控制MATLAB的搜索路径管理变量和工作空间Who 列出当前变量Whos 列出当前变量(长表) Load 从磁盘文件中恢复变量 Save 保存工作空间变量Clear 从内存中清除变量和函数Pack 整理工作空间内存Size 矩阵的尺寸Length 向量的长度disp 显示矩阵或与文件和*作系统有关的命令cd 改变当前工作目录Dir 目录列表Delete 删除文件Getenv 获取环境变量值! 执行DOS*作系统命令Unix 执行UNIX*作系统命令并返回结果Diary 保存MA TLAB任务控制命令窗口Cedit 设置命令行编辑 Clc 清命令窗口Home 光标置左上角Format 设置输出格式Echo 底稿文件内使用的回显命令more 在命令窗口中控制分页输出启动和退出MA TLABQuit 退出MATLABStartup 引用MA TLAB时所执行的M文件Matlabrc 主启动M文件*作符和特殊字符+ 加—减* 矩阵乘法.* 数组乘法^ 矩阵幂.^ 数组幂\ 左除或反斜杠/ 右除或斜杠./ 数组除Kron Kronecker张量积: 冒号( ) 圆括号[ ] 方括号. 小数点.. 父目录…继续, 逗号; 分号% 注释! 感叹号‘转置或引用= 赋值= = 相等< > 关系*作符& 逻辑与| 逻辑或~ 逻辑非xor 逻辑异或逻辑函数Exist 检查变量或函数是否存在Any 向量的任一元为真,则其值为真All 向量的所有元为真,则其值为真Find 找出非零元素的索引号三角函数Sin 正弦Sinh 双曲正弦Asin 反正弦Asinh 反双曲正弦Cos 余弦Cosh 双曲余弦Acos 反余弦Acosh 反双曲余弦Tan 正切Tanh 双曲正切Atan 反正切Atan2 四象限反正切Atanh 反双曲正切Sec 正割Sech 双曲正割Asech 反双曲正割Csc 余割Csch 双曲余割Acsc 反余割Acsch 反双曲余割Cot 余切Coth 双曲余切Acot 反余切Acoth 反双曲余切指数函数Exp 指数Log 自然对数Log10 常用对数Sqrt 平方根复数函数Abs 绝对值Argle 相角Conj 复共轭Image 复数虚部Real 复数实部特殊变量和常数Ans 当前的答案Eps 相对浮点精度Realmax 最大浮点数Realmin 最小浮点数Pi 圆周率I,j 虚数单位Inf 无穷大Nan 非数值Flops 浮点运算次数Nargin 函数输入变量数Nargout 函数输出变量数基本X—Y图形Plot 线性图形Loglog 对数坐标图形Semilogx半对数坐标图形(X轴为对数坐标) Semilogy 半对数坐标图形(Y轴为对数坐标) Fill 绘制二维多边形填充图特殊X—Y图形Polar 极坐标图Bar 条形图Stem 离散序列图或杆图Stairs 阶梯图Errorbar 误差条图Hist 直方图Rose 角度直方图Compass 区域图Feather 箭头图Fplot 绘图函数Comet 星点图图形注释Title 图形标题Xlabel X轴标记Ylabel Y轴标记Text 文本注释Gtext 用鼠标放置文本Grid 网格线MATLAB不但擅长於矩阵相关的数值运算,也适合用在各种科学目视表示(Scientific visualization)。
三维绘图2 基本XYZ立体绘图命令●mesh和plot是三度空间立体绘图的基本命令,mesh可画出立体网状图,plot则可画出立体曲面图,两者产生的图形都会依高度而有不同颜色。
下列命令可画出由函数形成的立体网状图:x=linspace(-2, 2, 25); % 在x轴上取25点y=linspace(-2, 2, 25); % 在y轴上取25点[xx,yy]=meshgrid(x, y); % xx和yy都是25x25的矩阵zz=xx.*exp(-xx.^2-yy.^2); % 计算函数值,zz也是21x21的矩阵mesh(xx, yy, zz); % 画出立体网状图●surf和mesh的用法类似:x=linspace(-2, 2, 25); % 在x轴上取25点y=linspace(-2, 2, 25); % 在y轴上取25点[xx,yy]=meshgrid(x, y); % xx和yy都是25x25的矩阵zz=xx.*exp(-xx.^2-yy.^2); % 计算函数值,zz也是25x25的矩阵surf(xx, yy, zz); % 画出立体曲面图●peaks为了方便测试立体绘图,MATLAB提供了一个peaks函数,可产生一个凹凸有致的曲面,包含了三个局部极大点及三个局部极小点,其方程式为:要画出此函数的最快方法即是直接键入peaks:peaksz = 3*(1-x).^2.*exp(-(x.^2) - (y+1).^2) - 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2) - 1/3*exp(-(x+1).^2 - y.^2)●我们亦可对peaks函数取点,再以各种不同方法进行绘图。
meshz可将曲面加上围裙:[x,y,z]=peaks;meshz(x,y,z);●waterfall可在x方向或y方向产生水流效果:[x,y,z]=peaks;waterfall(x,y,z);●下列命令产生在y方向的水流效果:[x,y,z]=peaks;waterfall(x',y',z');●meshc同时画出网状图与等高线:[x,y,z]=peaks;meshc(x,y,z);●surfc同时画出曲面图与等高线:[x,y,z]=peaks;surfc(x,y,z);●contour3画出曲面在三度空间中的等高线:contour3(peaks, 20);●contour画出曲面等高线在XY平面的投影:contour(peaks, 20);plot3可画出三度空间中的曲线:t=linspace(0,20*pi, 501);plot3(t.*sin(t), t.*cos(t), t);亦可同时画出两条三度空间中的曲线:t=linspace(0, 10*pi, 501);plot3(t.*sin(t), t.*cos(t), t, t.*sin(t), t.*cos(t), -t);三维绘图的主要功能:绘制三维线图绘制等高线图绘制伪彩色图绘制三维网线图绘制三维曲面图、柱面图和球面图绘制三维多面体并填充颜色(一)三维线图plot3 ——基本的三维图形指令调用格式:plot3(x,y,z) —— x,y,z是长度相同的向量plot3(X,Y,Z) —— X,Y,Z是维数相同的矩阵plot3(x,y,z,s) ——带开关量plot3(x1,y1,z1,’s1’,x2,y2,z2,’s2’,…)二维图形的所有基本特性对三维图形全都适用。
下面程序是12月26日课题组例会报告的附件内容,叙述了matlab中Cell的概念,列举了画图要用到的一些常用命令。
function [ output_args ] = plot_operation( input_args )%PLOT_OPERATION Summary of this function goes here% 本程序同时写明大规模程序的习惯,就是像学位论文那样分章节,有标题,即Cell的写法。
% 各Cell类似于Word中的各节。
% Cell的定义:在题头使用'%% XXXXXXX'的形式,即“双百分号+空格+Cell名称”的写法。
% 这种写法用工具条中的“双百分号”按钮可以直接查看各节标题,类似于word中的大纲视图。
% 采用'cell'的格式的程序,运行时请点击Ctrl+Shift+Enter,就可以一个Cell一个Cell的运行了。
%% 1. 初始定义clcclearclose all % 关闭所有绘图窗口t=0:0.01:2*pi;y1=sin(t);y2=cos(t);% 设置端点,调试时需要一段一段调试。
%% 2.画图,句柄操作f1=figure(1)l1=plot(t,y1,'.r');hold onl2=plot(t,y2,'b*');handle_legend=legend([l2],'余弦') %只显示第2条线的图例handle_title=title('正弦曲线'); %标题handlex=xlabel('时间/s') %横坐标名称handley=ylabel('alpha_\alpha^\beta') %上标,下标,希腊字母表达式set(gca,'XTicklabel','第1列|第3列| ') %设置坐标轴上的标识handle_t=text(2,-1.1,'sine curve y_1') %在图上显示文字set(gca,'fontsize',30) %设置坐标轴文字大小set(handley,'fontsize',30) %设置坐标轴名称文字大小hold off%% 3. 分图表示f2=figure(2) %上一张图不关掉,重新开一个新图%3.1 开一个6张小图的图片%在第1图位绘制正弦曲线subplot(321)plot(t,y1)title('正弦')%3.2 在第2图位绘制余弦曲线subplot(322)plot(t,y2)title('余弦');%3.3 在第4图位绘制限定坐标范围的余弦曲线subplot(324) %断点绘图plot(t,y1);axis([0,3,-1,1])title('断点绘图,指定坐标范围')%% 4.保存结果%4.1 保存第1幅图,然后关闭saveas(f1,'正余弦曲线.emf')close(f1)%4.2 保存第2幅图,然后关闭figure(2) %打开图片为当前图片% 控制图片大小,输出为emf格式style = hgexport('factorystyle');hgexport(gcf,'size_test',style,'format','meta','width','800','height','600',...'units','points','bounds','loose');% 3个点表示代码换行,对于很长的语句,这中写法便于阅读saveas(gcf,'综合绘图.emf')% gcf 为当前图片的句柄。
MATLAB中绘图命令介绍本节将介绍MATLAB基本xy平面及xyz空间的各项绘图命令,包含一维曲线及二维曲面的绘制。
plot是绘制一维曲线的基本函数,但在使用此函数之前,我们需先定义曲线上每一点的x 及y座标。
下例可画出一条正弦曲线:close all;x=linspace(0, 2*pi, 100); % 100个点的x坐标y=sin(x); % 对应的y坐标plot(x,y);小整理:MATLAB基本绘图函数plot: x轴与y轴均为线性刻度(Linear scale)loglog: x轴与y轴均为对数刻度(Logarithmic scale)semilogx: x轴为对数刻度,y轴为线性刻度semilogy: x轴为线性刻度,y轴为对数刻度若要画出多条曲线,只需将座标对依次放入plot函数即可:hold on 保持当前图形,以便继续画图到当前坐标窗口hold off 释放当前图形窗口title(’图形名称’)(都放在单引号内)xlabel(’x轴说明’)ylabel(’y轴说明’)text(x,y,’图形说明’)legend(’图例1’,’图例2’,…)plot(x, sin(x), x, cos(x));若要改变颜色,在座标对後面加上相关字串即可:plot(x, sin(x), 'c', x, cos(x), 'g');若要同时改变颜色及图线型态,也是在座标对後面加上相关字串即可:plot(x, sin(x), 'co', x, cos(x), 'g*');小整理:plot绘图函数的叁数字元、颜色元、图线型态,y 黄色 .点k 黑色o 圆w 白色x xb 蓝色++g 绿色* *r 红色- 实线c 亮青色: 点线m锰紫色-. 点虚线-- 虚线plot3 三维曲线作图图形完成后,我们可用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围: axis([0, 6, -1.2, 1.2]);axis函数的功能丰富,其常用的用法有:axis equal :纵横坐标轴采用等长刻度axis square:产生正方形坐标系(默认为矩形)axis auto:使用默认设置axis off:取消坐标轴axis on :显示坐标轴此外,MATLAB也可对图形加上各种注解与处理:xlabel('Input Value'); % x轴注解ylabel('Function Value'); % y轴注解title('Two Trigonometric Functions'); % 图形标题legend('y = sin(x)','y = cos(x)'); % 图形注解grid on; % 显示格线我们可用subplot来同时画出数个小图形於同一个视窗之中:subplot(2,2,1); plot(x, sin(x));subplot(2,2,2); plot(x, cos(x));subplot(2,2,3); plot(x, sinh(x));subplot(2,2,4); plot(x, cosh(x));MATLAB还有其他各种二维绘图函数,以适合不同的应用,详见下表。
MATLAB绘制空间曲面的命令一、介绍在科学计算和工程领域中,经常需要对空间曲面进行可视化分析。
MATLAB作为一种强大的数值计算软件,提供了丰富的绘图函数和命令,可以轻松地绘制空间曲面。
本文将介绍MATLAB中常用的绘制空间曲面的命令,并通过示例代码详细说明其使用方法和效果。
二、绘制三维曲面1. meshgrid函数在绘制空间曲面之前,我们需要先创建一个网格。
MATLAB中的meshgrid函数可以用来生成二维或三维的坐标网格,为后续的曲面绘制提供基础。
示例代码:[x, y] = meshgrid(-10:0.5:10, -10:0.5:10);2. 曲面方程在绘制空间曲面时,我们需要定义一个曲面方程。
曲面方程可以是一个显式方程,也可以是一个隐式方程。
示例代码:z = x.^2 + y.^2;3. surf函数一旦我们定义了曲面方程和坐标网格,就可以使用surf函数来绘制空间曲面了。
surf函数可以根据给定的坐标网格和曲面方程,生成一个三维曲面图。
示例代码:surf(x, y, z);4. 其他参数设置除了基本的绘图命令外,还可以通过设置其他参数来调整曲面的外观和显示效果。
例如,可以设置曲面的颜色、透明度、光照等。
示例代码:surf(x, y, z, 'FaceColor', 'interp', 'EdgeColor', 'none');light('Position', [1 1 1]);三、绘制特殊曲面1. 球面球面是一种常见的空间曲面,可以通过给定球心和半径来定义。
MATLAB中的sphere函数可以绘制球面。
示例代码:[x, y, z] = sphere;surf(x, y, z);2. 圆柱面圆柱面是由平行于一个固定直线的所有直线和平行于一个固定平面的所有直线形成的曲面。
MATLAB中的cylinder函数可以绘制圆柱面。
MATLAB中绘图命令介绍本节将介绍MATLAB基本xy平面及xyz空间的各项绘图命令,包含一维曲线及二维曲面的绘制。
plot是绘制一维曲线的基本函数,但在使用此函数之前,我们需先定义曲线上每一点的x 及y座标。
下例可画出一条正弦曲线:close all;x=linspace(0, 2*pi, 100); % 100个点的x坐标y=sin(x); % 对应的y坐标plot(x,y);小整理:MATLAB基本绘图函数plot: x轴和y轴均为线性刻度(Linear scale)loglog: x轴和y轴均为对数刻度(Logarithmic scale)semilogx: x轴为对数刻度,y轴为线性刻度semilogy: x轴为线性刻度,y轴为对数刻度若要画出多条曲线,只需将座标对依次放入plot函数即可:hold on 保持当前图形,以便继续画图到当前坐标窗口hold off 释放当前图形窗口title(’图形名称’)(都放在单引号内)xlabel(’x轴说明’)ylabel(’y轴说明’)text(x,y,’图形说明’)legend(’图例1’,’图例2’,…)plot(x, sin(x), x, cos(x));若要改变颜色,在座标对後面加上相关字串即可:plot(x, sin(x), 'c', x, cos(x), 'g');若要同时改变颜色及图线型态,也是在座标对後面加上相关字串即可:plot(x, sin(x), 'co', x, cos(x), 'g*');小整理:plot绘图函数的叁数字元、颜色元、图线型态,y 黄色.点k 黑色o 圆w 白色x xb 蓝色+ +g 绿色* *r 红色- 实线c 亮青色: 点线m 锰紫色-. 点虚线-- 虚线plot3 三维曲线作图图形完成后,我们可用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围: axis([0, 6, -1.2, 1.2]);axis函数的功能丰富,其常用的用法有:axis equal :纵横坐标轴采用等长刻度axis square:产生正方形坐标系(默认为矩形)axis auto:使用默认设置axis off:取消坐标轴axis on :显示坐标轴此外,MATLAB也可对图形加上各种注解与处理:xlabel('Input Value'); % x轴注解ylabel('Function Value'); % y轴注解title('Two Trigonometric Functions'); % 图形标题legend('y = sin(x)','y = cos(x)'); % 图形注解grid on; % 显示格线我们可用subplot来同时画出数个小图形於同一个视窗之中:subplot(2,2,1); plot(x,sin(x));subplot(2,2,2); plot(x,cos(x));subplot(2,2,3); plot(x,sinh(x));subplot(2,2,4); plot(x, cosh(x));MATLAB还有其他各种二维绘图函数,以适合不同的应用,详见下表。
小整理:其他各种二维绘图函数bar 长条图 errorbar 图形加上误差范围fplot 较精确的函数图形 polar 极座标图hist 累计图 rose 极座标累计图stairs 阶梯图 stem 针状图fill 实心图 feather 羽毛图compass 罗盘图 quiver 向量场图pie,pie3 饼图技巧:对于变化剧烈的函数,可用fplot来进行较精确的绘图对符号函数作图可用ezplot以下我们针对每个函数举例。
当资料点数量不多时,长条图是很适合的表示方式:close all; % 关闭所有的图形视窗x=1:10;y=rand(size(x));bar(x,y);如果已知资料的误差量,就可用errorbar来表示。
下例以单位标准差来做资的误差量:x = linspace(0,2*pi,30);y = sin(x);e = std(y)*ones(size(x));errorbar(x,y,e)对于变化剧烈的函数,可用fplot来进行较精确的绘图,会对剧烈变化处进行较密集的取样,如下例:fplot('sin(1/x)', [0.02 0.2]);% [0.02 0.2]是绘图范围若要产生极座标图形,可用polar:theta=linspace(0, 2*pi);r=cos(4*theta);polar(theta, r);对於大量的资料,我们可用hist来显示资料的分情况和统计特性。
下面几个命令可用来验证randn产生的高斯乱数分:x=randn(5000, 1);% 产生5000个 m=0,s=1 的高斯乱数hist(x,20); % 20代表长条的个数rose和hist很接近,只不过是将资料大小视为角度,资料个数视为距离,并用极座标绘制表示:x=randn(1000, 1);rose(x);stairs 可画出阶梯图: x=linspace(0,10,50); y=sin(x).*exp(-x/3); stairs(x,y);stems 可产生针状图,常被用来绘制数位讯号: x=linspace(0,10,50); y=sin(x).*exp(-x/3); stem(x,y);stairs 将资料点视为多边行顶点,并将此多边行涂上颜色: x=linspace(0,10,50); y=sin(x).*exp(-x/3); fill(x,y,'b'); % 'b'为蓝色 feather 将每一个资料点视复数,并以箭号画出:theta=linspace(0, 2*pi, 20); z = cos(theta)+i*sin(theta); feather(z);compass和feather很接近,只是每个箭号的起点都在圆点:theta=linspace(0, 2*pi, 20);z = cos(theta)+i*sin(theta);compass(z);1. 消隐处理例.比较网图消隐前后的图形z=peaks(50);subplot(2,1,1);mesh(z);title('消隐前的网图')hidden offsubplot(2,1,2)mesh(z);title('消隐后的网图')hidden oncolormap([0 0 1])2. 裁剪处理利用不定数NaN的特点,可以对网图进行裁剪处理例.图形裁剪处理P=peaks(30);subplot(2,1,1);mesh(P);title('裁剪前的网图')subplot(2,1,2);P(20:23,9:15)=NaN*ones(4,7); %剪孔meshz(P) %垂帘网线图title('裁剪后的网图')colormap([0 0 1]) %蓝色网线注意裁剪时矩阵的对应关系,即大小一定要相同.3. 三维旋转体的绘制为了一些专业用户可以更方便地绘制出三维旋转体,MATLAB专门提供了2个函数:柱面函数cylinder和球面函数sphere(1) 柱面图柱面图绘制由函数cylinder实现.[X,Y,Z]=cylinder(R,N) 此函数以母线向量R生成单位柱面.母线向量R是在单位高度里等分刻度上定义的半径向量.N为旋转圆周上的分格线的条数. 可以用surf(X,Y,Z)来表示此柱面.[X,Y,Z]=cylinder(R)或[X,Y,Z]=cylinder此形式为默认N=20且R=[1 1] 例.柱面函数演示举例zyxx=0:pi/20:pi*3; r=5+cos(x);[a,b,c]=cylinder(r,30); mesh(a,b,c)例.旋转柱面图. t=0:pi/12:3*pi;r=abs(exp(-0.25*t).*sin(t)); [a,b,c]=cylinder(r,30); mesh(a,b,c) colormap([1 0 0]) 椭球y=-1:0.1:1;x=2*cos(asin(y)); % 旋转曲面的母线 [X,Y,Z]=cylinder(x,20); % 形成旋转曲面 surf(Z,Y,X); % 画曲面 xlabel('z') ylabel('y') zlabel('x')(2).球面图球面图绘制由函数sphere 来实现[X,Y,Z]=sphere(N) 此函数生成3个(N+1)*(N+1)的矩阵,利用函数 surf(X,Y,Z) 可产生单位球面.[X,Y,Z]=sphere 此形式使用了默认值N=20. Sphere(N) 只是绘制了球面图而不返回任何值. 例.绘制地球表面的气温分布示意图. [a,b,c]=sphere(40); surf(a,b,c);axis('equal') %此两句控制坐标轴的大小相同. axis('square') colormap('hot') 椭球[a,b,c]=sphere(40); surf(c,5*b,2*a); colormap('cool')简单动画a=rand(1,10);b=sin(a);for i=1:10axis([0 1 -1 1])plot(a(i),b(i),'*')pause(0.3);hold onend有meshc(), meshz(), surfc(), surfl(), contourf(), waterfall()等。