常用Matlab作图命令
matlab常用指令
matlab常用指令MATLAB是一款非常实用的科学计算软件，在使用过程中，一些常用的指令是非常必要的。
在本篇文章中，我们将会介绍MATLAB常用指令，以使你更加熟练掌握MATLAB的使用。
一、基本数学运算+ 加- 减* 乘/ 除^ 幂（指数）sqrt 平方根exp 取指数log 取自然对数log10 取以10为底的对数sin 正弦cos 余弦tan 正切asin 反正弦acos 反余弦atan 反正切abs 绝对值rem 模运算fix 向零取整floor 向负无穷取整ceil 向正无穷取整round 四舍五入mod 取摸余数二、变量与矩阵1、赋值：通过等号将数值赋给变量，如：a=3;b=2.1;c=2+3i;2、数列：建立一个等差数组，例如：d=1:10; %1到10的等差数列e=linspace(0,2*pi,100); %0到2*pi之间的100个等间距点 a=[1 2 3;4 5 6;7 8 9];b=zeros(2,3);c=ones(3,2);d=rand(3,3);e=eye(4);4、矩阵元素操作：通过下标访问矩阵中的元素，例如：a(1,2) %输出a矩阵第一行第二列的元素b(2,3)=7 %将b矩阵第二行第三列的元素赋为75、矩阵运算：矩阵加减乘除，如：a+b %对应元素相加a-b %对应元素相减a*b %矩阵乘法a/b %矩阵除法a' %矩阵转置6、矩阵函数：除了使用基本操作外，还能使用各种矩阵相关函数完成矩阵计算，例如：inv(a) %矩阵求逆det(a) %矩阵求行列式trace(a) %矩阵求迹eig(a) %求特征值rank(a) %矩阵的秩size(a) %返回矩阵的大小max(a) %求矩阵元素最大值min(a) %求矩阵元素最小值sum(a) %求矩阵元素的和prod(a) %求矩阵所有元素的乘积mean(a) %求矩阵元素的平均值三、绘图1、二维绘图：绘制二维函数的曲线、散点图等，例如：x=linspace(-3,3,100); %生成-3到3之间的100个等间距点y=sin(x);plot(x,y); %绘制正弦函数曲线plot(x,y,'r--'); %绘制红色的正弦函数曲线，形状为虚线xlabel('x values');ylabel('y values');title('sine function');grid on;四、数据处理1、数据导入：在MATLAB中，可以通过各种方式将数据导入，如：a=load('filename.txt'); %从文件中载入数据b=xlsread('filename.xls'); %从Excel文件中载入数据五、编程1、条件语句：通过条件语句实现程序的分支结构，例如：if(a<0)disp('a is negative');elseif(a==0)disp('a is zero');elsedisp('a is positive');endfor i=1:10disp(i);end3、函数：在MATLAB中，可以自定义函数，函数调用格式为：function [out1,out2,...]=function_name(in1,in2,...)%函数说明%计算过程end4、脚本：在MATLAB中，脚本是一些命令或函数的集合，可以将脚本保存到文件中执行，例如：%脚本说明a=1;b=2;c=a+b;disp(c);以上便是MATLAB一些常用指令的详细介绍。
matlab里面plot用法
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仿真语言基本命令
matlab仿真语言基本命令MATLAB（Matrix Laboratory）是一种用于数值计算、数据分析和可视化的高级编程语言和环境。
以下是一些MATLAB仿真中常用的基本命令：1. 基本操作：•clear: 清除工作区变量。
•clc: 清除命令窗口。
•close all: 关闭所有图形窗口。
•help 或doc: 获取帮助文档。
2. 变量和赋值：• a = 5;: 将值5赋给变量a。
• b = [1, 2, 3];: 创建一个包含元素1、2、3的数组b。
• c = rand(3, 3);: 创建一个3x3的随机矩阵c。
3. 矩阵和数组操作：• A = [1, 2; 3, 4];: 创建一个2x2的矩阵A。
• B = eye(3);: 创建一个3x3的单位矩阵B。
• C = A * B;: 矩阵乘法。
• D = A .* B;: 元素相乘。
4. 绘图：•plot(x, y): 绘制二维图。
•surf(X, Y, Z): 绘制三维曲面图。
•imshow(image): 显示图像。
5. 控制流程：•if、else、elseif: 条件语句。
•for i = 1:5, disp(i); end: for循环。
•while condition, disp('Looping'); end: while循环。
6. 函数：•function result = myFunction(arg1, arg2): 定义函数。
•result = myFunction(2, 3);: 调用函数。
7. 文件和目录操作：•cd: 改变当前工作目录。
•ls 或dir: 列出当前目录的内容。
•save filename.mat variable: 将变量保存到.mat文件中。
•load filename.mat: 从.mat文件加载变量。
8. 统计和数学函数：•mean, std, min, max: 计算均值、标准差、最小值、最大值等。
MATLAB 常用函数3 绘图函数及命令
MATLAB 常用函数3 绘图函数及命令By D. J. Liu1绘图函数(1) plot (二维线图)plot (x, y, ‘r*-’, ‘linewidth’, 5, ‘markersize’, 5)linewidth 设置线条的宽度markersize 设置点的大小(2) plot3 (三维线图)plot3 (X, Y, Z, ‘r*-, ‘linewidth’, 5, ‘markersize’, 5)linewidth 设置线条的宽度markersize 设置点的大小(3) scatter (二维散点图)scatter (X, Y, ‘S’)S设置点的形式、大小及颜色等属性(4) scatter3 (三维散点图)scatter3 (X, Y, Z, ‘S’)S设置点的形式、大小及颜色等属性(5) subplot (子图绘制)subplot(m,n,p)m行, n列, p当前位置(将一个窗口分成m×n个小窗口)(6) mesh (三维网格图)mesh (X, Y, Z)注意：X和Y必须为向量，如果X和Y的长度分别为m和n，则Z必须为m×n的矩阵，即[m，n]=size（Z），在这种情况下网格线的顶点为（X（j），Y（i），Z（i，j））。
(7) surf (三维曲面图)surf (X, Y, Z)surf的调用方法与mesh命令类似，不同的是mesh函数绘制的图像是一个网格图，而surf命令绘制得到的是着色的三维曲面。
着色的方法是在得到相应的网格后，对每个网格依据该网格所代表的节点的色值来定义这一网格的颜色。
注意：第一，surf只支持笛卡尔坐标系（直角坐标系）。
第二，如果要让曲面圆滑，去掉网格，只需要在绘图命令后加入shading interp，这样matlab就会进行相应插值。
注意：X和Y必须为向量，如果X和Y的长度分别为m和n，则Z必须为m×n的矩阵，即[m，n]=size（Z），在这种情况下网格线的顶点为（X（j），Y（i），Z（i，j））(8) pie3(X)（三维饼图）pie3(X) 用X中的数据画一个三维饼形图(百分比例)，X中的每一个元素代表三维饼形图中的一部分。
matlab中figure(1)的用法
在Matlab中，figure(1)是一个常用的指令，它用于创建一个新的图形窗口，并将其标记为1号窗口。
Figure(1)的用法非常灵活，可以用于绘制各种类型的图形，包括折线图、散点图、柱状图等。
在本文中，我们将讨论Figure(1)的基本用法，并结合实际代码和示例，帮助读者更深入地了解如何在Matlab中使用Figure(1)来创建和定制图形。
1. figure(1)的基本语法在Matlab中，可以通过简单的指令figure(1)来创建一个新的图形窗口，该指令的基本语法如下：```matlabfigure(1)```执行这条指令后，Matlab将会创建一个新的图形窗口，并将其标记为1号窗口，以便区分其他图形窗口。
接下来，我们将具体讨论如何在这个新建的图形窗口中绘制各种图形。
2. 绘制折线图Figure(1)可以用于绘制折线图，下面是一个简单的示例代码：```matlabx = 0:0.1:2*pi;y = sin(x);figure(1)plot(x, y)title('Sine Wave')xlabel('x')ylabel('sin(x)')```执行上述代码后，Figure(1)窗口将会显示一个标有“Sine Wave”的折线图，横轴为x，纵轴为sin(x)。
通过这个示例，读者可以了解如何在Figure(1)窗口中绘制简单的折线图，并添加标题、横纵坐标标签。
3. 绘制散点图除了折线图，Figure(1)还可以用于绘制散点图，下面是一个简单的示例代码：```matlabx = rand(1, 100);y = rand(1, 100);figure(1)scatter(x, y)title('Scatter Plot')xlabel('x')ylabel('y')```执行上述代码后，Figure(1)窗口将会显示一个标有“Scatter Plot”的散点图，横轴为x，纵轴为y。
MATLAB曲线绘制大全
一、二维数据曲线图1.1绘制单根二维曲线plot函数的基本调用格式为：plot(x,y)其中x和y为长度相同的向量，分别用于存储x坐标和y坐标数据。
例1-1在0x2p区间内，绘制曲线y=2e-0.5xcos(4x)程序如下：x=0:pi/100:2*pi;y=2*exp(-0.5*x).*cos(4*pi*x);plot(x,y)例1-2绘制曲线。
程序如下：t=0:0.1:2*pi;x=t.*sin(3*t);y=t.*sin(t).*sin(t);plot(x,y);plot函数最简单的调用格式是只包含一个输入参数：plot(x)在这种情况下，当x是实向量时，以该向量元素的下标为横坐标，元素值为纵坐标画出一条连续曲线，这实际上是绘制折线图。
1.2绘制多根二维曲线1.plot函数的输入参数是矩阵形式(1)当x是向量，y是有一维与x同维的矩阵时，则绘制出多根不同颜色的曲线。
曲线条数等于y矩阵的另一维数，x被作为这些曲线共同的横坐标。
(2)当x,y是同维矩阵时，则以x,y对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。
(3)对只包含一个输入参数的plot函数，当输入参数是实矩阵时，则按列绘制每列元素值相对其下标的曲线，曲线条数等于输入参数矩阵的列数。
当输入参数是复数矩阵时，则按列分别以元素实部和虚部为横、纵坐标绘制多条曲线。
2含多个输入参数的plot函数调用格式为：plot(x1,y1,x2,y2,,xn,yn)(1)当输入参数都为向量时，x1和yl,x2和y2,,xn和yn分别组成一组向量对，每一组向量对的长度可以不同。
每一向量对可以绘制出一条曲线，这样可以在同一坐标内绘制出多条曲线。
(2)当输入参数有矩阵形式时，配对的x,y按对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。
例1-3分析下列程序绘制的曲线。
x1=linspace(0,2*pi,100);x2=linspace(0,3*pi,100);x3=linspace(0,4*pi,100);y1=sin(x1);y2=1+sin(x2);y3=2+sin(x3);x=[x1;x2;x3]';y=[y1;y2;y3]';plot(x,y,x1,y1-1)3.具有两个纵坐标标度的图形在MATLAB中，如果需要绘制出具有不同纵坐标标度的两个图形，可以使用plotyy绘图函数。
matlab画图常用命令
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) 绘制（三维）饼状图。
Matlab绘图命令
Matlab绘图命令：1.p lot(x,y):该命令用于描点法作图，给出x的范围，然后利用函数式计算出每一个x对应的y,注意，x间距取得越小，图像会画的越圆滑，这是因为plot描点法作图的原因.如：plot(x,y,’--’)2.f plot(y,[xstart xend]):该命令同样用于绘制图像，但是只需要给出x范围与y关于x的函数式即可，相比于plot 函数更加方便。
3.g rid on可以给图像加上网格4.h old on:用于连续绘图，如果想要在一个坐标系中画出多个图形的话，则需要在每次画完图形之前加入这个命令。
5.l egend(‘y’):在图上标出y曲线线型6.设置图线颜色：plot(x,y,’颜色’)7.a xis ( [xmin xmax ymin ymax] ):设置图像横纵坐标范围例：fplot('cos(x)',[-5 5],'r--'),axis ( [-5 5 -1 1] ),8.s ubplot(m, n, p):用于在一个窗口中绘制多个图像。
例：subplot(1, 2, 1);fplot('sin(x)',[-5 5],'b--');subplot(1, 2, 2);fplot('cos(x)',[-5 5],'r--');9.p lot(thea,r): 绘制极坐标图形，thea 是角变量，r是极轴长度例：a = 2;theta = [0:pi/90:2*pi];r = a*theta;polar(theta,r), title('阿基米德螺线')10.bar(x,y)：画直方图11.stem(t,f)：绘制针头图，这种图经常用来表示波动情况，或者各值和平均值的偏差，可以一目了然。
例：subplot(1, 2, 1)；plot(t,f),xlabel('时间(秒)'),ylabel('弹簧响应')；subplot(1, 2, 2)；stem(t,f),xlabel('时间(秒)'),ylabel('弹簧响应')；12.counter3(x,y,z,n):绘制三维等高图例：[x,y] = meshgrid(-2:0.1:2);z = y.*exp(-x.^2 - y.^2);contour3(x, y, z, 30)；surface(x,y,z,'EdgeColor',[.8 .8 .8],'FaceColor','none')；grid off；view(-15,20)；13.mesh(x,y,z):绘制普通三维图14.surf(x,y,z):绘制颜色渐变三维图15.。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2.3 在原图上继续作图 hold on
[x,y,z] = peaks; pcolor(x,y,z) shading interp hold on contour(x,y,z,20,'k') hold off
2.4 同一图中作多个图
1 0.8 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.2 -0.8 -1 -0.4 0.2 0.6 1
subplot(1, 2, 1); % 定义第一个子图区域。 subplot(1, 2, 2); % 定义第二个子图区域。
Default 10 8 6 4 2 10 8 6 4 2
X scale manipulated
0
2
4
6
8
1
3
7
Y scale manipulated 9 9
Both scales manipulated
data=normrnd (0,1,30,2);
0.3
p=capaplot(data,[-2,2]) p= 0.9199
0.25
0.2
0.15
0.1
0.05
0 -4
-3
-2
-1
0
1
2
3
4
1.10 附加有正态密度曲线的直方图
25
r = normrnd (10,1,100,1);
20
histfit(r)
0.3
0.25
Density
0.2
0.15
0.1
0.05
0
7
8
9
10
11 12 Critical Value
13
14
15
16
normspec([10 Inf],11.5,1.25)
1.12 二项分布的函数图
0.35
p = 0.2; % Probability of success for each trial n = 10; % Number of trials k = 0:n; % Outcomes
0.05
0.04
0.03
0.02
0.01
0 -40
-20
0
20
40
60
80
100
default','3*default') hold off
0.045 0.04 0.035 0.03 0.025 0.02 0.015 0.01 0.005 0 -10 normal epanechnikov box triangle
2.常用作图函数 2.1 普通双函数图
Two Y Axes 3 3
t = 0:pi/20:2*pi; y = exp(sin(t)); plotyy(t,y,t,y,'plot','stem') xlabel('X Axis') ylabel('Plot Y Axis') title('Two Y Axes')
2.5
2.5Biblioteka 22Plot Y Axis
1.5
1.5
1
1
0.5
0.5
0
0
1
2
3 X Axis
4
5
6
0 7
2.2 多数据集在同一图中
1 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.8 -1 sin(x) sin(x-.25) sin(x-.5)
0
1
2
3
4
5
6
7
x = 0:pi/100:2*pi; y = sin(x); y2 = sin(x-.25); y3 = sin(x-.5); plot(x,y,x,y2,x,y3) legend('sin(x)','sin(x-.25)','sin(x-.5)')
0
0.5
1
1.5
2
2.5
3
1.14 ksdensity 概率密度估计函数
Density estimate for MPG 0.06 default width 1/3 default 3*default
cars = load('carsmall','MPG','Origin'); MPG = cars.MPG; [f,x,u] = ksdensity(MPG); plot(x,f) title('Density estimate for MPG') hold on [f,x] = ksdensity(MPG,'width',u/3); plot(x,f,'r'); [f,x] = ksdensity(MPG,'width',u*3); plot(x,f,'g'); legend('default width','1/3
plot(x,y,x,z) gtext('N(0,1)') gtext('N(0,2)') title('正态分布密度曲线')
1.2 绘出 t-分布的密度函数曲线，并与标准正态密度曲线比较
0.4 0.35 0.3 0.25
概率密度 p
t分 布 标准正态密度
0.2 0.15 0.1 0.05 0 -5
1.3 绘制开方分布密度函数在 n 分别等于 1、5、15 的图
0.2 0.18 0.16 0.14 0.12 0.1 0.08 0.06 0.04 0.02 0
x=0:1:30;y1=chi2pdf(x,1); plot(x,y1,':') hold on y2=chi2pdf(x,5);plot(x,y2,'+') y3=chi2pdf(x,15);plot(x,y3,'O') Axis([0,30,0,0.2])
0.2
0.15
0.1
0.05
0
0
1
2
3
4
5
6
7
8
9
10
1.13 指数分布函数图
2 1.8 1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 0
lambda = 2; % Failure rate t = 0:0.01:3; % Outcomes f = exppdf(t,1/lambda); % Probability density vector plot(t,f) % Visualize the probability distribution grid on
6
6
2
0
2
4
6
8
2
1
3
7
clear; x = [1 3 7]; y = [6 9 2]; s1 = subplot(2,2,1); plot(x,y); grid;
title('Default'); s2 = subplot(2,2,2); plot(x,y); set(s2,'XTick',x); ％ 改变 X 轴标记 set(s2,'XGrid','on');％ 画 X 轴的格栅线 title('X scale manipulated'); s3 = subplot(2,2,3); plot(x,y); set(s3,'YTick',[2,6,9]);％ 改变 y 轴标记 set(s3,'YGrid','on');％ 画 y 轴的格栅线 set(s3,'GridLineStyle','-.');％ 使用虚线格栅 title('Y scale manipulated'); s4 = subplot(2,2,4); plot(x,y); set(s4,'XTick',x);％ 改变 xy 轴标记 set(s4,'YTick',[2 6 9]); grid; ％ 画 xy 轴的格栅线 title('Both scales manipulated'); clf reset %Clearing the
0.3
0.25
m = binopdf(k,n,p); % Probability mass vector bar(k,m) % Visualize the probability distribution set(get(gca,'Children'),'FaceColor',[.8 .8 1]) grid on
0
5
10
15
20
25
30
1.4 计算自由度是 50,10 的 F-分布的 0.9 的分位数，并给出概率与分 位数关系的图形
概率与分位数的关系 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 2.1171 2.5 p=0.9
x=finv(0.9,50,10) p=fcdf(x,50,10) t=0:0.1:4; y=fpdf(x,50,10); z=fpdf(t,50,10); plot(t,z,[x,x],[0,y]) text(x,0,'2.1171') gtext('p=0.9')
常用 Matlab 作图命令
1.概率统计作图
1.1 绘出正态分布的密度函数曲线
正态分布密度曲线 0.4 0.35 0.3 N(0,2) 0.25 0.2 0.15 0.1 0.05 0 -5 N(0,1)