matlab直仿真代码.docxVIP

1.在AWGN信道下对直扩通信系统进行性能仿真研究%D??￠???μ?êéè?a100kHz￡?2é?ù?μ?êéè?a40MHz￡??±???ú???μ?êéè?a5MHzcloseallclearall%D??′?￡?éSNR=[1:1:15]Rb=0.25; %???a?ù?ê0.25Mb/sTb=1/Rb; %???a????M=100; %D??￠????êyL=400; %???????a?úμ?2é?ùμ?êyN=M*L; %×ü2é?ùμ?êydt=Tb/L; %ê±óò2é?ù????df=1/(N*dt); %?μóò2é?ù????Bs=N*df/2; %?μí3′??íT=N*dt; %???ìê±??t=linspace(-T/2,T/2,N); %ê±óòoá×?±êx=sign(rand(1,M)-0.5); %2úéú????D?2?1éá???t0=linspace(0,T,M);fori=1:100;d((1+(i-1)*400):i*400)=x(i); %????D??￠???a?úo?400??2é?ùμ?endfigure(1)plot(t,d)xlabel(t(ms));ylabel(s(t)(V));title(D??′);axis([0,200,-1.5,1.5]);gridon%êy?Yμ÷???￡?éfs=4000; %??2¨?μ?êkHzz=reshape(d,2,N/2); %′?2￠±???Q=z(1,:); %sμ???êy??,?′Q(t)I=z(2,:); %sμ???êy??,?′I(t)tq=linspace(-T/2,T/2,length(I));QPSK=(I.*cos(2*fs*pi*tq)-Q.*sin(2*fs*pi*tq))/sqrt(2); %QPSKμ÷??figure(2)plot(tq,QPSK)xlabel(t(ms));ylabel(QPSK(t)(V));title(QPSKμ÷??D?o?);axis([0,200,-1.5,1.5]);gridon%à??μ??2úéú?￡?é%à??μ??3??μéè?¨user=1; %ó??§êyseq=1; %1￡omDòáD 2￡oGoldDòáD 3￡o?y??GoldDòáDstage=3; %DòáD?×êyptap1=[1 3]; %μúò?????D?ò?????′??÷μ??μêy￡o[1 3]ptap2=[2 3]; %μú?t????D?ò?????′??÷μ??μêy￡o[2 3]regi1=[1 1 1]; %μúò?????D?ò?????′??÷μ?3?ê??ˉ￡o[1 1 1]regi2=[1 1 1]; %μú?t????D?ò?????′??÷μ?3?ê??ˉ￡o[1 1 1]switchseqcase 1 %mDòáDcode=mseq(stage,ptap1,regi1,user); case 2 %GoldDòáD m1= mseq(stage,ptap1,regi1); m2= mseq(stage,ptap2,regi2); code=goldseq(m1,m2,user); case 3 %?y??GoldDòáD m1=mseq(stage,ptap1,regi1); m2=mseq(stage,ptap2,regi2); code=[goldseq(m1,m2,user),zeros(user,1)]; endxxx=codeclen=length(code);xlength=100*50; %????D??￠μ÷?????úo?50???±??x_code=zeros(1,xlength);fori=1:xlength%°?0,1DòáD??±????a-1,1μ÷????if mod(i,7)~=0 %mod(a,b)?íê??óμ?ê?a3yò?bμ?óàêyx_code(i)=sign(code(mod(i,7))-0.5); elseif mod(i,7)==0x_code(i)=sign(code(7)-0.5);endendk_code=zeros(1,length(QPSK));fori=1:5000k_code((1+(i-1)*4):i*4)=x_code(i); %?????±?????a?úo?4??2é?ùμ?endfigure(3)plot(tq,k_code)xlabel(t(ms));ylabel(k_code(t)(V));title(à??μ??);axis([20,25,-1.5,1.