matlab数理统计

Matlab在《概率论与数理统计》教学中的应用
Matlab提供了丰富的概率分布函数，可以帮助学生更好地理解不同的概率分布。

学生可以使用Matlab生成正态分布、二项分布、泊松分布等不同的概率分布，并画出相应的概率密度函数、累积分布函数等图形。

通过实际的计算和绘图，学生可以更直观地看到不同概率分布的特点，加深对概率分布的理解。

Matlab提供了各种统计函数，可以方便地进行数据的描述性统计和推断性统计。

学生可以使用Matlab计算样本的平均值、方差等描述性统计量，还可以使用Matlab进行假设检验、置信区间估计等推断性统计。

通过实际的计算和分析，学生可以更好地掌握统计学中的概念和方法。

Matlab还可以进行模拟实验，帮助学生理解概率和统计的原理。

学生可以使用Matlab 模拟抛硬币的实验，验证概率的定义和性质。

学生还可以使用Matlab模拟中心极限定理，观察样本均值的分布趋于正态分布的情况。

通过实际的模拟实验，学生可以更深入地理解抽样分布和极限定理等重要概念。

Matlab还可以用于数据的可视化。

学生可以使用Matlab绘制直方图、散点图、箱线图等图形，展示数据的分布和变化。

通过可视化的方式，学生可以更好地理解数据的特点和规律，并能够更直观地展示和解释统计分析的结果。

Matlab在《概率论与数理统计》教学中具有广泛的应用价值。

通过利用Matlab进行计算、模拟和可视化等任务，可以帮助学生更好地理解概率和统计的概念和方法，提高学习效果。

在教学中合理地使用Matlab可以有效地促进学生对概率论与数理统计的学习和理解。

Matlab 概率论与数理统计、matlab 基本操作 1. 画图【例01.01】简单画图hold off; x=0:0.1:2*pi; y=sin (x);plot(x,y, '-r'); x1=0:0.1:pi/2; y1=s in( x1); hold on;fill([x1, pi/2],[y1,1/2],'b');【例01.02】填充，二维均匀随机数hold off ;x=[0,60];y0=[0,0];y60=[60,60]; x1=[0,30];y1=x1+30; x2=[30,60];y2=x2-30;plot(x,y0, 'r' ,y0,x, plot(x1,y1, 'r' ,x2,y2, yr=u nifrnd (0,60,2,100);plot(yr(1,:),yr(2,:), axis( 'on'); axis( 'square' ); axis([-20 80 -20 80 ]);xv=[0 0 30 60 60 30 0];yv=[0 30 60 60 30 0 0]; fill(xv,yv, 'b');hold on ;'r' ,x,y60, 'r' ,y60,x,'r')'r');'m.')2. 排列组合kC=nchoosek(n,k) : CC n ，例 nchoosek(5,2)=10, nchoosek(6,3)=20.prod(n1:n2):从 n1 至U n2 的连乘【例01.03】至少有两个人生日相同的概率365 364|||(365 rs 1)rs365365 364 365 rs 1 365 365365rs=[20,25,30,35,40,45,50]; %每班的人数p1= on es(1,le ngth(rs)); p2=on es(1,le ngth(rs));%用连乘公式计算for i=1:le ngth(rs) p1(i)=prod(365-rs(i)+1:365)/365A rs(i); end%用公式计算(改进) for i=1:le ngth(rs)for k=365-rs(i)+1:365p2(i)=p2(i)*(k/365); end ; end%用公式计算(取对数) for i=1:le ngth(rs)p1(i)=exp(sum(log(365-rs(i)+1:365))-rs(i)*log(365)); end公式计算P 1n!C NN nN!1 (N n)!1N nN (N 1) (N n 1)、随机数的生成3. 均匀分布随机数rand(m,n);产生m行n列的(0,1)均匀分布的随机数rand(n);产生n行n列的(0,1)均匀分布的随机数【练习】生成(a,b)上的均匀分布4. 正态分布随机数randn(m,n); 产生m行n列的标准正态分布的随机数【练习】生成N(nu,sigma42)上的正态分布5. 其它分布随机数三、一维随机变量的概率分布1. 离散型随机变量的分布率(1) 0-1分布(2) 均匀分布_ k k n k(3) 二项分布：binopdf(x,n,p)，若X ~ B(n, p)，则P{X k} C n p (1 p),x=0:9 ;n=9;p=0.3;y= bin opdf(x ,n, p);plot(x,y,'b-',x,y,'r*')y=[ 0.0404, 0.1556, 0.2668, 0.2668, 0.1715, 0.0735, 0.0210, 0.0039, 0.0004, 0.0000 ]当n较大时二项分布近似为正态分布x=0:100; n=100;p=0.3;y= bin opdf(x ,n, p);plot(x,y,'b-',x,y,'r*')ke⑷泊松分布：piosspdf(x, lambda)，若X ~ ()，贝U P{ X k}k!x=0:9; lambda = 3;y= poisspdf (x,lambda);plot(x,y,'b-',x,y,'r*')y=[ 0.0498, 0.1494, 0.2240, 0.2240, 0.1680, 0.1008, 0.0504, 0.0216, 0.0081,0.0027]k 1⑸几何分布：geopdf (x, p)，贝U P{X k} p(1 p)x=0:9;p=0.3y= geopdf(x,p);plot(x,y,'b-',x,y,'r*')y=[ 0.3000, 0.2100, 0.1470, 0.1029, 0.0720, 0.0504, 0.0353, 0.0247, 0.0173, 0.0121 ] x=0:10;N=20;M=8; n=4;y= hygepdf(x,N,M, n); plot(x,y,'b-',x,y,'r*')y=[ 0.1022, 0.3633, 0.3814, 0.1387, 0.0144, 0, 0, 0, 0, 0, 0 ]2. 概率密度函数(1)均匀分布：unifpdf(x,a,b) , f (x)其它a=0;b=1;x=a:0.1:b; y= uni fpdf (x,a,b);1 2 厂(x )2 ■厂ex=-10:0.1:12;mu=1;sigma=4;y= no rmpdf(x,mu,sigma);rn=10000;z= normrnd (mu,sigma,1,rn); % 产生 10000 个正态分布的随机数 d=0.5;a=-10:d:12;b=(hist(z,a)/rn)/d;%以a 为横轴，求出10000个正态分布的随机数的频率(6)超几何分布：hygepdf(x,N,M,n)，则 P{Xk}C k nM CNC N(2)正态分布:normpdf(x,mu,sigma) , f (x)plot(x,y,'b-',a,b,'r.')1 _x⑶指数分布：exppdf(x,mu)， f (x)其它x=0:0.1:10;mu=1/2;■ t京■I_ey= exppdf(x,mu); plot(x,y,'b-',x,y,'r*')1n i F⑷2分布：chi2pdf(x,n) , f (x; n) 2n ^( n 2) % e x 0hold onx=0:0.1:30;n=4;y= chi2pdf(x, n);plot(x,y,'b');%blue n=6;y= chi2pdf(x, n);plot(x,y,'r');%red n=8;y=chi2pdf(x ,n );plot(x,y,'c');%cya n n=10;y= chi2pdf(x, n);plot(x,y,'k');%black lege nd(' n=4', 'n=6', 'n=8', 'n=10');n 1((n 1) 2) x2 2⑸t 分布：tpdf(x,n) , f (x; n) ------------------ 1 -J n (n. 2) nhold onx=-10:0.1:10;n=2;y= tpdf(x, n);plot(x,y,'b');%bluen=6;y= tpdf(x, n);plot(x,y,'r');%redn=10;y= tpdf(x ,n );plot(x,y,'c');%cya nn=20;y= tpdf(x, n);plot(x,y,'k');%black lege nd(' n=2', 'n=6', 'n=10', 'n=20');((m山m 门2n2) 2)小2% 2 1 5 % 2(n2 2) n2n2x 0(6) F 分布：fpdf(x,n1,n2) , f (x; n「n2) (E 2)0 x 0hold onx=0:0.1:10;n1=2; n2=6;y= fpdf(x, n1, n2);plot(x,y,'b');%bluen1=6; n2=10;y= fpdf(x, n1, n2);plot(x,y,'r');%red n1=10; n2=6;y= fpdf(x, n1, n2);plot(x,y,'c');%cyann1=10; n2=10;y= fpdf(x, n1,n 2);plot(x,y,'k');%black legend(' n仁2; n2=6', ' n1= 6; n2=10', ' n仁10;n2=6', ' n仁10; n2=10');3.分布函数F(x) P{X x}【例03.01】求正态分布的累积概率值设X ~ N(3,22)，求 P{2 X 5}, P{ 4 X 10}, P{ X 2}, P{X 3},14.逆分布函数，临界值y F(x) P{X x} , x F (y) , x称之为临界值【例03.02】求标准正态分布的累积概率值y=0:0.01:1;x=normin v(y,0,1);【例03.03】求2(9)分布的累积概率值hold offy=[0.025,0.975];x=ch i2in v(y,9);n=9;x0=0:0.1:30;y0=chi2pdf(x0, n); plot(x0,y0, 'r'); x1=0:0.1:x(1);y1=chi2pdf(x1, n);x2=x(2):0.1:30;y2=chi2pdf(x2 ,n);hold onfill([x1, x(1)],[y1,0], 'b');fill([x(2),x2],[0,y2], 'b');【练习1.1】二项分布、泊松分布、正态分布(1)对n 10, p 0.2二项分布，画出b(n,p)的分布律点和折线；(2)对np，画出泊松分布()的分布律点和折线；(3)对np, 2叩(1 p)，画出正态分布N( , 2)的密度函数曲线;(4)调整n, p，观察折线与曲线的变化趋势。

