y=randn(1,3000); % 生成1-by-3000的標準正態分佈隨機數 ymin=min(y); ymax=max(y); x=linspace(ymin,ymax,20); %將最大最小區間分紅20個等分點(19等分),而後分別計算各個區間的個數 yy=hist(y,x); %計算各個區間的個數 yy=yy/length(y); %計算各個區間樣本數所佔總樣本數的比例(機率) bar(x,yy) %畫出機率密度分佈圖 s=0; for i=2:length(x) s=[s,trapz(x([1:i]),yy([1:i]))]; % 使用trapz模擬積分 end figure; plot(x,s,'*-')