仿真高斯白噪声信道下QPSK的EbN0与误比特率之间的关系

键入文字QPSK调制与解调在MATLAB平台上的实现QPSK即四进制移向键控（四进制相移键控)，它利用载波的四种不同相位来表示数字信息，由于每一种载波相位代表两个比特信息，因此每个四进制码元可以用两个二进制码元的组合来表示。两个二进制码元中的前一个码元用a表示，后一个码元用b表示。QPSK信号可以看作两个载波正交2PSK信号的合成，下图表示QPSK正交调制器。由QPSK信号的调制可知，对它的解调可以采用与2PSK信号类似的解调方法进行解调。解调原理图如下所示，同相支路和正交支路分别采用相干解调方式解调，得到和，经过抽样判决和并/串交换器，将上下支路得到的并行数据恢复成串行数据。%调相法清除所有关闭所有t=-1:0.01:7-0.01；tt=长度(t)；x1=1(1，800)；i=1:tt吨吨如果(t(I)=1t(I)=1)|(t(I)=5t(I)=7)；x1(I)=1；否则x1(I)=-1；目标目标t1=0:0.01:8-0.01；T2=0:0.01:7-0.01；T3=-1:0.01:7.1-0.01；T4=0:0.01:8.1-0.01；tt1=长度(t1)；x2=1(1，800)；对于i=1:tt1如果(t1(I)=0 t1(I)=2)|(t1(I)=4 t1(I)=8)；x2(1)=1；否则x2(1)=-1；目标目标f=0:0.1:1xrc=0.50.5 * cos(pi * f)；y1=conv(x1，xrc)/5.5；y2=conv(x2，xrc)/5.5；n0=随机(大小(T2)；f1=1；i=x1 .* cos(2 * pi * f1 * t)；q=x2 .* sin(2 * pi * f1 * t1)；I=i(101:800).Q=q(1:700).QPSK=sqrt(1/2).*一平方英尺(1/2).*问；QPSK_n=(sqrt(1/2).*一平方英尺(1/2).* Q)n0；n1=随机(大小(T2)；i_rc=y1 .* cos(2 * pi * f1 * T3)；q_rc=y2 .* sin(2 * pi * f1 * T4)；一级钢筋混凝土=一级钢筋混凝土(101:800)；Q_rc=q_rc(1:700).QPSK_rc=(sqrt(1/2).* 1/2 .*驻地协调员)；QPSK RC n1=QPSK RC n1；图(1)子情节(4，1，1)；曲线图(t3，I _ RC)；轴线(-1 8-1 1)；伊拉贝尔(a)序列)；子情节(4，1，2)；曲线图(t4，q _ RC)；轴线(-1 8-1 1)；伊拉贝尔(b)序列)；子情节(4，1，3)；图(t2，QPSK特区)；轴线(-1 8-1 1)；伊拉贝尔（合成序列)；子情节(4，1，4)；曲线图(t2，QPSK _ RC _ n1)；轴线(-1 8-1 1)；伊拉贝尔（加入噪声)；效果图：%设定T=1，加入高斯噪声清除所有关闭所有%调制bit_in=randint(1e3，1，0 1)；位_ 1=位_ in(1:2:1 E3)；bit _ Q=bit _ in(2:2:1 E3)；数据1=-2 *比特1；数据_ Q=-2 *比特_ Q1；数据_I1=repmat(数据_ 1，20，1)；数据_Q1=repmat(数据_Q，20，1)；对于i=1:1e4数据_ I2(一)=数据_ I1(一)；数据_ Q2(一)=数据_ Q1(一)；结束；f=0:0.1:1xrc=0.50.5 * cos(pi * f)；数据_I2_rc=conv(数据_I2，xrc)/5.5；数据_Q2_rc=conv(数据_Q2，xrc)/5.5；f1=1；t1=0:0.1:1e3 0.9n0=兰特(大小(t1)；I _ rc=数据_I2_rc .* cos(2 * pi * f1 * t1)；Q _ rc=数据_Q2特区.* sin(2 * pi * f1 * t1)；QPSK_rc=(sqrt(1/2).* 1/2 .*驻地协调员)；QPSK _ RC _ n0=QPSK _ RC _ n0；%解调I_demo=QPSK_rc_n0 .* cos(2 * pi * f1 * t1)；Q_demo=QPSK_rc_n0 .* sin(2 * pi * f1 * t1)；%低通滤波恢复=conv；Q_recover=conv(Q_demo，xrc)；一=一_恢复(11:10010)；Q=Q _恢复(11:10010)；T2=0:0.05:1 E3-0.05；T3=0:0.1:1 E3-0.1；%抽样判决数据恢复=；对于i=1:20:10000数据恢复=数据恢复一(:1:一19)问(:1:一19)；结束；bit _ recover=；对于i=1:20:20000如果总和(数据恢复(输入：输入19)0数据恢复=1；位恢复=位恢复1；其他数据恢复-1；位恢复=位恢复-1；目标目标错误=0。DD=-2 *位_ 1；DDD=DD；ddd1=repmat(ddd，20，1)；对于i=1:2e4DD D2(一)=DD D1(一)；目标对于i=1:1e3if bit_recover(i)=ddd(i)错误=错误1；目标目标p=误差/1000。图(1)子情节(2，1，1)；曲线图(t2，DD D2)；坐标轴(0 100-2 2)；标题（原序列)；子情节(2，1，2)；绘图(t2，data _ recover _ a)；坐标轴(0 100-2 2)；标题（解调后序列)；效果图：%设定T=1，不加噪声清除所有关闭所有%调制bit_in=randint(1e3，1，0 1)；位_ 1=位_ in(1:2:1 E3)；bit _ Q=bit _ in(2:2:1 E3)；数据1=-2 *比特1；数据_ Q=-2 *比特_ Q1；数据_I1=repmat(数据_ 1，20，1)；数据_Q1=repmat(数据_Q，20，1)；对于i=1:1e4数据_ I2(一)=数据_ I1(一)；数据_ Q2(一)=数据_ Q1(一)；结束；t=0:0.1:1 E3-0.1；f=0:0.1:1xrc=0.50.5 * cos(pi * f)；数据_I2_rc=conv(数据_I2，xrc)/5.5；数据_Q2_rc=conv(数据_Q2，xrc)/5.5；f1=1；t1=0:0.1:1e3 0.9I _ rc=数据_I2_rc .* cos(2 * pi * f1 * t1)；Q _ rc=数据_Q2特区.* sin(2 * pi * f1 * t1)；QPSK_rc=(sqrt(1/2).* 1/2 .*驻地协调员)；%解调演示=QPSK .* cos(2 * pi * f1 * t1)；问_演示=QPSK特区.* sin(2 * pi * f1 * t1)；恢复=conv；Q_recover=conv(Q_demo，xrc)；一=一_恢复(11:10010)；Q=Q _恢复(11:10010)；T2=0:0.05:1 E3-0.05；T3=0:0.1:1 E3-0.1；数据恢复=；对于i=1:20:10000数据恢复=数据恢复一(:1:一19)问(:1:一19)；结束；DDD=-2 *比特1；ddd1=repmat(ddd，10，1)；对于i=1:1e4DD D2(一)=DD D1(一)；目标图(1)子情节(4，1，1)；图(t3，I)；坐标轴(0 20-6 6)；子情节(4，1，2)；曲线图(t3，Q)；坐标轴(0 20-6 6)；子情节(4，1，3)；绘图(t2，data _ recover)；坐标轴(0 20-6 6)；子情节(4，1，4)；图(DDD 2)；坐标轴(0 20-6 6)；效果图：% QPSK误码率分析SNRindB1=0:2:10SNRindB2=0:0.1:10i=1:length(SNRindB1)pb，PS=cm _ sm32(Snrindb1(1)；(1)=铅；符号误差脉冲宽度(1)=脉冲宽度；结束；i=1:length(SNRindB2)信噪比=经验（信噪比2(I)*对数(10)/10)；theo _ err _ prb(1)=Qfunct(2 *信噪比)；结束；标题(QPSK误码率分析)；符号(Snrindb1，smld_bit_err_prb，*)；坐标轴(0 10 10e-8 1)；坚持住。% s logy(Snrindb1，smld_symbol_err_prb，o)；符号学(Snrindb 2，theo _ err _ prb)；图例（仿真比特误码率，理论比特误码率)；拖延；函数y=Qfunct(x)y=(1/2)* erfc(x/sqrt(2)；函数pb，ps=cm_sm32(SNRindB)N=10000 .E=1 .snr=10(snrindb/10)；sgma=sqrt(信噪比)/2；s00=10；s01=0 1；S11=-10；S10=0-1；对于i=1:N牛温度=兰特；如果(温度0.25)数据源1(I)=0；数据源2(I)=0；elseif(温度0.5)数据源1(I)=0；数据源2(I)=1；elseif(温度0.75)数据源1(I)=1；数据源2(I)=0；其他数据源1(I)=1；数据源2(I)=1；结束；结束；数字符号错误=0；numbieterror=0；对于i=1:N牛n=sgma*randn(大小(s00)；如果(d源1(I)=0)(d源2(I)=0)r=s00 n；否则，如果(数据源1(I)=0)(数据源2(I)=1)r=s01 n；否则，如果(数据