Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Предоплата всего
Подписываем
МИНИСТЕРСТВО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ
НИЖЕГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
Кафедра «Авиационные приборы и устройства»
Практическая работа
Моделирование процессов с заданным законом распределения.
Вариант №7
Работу выполнил
Студенты группы АСП 11-1
Кошелев С.И.
Работу проверил
Преподаватель
Волков В.Л.
г. Арзамас, 2013г.
Исходные данные:
|
a |
b |
M |
||
|
4 |
5.5 |
2.5 |
1.2 |
6 |
Ход работы:
1.1.Моделирование процесса с равномерным распределением
clc
clear
a=1
b=2
M=0.7
sigm=0.2
N=10
n=200
L=1.5e10-3
% построение графика процесса для равномерно распределенной величины
q=rand(200,1);
E=a+q*(b-a);
t=[1:1:200];
figure(1)
plot(t,E), grid
%расчет гистограммы
dE=(b-a)/N
Ei(1)=a-dE
for i=1:N+2 m(i)=0; Ei(i+1)=Ei(i)+dE
for j=1:n if E(j)>Ei(i)&& E(j)<Ei(i+1)
m(i)=m(i)+1
x(1)=Ei(1)
x(i)=Ei(i)
f(i)=m(i)/n
end;
end;end
figure(2)
plot(x,f),grid
1.2.Процесс с нормальным законом распределения
clc
clear
a=0
b=1
M=0.5
sigm=0.1
N=10
n=200
L=1e10-3
k=1;
%Построение график процесса для нормально распределенной величины
while k<n+1
e(k)=0; Q=rand(12,1); for j=1:12 e(k)=e(k)+Q(j); end; e(k)=e(k)-6; k=k+1;end;
y=e*sigm+M;
t=[1:1:200]'; figure(3);
plot(t,y), grid
%Построение гистограммы для нормально распределенной величины
min=M-5*sigm;
max=M+5*sigm;
dy=(max-min)/N;
yi(1)=min
for i=1:N+2 m(i)=0; yi(i+1)=yi(i)+dy
for j=1:n if y(j)>yi(i)&& y(j)<yi(i+1)
f(i)=m(i)/n
m(i)=m(i)+1
p(i)=yi(i)
end;
end;end
figure(4)
plot(p,f),grid
1.3. Процесс с экспоненциальным законом распределением.
clc
clear
a=0
b=1
M=0.5
sigm=0.1
N=10
n=200
L=1e10-3
%Построение графика процесса для экспоненциально распределенной величины
q=rand(200,1);
x=-log(q)/L
t=[1:1:200]'; figure(5);
plot(t,x), grid
%Построение гистограммы для экспоненциально распределенной величины
dx=10*(1/L)/N
xi(1)=0
for i=1:N+2 m(i)=0; xi(i+1)=xi(i)+dx
for j=1:n if x(j)>xi(i)&& x(j)<xi(i+1)
f(i)=m(i)/n
m(i)=m(i)+1
p(i)=xi(i)
end;
end;end
figure(6)
plot(p,f),grid