【文章內(nèi)容簡(jiǎn)介】 smld_err_prb2(i)=wubitlv(SNRindB(i),Lc,A2,W0)。 smld_err_prb3(i)=wubitlv(SNRindB(i),Lc,A3,W0)。echo offend。%計(jì)算在無(wú)正弦干擾信號(hào)情況下的誤碼率echo onSNRindB4=0:1:8。for i=1:length(SNRindB4)mld_err_prb4(i)=wubitlv(SNRindB(i),Lc,A4,W0)。echo offend。%繪制仿真結(jié)果曲線，誤碼率用對(duì)數(shù)表示semilogy(SNRindB,smld_err_prb1,SNRindB,smld_err_prb2,SNRindB,smld_err_prb3,SNRindB4,smld_err_prb4)。xlabel(39。信噪比Eb/N039。)。ylabel(39。誤碼率Pb39。)。title(39。DS系統(tǒng)在不同振幅的正弦干擾下的誤碼率曲線39。)子函數(shù)：%計(jì)算誤碼率的m文件function [p]=wubitlv(snr_in_dB,Lc,Ac,W0)snr=10^(snr_in_dB/10)。sgma=1。Eb=2*sgma^2*snr。 %每位信息比特的能量E_chip=Eb/Lc。 %每個(gè)碼片的能量N=10000。 %待傳輸?shù)男畔⒎?hào)數(shù)num_of_err=0。%生成一個(gè)待傳的隨機(jī)信息符號(hào)for i=1:N temp=rand。 if (temp) data=1。 else data=1。 end。 % 將一位信息比特重復(fù)LC次 for j=1:Lc repeat_data(j)=data。 end。 % 生成PN碼 for j=1:Lc temp=rand。 if (temp) pn_seq(j)=1。 else pn_seq(j)=1。 end。 end。 % 待傳輸?shù)臄U(kuò)頻信號(hào) trans_sig=sqrt(E_chip)*repeat_data.*pn_seq。 noise=sgma*randn(1,Lc)。 %高斯白噪聲 % 生成單頻率正弦干擾信號(hào) n=(i1)*Lc+1:i*Lc。 interference=Ac*sin(W0*n)。 %接收機(jī)接收的信號(hào) rec_sig=trans_sig+noise+interference。 temp=rec_sig.*pn_seq。 %進(jìn)行解擴(kuò) %進(jìn)行輸出判決 decision_variable=sum(temp)。 if (decision_variable0) decision=1。 else decision=1。 end。 if (decision~=data) num_of_err=num_of_err+1。 %誤碼計(jì)數(shù) end。end。p=num_of_err/N。 %誤碼率計(jì)數(shù)器  仿真結(jié)果及其分析圖32 DS系統(tǒng)在不同振幅的正弦干擾下的誤碼率曲線仿真結(jié)果給出了由3個(gè)不同的正弦干擾振幅干擾下的誤碼率,同時(shí)也繪出無(wú)干擾時(shí)的誤碼率。在誤碼率運(yùn)算中加性噪聲的方差是固定的,總共進(jìn)行了10000次運(yùn)算而且在每次運(yùn)算中,信號(hào)能量取值都是所需的信噪比要求。參考附錄[12]2，由仿真結(jié)果可知,在無(wú)干擾的情況下S/N=3dB時(shí)，系統(tǒng)就可達(dá)到語(yǔ)音質(zhì)量為3(即PE＜103)的要求,在S/N=4dB時(shí)誤碼率就可達(dá)到0；在有干擾的情況下，從仿真曲線上可看出，對(duì)于不同的干擾幅度，當(dāng)S/N達(dá)到一定值后，彎曲部分就比較彎，也就是說(shuō)信噪比只要增加一點(diǎn)誤碼率就急劇下降?？梢奃SSS系統(tǒng)對(duì)單頻正弦干擾有較強(qiáng)的抑制作用。 抗窄帶干擾仿真及結(jié)果分析 建立抗窄帶干擾仿真模型DS擴(kuò)頻系統(tǒng)的抗窄帶性能仿真模型和上面模型大體相同，只是干擾信號(hào)不同罷了。這里我們用的是窄帶高斯隨機(jī)過(guò)程，由于其特性為，包絡(luò)服從瑞利分布，相位服從均勻分布，股我們可以間接的產(chǎn)生窄帶高斯隨機(jī)過(guò)程。由于MATLAB沒有直接產(chǎn)生瑞利分布的函數(shù)，這里我們借助均勻隨機(jī)過(guò)程到瑞利分布隨機(jī)過(guò)程映射的辦法來(lái)產(chǎn)生服從瑞利分布的包絡(luò)，具體過(guò)程參考附錄[12]。仿真模型如圖33所示。仿真的結(jié)果同樣也是在不同信噪比下其誤碼率的曲線。同樣誤碼率用對(duì)數(shù)表示。仿真結(jié)果如圖34所示。 均勻RNG 重復(fù)Lc次檢測(cè)器PN碼比較差錯(cuò)計(jì)數(shù)器WGN窄帶信號(hào)圖33 DS擴(kuò)頻系統(tǒng)的抗窄帶性能仿真模型 編寫仿真代碼主函數(shù)：echo onLc=20。W0=1。 %窄帶干擾信號(hào)的中心頻率SNRindB=0:2:30。for i=1:length(SNRindB) smld_err_prb(i)=wubitlv(SNRindB(i),Lc,W0)。 echo offend。% 繪制仿真結(jié)果曲線，誤碼率用對(duì)數(shù)表示echo onsemilogy(SNRindB,smld_err_prb)。axis([0,10,0,10^1])。xlabel(39。信噪比Eb/N039。)。ylabel(39。誤碼率Pb39。)。title(39。DS系統(tǒng)在窄帶干擾下的誤碼率曲線39。)子函數(shù)：%計(jì)算誤碼率的m文件function [p]=wubitlv(snr_in_dB,Lc,W0)snr=10^(snr_in_dB/10)。sgma=1。Eb=2*sgma^2*snr。E_chip=Eb/Lc。N=10000。num_of_err=0。%生成一個(gè)待傳的隨機(jī)信息符號(hào)for i=1:N temp=rand。 if (temp) data=1。 else data=1。 end。% 將一位信息比特重復(fù)LC次 for j=1:Lc repeat_data(j)=data。 end。 % 生成PN碼 for j=1:Lc temp=rand。 if (temp) pn_seq(j)=1。 else pn_seq(j)=1。 end。 end。% 待傳輸?shù)臄U(kuò)頻信號(hào)trans_sig=sqrt(E_chip)*repeat_data.*pn_seq。%高斯白噪聲noise=sgma*randn(1,Lc)。% 生成窄帶高斯隨機(jī)過(guò)程zhongzi=20。ray_fudu=ray_sig(zhongzi)。 %調(diào)用由均勻隨機(jī)變量到瑞利分布 %變量的函數(shù)，具體過(guò)程見附錄[12] junyun_xiangwei=rand(1,zhongzi)。 %產(chǎn)生窄帶干擾信號(hào)的均勻相位%產(chǎn)生窄帶干擾信號(hào)n=(i1)*Lc+1:i*Lc。interference=ray_fudu.*cos(W0*n+junyun_xiangwei)。%接收機(jī)接收的信號(hào)rec_sig=trans_sig+noise+interference。%進(jìn)行解擴(kuò)temp=rec_sig.*pn_seq。%進(jìn)行輸出判決decision_variable=sum(temp)。if (decision_variable0) decision=1。 else decision=1。 end。 if (decision~=data) num_of_err=num_of_err+1。 %誤碼計(jì)數(shù) end。end。p=num_of_err/N。 %誤碼率計(jì)數(shù)器 仿真結(jié)果及其分析圖34 DS系統(tǒng)在窄帶干擾下的誤碼率曲線仿真結(jié)果曲線如圖34所示，同樣參考附錄[12]2，從圖可以看出當(dāng)S/N=6dB后，系統(tǒng)得語(yǔ)音質(zhì)量就可達(dá)到3。而且從曲線可以看出，S/N=2dB后曲線就比較彎曲，也就是說(shuō)，當(dāng)S/N=2dB后，信噪比增加一點(diǎn)就可以很大程度上減小誤碼率，提高系統(tǒng)的語(yǔ)音質(zhì)量。且在S/N大約為8dB時(shí)誤碼率大約為0?？梢娤到y(tǒng)對(duì)窄帶干擾具有較好的抑制作用。 抗多徑效應(yīng)仿真及結(jié)果分析 建立抗多徑效應(yīng)的仿真模型 多徑效應(yīng)的時(shí)域表示實(shí)際上就是不同路徑來(lái)的信號(hào)存在隨機(jī)的時(shí)延差。這樣我們就在仿真模型上加個(gè)時(shí)延單元，讓帶傳輸信號(hào)經(jīng)過(guò)它進(jìn)行時(shí)延，然后將其和傳輸信號(hào)一起送給接收端進(jìn)行解擴(kuò)檢測(cè)。同以上仿真模型部分相同，只是信道不分一樣?？苟鄰叫?yīng)仿真模型如圖35所示，仿真結(jié)果如圖36所示。它是一時(shí)延為橫坐標(biāo)的，在不同信噪比要求下到曲線，誤碼率不用對(duì)數(shù)。均勻RNG 重復(fù)Lc次檢測(cè)器PN碼比較差錯(cuò)計(jì)數(shù)器WGN延時(shí)圖35 抗多徑效應(yīng)仿真模型 編寫仿真代碼主函數(shù)：echo onLc=20。Tc=0:20。 %多徑信號(hào)延時(shí)時(shí)間SNRindB1=0 。 %不同信噪比下的干擾情況SNRindB2=5。 SNRindB3=10。 SNRindB4=12。SNRindB5=15。SNRindB6=17。for i=1:length(Tc) smld_err_prb1(i)=wubitlv(SNRindB1,Lc,Tc(i))。 smld_err_prb2(i)=wubitlv(SNRindB2,Lc,Tc(i))。 smld_err_prb3(i)=wubitlv(SNRindB3,Lc,Tc(i))。 smld_err_prb4(i)=wubitlv(SNRindB4,Lc,Tc(i))。 smld_err_prb5(i)=wubitlv(SNRindB5,Lc,Tc(i))。 smld_err_prb6(i)=wubitlv(SNRindB6,Lc,Tc(i))。 echo offend。subplot(6,1,1)。 %繪制仿真結(jié)果曲線圖plot(Tc,smld_err_prb1)。title(39。DS系統(tǒng)在不同多徑信號(hào)的干擾下的誤碼率曲線39。)。ylabel(39。誤碼率Pb39。)。subplot(6,1,2)。plot(Tc,smld_err_prb2)。subplot(6,1,3)。plot(Tc,smld_err_prb3)。subplot(6,1,4)。plot(Tc,smld_err_prb4)。subplot(6,1,5)。plot(Tc,smld_err_prb5)。subplot(6,1,6)。plot(Tc,smld_err_prb6)。xlabel(39。延遲Tc39。)。子函數(shù)：function [p]=wubitlv(snr_in_dB,Lc,Tc) %計(jì)算誤碼率的m文件snr=10^(snr_in_dB/10)。sgma=1。Eb=2*sgma^2*snr。E_chip=Eb/Lc。N=10000。num_of_err=0。%生成一個(gè)待傳的隨機(jī)信息符號(hào)for i=1:N temp=rand。 if (temp) data=1。 else data=1。 end。 for j=1:Lc repeat_data(j)=data。 end。 for j=1:Lc %生成PN碼 temp=rand。 if (temp) pn_seq(j)=1。 else pn_seq(j)=1。 end。 end。 trans_sig=sqrt(E_chip)*repeat_data.*pn_seq。 %待傳輸?shù)臄U(kuò)頻信號(hào) noise=sgma*randn(1,Lc)。 %白噪聲 % 生成延時(shí)信號(hào) a=zeros(1,Tc)。 delay_sig=[a,trans_sig]。 rec_sig=[trans_sig,a]+[noise,a]+delay_sig。 %接收機(jī)接收的信號(hào) temp=rec_sig.*[pn_seq,a]。 %進(jìn)行解擴(kuò) decision_variable=sum(temp)。 %進(jìn)行輸出判決 if (decision_variable0) decision=1。 else decision=1。 end。 if (decision~=data) num_of_err=num_of_err+1。 end。end。p=num_of_err/N。圖36 DS系統(tǒng)抗多經(jīng)仿真結(jié)果仿真結(jié)果如上圖所示，信噪比越大它的抗多徑效應(yīng)越好。同樣參考附錄2，從仿真圖上可知，此處誤碼率不是對(duì)數(shù)表示，當(dāng)S/N=0dB時(shí)誤碼率太大(＞)，即語(yǔ)音質(zhì)量小于3，使系統(tǒng)不能正常工作，同時(shí)可以看出，