直接序列扩频系统matlab仿真(bpsk调制)

《直接序列扩频系统matlab仿真(bpsk调制)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、1、生成m序列及m序列性质实验产生7位01序列，频率100Hz，模拟线性反馈移位寄存器序列，原理图如下：clearall;clc;Xl=0;X2=0;X3=l;m=350;%重复50遍的7位单极性m序列fori=l:mY3=X3;Y2=X2;Y1=X1;X3=Y2;X2=Y1;Xl=xor(Y3,Yl);L(i)=Yl;endfori=l:mM⑴=1-2*L⑴;%将单极性m序列变为双极性m序列endk=l:1:m;figure(1)subplot(3,1,1)%做m序列图stem(k-l,M);axis([0,7,-1,1]);xlabel(’k

2、’)；ylabel(’M序列’);titleC移位寄存器产生的双极性7位M序列’)；subplot(3,1,2)ym=fft(M,4096);magm=abs(ym);%求双极性m序列频潜fm=(1:2048)*200/2048;plot(fm,magm(1:2048)*2/4096);title(’双极性7位M序列的频谱’)axis([90,140,0,0.1]);[a,b]=xcorr(M,’unbiased’)；subplot(3,1,3)%求双极性m序列自相关函数plot(b,a);axis([-20,20,-0.5,1.2]);titl

3、cC双极性7位M序列的自相关函数’)；移位寄存器产生的双极性7位MJT•列10.5宏A°s•0.5•1111«111—«()(>1111双极性7位M關的自招紐数由上图可以看出，7位ni序列为1，-1，-1，-1，1,-1,1。另外，自相关函数的图形比较尖锐，最大值为1,最小值为-1/7,符合理论结果。2、生成50位随机待发送二进制比特序列，并进行扩频编码生成的信息码频率为100/7HZ,利用m序列编码后，频率变为100Hz。N=50;a:0;x_rand=rand(1,N);%产生50个0与1之间随机数fori=l:Nifxrand(i)>=0.

4、5%大于等于0.5的取1，小于0.5的取0x(i)=l;a=a+l;elsex(i)=0;endendt=0:N-l;figure(2)%做信息码图subplot(2,1,1)stcm(t,x);titleC扩频前待发送二进制信息序列’)；tt=0:349;subplot(2,1,2)1=1:7*N;y⑴=0;fori=l:Nk=7*i-6;y(k)=x(i);k=k+l;y(k)=x(i);k=k+l;y(k)=x(i);k=k+l;y(k)=x(i);k=k+l;y(k)=x(i);k=k+l;y(k)=x(i);k=k+l;y(k)=x(i

5、);end“1)=0;fori=l:350%扩频后，码率变为100/7*7=100Hzs(i)=xor(L(i),y(i));endtt=0:7*N-l;stcm(tt,s);axis([0,350,0,1]);title(’扩频后的待发送序列码’)；扩频前待发送二进制信息序列扩频后的待发送序列码3、对扩频前后信号进行BPSK调制，观察其时域波形BPSK调制釆用2kHz信号cos(2*2000*t)作为载波figure(3)subplot(2,1,2)fs=2000;ts=0:0.00001:3.5-0.00001;%为了使信号看起来更光滑，作图

6、时采样频率为100kHz%ps=cos(2*pi*fs*ts);s_b=rectpulse(s,1000);%将冲激信号补成矩形信号s_bpsk=(l-2.*s_b).*cos(2*pi*fs*ts);%扩频后信号BPSK调制时域波形，(卜2.*s_b)是1，-1序列plot(ts,s_bpsk);xlabel(’s’)；axis([O.055,0.085,-1.2,1.2])title(’扩频后bpsk信号吋域波形’)；subplot(2,1,1)sbb=rectpulse(x,7000);sbpskb=(l~2.*sbb).*cos(2*pi

7、*fs*ts);%无扩频信号BPSK调制吋域波形plot(ts,sbpskb);xlabel(’s’)；axis([0.055,0.085,-1.2,1.2]);title(’扩频前bpsk信号吋域波形，)s扩频后bpsk信号时域波形10.50-0.5-1HIIIIIIIII0.0550.060.0650.070.0750.080.085可以看出，100/7Hz的无扩频信号每0.07s吋由于序列极性变换产生相位变4、计算并观察扩频前后BPSK调制信号的频谱对信号采用400000点fft计算，得到频谱figure(4)N=400000;ybb=ff

8、t(s.bpskb,N);%无扩频信号BPSK调制频谱magb=abs(ybb);fbb=(l:N/2)*100000/N;subplo