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