%rake receiver
fm_rec=fm_ch;%signal on the chnnels
%fm_rec=fm_ch_noised;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[h w]=size(fm_rec);
%chnnel1
cr1_abs=c_abs(:,1);
cr1_abs=hor_copy(cr1_abs,Fs/Rc);
cr1_ang=c_ang(:,1)+1/W;
cr1_ang=hor_copy(cr1_ang,Fs/Rc);
tr1=1/Fs:1/Fs:1/Rc;
tr1=ver_copy(tr1,h);
fm_ch1_s1=1.*cr1_abs.*cos(2*pi*Fc*t1-cr1_ang);%relayted to signal 1
fm_ch1_s2=-1.*cr1_abs.*cos(2*pi*Fc*t1-cr1_ang);%relayted to signal -1
sum_ch1_s1=fm_ch1_s1.*fm_rec;
sum_ch1_s2=fm_ch1_s2.*fm_rec;
sum_ch1_s1=sum(sum_ch1_s1');
sum_ch1_s2=sum(sum_ch1_s2');
%chnnel2
cr2_abs=c_abs(:,2);
cr2_abs=hor_copy(cr2_abs,Fs/Rc);
cr2_ang=c_ang(:,2)+2/W;
cr2_ang=hor_copy(cr2_ang,Fs/Rc);
tr2=1/Fs:1/Fs:1/Rc;
tr2=ver_copy(tr2,h);
fm_ch2_s1=1.*cr2_abs.*cos(2*pi*Fc*t1-cr2_ang);%relayted to signal 1
fm_ch2_s2=-1.*cr2_abs.*cos(2*pi*Fc*t1-cr2_ang);%relayted to signal -1
sum_ch2_s1=fm_ch2_s1.*fm_rec;
sum_ch2_s2=fm_ch2_s2.*fm_rec;
sum_ch2_s1=sum(sum_ch2_s1');
sum_ch2_s2=sum(sum_ch2_s2');
%chnnel3
cr3_abs=c_abs(:,3);
cr3_abs=hor_copy(cr3_abs,Fs/Rc);
cr3_ang=c_ang(:,3)+3/W;
cr3_ang=hor_copy(cr3_ang,Fs/Rc);
tr3=1/Fs:1/Fs:1/Rc;
tr3=ver_copy(tr3,h);
fm_ch3_s1=1.*cr3_abs.*cos(2*pi*Fc*t1-cr3_ang);%relayted to signal 1
fm_ch3_s2=-1.*cr3_abs.*cos(2*pi*Fc*t1-cr3_ang);%relayted to signal -1
sum_ch3_s1=fm_ch3_s1.*fm_rec;
sum_ch3_s2=fm_ch3_s2.*fm_rec;
sum_ch3_s1=sum(sum_ch3_s1');
sum_ch3_s2=sum(sum_ch3_s2');
%sum
Um1=sum_ch1_s1+sum_ch2_s1+sum_ch3_s1;
Um2=sum_ch1_s2+sum_ch2_s2+sum_ch3_s2;
%判决
signal_recive=double(Um1>Um2);
signal_trans=(x_t'+1)./2;
%
x_chip_recive=reshape(signal_recive,15,N);
PN_mesh=hor_copy(PN',N);
PN_mesh_t=hor_copy(1-PN',N);
d1=sum(x_chip_recive~=PN_mesh);
d2=sum(x_chip_recive~=PN_mesh_t);
x_r=double(d1<d2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_bit=sum(x_r~=x);
%E_bit=n_bit/N