prob_tgamma.asv

matlab从入门到精通第6章书后源码对初学者很有帮助

% prob_tgamma.m
% 绘制两个相关变量的随机分布图形
subplot(111);
% 设置概率分布参数
% n数据点个数，nu表示自由度
n=1000;rho=0.7;nu=1;
% 产生多元t分布随机数值矩阵
T=mvtrnd([1,rho;rho 1],nu,n);
% 计算t分布数值的累计概率分布数值
U=tcdf(T,nu);
% 产生两个概率分布的数值
X=[gaminv(U(:,1),2,1) tinv(U(:,2),5)];
% 计算两个直方图的数值
[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 ')