MATLAB仿真AM调制解调无线通信实验报告

1、 无线通信实 验 报 告院系名称： 信息科学与工程学院 专业班级： 电信 班 学生姓名： 学 号： 授课教师： 2014 年 11 月 6 日实验一 高斯衰落信道建模 一、基本原理QPSK信号可以看成是对两个正交的载波进行多电平双边带调制后所得信号的叠加，因此可以用正交调制的方法得到QPSK信号。QPSK信号的星座如图4.1.1所示：图1.1 QPSK信号星座图从AWGN信道中，在一个信号区间内接收到的带宽信号可以表示为这里和是加性噪声的两个正交分量。 可以将这个接收信号与，给出的和作相关，两个相关器的输出产生受噪声污损的信号分量，它们可表示为式中和定义为 这两个正交噪声分量和是零均值，互不相

2、关的高斯随机过程。这样，和。和的方差是 最佳检测器将接收信号向量r投射到M个可能的传输信号向量之一上去，并选取对应于最大投影的向量。据此，得到相关准则为 ，m=0,1,，M-1由于全部信号都具有相等的能量，因此，对数字相位调制一种等效的检测器标准是计算接收信号向量r=(,)的相位为 并从信号集中选取其相位最接近的信号。在AWGN信道中，因为二相相位调制与二进制PAM是相同的，所以差错概率为，式中是每比特的能量。四相相位调制可以看作两个正交载波上的二相相位调制系统，所以1个比特的差错概率与二相相位调制是一样的。对于M4的符号差错概率不存在简单的闭式表达式。对的一种好的近似式是 式中比特/符号。Q

3、PSK调制解调原理四相相位键控（QPSK）也称之为正交PSK，其调制原理如图4.1.2所示。 图1.2 QPSK调制原理图如果输入的二进制信息码流（假设+1V为逻辑1，-1V为逻辑0）串行进入比特分离器，产生2个码流以并行方式输出，分别被送入I（正交支路）通道及Q（同相支路）通道，又各自经过一个平衡调制器，与一个和参考振荡器同频的正交的载波（和）调制形成了四相相移键控信号即得到平衡器的输出信号后，经过一个带通滤波器，然后再进入行信号叠加，可以得到已经调制的QPSK信号。QPSK的4种（I，Q星座图合为4种0 0，0 1，1 0，1 1）输出相位有相等的幅度，而且2个相邻的相位相差值为90度，但

4、是输出相位并不满足（m=0,1,M-1）,信号相位移可以偏移45度和-45度，接受端仍可以得到正确的解码，实际中数字输入电压必须比峰值载波电压高出很多，以确保平衡器的正常工作。经过调制的信号通过信道传输到达用户端，需要进行解调，这样一过程是与调制相类似的逆过程。首先，QPSK信号经过功率分离器形成两路相同的信号，进入乘积检验波，用两个正交的载波信号（和）实现相干解调，然后各自通过一个低通滤波器滤波得到低频和直流的成分，再经过一个并行-串行变换器，得到解调信号。QPSK的解调原理如图4.1.3所示。 图1.3 QPSK解调原理图二、结果分析 由图可见QPSK仿真误码率曲线和理论误码率曲线重合在一

