m序列扩频

下面是直扩系统仿真程序，第一段是扩频和调制，第二段是m序列生成，求第一段程序接着编写解扩和解调的程序~（1）clcclear%.............................生成待传输信息，码元宽度为切普宽度的10倍...................message=randint(1,4);%随机生成4个01码code=[];fori=1:length(message)ifmessage(1,i)==0code1=zeros(1,10);elsecode1=ones(1,10);endcode=[code,code1];end%...............................生成m序列......................................a=[0100111];%m序列特征多项式b=[0000001];%移位寄存器初始状态c=length(code);%生成m序列长度m_sequence=mseq(a,b,c);dsss=code.*m_sequence;%扩频%.............................BPSK调制.......................................f=20;%载波频率，为码元速率二倍t=0:1/(f-1):1;cp=[];mod=[];forj=1:length(dsss)ifdsss(j)==0cp1=-ones(1,20);elsecp1=ones(1,20);endc=cos(2*pi*f*t);cp=[cp,cp1];mod=[mod,c];endbpsk=cp.*mod;plot(bpsk,'linewidth',1.5);gridonylabel('bpskmodulation');axis([020*length(dsss)-22]);（2）function[seq]=mseq(connections,registers,len)m=length(connections);L=2^m-1;%周期长度iflen==0len=L;endfan=0;fori=1:lenseq(i)=registers(m);forj=1:mfan=fan+connections(j)*registers(j);fan=(mod(fan,2));endfort=m:-1:2%寄存器移位registers(t)=registers(t-1);endregisters(1)=fan;fan=0;end