b.m

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

clc
clear
gama=[0.1,0.4,0.5,0.7,1,4/3,1.5,2.0,3.5,10];
for k=1:5
for j=1:9

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)+25*gama(1,j))
r(1,i)=1;
d1=d1+1;
else  r(1,i)=0;
        d0=d0+1;
end 
end
pd(k,j)=d1/N;
pm(k,j)=d0/N;
end
end