《概率论与数理统计》MATLAB上机实验实验报告一、实验目的1、熟悉matlab的操作。

了解用matlab解决概率相关问题的方法。

2、增强动手能力，通过完成实验内容增强自己动手能力。

二、实验内容1、列出常见分布的概率密度及分布函数的命令，并操作。

概率密度函数分布函数(累积分布函数) 正态分布normpdf(x,mu,sigma) cd f(‘Normal’,x, mu,sigma);均匀分布（连续）unifpdf(x,a,b) cdf(‘Uniform’,x,a,b);均匀分布（离散）unidpdf(x,n) cdf(‘Discrete Uniform’,x,n);指数分布exppdf(x,a) cdf(‘Exponential’,x,a);几何分布geopdf(x,p) cdf(‘Geometric’,x,p);二项分布binopdf(x,n,p) cdf(‘Binomial’,x,n,p);泊松分布poisspdf(x,n) cdf(‘Poisson’,x,n);2、掷硬币150次，其中正面出现的概率为0.5，这150次中正面出现的次数记为X(1) 试计算X=45的概率和X≤45 的概率；(2) 绘制分布函数图形和概率分布律图形。

答:(1)P(x=45)=pd =3.0945e-07P(x<=45)=cd =5.2943e-07(2)3、用Matlab软件生成服从二项分布的随机数，并验证泊松定理。

用matlab依次生成(n=300,p=0.5),(n=3000,p=0.05),(n=30000,p=0.005)的二项分布随机数，以及参数λ=150的泊松分布，并作出图线如下。

由此可以见得，随着n的增大，二项分布与泊松分布的概率密度函数几乎重合。

因此当n足够大时，可以认为泊松分布与二项分布一致。

