调制解调作业文档格式.docx

1、解调器根据星座图及接收到的载波信号的相位来判断发送端发送的信息比特。（b）误比特率为1e-2、1e-3及1e-4时的Eb/N0分别是多少？BPSKQPSK1e-2 4.35 4.351e-36.791e-48.40 （c）从物理意义上说明为什么bpsk、qpsk误比特率曲线是重叠的因为bpsk是双极性信号，成对出错概率为= 对等概率发送的二进制信号，成对出错率等于符号错误概率，也等于误比特率。Qpsk是那双正交信号，用相干接收机接收，成对错误概率为= 由联合边界计算的误码率是成对错误概率的两倍，误比特率是误码率的一半（d）当samples减少为100000，10000，1000时观察误比特率曲

2、线的变化，你得出什么结论。1e6个点1e5个点1e4个点1e3个点结论，取样点越少，在同等信噪比下误码率越大。取样点越多，曲线越光滑。2.用matlab运行“bpsk.m”、“qpsk.m” （a）在各程序中标注“注释”处加上注释（英文或中文）BPSK.mdata1=rand（1,nd）0.5; % 注释返回一组行向量，nd个码元%* Attenuation Calculation * code_rate=1; E=1; sigma=E/sqrt（2*SNR*code_rate）; %注释：计算编码后高斯分布随机变量标准差for i=1:nd if （data1（i）=0）, data4（i）

3、=-E+Gngauss（sigma）; % 注释：传0码时加入高斯分布白噪声干扰，双极性码 else data4（i）=E+Gngauss（sigma）;传1码时加入高斯分布白噪声干扰，双极性码 end;%* BPSK Demodulation * demodata1=data4 0; % 注释：data4中每个元素和0比较，比较结果为真（1）和假（0），然后再将这些0,1赋值给demodata1%* Bit Error Rate （BER） * noe2=sum（abs（data1-demodata1）;计算误码个数 nod2=length（data1）;计算总码数 noe=noe+noe2

4、; nod=nod+nod2;end %* Output result *ber（snr_num） = noe/nod %注释：计算误码率，误码除以总传递码数QPSK.mSNR=exp（snr_in_dB（snr_num）*log（10）/10% 注释：计算信噪比for iii=1:nloop%* Data generation * data=rand（1,nd*ml）生成 同相分量 data1=2*data-1 % 注释：原始数据流分成两部分，此为正交分量%* QPSK Modulation * tout=qpskmod（data1,1,nd,ml）;%* Attenuation Calcu

5、lation *%* gngauss channel* E=1/sqrt（2）; sigma=E/sqrt（2*SNR*code_rate） % 注释： for i=1:nd, gsrv1,gsrv2=Gngauss（sigma）; tout（i）=tout（i）+gsrv1+j*gsrv2; %sigma=E/sqrt（2*SNR*code_rate）; %n = 1/sqrt（2）*randn（1,nd） + j*randn（1,nd）; %sigma.*n %注释：另一种高斯白噪声产生方法%* QPSK Demodulation * demodata=qpskdemod（tout,1,n

6、d,ml）; demodata1=demodata 与原始码元对比判决，得到正确的码元数 demodata=demodata1;%* Bit Error Rate （BER） * noe2=sum（abs（data-demodata）;计算错误码元数nod2=length（data1）;总码元数noe=noe+noe2;nod=nod+nod2;ber（snr_num） = noe/nod;计算误码率end;isi = isi + 2.（ m2 - ii ） .* paradata2（1:para）,ii+count2）;产生同相分量isq = isq + 2.（ m2 - ii ） .* p

7、aradata2（1:para）,m2+ii+count2）;产生正交分量 （b）说明加性高斯白噪声的产生方法，请再给出一种加性高斯白噪声的产生方法，并验证其正确性。原始方法为Gngauss.mif nargin = 0, m=0; sgma=1;elseif nargin = 1, sgma=m;u=rand; % a uniform random variable in （0,1） z=sgma*（sqrt（2*log（1/（1-u）; % a Rayleigh distributed random variable % another uniform random variable in

8、 （0,1） gsrv1=m+z*cos（2*pi*u）;gsrv2=m+z*sin（2*pi*u）;nargin为输入变量的个数，如果nargin为零，即空号，产生均值为0，标准差为1的高斯分布系列。如果传号，产生均值为0，方差为m的随机序列。另一种方法：y = awgn（m,n,p） 产生一个m行n列的高斯白噪声的矩阵，p以dBW为单位指定输出噪声的强度。但与实际差距较大。data2=2*data1-1; data4=awgn（data2,snr_in_dB（snr_num） （c）参考“bpskqpsk125.m”的画图功能，给出Eb/N0-误比特率曲线和高斯信道下的理论误比特率曲线。

9、BPSK误比特率曲线和高斯信道下的理论误比特率曲线 QPSK误比特率曲线和高斯信道下的理论误比特率曲线clear all% qpsk.m% Simulation program to realize QPSK transmission system%* Preparation part%*ml=2; sr=25600; % ml:Number of modulation levels （BPSK:ml=1, QPSK:ml=2, 16QAM:ml=4）br=sr .* ml; % Bit ratend =100; % Number of symbols that simulates in ea

10、ch loopsnr_in_dB=0:18 ber=zeros（1,length（snr_in_dB）;for snr_num=1:length（snr_in_dB）SNR=exp（snr_in_dB（snr_num）*log（10）/10）;将SNR（dB）转化为倍数%* START CALCULATION *nloop=100; % Number of simulation loopsnoe = 0; % Number of error datanod = 0; % Number of transmitted data随机产生相应数目的0-1信号将单极性0-1信号转化为双极性的-1,1信号 sigma=E