基于matlab的qpsk與bpsk信號(hào)性能比較仿真(編輯修改稿)
2024-09-06 16:23 本頁面
【文章內(nèi)容簡介】 aseError=PhaseBias。 % 載波相位均值StdCarrierPhaseError=PhaseJitter。 % 相位誤差標(biāo)準(zhǔn)差MeanSymbolSyncError=TimingBias。 % 符號(hào)同步誤差均值StdSymbolSyncError=TimingJitter。 % 符號(hào)同步誤差標(biāo)準(zhǔn)差ChanGain=10^(ChanAtt/20)。 % 信道增益TxBitClock=Ts/2。 % 發(fā)送機(jī)時(shí)鐘RxBitClock=Ts/2。 % 接收機(jī)時(shí)鐘 %接收機(jī)輸入端噪聲標(biāo)準(zhǔn)差和信號(hào)幅度 RxNoiseStd=sqrt((10^((N030)/10))*(fs/2))。TxSigAmp=sqrt(10^((Eb30)/10)*SymRate)。 %分配緩存 SampPerSym=fs/SymRate。probe1=zeros((SymToSend+1)*SampPerSym,1)。probe1counter=1。probe2=zeros((SymToSend+1)*SampPerSym,1)。probe2counter=1。 %計(jì)已傳輸信號(hào)的個(gè)數(shù) TxSymSent=1。RxSymDemod=0。 %發(fā)送和接收數(shù)據(jù)緩沖區(qū)[unused,SourceBitsI]=random_binary(SymToSend,1)。[unused,SourceBitsQ]=random_binary(SymToSend,1)。 %差分編碼TxBitsI=SourceBitsI*0。TxBitsQ=SourceBitsQ*0。for k=2:length(TxBitsI) TxBitsI(k)=or(and(not(xor(SourceBitsI(k),SourceBitsQ(k))),xor(SourceBitsI(k),TxBitsI(k1))),and(xor(SourceBitsI(k),SourceBitsQ(k)),xor(SourceBitsQ(k),TxBitsQ(k1))))。 TxBitsQ(k)=or(and(not(xor(SourceBitsI(k),SourceBitsQ(k))),xor(SourceBitsQ(k),TxBitsQ(k1))),and(xor(SourceBitsI(k),SourceBitsQ(k)),xor(SourceBitsI(k),TxBitsI(k1))))。end % 產(chǎn)生復(fù)信號(hào) TxBits=((TxBitsI*2)1)+(sqrt(1)*((TxBitsQ*2)1))。 RxIntegrator=0。 %初始化接收機(jī)積分器TxBitClock=2*TSym。 %初始化發(fā)送機(jī) %設(shè)計(jì)信道濾波器，產(chǎn)生濾波器參數(shù)序列[b,a]=butter(2,ChanBW/(fs/2))。b=[1]。a=[1]。[junk,FilterState]=filter(b,a,0)。 % 開始仿真循環(huán) while TxSymSentSymToSend TxBitClock=TxBitClock+Ts。 if TxBitClockTSym TxSymSent=TxSymSent+1。 TxBitClock=mod(TxBitClock,TSym)。 TxOutput=TxBits(TxSymSent)*TxSigAmp。 end %信號(hào)經(jīng)過信道濾波器 [Rx,FilterState]=filter(b,a,TxOutput,FilterState)。 %加高斯白噪聲 Rx=(ChanGain*Rx)+(RxNoiseStd*(randn(1,1)+sqrt(1)*randn(1,1)))。 %基于接收機(jī)載波同步誤差的相位旋轉(zhuǎn) PhaseRotation=exp(sqrt(1)*2*pi*(MeanCarrierPhaseError+(randn(1,1)*StdCarrierPhaseError))/360)。 Rx=Rx*PhaseRotation。 probe1(probe1counter)=Rx。 probe1counter= probe1counter+1。 %更新接收機(jī)積分清除器 RxIntegrator=RxIntegrator+Rx。 probe2(probe2counter)=RxIntegrator。 probe2counter= probe2counter+1。 %更新接收機(jī)時(shí)鐘，判斷是不是適合采樣 RxBitClock=RxBitClock+Ts。 RxTSym=TSym*(1+MeanSymbolSyncError+(StdSymbolSyncError*randn(1,1)))。 if RxBitClockRxTSym %解調(diào)信號(hào) RxSymDemod=RxSymDemod+1。 RxBitsI(RxSymDemod)=round(sign(real(RxIntegrator))+1)/2。 RxBitsQ(RxSymDemod)=round(sign(imag(RxIntegrator))+1)/2。 RxBitClock=RxBitClockTSym。 RxIntegrator=0。 endend% 差分解碼SinkBitsI=SourceBitsI*0。SinkBitsQ=SourceBitsQ*0。 for k=2:RxSymDemod SinkBitsI(k)=or(and(not(xor
職業(yè)教育相關(guān)推薦
鄂ICP備17016276號(hào)-1