Matlab Simulink (Квантование)

p=твое число; %вероятность наступления события
m=твое число; %число реализация события
k=300; %число испытаний
sum = 0;
sum2 = 0;
load_system('название файла')
for i=1:k
sim('название файла')
n=simout(end,1); %текущее число заходов
sum=sum+n;
sum2=sum2+n^2;
end;
sr=sum/5 %среднее число заходов
sigma2=1/(k-1)*(sum2-1/k*sum^2) %дисперсия числа заходов
% вычисление доверительного интервала
% пусть задана доверительность 0.95
% тогда уровень значимости равен 0.05
% для определения квантиля нормального распределения уровня
% 0.95
% используем обратную функцию нормального распределения norminv(p,m,s)
% где p=(1-0.95)/2; m=0: s=1
t095=norminv(0.025,0,1);
% нижняя граница доверительного интервала
nmin=sr-t095*sqrt(double(sigma2/k))
% верхняя граница доверительного интервала
nmax=sr + t095*sqrt(double(sigma2/k))
% nmin < E[x] < nmax