【導讀】二零一零年五月。目錄。1

　　

【正文】 fs=20e6。 f0=30e6。 AI=1。 dt=fs/f0。 n=0:dt/7:dt。 %一個載波周期內(nèi)采樣八個點 cI=AI*cos(2*pi*f0*n/fs)。 for i=1:3000 signal_h((1+(i1)*8):i*8)=signal((1+(i1)*8):i*8).*cI。 end %解擴 jk_code=signal_h.*w_code。 %低通濾波 (代通采樣 ) wn=1/10。 %截至頻率 wn=fn/(fs/2),這里的 fn 為信息碼 (擴頻碼 )的帶寬 1M b=fir1(16,wn)。 H=freqz(b,1,16000)。 signal_d=filter(b,1, jk_code)。 figure(5) subplot(2,1,1)。 plot(signal_d) %解調(diào)信號圖形 title(39。解調(diào)后信息碼 39。)。 xlabel(39。t39。)。 subplot(2,1,2)。 plot(s)。 37 title(39。信源信息碼 39。)。 xlabel(39。t39。)。 end 3. 加噪聲后信號未經(jīng)過調(diào)制的仿真程序 code_length=100。 %信息碼元個數(shù) N=1: code_length。 rand(39。seed39。,0)。 x=sign(rand(1, code_length))。 %信息碼 for i=1:100 s((1+(i1)*800):i*800)=x(i)。 %每個信息碼元內(nèi)含 fs/f=800 個采樣點 end %產(chǎn)生偽隨機碼 length=100*100。 %偽碼頻率 1MHz,每個信息碼內(nèi)含 1MHz/10kHz=100 個偽碼 x_code=sign(mgen (19,8,length))。 %偽隨機碼 for i=1:10000 w_code((1+(i1)*8):i*8)=x_code(i)。 %每個偽碼碼元內(nèi)含 8 個采樣點 end %擴頻 k_code=s.*w_code。 %k_code 為擴頻碼 %加入噪聲的程序如下： gs=10e6。 g0=30e6。 dt=gs/g0。 m=0:dt/10:dt*(80001/10)。 cIl=cos(2*pi*g0*m/gs)。 noise=cIl.*(1/4+wgn(1,80000,3))。 %噪聲干擾 tran= k_code +noise。 figure(6) plot(tran)。%加噪聲的擴頻碼信號 grid on。 title(39。帶有噪聲的未調(diào)制信號 39。)。 xlabel(39。t39。)。 %解擴 jk_code= tran.*w_code。 %低通濾波 (代通采樣 ) wn=1/10。 %截至頻率 wn=fn/(fs/2),這里的 fn 為信息碼 (擴頻碼 )的帶寬 1M b=fir1(16,wn)。 H=freqz(b,1,16000)。 signal_d=filter(b,1, jk_code)。 figure(7) subplot(2,1,1)。 38 plot(signal_d)。 %加有噪聲信號的解擴信號 grid on。 title(39。帶有噪聲的未調(diào)制信號解擴后的信息碼 39。)。 xlabel(39。t39。)。 subplot(2,1,2)。 plot(s)。 grid on。 title(39。信源信息碼 39。)。 xlabel(39。t39。)。 4. 加入噪聲后信號經(jīng)過調(diào)制的仿真程序 %產(chǎn)生信號 code_length=100。 %信息碼元個數(shù) N=1: code_length。 rand(39。seed39。,0)。 x=sign(rand(1, code_length))。 %信息碼 for i=1:100 p((1+(i1)*800):i*800)=x(i)。 %每個信息碼元內(nèi)含 fs/f=800 個采樣點 end %產(chǎn)生偽隨機碼 length=100*100。%偽碼頻率 1MHz,每個信息碼內(nèi)含 1MHz/10kHz=100 個偽碼 x_code=sign(mgen (19,8,length))。 %偽隨機碼 for i=1:10000 w_code((1+(i1)*8):i*8)=x_code(i)。 %每個偽碼碼元內(nèi)含 8 個采樣點 end %擴頻 k_code=p.*w_code。 %k_code 為擴頻碼 %調(diào)制 fs=20e6。 f0=30e6。 for i=1:10000 AI=2。 dt=fs/f0。 n=0:dt/7:dt。 %一個載波周期內(nèi)采樣八個點 cI=AI*cos(2*pi*f0*n/fs)。 signal((1+(i1)*8):i*8)=k_code((1+(i1)*8):i*8).*cI。 end %加入噪聲的程序如下： gs=10e6。 g0=30e6。 dt=gs/g0。 m=0:dt/10:dt*(80001/10)。 39 cIl=cos(2*pi*g0*m/gs)。 noise=cIl.*(1/4+wgn(1,80000,3))。 %噪聲干擾 tran=signal+noise。 figure(8) plot(tran) %帶有噪聲的調(diào)制信號波形 grid on。 title(39。帶有噪聲的調(diào)制信號波形 39。)。 xlabel(39。t39。)。 %解調(diào) fs=20e6。 f0=30e6。 AI=1。 dt=fs/f0。 n=0:dt/7:dt。 %一個載波周期內(nèi)采樣八個點 cI=AI*cos(2*pi*f0*n/fs)。 for i=1:10000 signal_h((1+(i1)*8):i*8)= tran ((1+(i1)*8):i*8).*cI。 figure(9) plot(signal_h)%帶有噪聲的解調(diào)信號波形 grid on。 title(39。帶有噪聲的解調(diào)信號波形 39。)。 xlabel(39。t39。)。 end %解擴 jk_code=signal_h.*w_code。 %低通濾波 (代通采樣 ) wn=1/10。 %截至頻率 wn=fn/(fs/2),這里的 fn 為信息碼 (擴頻碼 )的帶寬 1M b=fir1(16,wn)。 H=freqz(b,1,16000)。 signal_d=filter(b,1, jk_code)。 figure(10) subplot(2,1,1)。 plot(signal_d)。 %帶有噪聲信號解擴信號 grid on。 title(39。帶有噪聲的信號 — 解擴后信息碼 39。)。 xlabel(39。t39。)。 subplot(2,1,2)。 plot(p)。 grid on。 title(39。信源信 息碼 39。)。 xlabel(39。t39。)。 40 BPSK 信號加噪聲誤碼率程序如下 %產(chǎn)生信號 code_length =100。 %信息碼元個數(shù) N=1: code_length。 rand(39。seed39。,0)。 s= randint (1, code_length)。 g=sign(s )。 %信息碼 max_loop=1*10^2。 x=[0:10]。 for ebn0=0:10 totalerror=0。 for loop=1:max_loop % get Eb eb=sum(g.^2)/ code_length。 %% getting noise power noise_power=sqrt(*eb*10.^(ebn0/10))。 %creating AWGN noise=randn(1, code_length)*noise_power。 transignal=g+noise。 % 發(fā)送信號 demosignal=transignal0。 noe=sum(abs(sdemosignal))。 totalerror=totalerror+noe。 end ber1(ebn0+1)=totalerror/( code_length *max_loop)。 pb1(ebn0+1)=*erfc(sqrt(10.^(ebn0/10)))。 end figure。 %print ideal graph。 semilogy(x,pb1,39。r39。)。 grid on。 hold on。 semilogy(x,ber1,39。b39。)。 title(39。信號加噪聲誤碼率 39。)。 legend(39。理論誤碼率 39。,39。實際誤碼率 39。)。 xlabel(39。eb/n0(dB)39。)。 ylabel(39。Pb39。)。 41 Lc=100。%碼元抽樣次數(shù) max_loop=1*10。 for snr_in_dB=0:10 totalerror=0。 for loop=1:max_loop snr=10^(snr_in_dB/10)。 sgma=1。 %撒函數(shù) . Eb=2*sgma^2*snr。 %碼元能量 E_chip=Eb/Lc。 % 每個碼片的能量 N=1000。 % 信息碼個數(shù) num_of_err=0。 %產(chǎn)生信息碼 for i=1:N, temp=rand。 if (temp), data=1。 else data=1。 end。 %每個信息碼元內(nèi)含 Lc個采樣點 for j=1:Lc, repeated_data(j)=data。 end。 %產(chǎn)生 PN碼 pn_seq=sign(mgen (19,8,Lc))。 %偽隨機碼 %擴頻 trans_sig=sqrt(E_chip)*repeated_data.*pn_seq。 %產(chǎn)出噪聲 noise=sgma*randn(1,Lc)。 % 接受信號 rec_sig=trans_sig+noise。 % 解擴 temp=rec_sig.*pn_seq。 decision_variable=sum(temp)。 % 判決 if (decision_variable0), decision=1。 else decision=1。 end。 % 判決信號和原信號對比 42 if (decision~=data), num_of_err=num_of_err+1。 end。 end。 totalerror=totalerror+num_of_err。 end。 ber1(snr_in_dB+1)=totalerror/(N*max_loop)。 %實際誤碼率 pb1(snr_in_dB+1)=*erfc(sqrt(10.^(snr_in_dB/10)))。 %理論誤碼率 end x=[0:10]。 figure。 %畫圖 。 semilogy(x,pb1,39。r39。)。 grid on。 hold on。 semilogy(x,ber1,39。b39。)。 title(39。信號 擴頻后加噪聲誤碼率 39。)。 legend(39。理論誤碼率 39。,39。實際誤碼率 39。)。 xlabel(39。snr_in_dB39。)。 ylabel(39。Pb39。)。 43 主要參考文獻 [1] 曾興 文 .《擴頻通信及多址技術》 西安電子科技大學出版社 。 [2] 曹麗娜、樊昌信 .《通信原理》 (第六版 ) 國防工業(yè)出版社 。 [3] 查光明、熊賢祚．《擴頻通信》 西安電子科技大學出版社 。 [4] 韋惠民 .《 擴頻通信技術及應用 》 西安電子科技大學出版社 。 [5] 樊昌信 . 《通信原理教程》 電 子工業(yè)出版社。 [6] 求是科技 .《 MATLAB 》 人民郵電 出版社 。 [7]