% prob_tgamma.m
% 繪制兩個相關變量的隨機分布圖形
% 設置概率分布參數(shù)
% n數(shù)據(jù)點個數(shù)，nu表示自由度
n=1000;rho=0.9;nu=1;
% 產(chǎn)生多元t分布隨機數(shù)值矩陣
T=mvtrnd([1,rho;rho 1],nu,n);
% 計算t分布數(shù)值的累計概率分布數(shù)值
U=tcdf(T,nu);
% 產(chǎn)生兩個概率分布的數(shù)值
% gaminv為Gamma累計概率數(shù)值取反
% tinv為t分布累計概率數(shù)值取反
X=[gaminv(U(:,1),3,1) tinv(U(:,2),5)];
% 計算兩個直方圖的數(shù)值
[n1,ctr1]=hist(X(:,1),20);
[n2,ctr2]=hist(X(:,2),20);
% 繪制概率分布圖形
subplot(2,2,2);
plot(X(:,1),X(:,2),'.');
axis([0 15 -10 10]);
h1=gca;
title('1000 simulated dependent t and gamma');
xlabel('X1：Gamma(3,1)');
ylabel('X2： t(5)');
% 繪制gamma自變量分布的直方圖
subplot(2,2,4);bar(ctr1,-n1,1);
axis([0 15 -max(n1)*1.1 0]);
axis('off');
h2=gca;
% 繪制t分量分布的直方圖
subplot(2,2,1);barh(ctr2,-n2,1);
axis([-max(n2)*1.1 0 -10 10]);
axis('off');
h3=gca;
% 設置三個子圖的位置
set(h1,'Position',[0.35 0.35 0.55 0.55]);
set(h2,'Position',[0.35 0.1  0.55 0.15]);
set(h3,'Position',[ 0.1 0.35 0.15 0.55]);
colormap([0.8 0.8 1]);