MATLAB函数和命令的用法Word文档下载推荐.docx
《MATLAB函数和命令的用法Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《MATLAB函数和命令的用法Word文档下载推荐.docx(21页珍藏版)》请在冰豆网上搜索。
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.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)为每次试验成功的概率。
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)
0.0176
质量检查技术员一天中最大可能检测出有缺陷的电路板是多少块?
defects=0:
200;
y=binopdf(defects,200,.02);
[x,i]=max(y);
defects(i)
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)
读取有分隔符的ASCII数值数据文件filename,并把数据给到矩阵M中。
文件名(filename)以字符串形式用单引号括起来。
dlmread从文件的格式中推断分隔符。
指定分隔符delimiter。
使用'
\t'
代表制表符tab作为分隔。
R和C指定了数据
读取数据,其左上角的值位于文件中的第R行,第C列。
R和C从0开始,所以R=0,C=0指向文件中的第一个值。
读取由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在导入任何复数时,将其作为一个整体导入到一个复数单元中。
下面是有效的复数格式:
±
<
real>
imag>
i|jExample:
5.7-3.1i
-7j
嵌入了空格的复数是不正确的格式,空格将被认为是分隔符。
例1
Exporta5-by-8testmatrixMtoafile,andreaditwithdlmread,firstwithnoargumentsotherthanthefilename:
M=gallery('
integerdata'
100,[58],0);
dlmwrite('
myfile.txt'
M,'
delimiter'
)
dlmread('
96776241621242
2446809436207585
612939282614553
498374421289421
9045189014204768
Nowreadaportionofthematrixbyspecifyingtherowandcolumnoftheupperleftcorner:
2,3)
9282614553
421289421
9014204768
Thistime,readadifferentpartofthematrixusingarangespecifier:
C1..G4'
62416212
8094362075
9392826145
744212894
例2
ExportmatrixMtoafile,andthenappendanadditionalmatrixtothefilethatisoffsetonerowbelowthefirst:
M=magic(3);
[M*5M/5],'
'
[M/3],'
-append'
...
roffset'
1,'
typemyfile.txt
405301.60.21.2
1525350.611.4
2045100.81.80.4
2.66670.333332
11.66672.3333
1.333330.66667
Whendlmreadimportsthesetwomatricesfromthefile,itpadsthesmallermatrixwithzeros:
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.
%6.4g'
3.142'
Canbeanumber,oranasterisk(*)torefertoanargumentintheinputlist.Forexample,theinputlist('
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