资源描述：
《基于最小均方误差的-mmse程序2》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、MMSE检测程序%m_sequence程序一functionmseq=m_sequence(fbconnection)n=length(fbconnection);N=2^n-1;register=[zeros(1,n-1)1];mseq(1)=register(n);fori=2:Nnewregister(1)=mod(sum(fbconnection.*register),2);forj=2:nnewregister(j)=register(j-1);endregister=newregister;mseq(i)
2、=register(n);end%%gold_seq.m(程序二)functiongoldseq=gold_seq(fbconnection1,fbconnection2)mseq1=m_sequence(fbconnection1);mseq2=m_sequence(fbconnection2);N=2^length(fbconnection1)-1;forshift_amount=0:N-1shift_mseq2=[mseq2(shift_amount+1:N)mseq2(1:shift_amount)];gol
3、dseq(shift_amount+1,:)=mod(mseq1+shift_mseq2,2);end%mmse_main(程序三)%MMSE解相关多用户检测器和CD传统多用户检测器的误码率比较（八个用户）clearall;snr_indb=1:8;fork=1:length(snr_indb)snr=10^(snr_indb(k)/10);sgma=1;eb=2*(sgma^2)*snr;LC=31;echip=eb/LC;N=10000;%numberofbitstransmitted%creatPNcodesf
4、bconnection=[01001];mseq=m_sequence(fbconnection);fbconnection1=[00101];fbconnection2=[01111];goldseq=gold_seq(fbconnection1,fbconnection2);%N=2^length(fbconnection)-1;ind1=find(mseq==0);mseq(ind1)=-1;ind2=find(goldseq==0);goldseq(ind2)=-1;%creat31*31goldsequen
5、cetemp=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,:);%computmatrixRpp=[pn_seq1;pn_seq2;pn_seq3;pn_seq4;pn_seq5;pn_seq6;pn_seq7;pn_seq8];RR=(1/LC)*pp*p
6、p';RR1=inv(eb*eye(8));R=inv(RR+RR1);%检测器输出%generatecodesfori=1:Ntemp1=rand;if(temp1<0.5),dsource1(i)=-1;elsedsource1(i)=1;end;temp2=rand;if(temp2<0.5),dsource2(i)=-1;elsedsource2(i)=1;end;temp3=rand;if(temp3<0.5),dsource3(i)=-1;elsedsource3(i)=1;end;temp4=rand;
7、if(temp4<0.5),dsource4(i)=-1;elsedsource4(i)=1;end;temp5=rand;if(temp5<0.5),dsource5(i)=-1;elsedsource5(i)=1;end;temp6=rand;if(temp6<0.5),dsource6(i)=-1;elsedsource6(i)=1;end;temp7=rand;if(temp7<0.5),dsource7(i)=-1;elsedsource7(i)=1;end;temp8=rand;if(temp8<0.5)
8、,dsource8(i)=-1;elsedsource8(i)=1;end;end;number_of_err=0;number_of_err1=0;number_of_err2=0;fori=1:N%将每个bitrepeatLC=31timesforj=1:LCrepeatdata1(j)=dsource1(i);repeatdata2(j)
