1.直方圖

　　clear all; %清除數據

　　clc; %清屏

　　mu1 = ［0］; %均值向量

　　sigma1=［1］; %協(xié)方差向量

　　a=input（‘請輸入數據個數：’） %a為數據個數

　　w1 = mvnrnd（mu1，sigma1，a）; %創(chuàng)建a個數據，均值與協(xié)方差分別為mu1和sigma1，服從正態(tài)分布

　　b=input（‘請輸入分組個數：’） %b為分組個數

　　max = max（w1）; %max為a個數據中的最大值

　　min = min（w1）; %min為a個數據中的最小值

　　h = （max-min）/b; %h為每個分組的寬度

　　sum = zeros（1，b）; %創(chuàng)建一個1*b的向量，用來計每組的數據數

　　for i=1:a %循環(huán)a個數據

　　for j=1:b %循環(huán)，尋找數據所在的分組

　　if（w1（i）《min+j*h） %如果在第j組

　　sum（j） = sum（j）+1; %第j組的數據數加1

　　break; %跳出循環(huán)

　　end %與if對應的結束

　　end %結束

　　end %結束

　　x=zeros（1，b）; %創(chuàng)建一個1*b的向量x，用來裝每組的中點

　　for i=1:b %循環(huán)b次，將b組的中點放進x

　　x（i） = min+（2*i-1）/2*h; %第i組的中點

　　end %標志for的結束

　　axis（［min max 0 1］） %控制x，y坐標軸的范圍

　　sum = sum/a/h; %將縱坐標變?yōu)楦琶?/p>

　　bar（x，sum） %畫出直方圖

　　100個數據，10組：