概率统计与MATLAB精品PPT课件
10MATLAB符号运算与概率统计.ppt
26.03.2019
MATLAB符号运算
dsolve('D2y = -a^2*y', 'y(0) = 1', 'Dy(pi/a) = 0') ans= cos(a*t) dsolve('Dx = y', 'Dy = -x') x= cos(t)*C1+sin(t)*C2 y= -sin(t)*C1+cos(t)*C2
MATLAB符号运算与概率统计
26.03.2019
1
MATLAB符号运算
1.设置符号变量（sym, syms） S = sym(A)：将已知量声明成符号变量。 x = sym(‘x’)：声明x为符号型变量。 x = sym('x','real') x = sym('x','unreal') S = sym(A,flag) flag 可以是 ‘r’, ‘d’, ‘e’, 或 'f'.
12
26.03.2019
MATLAB符号运算
A = solve('a*u^2 + v^2', 'u - v = 1', 'a^2 - 5*a + 6') A= a: [1x4 sym] u: [1x4 sym] v: [1x4 sym] 这里A.a = [ 2, 2, 3, 3]
26.03.2019
8
MATLAB符号运算
v = [(1 + a/x)^x, exp(-x)]; limit(v,x,inf,'left') => [exp(a), 0] 2.4.级数和 Symsum(s,t,a,b) s中t从a到b求和。
10MATLAB符号运算与概率统计.ppt
26.03.2019
8
MATLAB符号运算
v = [(1 + a/x)^x, exp(-x)]; limit(v,x,inf,'left') => [exp(a), 0] 2.4.级数和 Symsum(s,t,a,b) s中t从a到b求和。
26.03.2019
9
MATLAB符号运算
int([exp(t),exp(alpha*t)]) ans= [exp(t), 1/alpha*exp(alpha*t)] 2.3.极限 limit(F,x,a) limit(F,a)（x->a时F的极限） limit(F) （x->0时F的极限) limit(F,x,a,'right') limit(F,x,a,'left')
26.03.2019
4
MATLAB符号运算
2.2积分 R = int(S) R = int(S,v) R = int(S,a,b) R = int(S,v,a,b) 例如： int(-2*x/(1+x^2)^2) ans= 1/(1+x^2)
26.03.2019 5
26.03.2019
11
MATLAB符号运算
如： solve('a*x^2 + b*x + c') ans= [ 1/2/a*(-b+(b^2-4*a*c)^(1/2)), 1/2/a*(-b-(b^2-4*a*c)^(1/2))] solve('a*x^2 + b*x + c','b') ans= -(a*x^2+c)/x solve('x + y = 1','x - 11*y = 5') y = -1/3, x = 4/3
概率统计的数值实验MATLAB在概率统计教学中的应用-PPT精选
P A iA jA k P A k P A iA jA k n n 1 1 n 2 1 i j k n
PA1A2A3Ann1!
于是
n
P Ai i1
n
1, P Ai A j
1i jn
n 2 15
k 1
模拟Galton钉板试验的步骤： (1) 确定钉子的位置：将钉子的横、纵坐标存储在两个矩阵X和
Y中。 (2) 在Galton钉板试验中，小球每碰到钉子下落时都具有两种
4/5 0.5134
4/5 0.5086
4/5 0.5093
4/5 0.5093
π的近似值 3.1116 3.1165 3.1460 3.1418 3.1418
试验次数n
5千
1万
10万 100万 1000万
针长l/平行间 距d
相交频率
17/20 0.5432
17/20 0.5452
17/20 0.5420
概率
生日各不相同的概率 至少两人生日相同的概率 1
0.8
0.6
0.4
0.2
0 0 10 20 30 40 50 60 70 80 90 100 人数
• p1(30)=0.7063, p1(60)= 0.9941
分析：在30名学生中至少两人生日相同的概率为70.63％。 下面进行计算机仿真。
随机产生30个正整数，代表一个班30名学生的生日，然后观
解 记事件 A i 为第i个人拿到自已枪，事件 A i 为第i个人 没拿到自己枪，易知：
PAi
1 n
P Ai
n1 n
i1,2, ,n
又记 p 0 为没有一个人拿到自己枪的概率。
《概率统计》学习课件
134 第4章概率统计本章介绍MATLAB在概率统计中的若干命令和使用格式，这些命令存放于MatlabR12\Toolbox\Stats中。
4.1 随机数的产生4.1.1 二项分布的随机数据的产生命令参数为N，P的二项随机数据函数binornd格式R = binornd(N,P) %N、P为二项分布的两个参数，返回服从参数为N、P的二项分布的随机数，N、P大小相同。
R = binornd(N,P,m) %m指定随机数的个数，与R同维数。
R = binornd(N,P,m,n) %m,n分别表示R的行数和列数例4-1>> R=binornd(10,0.5)R =3>> R=binornd(10,0.5,1,6)R =8 1 3 7 6 4>> R=binornd(10,0.5,[1,10])R =6 8 4 67 5 3 5 6 2>> R=binornd(10,0.5,[2,3])R =7 5 86 5 6>>n = 10:10:60;>>r1 = binornd(n,1./n)r1 =2 1 0 1 1 2>>r2 = binornd(n,1./n,[1 6])r2 =0 1 2 1 3 14.1.2 正态分布的随机数据的产生命令参数为μ、σ的正态分布的随机数据函数normrnd格式R = normrnd(MU,SIGMA) %返回均值为MU，标准差为SIGMA的正态分布的随机数据，R可以是向量或矩阵。
R = normrnd(MU,SIGMA,m) %m指定随机数的个数，与R同维数。
R = normrnd(MU,SIGMA,m,n) %m,n分别表示R的行数和列数例4-2>>n1 = normrnd(1:6,1./(1:6))n1 =2.1650 2.31343.02504.0879 4.8607 6.2827>>n2 = normrnd(0,1,[1 5])n2 =0.0591 1.7971 0.2641 0.8717 -1.4462>>n3 = normrnd([1 2 3;4 5 6],0.1,2,3) %mu为均值矩阵n3 =0.9299 1.9361 2.96404.12465.0577 5.9864>> R=normrnd(10,0.5,[2,3]) %mu为10，sigma为0.5的2行3列个正态随机数R =9.7837 10.0627 9.42689.1672 10.1438 10.59554.1.3 常见分布的随机数产生常见分布的随机数的使用格式与上面相同表4-1 随机数产生函数表函数名调用形式注释Unifrnd unifrnd ( A,B,m,n) [A,B]上均匀分布(连续) 随机数Unidrnd unidrnd(N,m,n) 均匀分布（离散）随机数Exprnd exprnd(Lambda,m,n) 参数为Lambda的指数分布随机数Normrnd normrnd(MU,SIGMA,m,n) 参数为MU，SIGMA的正态分布随机数chi2rnd chi2rnd(N,m,n) 自由度为N的卡方分布随机数Trnd trnd(N,m,n) 自由度为N的t分布随机数Frnd frnd(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的韦伯分布随机数binornd binornd(N,P,m,n) 参数为N, p的二项分布随机数geornd geornd(P,m,n) 参数为p的几何分布随机数hygernd hygernd(M,K,N,m,n) 参数为M，K，N的超几何分布随机数Poissrnd poissrnd(Lambda,m,n) 参数为Lambda的泊松分布随机数4.1.4 通用函数求各分布的随机数据命令求指定分布的随机数135函数random格式y = random('name',A1,A2,A3,m,n) %name的取值见表4-2；A1，A2，A3为分布的参数；m，n指定随机数的行和列例4-3 产生12（3行4列）个均值为2，标准差为0.3的正态分布随机数>> y=random('norm',2,0.3,3,4)y =2.3567 2.0524 1.8235 2.03421.9887 1.94402.6550 2.32002.0982 2.2177 1.9591 2.01784.2 随机变量的概率密度计算4.2.1 通用函数计算概率密度函数值命令通用函数计算概率密度函数值函数pdf格式Y=pdf(name，K，A)Y=pdf(name，K，A，B)Y=pdf(name，K，A，B，C)说明返回在X=K处、参数为A、B、C的概率密度值，对于不同的分布，参数个数是不同；name为分布函数名，其取值如表4-2。
matlab概率统计第2讲ppt课件
A(i:k , j:l) 返回由二维矩阵A中的第i到k行行向量和
第j到l列列向量组成的子阵。
17
e. 特殊矩阵
命令 hilb(n) toeplitz(k,r) pascal(n) rosser gallery(n) wilkinson(n) magic(n)
……
运行结果 生成n×n的希尔波特矩阵 生成非对称的托普利兹矩阵 pascal矩阵(帕斯卡矩阵) rosser矩阵 数字分析中有名的n×n试验矩阵 返回wilkinson特征值测试矩阵 魔方矩阵
阵，B为n×p矩阵，则C=A*B为m×p矩阵。
39
(3) 矩阵除法 在MATLAB中,有两种矩阵除法运算：\和/,分别表
示左除和右除.如果A矩阵是非奇异方阵,则A\B和B/A 运算可以实现.A\B等效于A的逆左乘B矩阵,也就是 inv(A)*B,而B/A等效于A矩阵的逆右乘B矩阵,也就是 B*inv(A).
1 4
7
%A的2,3行 %A的第一列
33
键入： A=[1 2 3;4 5 6;7 8 9]; C=A(1:2, [1 3]) %A的第1，2行， %A的第1，3列
输出：C= 13 46
还有A(1:2:3, 3:-1:1)
34
5) 矩阵的拼接
将几个矩阵接在一起称为拼接，左右拼 接行数要相同，上下拼接列数要相同。
15
c. 三角矩阵
命令
triu(A) triu(A,k) tril(A) tril(A,k)
运行结果
生成一个和A维数相同的上三角矩阵。该矩 阵主对角线及以上元素取自A中相应元素。 其余元素为0。
生成一个和A维数相同的上三角矩阵。该矩 阵第k条对角线及以上元素取自A中相应元 素。其余元素为0。
matlab课件--第5讲 概率统计实验
2. exprnd函数
例7: 产生4行5列的指数分布的随机数. 程序如下： y=exprnd(3,4,5) %参数=3
Matlab 软件实习
三、随机变量与概率分布密度
1. 几个常用的离散型分布密度函数(…pdf )
(1)均匀分布 P(X=xn)=1/n 密度函数调用格式：y=unidpdf(X,N) 例8：求X取值为1,2,3,4,5,6,7,8时服从均匀分布的概率值. 程序如下: X=1:8,N=8; Y=unidpdf(X,N)
% 参数SIGMA为正数
Matlab 软件实习
四、随机变量与概率分布函数
累积分布函数(…cdf )—在工具箱中分布函数亦称累积分布函 数,即表示事件的概率P{Xx}。
累积分布函数表 分布类型名称 函数名称 函数调用格式
离散均匀分布
二项分布 泊松分布 几何分布
unidcdf
binocdf poisscdf geocdf
程序如下: Y=[1500 2000 2500 3000]; P=[0.0952 0.0861 0.0779 0.7408]; EX=Y*P’
Matlab 软件实习 (2) 连续型 EX=int(x*f(x),-inf,inf)
例2：
1 , f (v ) a 0,
0va 其它
EV=int(v*1/a,0,a)
DV=int(v^2*1/a,0,a)-EV^2
Matlab 软件实习
3. 常见分布的期望与方差函数
分布类型名称 离散均匀分布 二项分布 几何分布 超几何分布 函数名称 unidstat binostat geostat hygestat 函数调用格式 [E,D]=unidstat(N) [E,D]=binostat(N,P) [E,D]=geostat(P) [E,D]=hygestat(M,K,N)
基于MATLAB的概率统计数值实验ppt课件
快捷的学习可借助MATLAB的系统帮助，通过指令doc 获得具体函数的详细信息，语法是 doc <函数名>
5/60
2. 二项分布实验
已知Y~b(20, 0.3)求Y分布率的值，并划出图形
在Matlab中输入以下命令：
binopdf(10,20,0.2) x=0:1:20; y=binopdf(x,20,0.2) plot(x,y,’r.’)
例9 某种重大疾病的医疗险种，每份每年需交保险费100元，若在 这一年中，投保人得了这种疾病，则每份可以得到索赔额10000元， 假设该地区这种疾病的患病率为0.0002，现该险种共有10000份保 单，问： (1)保险公司亏本的概率是多少? (2)保险公司获利不少于80万元的概率是多少?
15/60
(1) 每小时恰有4次呼叫的概率
(2) 一小时内呼叫不超过5次的概率 (3) 画出分布律图像
(1)
( 2)
P ( X 4)
4
4!
5
e
3k 3 P ( X 5) P ( X k ) e k 0 k 0 k!
5
34 3 e 4!
在Matlab中输入以下命令： (1)p1= poisspdf(4,3) (2)p2= poisscdf(5,3) (3)x=0:1:20;y=poisspdf(x,3);plot(x,y)
(2) σ=0.5, μ=1,2,3,4
(1)命令： x=-6:0.1:6; y1=normpdf(x,3,0.5); y2=normpdf(x,3,0.7); y3=normpdf(x,3,1); y4=normpdf(x,3,1.5); y5=normpdf(x,3,2); plot(x,y1,'.',x,y2,'+',x,y3,'*',x,y4,'d',x,y5)
matlab概率统计ppt课件
例7: 参数方程绘圆心在原点,半径为2的圆.
程序如下: t=linspace(-2,2,30)*pi; x=2*cos(t);y=2*sin(t); plot(x,y,’r’); text(-0.25,0,’x^2+y^2=4’); % 在图中标记曲线方程
2. 一元符号表达式函数绘图法
① 首先定义x是符号变量,再定义y是x的符号表达式函数; ② 用绘图命令ezplot或fplot绘图.
程序如下: t= -pi:pi/200:pi; comet(t,tan(sin(t))-sin(tan(t)))
(2) 三维动态轨线图
程序如下: t= 0:0.05:100; x=t;y=sin(t);z=sin(2*t); comet3(x,y,z)
4. 统计图形绘制
(1) 条形图
bar(x,y) ％竖直条形图,其中x是横坐标向量,y是向量或矩阵 barh(x,y) ％水平条形图 bar3(x,y) ％三维条形图 bar3h(x,y) ％三维水平条形图
例6 :在同一坐标内,分别用不同线型和颜色绘制曲线 y1=0.2e-0.5xcos(4x) 和y2=2e-0.5xcos(x),标记两曲线 交叉点.
程序如下：
x=linspace(0,2*pi,1000);
y1=0.2*exp(-0.5*x).*cos(4*pi*x);
y2=2*exp(-0.5*x).*cos(pi*x);
② ezplot(f,[a,b])表示在给定区间a<x<b上绘制函数f(x)的图形.
③ ezplot(f(x,y))表示在默认区域－2*pi<x<2*pi, －2*pi<y<2*pi上 绘制隐函数f(x,y)=0的函数图形.
MATLAB教程PPT课件
6
五、Matlab程序设计
（一）选择结构
选择结构的语句有if语句和switch语句。
1． if语句
格式一： if 条件 语句组 end 格式二： if 条件 语句组1 else 语句组2 end
格式三： if 条件1 语句组1 elseif 条件2 语句组2 …… elseif 条件m 语句组m else 语句组m+1 end
4
三、Matlab 的输出
若Matlab语句结束后加“;”，则不会输出该语句的 执行结果，否则，执行结果将会输出。
Matlab 以双精度执行所有的运算，运算结果可以在 屏幕上输出，同时赋给指定变量；若无指定变量，则系 统会自动将结果赋给变量 “ans”
，a+b 减法，a-b 乘法，a×b 除法，a÷b 乘幂，a*a*…*a 符号 + * / or \ ^ 举例 5+3 5-3 5*3 48/4=4\48=12 5^2=25
function bool = isdiff( vec ) %检验一个向量中是否有相同的元素 %@param：vec 目标向量 [row,col] = size( vec ); new_vec=unique(vec); %去除vec中相同元素 [row_new,col_new] = size( new_vec ); if(col == col_new) bool = 0; return; end bool = 1; end
MATLAB
《概率与统计》 MATLAB上机指导
1
Matlab 的工作界面
当前工 作目录
命令 提示符
当前工 作空间 命令窗口
输入命令的 历史记录
2
一、变量及矩阵输入
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
功能：产生M lambda)
功能：计算分布密度p(x)在x的值
21.10.2020
x0 x0
7
§1 随机变量及其分布
均匀分布X~U(a,b) 命令1：Fx=unifcdf(x, a,b) 功能：计算累积概率Fx=P{X≤x}=F(x) 命令2：x=unifinv(p, a,b) 功能：计算随机量x，使得p=P{X≤x} 命令3：X=unifrnd(a,b,M,N) 功能：产生M*N维随机数矩阵X 命令4：Px=unifpdf(x, a,b) 功能：计算分布密度p(x)在x的值 补充：rand()---(0,1)均匀分布随机数
21.10.2020
12
§1 随机变量及其分布
例1某人向空中抛硬币100次，落下为正面的概率 为0.5。这100次中正面向上的次数记为X： (1)试计算x=45的概率和x≤45的概率； (2)绘制分布函数图象和分布列图象。
程序：》clear;
px=binopdf(45,100,0.5) % 计算x=45的概率
命令2：x=hygeinv(p,M, N,K)
功能：在已知参数M、N 、 K和p的情况下计算随 机量x，使得p=P{0≤次品数X≤x}
命令3：X=hygernd(M,N,K,m,n)
功能：在已知参数M,N ,K的情况下产生m*n维符合
超几何分布的随机数矩阵X
21.10.2020
2
§1 随机变量及其分布
21.10.2020
6
§1 随机变量及其分布
指数分布X~exp(λ)
1ex
P{Xx}
0
命令1：Fx=expcdf(x, lambda)
功能：计算累积概率Fx=P{X≤x}=F(x)
命令2：x=expinv(p, lambda)
功能：计算随机量x，使得p=P{X≤x}
命令3：X=exprnd(lambda,M,N)
21.10.2020
9
§1 随机变量及其分布
2分布
密度函f数 2(x)： 2k21(k)xk21e2x
2
0
x0 x0
命令:chi2cdf(x, k), chi2inv(p, k),chi2pdf(x, k) chi2rnd(k,m,n)
21.10.2020
10
§1 随机变量及其分布
T分布
密度函fT(数 x)： k(k 2(1k))(1xk2)k21
命令1：Fx=binocdf(x,n,p) 功能：计算二项分布的累积概率Fx=P{X≤x}=F(x) 命令2：x=binoinv(y, n,p) 功能：计算随机量x，使得y=P{X≤x} 命令3：X=binornd(n,p,M,N) 功能：产生M*N维符合二项分布的随机数矩阵X 命令4：Px=binopdf(x,n, p) 功能：计算试验中事件恰好发生x次的概率Px=P{X=x}
命令2：x=norminv(p, mu,sigma)
功能：计算随机量x，使得p=P{X≤x}
命令3：X=normrnd(mu,sigma,M,N)
功能：产生M*N维随机数矩阵X
命令4：Px=normpdf(x, mu,sigma)
功能：计算分布密度p(x)在x的值
补充：randn()---标准正态分布随机数
2
命令:tcdf(x, k), tinv(p, k),tpdf(x, k) trnd(k,m,n)
21.10.2020
11
§1 随机变量及其分布
F分布
密度函 fF(x)数 ： ((2pp) 2(qq2))p2pqq20x2p1(pq)xp2q
x0 x0
命令:fcdf(x, p,q), finv(F,p,q),fpdf(x, p,q) frnd(p,q,m,n)
21.10.2020
8
§1 随机变量及其分布
Γ分布
密度函数：f (x)
a
x
a1
(a)
e x
0
x0 x0
其中(a) xa1ex dx 0
满足：(a) a(a 1),(1) 1, (1)
2
命令:gamcdf(x, a, lambda), gaminv(p, a, lambda) gampdf(x, a,lambda), gamrnd(a, lambda,m,n)
功能：产生M*N维随机数矩阵X
命令4：Px=poisspdf(x,lambda)
功能：计算概率Px=P{X=x}
21.10.2020
5
§1 随机变量及其分布
正态分布X~N(μ,σ2)
x
P{Xx}
1
1 (x)2
e22 dx
命令1：Fx=normcdf(x,
2
mu,sigma)
功能：计算累积概率Fx=P{X≤x}=F(x)
21.10.2020
4
§1 随机变量及其分布
泊松分布X~P(λ)
P{Xk}e k
k!
命令1：Fx=poisscdf(x,lambda)
功能：计算累积概率Fx=P{X≤x}=F(x)
命令2：x=poissinv(p, lambda)
功能：计算随机量x，使得p=P{X≤x}
命令3：X=poissrnd(lambda,M,N)
命令4：Px=hygepdf(x,M, N, K)
功能：总共M件产品，其中次品N 件，抽取K件检查， 计算发现恰好x件次品的概率Px=P{X=x}
命令5：stairs(x,Px)
功能：绘制以 x为横坐标，Px为纵坐标的阶梯平面图； 当Px是分布列(或密度)时，绘制概率密度分布图； 当Px是累积分布时，绘制概率分布函数图
注：以后碰到命令末尾为：
rnd----产生随机数X； cdf----产生分布函数F(x)
pdf----产生密度函数p(x)或分布列Px=P{X=x}
inv----计算x=F-1(p)→ p=F (x)
21.10.2020
3
§1 随机变量及其分布
二项分布B(n,p) P { X k } C k np k(1 p )n k
概率统计与MATLAB
随机变量及其分布 随机变量函数的分布 随机变量的数字特征 参数估计 假设检验
回归分析
21.10.2020
1
§1 随机变量及其分布
超几何分布H(n,M,N)
P{X
命令1：Fx=hygecdf(x,M,N,K)
i}
CiNCKMiN CKM
功能：计算超几何分布的累积概率，总共M件产品， 其中次品N 件，抽取K件检查，计算发现次品不 多于x件的概率Fx=P{次品数X≤x}=F(x)