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

1、直接序列擴(kuò)頻通信系統(tǒng)仿真程序楊晶超 s20070416>> code_length=20; %信息碼元個(gè)數(shù)>> N=1:code_length;>> rand('seed',0);>> x=sign(rand(1,code_length)-0.5); %信息碼>> for i=1:20 s(1+(i-1)*800):i*800)=x(i); %每個(gè)信息碼元內(nèi)含fs/f=800個(gè)采樣點(diǎn) end生成的信息碼的波形圖如圖1所示。圖1 信源信息碼>> %產(chǎn)生偽隨機(jī)碼，調(diào)用的mgen函數(shù)見(jiàn)附錄>> len

2、gth=100*20; %偽碼頻率5MHz,每個(gè)信息碼內(nèi)含5MHz/50kHz=100個(gè)偽碼>> x_code=sign(mgen(19,8,length)-0.5); %把0,1序列碼變換為-1,1調(diào)制碼>> for i=1:2000 w_code(1+(i-1)*8):i*8)=x_code(i); %每個(gè)偽碼碼元內(nèi)含8個(gè)采樣點(diǎn) end生成的PN碼波形如圖2所示。>> %擴(kuò)頻>> k_code=s.*w_code; %k_code為擴(kuò)頻碼擴(kuò)頻碼如圖3所示。圖2 PN碼圖3 擴(kuò)頻碼>> %調(diào)制>> fs=20e6;>

3、;> f0=30e6;>> for i=1:2000 AI=2; dt=fs/f0; n=0:dt/7:dt; %一個(gè)載波周期內(nèi)采樣八個(gè)點(diǎn) cI=AI*cos(2*pi*f0*n/fs); signal(1+(i-1)*8):i*8)=k_code(1+(i-1)*8):i*8).*cI; endPSK調(diào)制后的波形如圖4所示。圖4 PSK調(diào)制后的波形>> %解調(diào)>> AI=1;>> dt=fs/f0;>> n=0:dt/7:dt; %一個(gè)載波周期內(nèi)采樣八個(gè)點(diǎn)>> cI=AI*cos(2*pi*f0*n/fs);>

4、;> for i=1:2000 signal_h(1+(i-1)*8):i*8)=signal(1+(i-1)*8):i*8).*cI; end解調(diào)后的波形如圖5所示。圖5 解調(diào)后的波形>> %解擴(kuò)>>jk_code=signal_h.*w_code;>> %低通濾波>> wn=5/10000000; %截止頻率wn=fn/(fs/2),這里的fn為信息碼(擴(kuò)頻碼)的帶寬5M>> b=fir1(16,wn);>> H=freqz(b,1,16000);>> signal_d=filter(b,1, jk_code);解擴(kuò)并濾波后的波形如圖6所示。圖6 解擴(kuò)并濾波后的波形從圖形整體看，解擴(kuò)出來(lái)的信息碼與信源信息碼基本相同。輸入與輸出對(duì)比如圖7所示。圖7 輸入與輸出對(duì)比附錄：%mgen.mfunctionout=mgen(g,state,N)gen=dec2bin(g)-48;M=length(gen);curState=dec2bin(state,M-1)-48;for k=1:N out(k)=curStat