5、起，QPSK仿真误比特率曲线和理论误比特率曲线也重合在一起，误码率约是误比特率的两倍，说明实验方法是正确可行的。QPSK信号的误码率：QPSK信号的误比特率：误码率是误比特率的两倍。三、源程序%2014.10.25%高斯衰落信道建模%clcclearSNR_DB=0:1:12;sum=1000000;data= randsrc(sum,2,0 1); %生成0、1等概出现的sum*2矩阵a1,b1=find(data(:,1)=0&data(:,2)=0); %第一列为0，第二列也为0时message(a1)=-1-j; %输出-1-ja2,b2=find(data(:,1)=0&am

6、p;data(:,2)=1); %第一列为0，第二列为1时message(a2)=-1+j; %输出-1+ja3,b3=find(data(:,1)=1&data(:,2)=0); %第一列为1，第二列为0时message(a3)=1-j; %输出1-ja4,b4=find(data(:,1)=1&data(:,2)=1); %第一列为1，第二列为1时message(a4)=1+j; %输出1+jA=1;Tb=1;Eb=A*A*Tb;P_signal=Eb/Tb;NO=Eb./(10.(SNR_DB/10);P_noise=P_signal*NO;sigma=sqrt(P_no

7、ise);for Eb_NO_id=1:length(sigma) noise1=sigma(Eb_NO_id)*randn(1,sum); noise2=sigma(Eb_NO_id)*randn(1,sum); receive=message+noise1+noise2*j; resum=0; total=0; m1=find(angle(receive)<=pi/2&angle(receive)>0);%收到信号角度大于0小于pi/2时 remessage(1,m1)=1+j; %信号恢复为1+j即（1,1） redata(m1,1)=1; redata(m1,2)=

8、1; m2= find( angle(receive)>pi/2&angle(receive)<=pi);%收到信号角度大于pi/2小于pi时 remessage(1,m2)=-1+j; %信号恢复为-1+j即（0,1） redata(m2,1)=0; redata(m2,2)=1; m3=find( angle(receive)>-pi&angle(receive)<=-pi/2);%收到信号角度大于-pi小于-pi/2时 remessage(1,m3)=-1-j; %信号恢复为-1-j即（0,0） redata(m3,1)=0; redata(m3,

9、2)=0; m4=find( angle(receive)>-pi/2&angle(receive)<=0);%收到信号角度大于-pi/2小于0时 remessage(1,m4)=1-j; %信号恢复为1+j即（1,0） redata(m4,1)=1; redata(m4,2)=0; resum,ratio1=symerr(data,redata); pbit(Eb_NO_id)=resum/(sum*2); total,ratio2=symerr(message,remessage); pe(Eb_NO_id)=total/sum;endPe=1-(1-1/2*erfc(

10、sqrt(10.(SNR_DB/10)/2).2; %计算输出误码率Pbit=1/2*erfc(sqrt(10.(SNR_DB/10)/2); %计算输出误比特率semilogy(SNR_DB,pe,':s',SNR_DB,pbit,'-o') %画出PSK仿真误码率', legend('QPSK仿真误码率','QPSK仿真误比特率',1) %'QPSK仿真误比特率的曲线xlabel('信噪比/dB')ylabel('概率P')实验二 AM调制信号的MATLAB实现一、实验原理 标准

11、调幅波（AM）产生原理：调制信号是来自信源的调制信号（基带信号），这些信号可以是模拟的，亦可以是数字的。为首调制的高频振荡信号可称为载波，它可以是正弦波，亦可以是非正弦波(如周期性脉冲序列)。载波由高频信号源直接产生即可，然后经过高频功率放大器进行放大，作为调幅波的载波，调制信号由低频信号源直接产生，二者经过乘法器后即可产生双边带的调幅波。 设载波信号的表达式为，调制信号的表达式为 ，则调幅信号的表达式为： 图1.1 标准调幅波示意图所谓非相干解调（包络检波）是在接收端解调信号时不需要本地载波，而是利用已调信号中的包络信号来恢复原基带信号。因此，非相干解调一般只适用幅度调制（AM）系统。忧郁包

12、络解调器电路简单，效率高，所以几乎所有的幅度调制（AM）接收机都采用这种电路。如下为串联型包络检波器的具体电路。图2.2 AM信号的非相干解调原理当RC满足条件时，包络检波器的输出基本与输入信号的包络变化呈线性关系，即 其中，。隔去直流后就得到原信号二、结果分析 包络检波是在的情况下才能实现的，由上图A=2和A=1两种情况下可以看出当A=2时满足条件，AM波的包络与调制信号m(t)的形状完全一样，可以通过隔直流的方法（包络解调）得到原信号；当A=1时，不满足条件出现欠调幅，AM波的包络失真无法通过隔直流的方法（包络解调）得到原信号，所以在这种情况下不能用包络检波的方法，只能通过相干解调。 三、

13、源代码%2014.11.1%AM调制信号的MATLAB实现 %clear all;close all;dt=0.0001; %时间采样频谱fc=10; %载波中心频率T=5; %信号时长N=T/dt;t=0:N-1*dt;mt=sqrt(2)*cos(2*pi*t); %信源A=2;s_am=(A+mt).*cos(2*pi*fc*t);subplot(211);plot(t,s_am);hold on; %画出AM信号波形plot(t,A+mt,'r-');title('AM调制信号及其包络 A=2');A=1;s_am=(A+mt).*cos(2*pi*fc

14、*t);subplot(212);plot(t,s_am);hold on; %画出AM信号波形plot(t,A+mt,'r-');title('AM调制信号及其包络 A=1');实验三 单极性归零码的matlab实现一、实验原理 单极性归零码（RZ）即是以高电平和零电平分别表示二进制码1 和0，而且在发送码1 时高电平在整个码元期间T 只持续一段时间，其余时间返回零电平在单极性归零码中，/T 称为占空比单极性归零码的主要优点是可以直接提取同步信号，因此单极性归零码常常用作其他码型提取同步信号时的过渡码型也就是说其他适合信道传输但不能直接提取同步信号的码型，可先

15、变换为单极性归零码，然后再提取同步信号。二、结果分析 由上图可知10111101 其中1前半个码元期间是高电平，后半个码元期间是低电平（即后半段归零），0是一个码元期间都是低电平。该码的占空比为1/2,由波形可以看出单极性归零码在matlab上得到了验证。三、源代码%2014.11.1%单极性归零码的matlab实现%clear all;close all;x=1 0 1 1 1 1 0 1;y=x;t0=200;t=0:1/t0:length(x);for i=1:length(x) if x(i)=1 for j=1:t0/2 y(2*i-2)*t0/2+j)=1; y(2*i-1)*t0

16、/2+j)=0; end else for j=1:t0 y(i-1)*t0+j)=0; end endendy=y,x(i);plot(t,y);title('单极性归零码:10111101');grid on;axis(0,i,-0.1,1.1);xlabel('t/s');ylabel('幅 度');实验四 线性分组码的matlab实现一、实验原理 线性分组码是一类奇偶校验码，它可以由（n，k）形式表示。编码器将一个k比特信息分组（信息矢量）转变为一个更长的由给定元素符号集组成的n比特编码分组。当这个符号集包含两个元素（0和1），与二进制相对，称为二进制编码。 将新的信源序列与旧的信源序列之间的关系用矩阵表达出来，然后通过矩阵的运算生成新的序列。本实验的旧信源序列为M=(0 0 0)(1 1 1)共八个，通过给出的关系求出生成矩阵为G= 1 0 0 1 1 1; 0 1 0 1 1 0;0 0 1 0 1 1 ;C=M*G求出新的序列。二、结果分析C = c5 c4 c3 c2 c1 c0 0 0 0 0 0 0 0 0 1 0 1 1 0 1 0 1 1 0