?? kuopin.m
字號:
sync_dl=[1 0 1 1 0 0 1 1 1 0 1 0 0 1 1 1 1 1 0 0 1 1 0 0 0 0 0 0 0 1 0 1 1 0 1 0 1 0 0 1 1 0 0 0 0 1 1 0 1 0 0 0 1 0 0 0 1 1 1 0 0 1 0 0];
sync_ul=randsrc(1,128);
gp1=zeros(1,32);
gp2=zeros(1,96);
m=1;
%%%%%%%%%%%%%%%%%%%
for i=1:64;
sync_dl(i)= sync_dl(i)+(sync_dl(i)==0)*(-1);
end
%%%%%%%%%%%%%%%%%%%%%%%%
d=[gp1,sync_dl,gp2,sync_ul,gp1];
q=randsrc(1,32000);
k=100;
%%%%%%%%%%%%%
for i=0:4
a=k+i*6400;
for j=1:352;
q(a+j)=d(j);
end
end
%%%%%%%%%%%%%%
for snr=-0.8:0.01:0.8
%%%%%%%%
b=0;
for n=1:500;
q2=awgn(q,snr);
q3=reshape(q2,6400,5);
y=zeros(1,6400);
%%%%%%
for i=1:5
y=y+q3(:,i)';
end
%%%%%%
for i=1:6273;
p1=0;
p2=0;
p3=0;
for j=0:31;%
p1=p1+y(i+j)^2;%
p2=p2+y(i+96+j)^2;%
end%
for t=32:95;
p3=p3+y(i+t)^2;
end
z(i)=(p1+p2)/p3;
end
%%%%%%%%%%%%%
[c,d]=min(z);
if (d==k+1)
b=b+1;
else
b=b;
end
end
p(m)=b/500;
m=m+1;
end
snr=-0.8:0.01:0.8;
plot(snr,p);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -