蒙特卡洛MATLAB程序

%基于Monte Carlo方法的试题难度分布研究
% 作者：姚进等
% 计算机应用与软件 2004.9

%基本思路：
%1.定义有效的难度值范围[0，1]，并保证正态分布下的该区间面积不小于0.995。找到最大sigma值。
%2.作者认为，难度一般分布在[0.3，0.7]。
%3.固定u值为0.5，0.6，0.7，0.8，求最大sigma值，公式为P{（x‘-Z*u）/sigma=0.995 
%4.得到难度的pdf后，进行模拟。
%5.统计2000次的模拟值，包括：各次的平均难度、各难度区间的频率。
%6.求模拟统计值的期望。得到最终结果。

%几个基本的变量
% 难度等级A1=‘易’；A2=‘较易’；A3=‘中等’；A4=‘较难’；A5=‘难’
% 题目数：itemNUM
% 模拟次数：simuNUM
function main（）
itemNUM=2000;
simuNUM=2000;
cols=itemNUM;
A=zeros（simuNUM5）;
for simuLoop=1:simuNUM
    nandu=myrandn（0.60.1551cols）;
    expectation（simuLoop）=mean（nandu）;
    frequency=getndfreq（nandu）;
    frequency2（simuLoop:）=frequency;
end
nandu_Mean=mean（expectation）;
nandu_Sigma=std（expectation1）;
nandu_Distribution=mean（frequency2）;

% fprintf（‘    组卷参数                                                 计算结果