1.簡単な説明
確率分布則と密度関数値の計算
Matlab は、確率密度関数の値を計算する一般関数を直接提供します。これらは pdf 関数と namepdf 関数であり、次のように使用されます: Y=pdf
('name', K, A, B) または: namepdf (K, A , B)
上記 この関数は、パラメーター A、B、および C を使用して X=K における確率値または密度値を返すことを意味します。分布が異なると、パラメーターの数も異なります。名前は分布関数の名前です。使用時に対応するディストリビューションに応じて変更されます。関数名を次の表にまとめます。
名前の値 関数 説明
'beta' または 'Beta' ベータ分布
'bino' または 'Binomial' 二項分布
'chi2' または 'Chisquare' カイ二乗分布
'exp' または 'Exponential' 指数分布
'f' または 'F ' F 分布
'gam' または 'Gamma' ガンマ分布
'geo' または 'Geometric' 幾何分布
'hyge' または 'Hypergeometric' 超幾何分布 '
logn' または 'Lognormal' 対数正規分布
'nbin' または 'Negative Binomial' ' 負の二項分布
'ncf' または '非心 F' 非心 F 分布
'nct' または '非心 t' 非心 t 分布
'ncx2' または '非心カイ二乗' 非心カイ二乗分布
'norm' または 'Normal' 正規分布
'poiss' または'Poisson' ポアソン分布
'rayl' または 'Rayleigh' レイリー分布
't' または 'T' T 分布
'unif' または 'Uniform' 連続一様分布
'unid' または 'Discrete' Unique' 離散一様分布
「weib」または「Weibull」ワイブル分布
2. コードと実行結果
%% 二項分布の密度関数
すべてクリア;
x=1:20;
y=binopdf(x,200,0.06);
図;
プロット(x,y,'r*');
title('二項分布 (n =200) 、p=0.06)');
%% ポアソン分布密度関数
すべてクリア;
x=1:20;
y=poisspdf(x,20); %ポアソン分布
図;
plot(x,y,'r+');
title('ポアソン分布' );
%% 幾何分布
すべてクリア;
x=1:10;
y=geopdf(x,0.4); % 幾何分布
図;
lot(x,y,'rx');
title('幾何分布');
%% 一様分布 (離散)
すべてクリア;
n=10;
x=1:n;
y=unidpdf(x,n); % 一様分布 (離散)
図;
plot(x,y,'ro');
title( '一様分布 (離散)');
%% 一様分布 (連続)
すべてクリア;
x=-2:0.1:15;
y=unifpdf(x,0,6); % 0 と 6 の間の一様分布 (連続)
Figure;
plot(x,y,' r :');
title('均一分布(連続)');
%% 指数分布
をすべてクリア;
x=0:0.1:10;
y=exppdf(x,2); % 指数分布の
図;
Lot(x,y,'r:');
title('指数分布');
%% 正規分布
すべてクリア;
x=-5:0.1:5;
y1=normpdf(x,0,1); % 標準正規分布
y2=normpdf(x,3,3); % 非標準正規分布
図;
プロット(x,y1,x,y2,':');
legend('標準正規分布','非標準正規分布');
x1=-5:0.1:5;
y3=normpdf(x1, 3,1 ); %SIGMA=1
y4=normpdf(x1,3,2); %SIGMA=2
y5=normpdf(x1,3,3); %SIGMA=3
形
プロット(x1,y3,'r-',x1,y4,'b:',x1,y5,'k--');
legend('SIGMA=1','SIGMA=2','SIGMA=3');
y6=normpdf(x1,0,2); %MU=0
y7=normpdf(x1,2,2); %MU=2
y8=normpdf(x1,4,2); %MU=4
形
プロット(x1,y6,'r-',x1,y7,'b:',x1,y8,'k--');
legend('MU=0','MU=2','MU=4');
%% 3 つの主要な標本分布の確率密度関数
%% カイ二乗分布
すべてクリア;
x=0:0.1:15;
y1=chi2pdf(x,2); %カイ二乗分布 n=2
y2=chi2pdf( x,3) ; % カイ二乗分布 n=3
図;
ちょっと待って;
プロット(x,y1);
プロット(x,y2,':');
凡例('n=2','n=3') ;
title('カード広場配布');
%% t 分布
をすべてクリア;
x=-5:0.1:5;
y1=tpdf(x,2); %t 分布 (n=2)
y2=tpdf(x,10); %t 分布 (n=10)
図;
プロット(x,y1,'r:',x,y2,'b-');
凡例('n=2','n=10');
タイトル('t分布');
%% F 分布
をすべてクリア;
x=0.1:0.1:5;
y=fpdf(x,2,5); %F 分布
図;
plot(x,y,'r:');
title('F 分布 (m =2,n=5)');