4、 设22221),(y x e y x f +−=π是一个二维随机变量的联合概率密度函数，画出这一函数的联合概率密度图像。

Matlab数理统计工具箱应用简介1．概述Matlab的数理统计工具箱是Matlab工具箱中较为简单的一个，其牵扯的数学知识是大家都很熟悉的数理统计，因此在本文中，我们将不再对数理统计的知识进行重复，仅仅列出数理统计工具箱的一些函数，这些函数的意义都很明确，使用也很简单，为了进一步简明，本文也仅仅给出了函数的名称，没有列出函数的参数以及使用方法，大家只需简单的在Matlab工作空间中输入“help 函数名”，便可以得到这些函数详细的使用方法。

2．参数估计betafit 区间3．累积分布函数betacdf β累积分布函数binocdf 二项累积分布函数cdf 计算选定的累积分布函数chi2cdf 累积分布函数2χexpcdf 指数累积分布函数fcdf F累积分布函数gamcdf γ累积分布函数geocdf 几何累积分布函数hygecdf 超几何累积分布函数logncdf 对数正态累积分布函数nbincdf 负二项累积分布函数ncfcdf 偏F累积分布函数nctcdf 偏t累积分布函数ncx2cdf 偏累积分布函数2χnormcdf 正态累积分布函数poisscdf 泊松累积分布函数raylcdf Reyleigh累积分布函数tcdf t 累积分布函数unidcdf 离散均匀分布累积分布函数unifcdf 连续均匀分布累积分布函数weibcdf Weibull累积分布函数4．概率密度函数betapdf β概率密度函数binopdf 二项概率密度函数chi2pdf 概率密度函数2χexppdf 指数概率密度函数fpdf F概率密度函数gampdf γ概率密度函数geopdf 几何概率密度函数hygepdf 超几何概率密度函数lognpdf 对数正态概率密度函数nbinpdf 负二项概率密度函数ncfpdf 偏F概率密度函数nctpdf 偏t概率密度函数ncx2pdf 偏概率密度函数2χnormpdf 正态分布概率密度函数pdf 指定分布的概率密度函数poisspdf 泊松分布的概率密度函数raylpdf Rayleigh概率密度函数tpdf t概率密度函数unidpdf 离散均匀分布概率密度函数unifpdf 连续均匀分布概率密度函数weibpdf Weibull概率密度函数5．逆累积分布函数Betainv 逆β累积分布函数binoinv 逆二项累积分布函数chi2inv 逆累积分布函数2χexpinv 逆指数累积分布函数finv 逆F累积分布函数gaminv 逆γ累积分布函数geoinv 逆几何累积分布函数hygeinv 逆超几何累积分布函数logninv 逆对数正态累积分布函数nbininv 逆负二项累积分布函数ncfinv 逆偏F累积分布函数nctinv 逆偏t累积分布函数ncx2inv 逆偏累积分布函数2χnorminv 逆正态累积分布函数possinv 逆正态累积分布函数raylinv 逆Rayleigh累积分布函数tinv 逆t累积分布函数unidinv 逆离散均匀累积分布函数unifinv 逆连续均匀累积分布函数weibinv 逆Weibull累积分布函数6．分布矩函数betastat 计算β分布的均值和方差binostat 二项分布的均值和方差chi2stat 计算分布的均值和方差2χexpstat 计算指数分布的均值和方差fstat 计算F分布的均值和方差gemstat 计算γ分布的均值和方差geostat 计算几何分布的均值和方差hygestat 计算超几何分布的均值和方差lognstat 计算对数正态分布的均值和方差nbinstat 计算负二项分布的均值和方差ncfstat 计算偏F分布的均值和方差nctstat 计算偏t分布的均值和方差ncx2stat 计算偏分布的均值和方差2χnormstat 计算正态分布的均值和方差poissstat 计算泊松分布的均值和方差raylstat 计算Rayleigh分布的均值和方差tstat 计算t分布的均值和方差unidstat 计算离散均匀分布的均值和方差unifstat 计算连续均匀分布的均值和方差weibstat 计算Weibull分布的均值和方差7．统计特征函数corrcoef 计算互相关系数cov 计算协方差矩阵geomean 计算样本的几何平均值harmmean 计算样本数据的调和平均值iqr 计算样本的四分位差kurtosis 计算样本的峭度mad 计算样本数据平均绝对偏差mean 计算样本的均值median 计算样本的中位数moment 计算任意阶的中心矩prctile 计算样本的百份位数range 样本的范围skewness 计算样本的歪度std 计算样本的标准差trimmean 计算包含极限值的样本数据的均值var 计算样本的方差8．统计绘图函数boxplot 在矩形框内画样本数据errorbar 在曲线上画误差条fsurfht 画函数的交互轮廓线gline 在图中交互式画线gname 用指定的标志画点lsline 画最小二乘拟合线normplot 画正态检验的正态概率图pareto 画统计过程控制的Pareto图qqplot 画两样本的分位数-分位数图refcurve 在当前图中加一多项式曲线refline 在当前坐标中画参考线surfht 画交互轮廓线weibplot 画Weibull概率图9．统计处理控制capable 处理能力索引capaplot 画处理能力图ewmaplot 画指数加权移动平均图histfit 叠加正态密度直方图normspec 在规定的极限内画正态密度图schart 画标准偏差图xbarplot 画水平条图10．假设检验Ranksum 计算母体产生的两独立样本的显著性概率和假设检验的结果signrank 计算两匹配样本中位数相等的显著性概率和假设检验的结果signtest 计算两匹配样本的显著性概率和假设检验的结果ttest 对单个样本均值进行t检验ttest2 对两样本均值差进行t检验ztest 对已知方差的单个样本均值进行z检验11．试验设计cordexch 配位交叉算法D-优化试验设计daugment D-优化增强试验设计dcovary 使用指定协变数的D-优化试验设计ff2n 两水平全因素试验设计fullfact 全因素试验设计hadamard Hadamard正交试验rowexch 行交换算法D-优化试验设计。

Matlab在数理统计中的运用摘要：概率论与数理统计是现代数学的重要分支，近年来随着计算机的普及，概率论在经济，管理，金融，保险，生物，医学等方面都发挥着越来越大的作用。

使得概率统计成为今天各类各专业大学生最重要的数学必修课之一。

然而,传统的概率统计教学过于偏重理论的阐述、公式的推导、繁琐的初等运算;同时,缺乏与计算机的结合,给学生的学习带来很多困难。

本文介绍概率统计中的主要问题在Matlab中的实现,让我们从繁琐的计算中解放出来,把更多的时间和精力用于基本概念和基本理论的思考和方法的创新,从而提高教师的教学效率和学生的学习效率。

关键词：区间估计，matlab，概率统计一、常用概率密度的计算Matlab中计算某种概率分布在指定点的概率密度的函数,都以代表特定概率分布的字母开头,以pdf (probability density function)结尾,例如:unidpdf(X, N):计算1到N上的离散均匀分布在X每一点处的概率密度;poisspdf(X, Lambda):计算参数为Lambda的泊松分布在X每一点处的概率密度;exppdf(X, mu):计算参数为mu的指数分布在X每一点处的概率密度;normpdf(X, mu, sigma):计算参数为mu, sigma的正态分布在X每一点处的概率密度。

其他如连续均匀分布、二项分布、超几何分布等也都有相应的计算概率密度的函数。

除计算概率密度的函数外,Matlab中还有计算累积概率密度、逆概率分布函数及产生服从某分布的随机数的函数,分别以cdf,inv和rnd结尾。

下面我们来用一个具体的例子说明一下：例1：计算正态分布N（0，1）的随机变量X在点0.6578的密度函数值。

解：>> pdf('norm',0.6578,0,1)ans =0.3213例2：自由度为8的卡方分布，在点2.18处的密度函数值。

解：>> pdf('chi2',2.18,8)ans = 0.0363二、随机变量数字特征的计算(一)数学期望与方差对离散型随机变量,可利用Matlab矩阵运算计算出其数学期望和方差;而对于连续型随机变量,则可以利用Matlab符号运行计算。

MATLAB 在概率论和数理统计中的应用一、 引言概率论与数理统计作为现代数学的重要分支，在自然科学、社会科学和工程技术等领域都具有极为广泛的应用。

概率论和数理统计是研究随机现象的客观规律并付诸应用的数学学科。

用概率论和数理统计的知识来解决实际问题时，大致遵循以下流程图。

实际问题数学表达概率论与数理统计模型符合实际结束分析求解 检验是否随着计算机技术的普及和开展，我们可以用计算机语言轻松的完成以上过程中的求解和建立模型过程。

可以大大提高准确率和使用者的效率。

二、 MATLAB 软件介绍及其特点1984年美国MathWorks 公司推出了MATLAB 软件。

MATLAB 是以矩阵作为数据操作的根本单位的程序设计语言,是主要面对科学计算、数据可视化、系统仿真，以及交互式程序设计的高科技计算环境。

为科学研究、工程设计以及必须进展有效数值计算的众多科学领域提供了一种全面的解决方案。

MATLAB 软件具有以下特点[1]：I,具有强大的数值计算和符号运算功能II，操作界面简单，编程语言自然III，具有先进的数据可视化功能IV，具有强大的开放性和可扩展性总之，MATLAB是工程师和科研者使用最广泛的软件之一。

三、MATLAB中关于概率统计的命令和函数MATLAB中的Statistics Toolbox提供了丰富的关于概率统计的命令和函数，用于解决概率论和数理统计中的常见问题。

下表将列举常用的概率统计中的命令和函数。

利用上述函数产生一个44矩阵的标准正态随机数，MATLAB代码如下：>> R=normrnd(0,1,4,4)R =-0.8095 -0.7549 -0.2414 -0.0301 -2.9443 1.3703 0.3192 -0.16491.4384 -1.7115 0.3129 0.62770.3252 -0.1022 -0.8649 1.0933 产生的随机数可以在工作窗口查看。

假设要想求参数为=1=2=2=3μσμσ，；，的正态分布的期望和方差，那么相应的MATLAB 的代码为： >> clear >> a=[1 2]; >> b=[2 3];>> [m v]=normstat(a,b) m =1 2 v =4 9在MATLAB 的统计工具箱中提供了一个演示程序disttool ，可以直观的演示常见分布的分布函数和概率密度函数。

MATLAB基础知识一、MATLAB软件简介1967年美国Mathwork公司推出了、基于矩阵运算的“Matrix Laboratory”(缩写为MATLAB) 的交互式软件包. MATLAB既是一种直观、高效的计算机语言, 同时又是一个科学计算平台. 它为数据分析和数据可视化、算法和应用程序开发提供了最核心的数学和高级图形工具. 根据它提供的500多个数学和工程函数, 工程技术人员和科学工作者可以在它的集成环境中交互或编程以完成各自的计算. MATLA-B一般用于线性代数、概率统计、图像处理、样条分析、信号处理、小波分析、振动理论、神经网络、自动控制、系统识别、算法优化和财政金融等各个方面.不过, MATLAB作为一种新的计算机语言, 要想运用自如, 充分发挥它的威力, 也需要系统的学习. 但由于使用MATLAB编程运算与人进行科学计算的思路和表达方式完全一致, 所以不像学习其他高级语言如Basic、Fortan和C语言等那样难于掌握. 下面的内容均是基于MATLAB7.5版本.1、MATLAB的主要功能(1) 数值计算功能(Numeric)(2) 符号计算功能(Symblic)(3) 图形和可视化功能(Graphic)(4) MATLAB的活笔记本功能(Notebook)(5) 可视化建模和仿真功能(Simulink)2、MATLAB的工作环境MATLAB的工作环境主要包括:·【Command Window】命令窗口;·【File Editor】文本编辑窗口;·【Figure Window】图形窗口.图0-1MATLAB 6.x的命令窗、文本编辑窗、图形窗、菜单栏和工具栏MATLAB 7.5还包含几个辅助视窗, 组成其“桌面系统”. 它们分别为:·【Workspace】工作台窗口;·【Command History】指令历史纪录窗口;·【Current Directory】当前目录选择窗口.图0-2MATLAB 7.5的桌面系统和命令窗口3、MATLAB的工作原理(1) 语言结构:MATLAB语言= 窗口命令+ M文件(2) 窗口命令:在MATLAB命令窗口中输入的MATLAB语句, 并直接执行它们完成相应的运算、绘图等.(3) M文件:在MATLAB文本编辑窗口中用MATLAB语句编写的磁盘文件, 扩展名为“.M”.二、MATLAB入门1、数学运算符及特殊字符数组的算术运算符: + - .* ./ .\ .^矩阵的算术运算符: + - * / \ ^关系运算符: < <= > >= = = ∽=逻辑运算符: & 与; | 或; ~ 非三种运算的顺序依次为: 算术运算、关系运算、逻辑运算.pi 数学常数, 即3.1415926535897....2 eps 系统的浮点(Floating-ponit) 精确度. 在PC机上, 它等于52Inf 正无穷大, 定义为1 0ans 计算结果的默认变量名NaN 不定值, 由Inf/Inf或0/0等运算产生2、基本库函数(1) 常用三角函数:sin, cos, tan, cot, sec, csc, asin, acos, atan, acot, asec, acsc等(2) 常用基本函数:sqrt(x)—开平方abs(x)—取绝对值exp(x)—以e为底的指数log(x)—自然对数log10(x)—以10为底的对数log2(x)—以2为底的对数sum(x)—求和prod(x)－求积max(x)—最大值min(x)—最小值fix(x)—对称取整sign(x)—符号函数length(x)—矩阵行数与列数中的最大值size(x)—矩阵的行数与列数注意: (1) 由于MATLAB是基于矩阵的运算,所以上面的x均表示矩阵, 数可看作是1×1的矩阵.(2) 对非向量型矩阵, 如不作特殊说明, 都是列优先.3、命令行的编写随时输入指令并按回车键, 即时给出结果;在指令最后不用任何符号并按回车键, 将显示最后结果;在指令最后用“; ”并按回车键, 将只计算但不显示最后结果.同时输入几条指令时, 用“, ”或“; ”隔开.【例0-1】数学运算符、特殊字符与基本库函数的应用>>3*(-5), 2/5, [1 2 3].*[2 4 5], [1 2 3]./[2 4 5], [2,4,5].^2ans = -15ans = 0.4000ans = 2 8 15ans = 0.5000 0.5000 0.6000ans = 4 16 25>> sin(pi/4), log(exp(1))ans = 0.7071ans = 14、变量与表达式在MATLAB中, 把由下标表示次序的标量的集合称为矩阵或数组. MATLAB是基于矩阵运算的, 因此其基本数据结构只有一个: 矩阵. 一个数也是矩阵, 只不过它是1行×1列的矩阵. MATLAB中的变量可用来存放数据, 也可用来存放向量或矩阵, 并进行各种运算.变量命名的规则为:·变量名、函数名是要区分大小写字母的;·第一个字符必须是英文字母;·字符间不可留空格;·最多只能有31个字符(只能有英文字母、数字和下连字符) .表达式由变量名、运算符和函数名等组成. 如x/sin(x), 其中x为变量名, /为运算符, sin为函数名.MATLAB语句有两种最常见形式: 1) 表达式; 2) 赋值语句: 变量= 表达式.【例0-2】赋值语句的使用>> x=1; y=x/sin(x)y = 1.1884>> x=[pi/6,pi/4,pi/3,pi/2]; sin(x)ans =0.5000 0.7071 0.8660 1.0000>> x=0:0.1:2*pi; y=sin(x); plot(x,y)图0-3y=sin(x)的曲线图5、M文件的建立、编写、保存与调用(1) 进入文本编辑窗口的方式: 在菜单栏“File”下直接点击“新建…”进入文本编辑窗口.(2) M文件的分类与格式:①命令文件: 由一系列MATLAB语句组成, 运行时将自动执行一系列命令直至给出最后结果, 而不交互地等待键盘输入. 命令文件定义的变量为全局变量, 存放于内存.②函数文件: 第一行必须包含“function”, 主要功能是建立一个函数. 函数文件定义的变量为局部变量.function 因变量名= 函数名(自变量名)注意: 函数文件要求函数名和文件名相同, 且函数名、文件名与变量名的命名规则一样.(3) 退出文本编辑窗口: 录入完毕, 存盘退出文本编辑窗口则可.【例0-3】已知1232,3,1x x x =-==, 而2112112122233,y z z z x y z z z x x =+⎧=⎧⎨⎨=-=+⎩⎩, 试求12,y y 的值.·在文本编辑窗口中编写命令文件f0_3.m:x1=-2;x2=3;x3=1; z1=3*x1^2; z2=x2+x3; y1=z1+z2 y2=z1-z2·在命令窗口中运行命令文件f0_3.m:>> f0_3y1 = 16 y2 = 8【例0-4】求()f x =分别在0,5,10x x x ===处的函数值. ·在文本编辑窗口中编写函数文件f0_4.m:function y=f0_4(x)y=log10(sqrt((x-5).^2+(x-100).^2)); ·在命令窗口中调用函数文件f0_4.m:>>x=[0,5,10]; y=f0_4(x);y = 2.0005 1.9777 1.95496、MATLAB 的在线帮助 (1) 从菜单栏上的“help ”进入 (2) 其它命令窗口帮助clc——清除显示屏上的内容clear —— 清除内存变量和函数what —— 列出当前目录下的M 、MAT 、MEX 文件 who——列出当前工作空间 (Workspace) 的变量名7、路径的设置在保存M 文件时, MATLAB 的默认位置是C:\MATLAB6p5\work. 如果用户将编写的M 文件保存在E:\experiment 目录下, 则从MATLAB 窗口的“File ”菜单中单击子菜单“Save As …”, 选择E:\experiment, 再输入本M 文件的文件名, 按“保存”键返回则可.第一章 数理统计的基本概念一、直方图与经验分布函数图的绘制hist(A,n) ——对矩阵A 按列作统计频数直方图, n 为条形图的条数ni=hist(A,n)—— 对矩阵A 按列得各划分区间内的统计频数注意: 当A 为向量时, 上述所有命令直接作用在向量上, 而不是列优先.[Fn,x0]=ecdf(x) —— 得到样本x 的经验分布函数值Fn, 当x 中有m 个不同的数 (记为向量x0) 时, 则Fn 的个数为m+1个ecdfhist(Fn,x0, m) —— 绘制数据x 的频率(密度)直方图, 其中Fn 与x0是由ecdf 函数得到的样本x 的经验分布函数值Fn 与分段点x0, m 为条形的个数, m 的默认值为10cdfplot(x) —— 绘制样本x 的经验分布函数图例如:>> x = [6 4 5 3 6 8 6 7 3 4]; >> [Fn,x0]=ecdf(x)Fn = 0 0.2000 0.4000 0.5000 0.8000 0.9000 1.0000 x0 = 3 3 4 5 67 8>> cdfplot(x)图1-1 经验分布函数图【例1】在齿轮加工中, 齿轮的径向综合误差i F ''∆是个随机变量, 今对200件同样的齿轮进行测量, 测得i F ''∆的数值 (mm) 如下, 求作i F ''∆的频率密度直方图, 并作出i F ''∆的经验分布函数图形.16 25 19 20 25 33 24 23 20 24 25 17 15 21 22 26 15 23 22 24 20 14 16 11 14 28 18 13 27 31 25 24 16 19 23 26 17 14 30 2118 16 18 19 20 22 19 22 18 26 26 13 21 13 11 19 23 18 24 2813 11 25 15 17 18 22 16 13 12 13 11 09 15 18 21 15 12 17 1314 12 16 10 08 23 18 11 16 28 13 21 22 12 08 15 21 18 16 1619 28 19 12 14 19 28 28 28 13 21 28 19 11 15 18 24 18 16 2819 15 13 22 14 16 24 20 28 18 18 28 14 13 28 29 24 28 14 1818 18 08 21 16 24 32 16 28 19 15 18 18 10 12 16 26 18 19 3308 11 18 27 23 11 22 22 13 28 14 22 18 26 18 16 32 27 25 2417 17 28 33 16 20 28 32 19 23 18 28 15 24 28 29 16 17 19 18 ·编写命令文件example1_6.m:F=[16 25 19 20 25 33 24 23 20 24 25 17 15 21 22 26 15 23 22 24....20 14 16 11 14 28 18 13 27 31 25 24 16 19 23 26 17 14 30 21....18 16 18 19 20 22 19 22 18 26 26 13 21 13 11 19 23 18 24 28....13 11 25 15 17 18 22 16 13 12 13 11 09 15 18 21 15 12 17 13....14 12 16 10 08 23 18 11 16 28 13 21 22 12 08 15 21 18 16 16....19 28 19 12 14 19 28 28 28 13 21 28 19 11 15 18 24 18 16 28....19 15 13 22 14 16 24 20 28 18 18 28 14 13 28 29 24 28 14 18....18 18 08 21 16 24 32 16 28 19 15 18 18 10 12 16 26 18 19 33....08 11 18 27 23 11 22 22 13 28 14 22 18 26 18 16 32 27 25 24....17 17 28 33 16 20 28 32 19 23 18 28 15 24 28 29 16 17 19 18];%(1)下面作频数直方图figure(1)hist(F,8)title('频数直方图');xlabel('齿轮的径向综合误差(mm)');%(2)下面作频率(密度)直方图[Fn,x0]=ecdf(F);figure(2)ecdfhist(Fn,x0,8);title('频率(密度)直方图');xlabel('齿轮的径向综合误差(mm)');%(3)下面作经验分布函数图figure(3)cdfplot(F)title('经验分布函数图');xlabel('齿轮的径向综合误差(mm)');·运行命令文件example1_6.m:>> example1_6图1-2二、常见的概率分布表1-1 常用概率分布及代码频数直方图齿轮的径向综合误差(mm)频率(密度)直方图齿轮的径向综合误差(mm)齿轮的径向综合误差(mm)F (x )经验分布函数图三、MATLAB 为常见分布提供的五类函数1) 概率密度函数(分布名+pdf) 2) (累积)分布函数(分布名+cdf) 3) 逆(累积)分布函数(分布名+inv) 4) 随机数发生器(分布名+rnd) 5) 均值和方差(分布名+stat) 1、概率密度函数表1-2 概率密度函数(pdf)注意: Y=normpdf (X, mu, sigma)的sigma 是指标准差σ, 而非2σ. 【例2】 绘制标准正态分布(0,1)N 的概率密度图.x=-4:0.1:4; y=normpdf(x,0,1); plot(x,y)title('N(0,1)的概率密度曲线图')图1-3 标准正态分布的概率密度图2、累积分布函数表1-3 累积分布函数(cdf)【例3】求服从标准正态分布的随机变量落在区间[-2, 2]上的概率.>> P=normcdf (2,0,1)-normcdf(-2,0,1)ans = 0.95453、逆累积分布函数 (用于求分位点)表1-4 逆累积分布函数(inv)【例4】 求下列分位数: (i) 0.9u ;(ii) 0.25(4)t ;(iii) 0.1(14,10)F ;(iv) 20.025(50)χ.>> u_alpha=norminv(0.9,0,1)u_alpha = 1.2816 >> t_alpha=tinv(0.25,4)t_alpha = -0.7407 >> F_alpha=finv(0.1,14,10)F_alpha = 0.4772 >> X2_alpha=chi2inv(0.025,50)X2_alpha = 32.35744、随机数发生函数表1-5 随机数发生函数(rnd)5、均值和方差表1-6 常见分布的均值和方差函数(stat)注意: (1) MATLAB 中的指数分布的概率密度函数是1,0()0,0xue xf x u x -⎧>⎪=⎨⎪≤⎩.(2) 如果省略调用格式左边的[M, V], 则只计算出均值.四、常用的统计量表1-7 常用统计量说明:(1) y=var(X) ——计算X 中数据的方差, 其中211var()()1ni i X x x n ==--∑. y=var(X, 1) ——211var(,1)()n i i X x x n ==-∑, 得到样本的二阶中心矩 (转动惯量).(2) C ＝cov(X) ——返回一个协方差矩阵, 其中输入矩阵X 的每列元素代表着一个随机变量的观测值. 如果X 为n ×m 的矩阵, 则C 为m ×m 的矩阵.(3) var(X)=diag(cov(X)), std(X)=sqrt(diag(cov(X))).。

