Поможем написать учебную работу
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
НИЖЕГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ
УНИВЕРСИТЕТ им Р.Е. АЛЕКСЕЕВА
АРЗАМАССКИЙ ПОЛИТЕХНИЧЕСКИЙ ИНСТИТУТ
Кафедра «Авиационные приборы и устройства»
Отчет о выполнении лабораторной работы
«Оценка энтропии ансамбля случайных сообщений»
Вариант 14
Работу выполнил студент:
группы АСП 11-1
Перфилова И.
Проверил преподаватель:
Волков В.Л.
Арзамас 2014
clc;
clear;
F=load('Fx.txt');
var=14;
dt=0.01;
n=3000;
K=15;
t=0:1:2999;
for i=1:n X(i)=F(i,var); end
minX=11; maxX=max(X);
dX=(maxX-minX)/15;
Xe(1)=minX+dX/2; m(1)=0;
for i=2:15 Xe(i)=Xe(i-1)+dX;m(i)=0;end
for i=1:n for j=1:K
if (X(i)>Xe(j)-dX/2)&&(X(i)<=Xe(j)+dX/2) m(j)=m(j)+1; end
end; end
for j=1:K f(j)=m(j)/n; end
for j=1:K P(j)=f(j); end
P=P(P>0);
Hx=-sum(P.*log2(P))
figure(1)
plot(t,X);grid
figure(2)
plot(Xe,f);grid
for i=1:13
DHx(i)=P(i)*(-log(P(i))-Hx)^2;
end
Dhx=sum(DHx)
Dx=(dX^2)/12
Hx =3.2532
Dhx =1.3818
Dx =0.0088
F=load('Fy.txt');
var=14;
dt=0.01;
n=3000;
K=15;
t=0:1:2999;
for i=1:n X(i)=F(i,var); end
minX=min(X); maxX=max(X);
dX=(maxX-minX)/15;
Xe(1)=minX+dX/2; m(1)=0;
for i=2:15 Xe(i)=Xe(i-1)+dX;m(i)=0;end
for i=1:n for j=1:K
if (X(i)>Xe(j)-dX/2)&&(X(i)<=Xe(j)+dX/2) m(j)=m(j)+1; end
end; end
for j=1:K f(j)=m(j)/n; end
for j=1:K P(j)=f(j); end
Hy=-sum(P.*log2(P))
figure(3)
plot(t,X);grid
figure(4)
plot(Xe,f);grid
for i=1:15
DHy(i)=P(i)*(-log(P(i))-Hy)^2;
end
Dhy=sum(DHy)
Dy=(dX^2)/12
Hy =3.9013
Dhy =1.4372
Dy =0.0089
F=load('Fz.txt');
var=14;
dt=0.01;
n=3000;
K=15;
t=0:1:2999;
for i=1:n X(i)=F(i,var); end
minX=min(X); maxX=max(X);
dX=(maxX-minX)/15;
Xe(1)=minX+dX/2; m(1)=0;
for i=2:15 Xe(i)=Xe(i-1)+dX;m(i)=0;end
for i=1:n for j=1:K
if (X(i)>Xe(j)-dX/2)&&(X(i)<=Xe(j)+dX/2) m(j)=m(j)+1; end
end; end
for j=1:K f(j)=m(j)/n; end
for j=1:K P(j)=f(j); end
P=P(P>0);
Hz=-sum(P.*log2(P))
figure(5)
plot(t,X);grid
figure(6)
plot(Xe,f);grid
for i=1:14
DHz(i)=P(i)*(-log(P(i))-Hz)^2;
end
Dhz=sum(DHz)
Dz=(dX^2)/12
Hz =2.1187
Dhz =1.3328
Dz =6.5093e-004
a=[1 2 3];
H=[Hx Hy Hz];
figure(7)
plot(a,H);grid
Контрольный вопрос
Понятие гауссовского процесса, равномерно распределенного и процесса с экспоненциальным распределением.