基于matlab的QPSK与BPSK信号性能比较仿真(3)

┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ 装 ┊ ┊ ┊ ┊ ┊ 订 ┊ ┊ ┊ ┊ ┊ 线 ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊

瑞利频率选择性衰落

4.1.2 横向比较分析

相同信道下，BPSK调制的系统误码率小于QPSK调制。本次仿真过程中，莱斯平坦衰落信道下两者误码率差别显著。

4.2 仿真结果分析

4.2.1 误码率分析

由前面介绍的误码率内容可以看出，QPSK判决门限为90，BPSK的判决门限为180。因此相同系统情况下的误码率BPSK优于QPSK。

??4.2.2 频带利用率比较

在传码率相同的情况下，四进制数字调制系统的信息速率是二进制系统的2倍。频带利用率公式

??Rb B在相同信号速率的情况下，QPSK和BPSK系统的带宽是相同的，但是由于QPSK每个信号都是四进制的，QPSK每个信号包含2bit信息，所以比特率就是BPSK的两倍，因而其频带利用率即为BPSK的两倍。BPSK系统理论的频带利用率最大为1，但是在实际的实现中不能达到1，而在QPSK系统中，频带利用率可以超过1。

┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ 装 ┊ ┊ ┊ ┊ ┊ 订 ┊ ┊ ┊ ┊ ┊ 线 ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊

附 录

代码1

function qpskconstellationA(M) M=4;

x=[0:M-1];

scatterplot(pskmod(x,M)); % A方式QPSK信号 grid on;

return

function qpskconstellationB(M) M=4;

x=[0:M-1];

scatterplot(pskmod(x,M,pi/4)); % B方式QPSK信号 grid on;

return

代码2

function [x,bits]=random_binary(nbits,nsamples)

x=zeros(1,nbits*nsamples); bits=round(rand(1,nbits)); for m=1:nbits

for n=1:nsamples

index=(m-1)*nsamples+n; x(1,index)=(-1)^bits(m); end end return

function [c,lags]=vxcorr(a,b)

%计算a b 的互相关系数，返回C是列向量，长度为ab的长度和减1 a=a(:); b=b(:);

M=length(a); maxlag=M-1;

lags=[-maxlag:maxlag]';

A=fft(a,2^nextpow2(2*M-1));%快速傅立叶变换 B=fft(b,2^nextpow2(2*M-1));

c=ifft(A.*conj(B));%快速傅立叶反变换

c=[c(end-maxlag+1:end,1);c(1:maxlag+1,1)];

[nr nc]=size(a); if(nr>nc) c=c';

lags=lags.'; end return

┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ 装 ┊ ┊ ┊ ┊ ┊ 订 ┊ ┊ ┊ ┊ ┊ 线 ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊

function BER_MC=MCQPSKrun(N,Eb,N0,ChanAtt,TimingBias,TimingJitter,PhaseBias,PhaseJitter)

fs=1e+6; % 采样速率 SymRate=1e+5; % 信号速率 Ts=1/fs; % 采样周期 TSym=1/SymRate; % 信号周期 SymToSend=N; % 发送信号数 ChanBW=4.99e+5; % 带宽

MeanCarrierPhaseError=PhaseBias; % 载波相位均值 StdCarrierPhaseError=PhaseJitter; % 相位误差标准差

MeanSymbolSyncError=TimingBias; % 符号同步误差均值 StdSymbolSyncError=TimingJitter; % 符号同步误差标准差 ChanGain=10^(-ChanAtt/20); % 信道增益 TxBitClock=Ts/2; % 发送机时钟 RxBitClock=Ts/2; % 接收机时钟

%接收机输入端噪声标准差和信号幅度

RxNoiseStd=sqrt((10^((N0-30)/10))*(fs/2)); TxSigAmp=sqrt(10^((Eb-30)/10)*SymRate);

%分配缓存

SampPerSym=fs/SymRate;

probe1=zeros((SymToSend+1)*SampPerSym,1); probe1counter=1;

probe2=zeros((SymToSend+1)*SampPerSym,1); probe2counter=1;

%计已传输信号的个数

TxSymSent=1; RxSymDemod=0;

%发送和接收数据缓冲区

[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(k-1))),and(xor(SourceBitsI(k),SourceBitsQ(k)),xor(SourceBitsQ(k),TxBitsQ(k-1))));

TxBitsQ(k)=or(and(not(xor(SourceBitsI(k),SourceBitsQ(k))),xor(SourceBitsQ(k),TxBitsQ(k-1))),and(xor(SourceBitsI(k),SourceBitsQ(k)),xor(SourceBitsI(k),TxBitsI(k-1))));

end

% 产生复信号

TxBits=((TxBitsI*2)-1)+(sqrt(-1)*((TxBitsQ*2)-1));

RxIntegrator=0; %初始化接收机积分器 TxBitClock=2*TSym; %初始化发送机

%设计信道滤波器，产生滤波器参数序列 [b,a]=butter(2,ChanBW/(fs/2)); b=[1]; a=[1];

[junk,FilterState]=filter(b,a,0);

% 开始仿真循环

┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ 装 ┊ ┊ ┊ ┊ ┊ 订 ┊ ┊ ┊ ┊ ┊ 线 ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊ ┊

while TxSymSent

TxBitClock=TxBitClock+Ts; if TxBitClock>TSym

TxSymSent=TxSymSent+1;

TxBitClock=mod(TxBitClock,TSym);

TxOutput=TxBits(TxSymSent)*TxSigAmp; end

%信号经过信道滤波器

[Rx,FilterState]=filter(b,a,TxOutput,FilterState); %加高斯白噪声

Rx=(ChanGain*Rx)+(RxNoiseStd*(randn(1,1)+sqrt(-1)*randn(1,1))); %基于接收机载波同步误差的相位旋转

PhaseRotation=exp(sqrt(-1)*2*pi*(MeanCarrierPhaseError+(randn(1,1)*StdCarrierPhaseError))/360); Rx=Rx*PhaseRotation; probe1(probe1counter)=Rx;

probe1counter= probe1counter+1; %更新接收机积分清除器 RxIntegrator=RxIntegrator+Rx;

probe2(probe2counter)=RxIntegrator; probe2counter= probe2counter+1;

%更新接收机时钟，判断是不是适合采样 RxBitClock=RxBitClock+Ts;

RxTSym=TSym*(1+MeanSymbolSyncError+(StdSymbolSyncError*randn(1,1))); if RxBitClock>RxTSym %解调信号 RxSymDemod=RxSymDemod+1;

RxBitsI(RxSymDemod)=round(sign(real(RxIntegrator))+1)/2; RxBitsQ(RxSymDemod)=round(sign(imag(RxIntegrator))+1)/2; RxBitClock=RxBitClock-TSym; RxIntegrator=0; end end

% 差分解码

SinkBitsI=SourceBitsI*0; SinkBitsQ=SourceBitsQ*0;

for …… 此处隐藏：4522字，全部文档内容请下载后查看。喜欢就下载吧 ……