Matlab在《概率论与数理统计》教学中的应用1. 引言1.1 研究背景概率论与数理统计作为现代科学研究的基础，广泛应用于物理、生物、经济、工程等各个领域。

在教学中，传统的概率论与数理统计教学往往通过纸笔计算和手工绘图进行，这样的方式在一定程度上限制了学生对概念的理解和实际应用能力的培养。

而引入Matlab这样的数学计算软件，可以极大地提高教学效率，使学生更直观地理解抽象的数学概念，提高他们的学习兴趣和动手能力。

通过将Matlab与概率论与数理统计相结合，可以更好地展示概率分布、统计分析、随机模拟等概念，加深学生对这些内容的理解和掌握。

研究Matlab在概率论与数理统计教学中的应用具有重要意义。

本文将探讨Matlab在概率论与数理统计教学中的具体应用，分析其在教学中的优势和未来发展方向。

1.2 研究意义概率论与数理统计作为数学学科中重要的分支，旨在研究事件的发生规律以及数据的分布特征，对现代科学、技术和社会管理等领域都具有重要的应用价值。

在教学中，采用Matlab作为工具可以加深学生对概率与统计理论的理解，提高他们的计算和分析能力，培养他们解决实际问题的能力。

通过引入Matlab，学生可以更加直观地掌握数学模型的建立和计算方法，提高他们对概率与统计学习的兴趣和积极性，进一步激发他们学习的潜力。

