我写了两天了，现在头晕，不知道怎么错了！请高手帮一下忙这个仿真程序用于实现DS-CDMA通信系统仿真 clear all; clc sr=256000.0; ml=2; br=sr*ml; nd=100; SNR=-5:1:10; irfn=21; IPOINT=8;alfs=0.5;[xh]=hrollfcoef(irfn,IPOINT,sr,alfs,1); [xh2]=hrollfcoef(irfn,IPOINT,sr,alfs,0); user=1; seq=1; stage=3; ptap1=[1 3];ptap2=[2 3]; regi1=[1 1 1];regi2=[1 1 1]; disp('--------------start-------------------');for ebn0=-5:1:10 switch seq case 1 code=mseq(stage,ptap1,regi1,user); case 2 m1=mseq(stage,ptap1,regi1); m2=mseq(stage,ptap2,regi2); code=goldseq(m1,m2,user); case 3 m1=mseq(stage,ptap1,regi1); m2=mseq(stage,ptap2,regi2); code=[goldseq(m1,m2,user),zeros(user,1)]; end codecode=code*2-1; clen=length(code); rfade=0; itau=[0,8]; dlvl1=[0.0,40.0]; n0=[6,7]; th1=[0.0,0.0];itnd1=[3001,4004]; now1=2; tstp=1/sr/IPOINT/clen; fd=160;flat=1; itndel=nd*IPOINT*clen*30; nloop=1000; noe=0; nod=0; for ii=1:nloop data=rand(user,nd*ml)>0.5; [ich, qch]=qpskmod(data,user,nd,ml); [ich1,qch1]=spread(ich,qch,code); [ich2,qch2]=compoversamp2(ich1,qch1,IPOINT); [ich3,qch3]=compconv2(ich2,qch2,xh); if user==1 ich4=ich3; qch4=qch3; else ich4=sum(ich3); qch4=sum(qch3); end if rfade==0 ich5=ich4; qch5=qch4; else [ich5,qch5]=sefade(ich4,qch4,itau,dlvl1,th1,n0,itnd1,… now1,length(ich4),tstp,fd,flat); itnd1itnd1=itnd1+itndel; end spow=sum(rot90(ich3.^2+qch3.^2))/nd; attn=sqrt(0.5*spow*sr/br*10^(-ebn0/10)); [ich6,qch6]=comb2(ich5,qch5,attn); [ich7,qch7]=compconv2(ich6,qch6,xh2); sampl=irfn*IPOINT+1; ich8=ich7(:,sampl:IPOINT:IPOINT*nd*clen+samp1-1); qch8=qch7(:,sampl:IPOINT:IPOINT*nd*clen+sampl-1); [ich9 qch9]=despread(ich8,qch8,code); demodata=qpskdemod(ich9,qch9,user,nd,ml); noe2=sum(sum(abs(data-demodata))); nod2=user*nd*ml; noenoe=noe+noe2; nodnod=nod+nod2; end ber=noe/nod; fprintf('%d/t%d/t%d/t%e/n',ebn0,noe,nod,noe/nod); fid=fopen('BER.dat','a'); fprintf(fid,'%d/t%e/t%f/t%f/t/n',ebn0,noe/nod,noe,nod); fclose(fid); err_rate_final(ebn0+6)=ber; end figure semilogy(SNR,err_rate_final,'b-*'); xlabel('信噪比/dB') ylabel('误码率') axis([-5,10,0,1]) grid on disp('--------------end-------------------');
|
