扩频MATLAB.docx

clcclearCodeNum=1000;a=ones(1,10);fback=xor(xor(a(10),a(9)),xor(xor(a(6),a(5)),xor(a(4),a(2))));for?i=1:CodeNumshifttemp(i,1:10)=a;b(i)=a(10);a=[fback?a(1,1:9)];fback=?xor(xor(a(10),a(9)),xor(xor(a(6),a(5)),xor(a(4),(2))));end?for?i=1:1:CodeNumb(1,i)=1-2*b(1,i);endfor?i=1:1:CodeNumsampleb(1,(i-1)*320+1:i*320)=repmat(b(1,i),1,320);endd=[sampleb?sampleb];??N=2;Fd=randint(1,N);for?i=1:1:Nsampleb(1,(i-1)*320000+1:i*320000)=repmat(b(1,i),1,320000);endc=d.*sampleb;%figure(1);%plot(c);?%调制fs=160e6;t=0:1/fs:(4e-3)-1/fs;y1=1*sin(2*pi*5e6*t);BPSKsignal=y1.*c;figure(2);plot(BPSKsignal);figure(3);plot(abs(fft(BPSKsignal)));?%上下变频y2=1*sin(2*pi*15e6*t);figure(4);plot(abs(fft(y2)));y3=BPSKsignal.*y2;figure(5);plot(abs(fft(y3)));Wn=0.25;filterparemater=fir1(64,Wn,high);filterResult=fftfilt(filterparemater,y3);figure(6);plot(filterResult);?y4=2*y3;Wn2=[0.24375,0.25625];filterparemater2=fir1(63,Wn2,band);filterResult2=fftfilt(filterparemater2,y4);y2=1*sin(2*pi*15e6*t);y5=y4.*y2;figure(7);plot(abs(fft(y5)));?Wn3=0.0625;filterparemater3=fir1(31,Wn3,low);filterResult3=fftfilt(filterparemater3,y5);figure(8);plot(filterResult3);?%解调DemodulatSignal=filterResult3.*y1;wn2=5e6/80e6;f1=fir1(31,wn2);demodulatPcode=fftfilt(f1,DemodulatSignal);figure(9);plot(demodulatPcode);demodulatPcode1=demodulatPcode/abs(demodulatPcode)；figure(10);plot(demodulatPcode);subplot(2,2,1)plot(c);axis([160000,170000,-1.1,1.1])subplot(2,2,2)plot(BPSKsignal);subplot(2,2,3)plot(fs,abs(fft(BPSKsignal)));subplot(2,2,4)plot(demodulatPcode1);axis([160000,170000,-1.1,1.1])for?k=1:1:32e4/32?RecCode(k)=demodulatPcode(k*32+32e4);?%滑动相关捕获??endLocCode=rectpulse(b,10);CodeNum=k;RecCode=[RecCode(1,9001:10000)?RecCode(1,1:9000)];%RecCode接收扩频码?LocCode本地扩频码CorResultLoc(1)=sum(LocCode.*RecCode)/CodeNum;CorResultRec(1)=sum(LocCode.*RecCode)/CodeNum;for?i=2:1:CodeNum?BackCode=[LocCode(1,i:CodeNum)?LocCode(1,1:i-1)];?CorResultLoc(i)=sum(RecCode.*BackCode)/CodeNum;??FrontCode=[LocCode(1,i+1:CodeNum)?LocCode(1,1:i)];%前后码片差1/10码周期?CorResult