Matlab在教学中的应用也有助于培养学生的动手能力和实际解决问题的能力，提高他们的实践能力和创新思维。

教师可以结合具体案例，引导学生运用Matlab工具分析问题，并进行模拟实验和数据处理，使学生在实践中不断探索、思考和总结，从而提高他们的学习效果和实际应用能力。

Matlab在概率论与数理统计教学中的应用具有重要的意义和价值。

2. 正文2.1 Matlab在概率论教学中的基本概念应用Matlab可以用来计算概率。

通过编写简单的代码，可以计算各种随机事件发生的概率，例如掷硬币、抛骰子等。

这样的实践可以帮助学生深入理解概率的概念，同时提高他们的计算能力。

Matlab概率论与数理统计一、 matlab 基本操作1.画图【例】简单画图hold off;x=0:0.1:2*pi;y=sin(x);plot(x,y,'-r');x1=0:0.1:pi/2;y1=sin(x1);hold on;fill([x1, pi/2],[y1,1/2],'b' );【例】填充，二维平均随机数hold off;x=[0,60];y0=[0,0];y60=[60,60];x1=[0,30];y1=x1+30;x2=[30,60];y2=x2-30;xv=[0 0 30 60 60 30 0];yv=[0 30 60 60 30 0 0];fill(xv,yv,'b');hold on ;plot(x,y0,'r',y0,x,'r',x,y60,'r' ,y60,x,'r' );plot(x1,y1,'r',x2,y2,'r');yr=unifrnd (0,60,2,100);plot(yr(1,:),yr(2,:),'m.')axis('on');axis('square');2.排列组合C=nchoosek(n,k) ：C C n k，例 nchoosek(5,2)=10, nchoosek(6,3)=20.prod(n1:n2) ：从 n1 到 n2 的连乘【例】最少有两个人寿辰相同的概率n!C N nN !( N n)!N(N1)(N n1)公式计算 p 111N nN n N n365 364 (365rs1)365364365rs 1 1365rs1365365365rs=[20,25,30,35,40,45,50];%每班的人数p1=ones(1,length(rs));p2=ones(1,length(rs));%用连乘公式计算for i=1:length(rs)p1(i)=prod(365-rs(i)+1:365)/365^rs(i);end%用公式计算（改进）for i=1:length(rs)for k=365-rs(i)+1:365p2(i)=p2(i)*(k/365);end ;endp1(i)=exp(sum(log(365-rs(i)+1:365))-rs(i)*log(365));endp_r1=1-p1;p_r2=1-p2;Rs =[20253035404550 ]P_r=[0.4114 0.5687 0.7063 0.8144 0.8912 0.9410 0.9704]二、随机数的生成3.平均分布随机数rand(m,n); 产生 m 行 n 列的 (0,1) 平均分布的随机数rand(n); 产生 n 行 n 列的 (0,1)平均分布的随机数【练习】生成(a,b)上的平均分布4.正态分布随机数randn(m,n); 产生 m 行 n 列的标准正态分布的随机数【练习】生成N(nu,sigma.^2) 上的正态分布5.其他分布随机数函数名调用形式注释Unidrnd unid rnd (N,m,n)平均分布（失散）随机数binornd bino rnd (N,P,m,n)参数为 N, p的二项分布随机数Poissrnd poiss rnd (Lambda,m,n)参数为 Lambda的泊松分布随机数geornd geornd (P,m,n)参数为 p 的几何分布随机数hygernd hygernd (M,K,N,m,n)参数为 M， K， N 的超几何分布随机数Normrnd normrnd (MU,SIGMA,m,n)参数为 MU， SIGMA的正态分布随机数，SIGMA是标准差Unifrnd unif rnd ( A,B,m,n)[A,B] 上平均分布 ( 连续 ) 随机数Exprnd exprnd (MU,m,n)参数为 MU的指数分布随机数chi2rnd chi2 rnd(N,m,n)自由度为 N 的卡方分布随机数Trnd t rnd(N,m,n)自由度为 N 的 t分布随机数Frnd f rnd(N1, N2,m,n)第一自由度为N1, 第二自由度为 N2 的 F 分布随机数gamrnd gamrnd(A, B,m,n)参数为 A, B的分布随机数betarnd betarnd(A, B,m,n)参数为 A, B的分布随机数lognrnd lognrnd(MU, SIGMA,m,n)参数为 MU, SIGMA的对数正态分布随机数nbinrnd nbinrnd(R, P,m,n)参数为 R，P 的负二项式分布随机数ncfrnd ncfrnd(N1, N2, delta,m,n)参数为 N1， N2， delta 的非中心 F 分布随机数nctrnd nctrnd(N, delta,m,n)参数为 N，delta的非中心 t 分布随机数ncx2rnd ncx2rnd(N, delta,m,n)参数为 N，delta的非中心卡方分布随机数raylrnd raylrnd(B,m,n)参数为 B 的瑞利分布随机数weibrnd weibrnd(A, B,m,n)参数为 A, B的韦伯分布随机数三、一维随机变量的概率分布1.失散型随机变量的分布率(1)0-1 分布(2)平均分布(3) 二项分布： binopdf(x,n,p) ，若X ~ B(n, p)，则P{ X k} C n k p k (1p) n k，x=0:9;n=9;p=0.3;y= binopdf(x,n,p);plot(x,y,'b-',x,y,'r*')y=[ 0.0404, 0.1556, 0.2668, 0.2668, 0.1715, 0.0735, 0.0210, 0.0039, 0.0004, 0.0000 ]‘当 n 较大时二项分布近似为正态分布x=0:100;n=100;p=0.3;y= binopdf(x,n,p);plot(x,y,'b-',x,y,'r*')(4) 泊松分布： piosspdf(x, lambda) ，若X ~k e ( ) ，则 P{ X k}k !x=0:9; lambda = 3;y= poisspdf (x,lambda);plot(x,y,'b-',x,y,'r*')y=[ 0.0498, 0.1494, 0.2240, 0.2240, 0.1680, 0.1008, 0.0504, 0.0216, 0.0081, 0.0027 ](5) 几何分布： geopdf (x, p)，则P{ X k} p(1p) k 1y= geopdf(x,p);plot(x,y,'b-',x,y,'r*')y=[ 0.3000, 0.2100, 0.1470, 0.1029, 0.0720, 0.0504, 0.0353, 0.0247, 0.0173, 0.0121 ]C M k C N n k Mx=0:10;N=20;M=8;n=4;y= hygepdf(x,N,M,n);plot(x,y,'b-',x,y,'r*')y=[ 0.1022, 0.3633, 0.3814, 0.1387, 0.0144, 0, 0, 0, 0, 0, 0 ]2.概率密度函数1a x b(1)平均分布： unifpdf(x,a,b) ，f ( x)b a0其他a=0;b=1;x=a:0.1:b;y= unifpdf (x,a,b);112(2)正态分布： normpdf(x,mu,sigma) ，f ( x)e2 2 ( x)2x=-10:0.1:12;mu=1;sigma=4;y= normpdf(x,mu,sigma);rn=10000;z= normrnd (mu,sigma,1,rn); %产生 10000 个正态分布的随机数d=0.5;a=-10:d:12;b=(hist(z,a)/rn)/d;% 以 a 为横轴，求出 10000 个正态分布的随机数的频率plot(x,y,'b-',a,b,'r.')(3) 指数分布： exppdf(x,mu) ，f (x)1 e1xa x by= exppdf(x,mu); plot(x,y,'b-',x,y,'r*')1n1(4)2分布： chi2pdf(x,n) ， f (x; n)2n 2x2( n 2)hold on x=0:0.1:30;n=4;y= chi2pdf(x,n);plot(x,y,'b');%blue n=6;y= chi2pdf(x,n);plot(x,y,'r');%redn=8;y= chi2pdf(x,n);plot(x,y,'c');%cyan n=10;y= chi2pdf(x,n);plot(x,y,'k');%black legend('n=4', 'n=6', 'n=8', 'n=10');(( n 1) 2) x 2(5) t 分布： tpdf(x,n) ， f (x; n)(n 2)1nnhold on x=-10:0.1:10;n=2;y= tpdf(x,n);plot(x,y,'b');%blue e 2n 1 2x 0x 0n=20;y= tpdf(x,n);plot(x,y,'k');%black legend('n=2', 'n=6', 'n=10', 'n=20');n1n1 2n1n222(6) F 分布： fpdf(x,n1,n2) ，f ( x; n1, n2)(( n1n2 ) 2) n1x 21n1x x 0 (n1 2)(n2 2) n2n20x 0hold onx=0:0.1:10;n1=2; n2=6;y= fpdf(x,n1,n2);plot(x,y,'b');%bluen1=6; n2=10;y= fpdf(x,n1,n2);plot(x,y,'r');%redn1=10; n2=6;y= fpdf(x,n1,n2);plot(x,y,'c');%cyann1=10; n2=10;y= fpdf(x,n1,n2);plot(x,y,'k');%blacklegend(' n1=2; n2=6', ' n1=6; n2=10', ' n1=10; n2=6', ' n1=10; n2=10');3.分布函数 F (x) P{ X x}【例】求正态分布的累积概率值设 X ~ N(3,22)，求P{2X 5},P{ 4 X 10},P{ X 2}, P{X3} ，4.逆分布函数，临界值y F (x) P{ X x} ， x F 1 ( y) ， x 称之为临界值【例】求标准正态分布的累积概率值y=0:0.01:1;x=norminv(y,0,1);【例】求2 (9) 分布的累积概率值hold offy=[0.025,0.975];x=chi2inv(y,9);n=9;x0=0:0.1:30;y0=chi2pdf(x0,n);plot(x0,y0,'r');x1=0:0.1:x(1);y1=chi2pdf(x1,n);x2=x(2):0.1:30;y2=chi2pdf(x2,n);hold onfill([x1, x(1)],[y1,0],'b');fill([x(2),x2],[0,y2],'b');5. 数字特色函数名调用形式注释sort sort(x),sort(A)排序 ,x 是向量， A 是矩阵，按各列排序sortrows sortrows(A) A 是矩阵，按各行排序mean mean(x)向量 x 的样本均值var var(x)向量 x 的样本方差std std(x)向量 x 的样本标准差median median(x)向量 x 的样本中位数geomean geomean(x)向量 x 的样本几何平均值harmmean harmmean(x)向量 x 的样本调停平均值skewness skewness(x)向量 x 的样本偏度max max(x)向量 x 的最大值min min(x)向量 x 的最小值cov cov(x), cov(x,y)向量 x 的方差，向量x,y 的协方差矩阵corrcoef corrcoef(x,y)向量 x,y 的相关系数矩阵【练习】二项分布、泊松分布、正态分布（ 1）对n10, p 0.2 二项分布，画出 b(n, p) 的分布律点和折线；（ 2）对np ，画出泊松分布( ) 的分布律点和折线；（ 3）对np,2np(1 p) ，画出正态分布N ( , 2 )的密度函数曲线；（ 4）调整 n, p ，观察折线与曲线的变化趋势。