a.asv

用matlab完成统计信号处理课程实验！

clc
clear
gama=[0.1,0.4,0.7,1,1.5,2.0];
for j=1:5

N=10000;
r(1,N)=0;
d1=0;
d0=0;
for i=1:N
t=0:199;
s=[(0:50)/50,(100-(51:150))/50,((151:199)-200)/50];
sigama=25;
mu=0;
x =mu+ randn(1,200)*sqrt(sigama);
sum_xs=sum(x.*s);
s2=sum(s.^2);
if sum_xs>(0.5*s2)
r(1,i)=1;
d1=d1+1;
else  r(1,i)=0;
        d0=d0+1;
end 
end
pd(1,j)=d1/N;
pm(1,j)=d0/N;
end