b.asv

来自「用matlab完成统计信号处理课程实验！」· ASV 代码 · 共 23 行

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=d1/N;
pm=d0/N;