ch8example1prog1.m

清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB和SIMULINK代码

% ch8example1prog1.m
n1=5; n2=4; % F分布的参数
n=10000;    % 随机数样本数量
x=frnd(n1,n2,n,1);  % 随机数样本产生
a=min(x); b=max(x); % 样本值域区间计算
m=200;              % 分组区间数 or m=500 等
Delta=(b-a)/m;      % 分组间隔
[r,xout]=hist(x,[a:Delta:b]);    % 计算直方图数据
pdf=r./(n*Delta);                % 计算统计频率密度
bar(xout,pdf);hold on;           % 作出频率密度直方图
h = findobj(gca,'Type','patch'); % 修改直方图样式
set(h,'FaceColor',[0.7,0.7,0.7],'EdgeColor','k');
X=0:0.01:10; % 计算并画出F分布的理论概率密度函数曲线
Y = fpdf(X,n1,n2); 
plot(X,Y,'k-');axis([0 10 0 1]);