直接序列擴(kuò)頻通信系統(tǒng)仿真程序

1、直接序列擴(kuò)頻通信系統(tǒng)仿真程序楊晶超 S20070416 code_le ngth=20;%信息碼元個(gè)數(shù) N=1:code_le ngth; ran d(seed,0); x=sig n(ran d(1,code_le ngth)-0.5);% 信息碼 for i=1:20s(1+(i-1)*800):i*800)=x(i);%每個(gè)信息碼元內(nèi)含fs/f=800 個(gè)采樣點(diǎn)end圖1信源信息碼 %產(chǎn)生偽隨機(jī)碼，調(diào)用的mgen函數(shù)見附錄 length=100*20;%偽碼頻率5MHz,每個(gè)信息碼內(nèi)含 5MHz/50kHz=100個(gè)偽碼 x_code=sign(mgen(19,8,length)-0.5

2、);%把 0,1 序列碼變換為-1,1 調(diào)制碼 for i=1:2000w_code(1+(i-1)*8):i*8)=x_code(i);%每個(gè)偽碼碼元內(nèi)含8個(gè)采樣點(diǎn)end生成的PN碼波形如圖2所示。%k_code為擴(kuò)頻碼 %擴(kuò)頻 k_code=s.*w_code;擴(kuò)頻碼如圖3所示。PF網(wǎng)圖2 PN碼圖3擴(kuò)頻碼 %調(diào)制 fs=20e6; f0=30e6; for i=1:2000AI=2;dt=fs/fO;n=0:dt/7:dt;% 一個(gè)載波周期內(nèi)采樣八個(gè)點(diǎn)cl=AI*cos(2*pi*fO* n/fs);sig nal(1+(i-1)*8):i*8)=k_code(1+(i-1)*8):i*

3、8).*cl;endPSK調(diào)制后的波形如圖4所示。pskVfl制陽(yáng)豹激飛圖4 PSK調(diào)制后的波形 %解調(diào) Al=1; dt=fs/fO; n=0:dt/7:dt;% 一個(gè)載波周期內(nèi)采樣八個(gè)點(diǎn) cl=AI*cos(2*pi*fO* n/fs); for i=1:2000sig nal_h(1+(i-1)*8):i*8)=sig nal(1+(i-1)*8):i*8).*cl; end解調(diào)后的波形如圖5所示。圖5解調(diào)后的波形 %解擴(kuò)jk_code=sig nal_h.*w_code; %低通濾波 wn=5/10000000;%截止頻率 wn=fn/(fs/2),這里的fn為信息碼(擴(kuò)頻碼)的帶寬5M b=fir1(16,w n); H=freqz(b,1,16000); sig nal_d=filter(b,1, jk_code);解擴(kuò)并濾波后的波形如圖 6所示。解護(hù)并譙彼后的彼形圖6解擴(kuò)并濾波后的波形7所示。圖7輸入與輸岀對(duì)比附錄：%mgen.m functionout=mgen(g,state,N) gen=dec2bin(g)-48;M=length(gen);curState=dec2bin(state,M-1)-48;for k=1:Nout(k)=curState(M-1);a=rem