MATLAB函数和命令的用法.docx

MATLAB函数和命令的用法.docx

《MATLAB函数和命令的用法.docx》由会员分享，可在线阅读，更多相关《MATLAB函数和命令的用法.docx（21页珍藏版）》请在冰豆网上搜索。

MATLAB函数和命令的用法

Binocdf二项式累积分布函数

语法格式

Y=binocdf（X,N,P）

函数功能

Y=binocdf（X,N,P）

计算X中每个X（i）的二项式累积分布函数，其中，N中对应的N（i）为试验数，P中对应的P（i）为每次试验成功的概率。

Y,N,和P的大小类型相同，可以是向量、矩阵或多维数组。

输入的标量将扩展成一个数组，使其大小类型与其它输入相一致。

ThevaluesinNmustallbepositiveintegers,thevaluesinXmustlieontheinterval[0,N],andthevaluesinPmustlieontheinterval[0,1].

Thebinomialcdfforagivenvaluexandagivenpairofparametersnandpis

Theresult,y,istheprobabilityofobservinguptoxsuccessesinnindependenttrials,wheretheprobabilityofsuccessinanygiventrialisp.TheindicatorfunctionI（0,1,...,n）（i）ensuresthatxonlyadoptsvaluesof0,1,...,n.

示例

若一个棒球队在一个赛季要比赛162场，每场比赛取胜的机会是50-50，则该队取胜超过100场的概率为：

>>1-binocdf（100,162,0.5）

ans=

0.0010433

相关函数

binofit|binoinv|binopdf|binornd|binostat|cdf

附：

二项式分布（binomialdistribution）

定义

二项分布的概率密度函数为

wherekisthenumberofsuccessesinntrialsofaBernoulliprocesswithprobabilityofsuccessp.

Thebinomialdistributionisdiscrete,definedforintegersk=0,1,2,...n,whereitisnonzero.

背景

Thebinomialdistributionmodelsthetotalnumberofsuccessesinrepeatedtrialsfromaninfinitepopulationunderthefollowingconditions:

Onlytwooutcomesarepossibleoneachofn trials.

Theprobabilityofsuccessforeachtrialisconstant.

Alltrialsareindependentofeachother.

ThebinomialdistributionisageneralizationoftheBernoullidistribution;itgeneralizestothemultinomialdistribution.

参数

Supposeyouarecollectingdatafromawidgetmanufacturingprocess,andyourecordthenumberofwidgetswithinspecificationineachbatchof 100.Youmightbeinterestedintheprobabilitythatanindividualwidgetiswithinspecification.Parameterestimationistheprocessofdeterminingtheparameter,p,ofthebinomialdistributionthatfitsthisdatabestinsomesense.

Onepopularcriterionofgoodnessistomaximizethelikelihoodfunction.Thelikelihoodhasthesameformasthebinomialpdfabove.Butforthepdf,theparameters（n and p）areknownconstantsandthevariableis x.Thelikelihoodfunctionreversestherolesofthevariables.Here,thesamplevalues（thex's）arealreadyobserved.Sotheyarethefixedconstants.Thevariablesaretheunknownparameters.MLEinvolvescalculatingthevalueofpthatgivethehighestlikelihoodgiventheparticularsetofdata.

ThefunctionbinofitreturnstheMLEsandconfidenceintervalsfortheparametersofthebinomialdistribution.Hereisanexampleusingrandomnumbersfromthebinomialdistributionwithn = 100andp = 0.9.

>>r=binornd（100,0.9）

r=

85

>>[phat,pci]=binofit（r,100）

phat=

0.85

pci=

0.764690.91355

TheMLEforparameter pis 0.8800,comparedtothetruevalueof 0.9.The95%confidenceintervalfor pgoesfrom0.7998to 0.9364,whichincludesthetruevalue.Inthismade-upexampleyouknowthe"truevalue"of p.Inexperimentationyoudonot.

示例

Thefollowingcommandsgenerateaplotofthebinomialpdfforn=10andp=1/2.

x=0:

10;

y=binopdf（x,10,0.5）;

plot（x,y,'+'）

相关内容

DiscreteDistributions

附：

二项式分布（网上）

定义

若某事件概率为p，现重复试验n次，该事件发生k次的概率为:

P=C（k,n）pk（1-p）（n-k）

C（k,n）表示组合数，即从n个事物中拿出k个的方法数。

二项分布的概念

考虑只有两种可能结果的随机试验，当成功的概率（π）是恒定的，且各次试验相互独立，这种试验在统计学上称为贝努里试验（Bernoullitrial）。

如果进行n次贝努里试验，取得成功次数为X（X=0,1,…,n）的概率可用下面的二项分布概率公式来描述：

P=C（X,n）πX（1-π）（n-X）

式中的n为独立的贝努里试验次数，π为成功的概率，（1-π）为失败的概率，X为在n次贝努里试验中出现成功的次数，C（X,n）表示在n次试验中出现X的各种组合情况，在此称为二项系数（binomialcoefficient）。

内容简介

二项分布，伯努里分布：

进行一系列试验,如果

