%基于Monte?Carlo方法的試題難度分布研究
%?作者：姚進(jìn)等
%?計(jì)算機(jī)應(yīng)用與軟件?2004.9

%基本思路：
%1.定義有效的難度值范圍[0，1]，并保證正態(tài)分布下的該區(qū)間面積不小于0.995。找到最大sigma值。
%2.作者認(rèn)為，難度一般分布在[0.3，0.7]。
%3.固定u值為0.5，0.6，0.7，0.8，求最大sigma值，公式為P{（x‘-Z*u）/sigma=0.995?
%4.得到難度的pdf后，進(jìn)行模擬。
%5.統(tǒng)計(jì)2000次的模擬值，包括：各次的平均難度、各難度區(qū)間的頻率。
%6.求模擬統(tǒng)計(jì)值的期望。得到最終結(jié)果。

%幾個(gè)基本的變量
%?難度等級(jí)A1=‘易’；A2=‘較易’；A3=‘中等’；A4=‘較難’；A5=‘難’
%?題目數(shù)：itemNUM
%?模擬次數(shù)：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（‘????組卷參數(shù)?????????????????????????????????????????????????計(jì)算結(jié)果?