前言:繪製連續型的分布函式很容易,直接根據分布函式計算函式值即可。但是對於離散型隨機變數而言,沒有已知的分布函式,只能使用經驗分布函式或者說是累積分布函式進行近似。下文以離散分布中比較經典的二項分布、泊松分布以及幾何分布為例繪製它們的經驗分布函式。
% 二項分布bin(n,p)
n =50
;% 二項分布引數n
p =0.5
;% 二項分布引數p
f_bin = bin(n,p,len);%
%% 二項分布bin(n,p)
function f_bin = bin(n,p,len)
u = rand(
1,len)
;bin
= binoinv(u,n,p)
;% 二項分布逆函式
% 畫分布函式
figure
xlim([0
round
(1.1
*n)]
)for i =1:
round
(1.1
*n) f_bin(i)
= length(find(
bin)/len;
% 二項分布的分布函式值
line(
[i-1 i]
,[f_bin(i) f_bin(i)],
'linewidth'
,1.1
)% 繪製分段的分布函式影象
hold on
endstr1 =
'二項分布分布函式影象'
;str2 =
['n = ' num2str(n)
' p = ' num2str(p)];
title(
)grid on
set(gca,
'linewidth'
,1.1
)end
經驗分布圖形為:
以下泊松分布理論部分摘自泊松分布部落格。
% 泊松分布poi(
lambda
)lambda=2
;x = poiss(
lambda
, len)
;% 泊松分布poi(
lambda
)function f_poiss = poiss(
lambda
,len)
u = rand(
1,len)
;poiss = poissinv(u,
lambda);
% 泊松逆函式
% 畫分布函式圖
figure
xlim([0
round
(max
(poiss))]
)for i =1:
round
(max
(poiss)
) f_poiss(i)
= length(find(poiss
)/len;
% 泊松分布的分布函式值
line(
[i-1 i]
,[f_poiss(i) f_poiss(i)],
'linewidth'
,1.1);
% 繪製分段的分布函式影象
hold on
endstr1 =
'泊松分布分布函式影象'
;str2 =
['lambda = ' num2str(
lambda)]
;title(
)grid on
set(gca,
'linewidth'
,1.1
)end
經驗分布函式影象為:
% 幾何分布
p =0.5
;f_geo = geo(p,len);%
%% 幾何分布geo(p)
function f_geo = geo(p,len)
u = rand(
1,len)
;geo = geoinv(u,p)
;% 幾何分布逆函式
% 畫分布函式
figure
xlim([0
max(geo)])
for i =1:
max(geo)
f_geo(i)
= length(find(geo
)/len;
% 二項分布的分布函式值
line(
[i-1 i]
,[f_geo(i) f_geo(i)],
'linewidth'
,1.1
)% 繪製分段的分布函式影象
hold on
endstr1 =
'二項分布分布函式影象'
;str2 =
['p = ' num2str(p)];
title(
)grid on
set(gca,
'linewidth'
,1.1
)end
經驗分布函式影象為:
概率論知識回顧(五) 隨機變數,離散隨機變數分布
重點 隨機變數,離散隨機變數分布 知識回顧用於鞏固知識和查漏補缺。知識回顧步驟 檢視知識回顧中的問題,嘗試自己解答 自己解答不出來的可以檢視下面的知識解答鞏固知識。對知識解答有疑問的,說明有關這一點的知識或者公式沒有理解透徹或者沒有記住,要重新翻看書籍。什麼是隨機變數?隨機變數的作用是什麼?為什麼要...
連續型隨機變數
1.對於乙個連續型隨機變數,它取任何固定值的概率都等於0。因此,對於連續隨機變數,下式成立 f a a f x dx p a f x dx df a da f a f a 可看作隨機變數取值於點a附近的可能性的乙個度量。3.連續型隨機變數的期望e x xf x dx,方差可根據var x e x2 ...
隨機變數函式的分布與聯合分布
以下文字均假設所求分布存在。設函式y f x 單調遞增,若隨機變數x變化範圍是 x,x dx 且由此引起的y變化範圍是 y,y dy 那麼 p x x x dx p y y y dy 若單調遞減,dy 0,上式就成為 p x x x dx p y dy y 也就是 x,x dxf x dx y,y ...