1、第四講 統(tǒng)計(jì)作圖 11、直方圖根據(jù)實(shí)驗(yàn)數(shù)據(jù)，會(huì)出直方圖，來(lái)顯示數(shù)據(jù)的分布特征，進(jìn)而觀察實(shí)驗(yàn)數(shù)據(jù)所反映的統(tǒng)計(jì)規(guī)律調(diào)用格式：n = hist(Y) 將Y中的元素分到10個(gè)間隔相同的條形中，并返回每個(gè)條形中的元素個(gè)數(shù)n = hist(Y,x) x為向量，返回Y的分布，如x為一個(gè)5個(gè)元素向量，則hist函數(shù)將Y中的元素分配到五組分組中直方圖的x軸反映Y中值得范圍，y軸顯示落到組中元素的個(gè)數(shù)。21、直方圖【例1】x = -4:0.1:4; y=normrnd(0,1,1,10000); hist(y,x) 32、經(jīng)驗(yàn)累積分布函數(shù)圖形 函數(shù) cdfplot格式 cdfplot(X) %作樣本X（向量）的累

2、積分布函數(shù)圖形h = cdfplot(X) %h表示曲線的環(huán)柄 h,stats = cdfplot(X) %stats表示樣本的一些特征 4【例2】 X=normrnd (0,1,50,1); h,stats=cdfplot(X)h = 3.0013stats = min: -1.8740 %樣本最小值 max: 1.6924 %最大值 mean: 0.0565 %平均值 median: 0.1032 %中間值 std: 0.7559 %樣本標(biāo)準(zhǔn)差 圖4-105【例3】比較經(jīng)驗(yàn)累積分布函數(shù)和理論分布函數(shù)。 x=normrnd(0,1,50,1); h,stats=cdfplot(x); hol

3、d on y=normcdf(x,0,1); plot(x,y,r+)63、最小二乘擬合直線73、最小二乘擬合直線 函數(shù) lsline格式 lsline %最小二乘擬合直線 h = lsline %h為直線的句柄8【例4】 X = 2 3.4 5.6 8 11 12.3 13.8 16 18.8 19.9; plot(X,+) lsline94、繪制正態(tài)分布概率圖形 函數(shù) normplot作用：判斷總體或者樣本是否服從正態(tài)分布格式 normplot(X) %若X為向量，則顯示正態(tài)分布概率圖形，若X為矩陣，則顯示每一列的正態(tài)分布概率圖形。h = normplot(X) %返回繪圖直線的句柄說(shuō)明

4、樣本數(shù)據(jù)在圖中用“+”顯示；如果數(shù)據(jù)來(lái)自正態(tài)分布，則圖形顯示為直線，而其它分布可能在圖中產(chǎn)生彎曲。 10【例5】產(chǎn)生50個(gè)標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)數(shù)和指數(shù)分布的隨機(jī)數(shù)，并畫(huà)出它們的正態(tài)分布概率圖形 X=normrnd(0,1,50,1); y=exprnd(1,1,50) normplot(X) normplot(y)115、樣本數(shù)據(jù)的盒圖 盒圖或箱型圖，用來(lái)比較數(shù)據(jù)的分布特征函數(shù) boxplot格式 boxplot(X) %產(chǎn)生矩陣X的每一列的盒圖和“須”圖，“須”是從盒的尾部延伸出來(lái)，并表示盒外數(shù)據(jù)長(zhǎng)度的線，如果“須”的外面沒(méi)有數(shù)據(jù)，則在“須”的底部有一個(gè)點(diǎn)。boxplot(X,notch) %

5、當(dāng)notch=1時(shí)，產(chǎn)生一凹盒圖，notch=0時(shí)產(chǎn)生一矩箱圖。boxplot(X,notch,sym) %sym表示圖形符號(hào)，默認(rèn)值為“+”。boxplot(X,notch,sym,vert) %當(dāng)vert=0時(shí)，生成水平盒圖，vert=1時(shí)，生成豎直盒圖（默認(rèn)值vert=1）。boxplot(X,notch,sym,vert,whis) %whis定義“須”圖的長(zhǎng)度，默認(rèn)值為1.5，若whis=0則boxplot函數(shù)通過(guò)繪制sym符號(hào)圖來(lái)顯示盒外的所有數(shù)據(jù)值 12【例6】x1 = normrnd(5,1,100,1);x2 = normrnd(6,1,100,1);x = x1 x2; b

6、oxplot(x) boxplot(x,1) boxplot(x,1,g+,1,0) boxplot(x,1,g+,0,0)136、樣本的概率圖形 獲得隨機(jī)變量落在指定范圍內(nèi)的概率，獲得樣本分布的概率密度圖形函數(shù) capaplot格式 p = capaplot(data,specs) %data為所給樣本數(shù)據(jù)，specs指定范圍，p表示在指定范圍內(nèi)的概率。14【例7】 data=normrnd (0,1,30,1); p=capaplot(data,-2,2)p = 0.9199 157、 附加有正態(tài)密度曲線的直方圖 函數(shù) histfit格式 histfit(data) %data為向量，返回

7、直方圖和正態(tài)曲線histfit(data,nbins) % nbins指定bar的個(gè)數(shù)，缺省時(shí)為data中數(shù)據(jù)個(gè)數(shù)的平方根。 16【例8】產(chǎn)生100個(gè)均值為10，標(biāo)準(zhǔn)差為1的正態(tài)分布的隨機(jī)數(shù)，畫(huà)出它們的直方圖并附加正態(tài)密度曲線，觀察它們之間的擬合程度。r = normrnd (10,1,100,1);histfit(r)178、在指定的界線之間畫(huà)正態(tài)密度曲線 函數(shù) normspec格式 p = normspec(specs,mu,sigma) %specs指定界線，mu,sigma為正態(tài)分布的參數(shù)p 為樣本落在上、下界之間的概率。 18【例9】normspec(10 Inf,11.5,1.25

8、)199、散度圖散度圖是顯示一個(gè)變量對(duì)另一個(gè)變量的一種簡(jiǎn)單圖形，它可以用于確定兩個(gè)變量的值或者兩個(gè)變量之間的關(guān)系是否屬于同一組。格式：gscatter(x,y,group)gscatter(x,y,group,clr,sym,siz)gscatter(x,y,group,clr,sym,siz,doleg)gscatter(x,y,group,clr,sym,siz,doleg,xnam,ynam)h = gscatter(.) 其中：x，y是具有相同大小的向量；group是組標(biāo)記；clr,sym是繪圖的顏色和符號(hào)；siz是大小的向量；doleg是否顯示圖標(biāo)記，nam,ynam是x軸和y軸的名稱；h返回圖形中直線的句柄。20【例10】比較3種不同年代汽車的重量和里程數(shù) load carsmall gscatter(Weight,MPG,Model_Year,xos); xlabel(重量); ylabel(里程量);21【例11】比較城市氣候和住房密度的關(guān)系 load discrim gscatter(ratings(:,1),ratings(:,2),group,br,xo)2223作業(yè)1、產(chǎn)生100個(gè)標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)數(shù)，指出它們的分布特征，并畫(huà)出經(jīng)驗(yàn)累積分布函數(shù)圖。2、產(chǎn)生100個(gè)均勻分布的隨機(jī)數(shù)，對(duì)這100個(gè)數(shù)據(jù)的列向量，用“*”