%本程序用于對(duì)單用戶檢測(cè)，傳統(tǒng)多用戶檢測(cè)，解相關(guān)多用戶檢測(cè)
%的性能進(jìn)行比較
%本程序只考慮到噪聲和干擾的抵消，不針對(duì)編解碼的問(wèn)題
%%本例針對(duì)同步cdma系統(tǒng)
%僅考慮了遠(yuǎn)近效應(yīng)和功率控制，信道考慮了AWGN和正弦干擾的影響
%異步和多徑?jīng)]有考慮
%擴(kuò)頻碼是長(zhǎng)度為31的gold序列
  echo on
  Lc=31;%用于對(duì)符號(hào)擴(kuò)頻的長(zhǎng)度
  A2=7;%正弦干擾信號(hào)的幅度
  w0=1;

  %
  %gold 序列，參考<<MATLAB與通信仿真>>
  %也可以考慮使用一個(gè)函數(shù)來(lái)產(chǎn)生
  % PN1=[0  1  1  1 0 0 0 0  1 0 0 0 0  1  1 0 0  1 0 0  1 0  1  1  1  1 0 0 0 0 0]   %sum()=7
   PN1=[1 -1 -1 -1 1 1 1 1 -1 1 1 1 1 -1 -1 1 1 -1 1 1 -1 1 -1 -1 -1 -1 1 1 1 1 1];
  % PN2=[0  1 0 0 0  1  1  1 0  1 0  1 0  1  1  1 0  1  1 0 0 0 0 0  1  1 0 0  1  1  1]%sum()=-1
   PN2=[1 -1 1 1 1 -1 -1 -1 1 -1 1 -1 1 -1 -1 -1 1 -1 -1 1 1 1 1 1 -1 -1 1 1 -1 -1 -1];
  % PN3=[ 1  1  1 0 0  1 0  1 0 0 0 0 0 0 0  1 0  1  1  1  1  1 0 0 0 0 0  1 0 0 0 ]  %sum()=7
   PN3=[-1 -1 -1 1 1 -1 1 -1 1 1 1 1 1 1 1 -1 1 -1 -1 -1 -1 -1 1 1 1 1 1 -1 1 1 1];
  % PN4=[0 0  1  1 0  1 0 0 0 0  1 0  1 0  1 0 0  1  1  1 0 0  1 0 0  1  1  1  1  1  1]%sum()=-1
   PN4=[1 1 -1 -1 1 -1 1 1 1 1 -1 1 -1 1 -1 1 1 -1 -1 -1 1 1 -1 1 1 -1 -1 -1 -1 -1 -1];
  
  %
  %correlatting index
  R12=sum(PN1.*PN2)/Lc;
  R13=sum(PN1.*PN3)/Lc;
  R14=sum(PN1.*PN4)/Lc;
  R23=sum(PN2.*PN3)/Lc;
  R24=sum(PN2.*PN4)/Lc;
  R34=sum(PN3.*PN4)/Lc;
  %
  %  correlating matrix
  R=[1 R12 R13 R14;R12 1 R23 R24;R13 R23 1 R34;R14 R24 R34 1];
  %
  %the inverse matrix of R
  %S=inv(R);
  %
  %在不同的信噪比條件下誤碼測(cè)試
  snrindb=-4:2:8;
  for i=1:length(snrindb),
    
   [err_prbmul(i),err_prbmuldecorr(i),err_mmse(i)]=ss_pemmse(snrindb(i),R,PN1,PN2,PN3,PN4);
  
  end;
  %
  %畫(huà)誤碼率曲線
  
  semilogy(snrindb,err_prbmul,'-b');%傳統(tǒng)多用戶的誤碼率曲線
  hold on
  semilogy(snrindb,err_prbmuldecorr,'-r');%解相關(guān)多用戶的誤碼率曲線
  hold on
  semilogy(snrindb,err_mmse,'-g');%mmse多用戶的誤碼率曲線