1.在每次试验中只有两种可能的结果,而且是互相对立的；

2.每次实验是独立的，与其它各次试验结果无关；

3.结果事件发生的概率在整个系列试验中保持不变，则这一系列试验称为伯努力试验。

在这试验中，事件发生的次数为一随机事件,它服从二次分布。

二项分布可以用于可靠性试验。

可靠性试验常常是投入n个相同的式样进行试验T小时，而只允许k个式样失败，应用二项分布可以得到通过试验的概率。

一个事件必然出现，就说它100%要出现。

100%=1，所以100%出现的含义就是出现的概率P=1。

即必然事件的出现概率为1。

若掷一枚硬币，正面向上的结果的概率为0.5。

反面向上的结果的概率也是0.5。

则出现正面向上事件或者反面向上事件的概率就是0.5+0.5=1，即二者必居其一。

若掷两次硬币，由独立事件的概率乘法定理那么两次都是正面（反面）向上的概率是0.5×0.5=0.25。

另外第一个是正第二个是反的出现概率也是0.5×0.5=0.25。

同理第一个反第二个正的出现概率也是0.5×0.5=0.25。

于是一正一反的概率是前面两个情况的和，即0.25+0.25=2×0.25=0.5。

它们的合计值仍然是1。

binopdf二项分布的概率密度函数

语法格式

Y=binopdf（X,N,P）

函数功能

Y=binopdf（X,N,P）

计算X中每个X（i）的概率密度函数，其中，N中对应的N（i）为试验数，P中对应的P（i）为每次试验成功的概率。

Y,N,和P的大小类型相同，可以是向量、矩阵或多维数组。

输入的标量将扩展成一个数组，使其大小类型与其它输入相一致。

N中的值必须是正整数，0≤P≤1。

对于给出的x和两个参数n和p，二项分布概率密度函数为

其中q=1–p。

y为n次独立试验中成功x次的概率，其中，每次试验成功的概率为p。

指示器函数I（0,1,...,n）（x）确保x取值为0,1,...,n。

示例

一个质量检查技术员一天能测试200块电路板。

假设有2%的电路板有缺陷，该技术员在一天的测试中没有发现有缺陷的电路板的概率是多少？

>>binopdf（0,200,0.02）

ans=

0.0176

质量检查技术员一天中最大可能检测出有缺陷的电路板是多少块？

>>defects=0:

200;

>>y=binopdf（defects,200,.02）;

>>[x,i]=max（y）;

>>defects（i）

ans=

4

相关函数

binocdf|binofit|binoinv|binornd|binostat|pdf

dlmread将以ASCII码分隔的数值数据文件读入到矩阵中

语法格式

M=dlmread（filename）

M=dlmread（filename,delimiter）

M=dlmread（filename,delimiter,R,C）

M=dlmread（filename,delimiter,range）

函数功能

M=dlmread（filename）

读取有分隔符的ASCII数值数据文件filename，并把数据给到矩阵M中。

文件名（filename）以字符串形式用单引号括起来。

dlmread从文件的格式中推断分隔符。

M=dlmread（filename,delimiter）

指定分隔符delimiter。

使用'\t'代表制表符tab作为分隔。

M=dlmread（filename,delimiter,R,C）

R和C指定了数据

读取数据，其左上角的值位于文件中的第R行，第C列。

R和C从0开始，所以R=0,C=0指向文件中的第一个值。

M=dlmread（filename,delimiter,range）

读取由range=[R1C1R2C2]指定的区域块，（R1,C1）是左上角，（R2,C2）是右下角。

也可以使用电子表格表示法来指定，如range='A1..B7'。

备注

•Alldataintheinputfilemustbenumeric.dlmreaddoesnotreadfilesthatcontainnonnumericdata,evenifthespecifiedrowsandcolumnscontainonlynumericdata.

•若没有指定分隔符，当从文件格式中推断分隔符时，连续的空格符当作一个分隔符对待。

若指定了分隔符，则重复的分隔符将分别作为单独的分隔符对待。

•IfyouwanttospecifyanR,C,orrangeinput,butnotadelimiter,setthedelimiterargumenttotheemptystring,（twoconsecutivesinglequoteswithnospacesinbetween,''）.Forexample,

M=dlmread（'myfile.dat','',5,2）

Inthiscase,dlmreadtreatsrepeatedwhitespacesasasingledelimiter.

•Dlmread将用0填充没有边界的区域。

有多行的数据文件，若以非空格分隔符结束，例如分号，则导入后会多产生全0的一列在最后。

•Dlmread在导入任何复数时，将其作为一个整体导入到一个复数单元中。

下面是有效的复数格式：

±±i|jExample:

5.7-3.1i

±i|jExample:

-7j

嵌入了空格的复数是不正确的格式，空格将被认为是分隔符。

示例

例1

Exporta5-by-8testmatrixMtoafile,andreaditwithdlmread,firstwithnoargumentsotherthanthefilename:

M=gallery（'integerdata',100,[58],0）;

dlmwrite（'myfile.txt',M,'delimiter','\t'）

dlmread（'myfile.txt'）

ans=

96776241621242

2446809436207585

