基于最小均方误差的-MMSE程序2.docVIP

MMSE检测程序 % m_sequence 程序一 function mseq=m_sequence(fbconnection) n=length(fbconnection); N=2^n-1; register=[zeros(1,n-1) 1]; mseq(1)=register(n); for i=2:N newregister(1)=mod(sum(fbconnection.*register),2); for j=2:n newregister(j)=register(j-1); end register=newregister; mseq(i)=register(n); end % %gold_seq.m (程序二) function goldseq=gold_seq(fbconnection1,fbconnection2) mseq1=m_sequence(fbconnection1); mseq2=m_sequence(fbconnection2); N=2^length(fbconnection1)-1; for shift_amount=0:N-1 shift_mseq2=[mseq2(shift_amount+1:N) mseq2(1:shift_amount)]; goldseq(shift_amount+1,:)=mod(mseq1+shift_mseq2,2); end % mmse_main (程序三) %MMSE解相关多用户检测器和CD传统多用户检测器的误码率比较（八个用户） clear all; snr_indb=1:8; for k=1:length(snr_indb) snr=10^(snr_indb(k)/10); sgma=1; eb=2*(sgma^2)*snr; LC=31; echip=eb/LC; N=10000;%number of bits transmitted %creat PN codes fbconnection=[0 1 0 0 1]; mseq=m_sequence(fbconnection); fbconnection1=[0 0 1 0 1]; fbconnection2=[0 1 1 1 1]; goldseq=gold_seq(fbconnection1,fbconnection2); %N=2^length(fbconnection)-1; ind1=find(mseq==0); mseq(ind1)=-1; ind2=find(goldseq==0); goldseq(ind2)=-1; %creat 31*31 gold sequence temp=goldseq; pn_seq1=temp(3,:); pn_seq2=temp(4,:); pn_seq3=temp(7,:); pn_seq4=temp(10,:); pn_seq5=temp(15,:); pn_seq6=temp(20,:); pn_seq7=temp(26,:); pn_seq8=temp(31,:); % comput matrix R pp=[pn_seq1;pn_seq2;pn_seq3;pn_seq4;pn_seq5;pn_seq6;pn_seq7;pn_seq8]; RR=(1/LC)*pp*pp; RR1=inv(eb*eye(8)); R=inv(RR+RR1); %检测器输出 %generate codes for i=1:N temp1=rand; if(temp10.5),dsource1(i)=-1; else dsource1(i)=1; end; temp2=rand; if(temp20.5),dsource2(i)=-1; else dsource2(i)=1; end; temp3=rand; if(temp30.5),dsource3(i)=-1; else dsource3(i)=1; end; temp4=rand; if(temp40.5),dsource4(i)=-1; else dsource4(i)=1; end; temp5=rand; if(temp50.5),dsource5(i)=-1; else dsource5(i)=1; end; temp6=rand; if(temp60.5),dsource6(i)=-1; else dsource6(i)=1; end; temp7=rand; if(temp70.5),dsource7(i)=-1; else dsource7(i)=1; end; temp8=rand; if(temp80.5),dsou