612939282614553

498374421289421

9045189014204768

Nowreadaportionofthematrixbyspecifyingtherowandcolumnoftheupperleftcorner:

dlmread（'myfile.txt','\t',2,3）

ans=

9282614553

421289421

9014204768

Thistime,readadifferentpartofthematrixusingarangespecifier:

dlmread（'myfile.txt','\t','C1..G4'）

ans=

62416212

8094362075

9392826145

744212894

例2

ExportmatrixMtoafile,andthenappendanadditionalmatrixtothefilethatisoffsetonerowbelowthefirst:

M=magic（3）;

dlmwrite（'myfile.txt',[M*5M/5],''）

dlmwrite（'myfile.txt',[M/3],'-append',...

'roffset',1,'delimiter',''）

typemyfile.txt

405301.60.21.2

1525350.611.4

2045100.81.80.4

2.66670.333332

11.66672.3333

1.333330.66667

Whendlmreadimportsthesetwomatricesfromthefile,itpadsthesmallermatrixwithzeros:

dlmread（'myfile.txt'）

40.00005.000030.00001.60000.20001.2000

15.000025.000035.00000.60001.00001.4000

20.000045.000010.00000.80001.80000.4000

2.66670.33332.0000000

1.00001.66672.3333000

1.33333.00000.6667000

替代

作为dlmread的替代，可使用导入向导。

选择菜单“File|ImportData”激活导入向导。

相关函数

dlmwrite|textscan

fprintf写数据到文本文件或输出到命令窗口

语法格式

fprintf（fileID,format,A,...）

fprintf（format,A,...）

count=fprintf（...）

函数功能

fprintf（fileID,format,A,...）

fileID为文件句柄，指定要写入的文件；format是用来控制所写数据格式的格式符；A是用来存放数据的矩阵。

appliestheformattoallelementsofarrayAandanyadditionalarrayargumentsincolumnorder,andwritesthedatatoatextfile.fprintfusestheencodingschemespecifiedinthecalltofopen.

fprintf（format,A,...）

将A中的数据按格式format输出到命令窗口。

count=fprintf（...）

返回fprintf写的字节数。

输入参数

fileID

取值为下面中的一种：

✧从文件打开时所得到的文件句柄，为一个整数值。

✧1为标准输出设备（屏幕）。

✧2forstandarderror.

fileID缺省时，输出到屏幕，即1为默认值。

format

单引号括起来的是“转换控制字符串”，指定输出的格式。

由下面给出的各部分的组合而成：

✧百分号后紧跟的格式字符，如'%s'。

✧输出区域的宽度，精度，和其它选项（附加格式说明符）。

✧要原样输出的文本（普通字符）。

✧转义字符，包括：

''单引号

%%百分号

\\反斜线

\a响铃

\b退格符

\fFormfeed

\n换行

\r回车

\t水平制表符

\v垂直制表符

\xN十六进制数,N

\N八进制数,N

转换字符和可选的操作符按以下顺序排列（包括空格）：

下表列出可用的转换字符和限定类型字符。

数据类型

转换

说明

Integer,signed

%d或%i

十进制整数

%ld或%li

64-bitbase10values

%hd或%hi

16-bitbase10values

Integer,unsigned

%u

十进制整数

%o

八进制整数（octal）

%x

十六进制整数（hexadecimal）,使用字母a–f

%X

十六进制整数（hexadecimal）,使用字母A–F

%lu

%lo

%lxor%lX

64-bitvalues,十、八、十六进制整数

%hu

%ho

%hxor%hX

16-bitvalues,十、八、十六进制整数

Floating-pointnumber

%f

定点表示法

%e

指数表示法，如3.141593e+00

%E

同%e，但用E，如3.141593E+00

%g

%e或%f的紧凑格式，没有尾部的0

%G

%E或%f的紧凑格式，没有尾部的0

%bxor%bX

%bo

%bu

双精度的十六、八、十进制数值

如：

%bxprintspias400921fb54442d18

%txor%tX

%to

%tu

单精度的十六、八、十进制数值

如：

%txprintspias40490fdb

Characters

%c

单个字符

%s

字符串

附加格式说明符，包括：

✧字段宽度

Minimumnumberofcharacterstoprint.Canbeanumber,oranasterisk（*）torefertoanargumentintheinputlist.Forexample,theinputlist（'%12d',intmax）isequivalentto（'%*d',12,intmax）.

✧精度

For%f,%e,or%E:

Numberofdigitstotherightofthedecimalpoint.

Example:

'%6.4f'printspias'3.1416'

For%gor%GNumberofsignificantdigits.

Example:

'%6.4g'printspias' 3.142'

Canbeanumber,oranasterisk（*）torefertoanargumentintheinputlist.Forexample,theinputlist（'%6.4f',pi）isequivalentto（'%*.*f',6,4,pi）.

✧Flags

Action

Flag

Example

Left-justify.

'–'

%-5.2f

Printsigncharacter（+or–）.

'+'

%+5.2f

Insertaspacebeforethevalue.

''

%5.2f

Padwithzeros.

'0'

%05.2f

Mo