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

1、®丙姑讹摩无线通信实验报告院系名称：信息科学与工程学院专业班级：电信班学生姓名：学 号：授课教师：2014年 11 月 6日实验一高斯衰落信道建模一、基本原理QPSK信号可以看成是对两个正交的载波进行多电平双边带调制后所得信号 的叠加，因此可以用正交调制的方法得到 QPSK言号。图1.1 QPSK信号星座图从AWG信道中，在一个信号区间内接收到的带宽信号可以表示为r t =Um(t) n(t)= Um(t) nc(t)cos(2二 fct) -ns(t)sin(2 二 f£)这里nc(t)和ns(t)是加性噪声的两个正交分量。可以将这个接收信号与 ft)二gT(t)cos(

2、2二fct), 'J =-gT(t)s in(2二fct)给出的H(t)和'- 2(t)作相关，两个相关器的输出产生受噪声污损的信号分量，它们可表示为r = sm n=(.一 s cosns)式中nc和ns定义为1 :nc；gT(t)nc(t)dt1n- _:.gT(t)ns(t)dt这两个正交噪声分量nc(t)和ns(t)是零均值，互不相关的高斯随机过程。这 样，E(nc) = E(Hs) = 0和 E(njs)=O。n和 n 的方差是E(n2)= E( nf)=叫2最佳检测器将接收信号向量r投射到M个可能的传输信号向量务之一上 去，并选取对应于最大投影的向量。据此，得到相关

3、准则为C(r,為)=r *sm， m=0,1,，M-1由于全部信号都具有相等的能量，因此，对数字相位调制一种等效的检测器标准 是计算接收信号向量r=( rc, rs)的相位为并从信号集 Sm中选取其相位最接近 片的信号。在AWG信道中，因为二相相位调制与二进制 PAM是相同的，所以差错概率 为巳=0|££1式中-b是每比特的能量。四相相位调制可以看作两个正交载_ No波上的二相相位调制系统，所以1个比特的差错概率与二相相位调制是一样的。对于M - 4的符号差错概率不存在简单的闭式表达式。 对PM的一种好的近似式是式中k =log2 M比特/符号QPS碉制解调原理四相相位键控

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

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

6、PSK仿真误码率QPSK仿真误比特率口-J-51006信噪比/dB1012o10o10由图可见QPSK仿真误码率曲线和理论误码率曲线重合在一起，QPSK仿真误比特率曲线和理论误比特率曲线也重合在一起，误码率约是误比特率的两倍，说明实验方法是正确可行的。QPSK信号的误码率:Peo10o101 QPSK信号的误比特 Pe = 2erfc “ / 2 率:误码率是误比特率的两倍三、源程序%2014.10.25%高斯衰落信道建模%clcclearSNR_DB=0:1:12;sum=1000000;%生成0、1等概出现的sum*2矩阵data= ran dsrc(sum,2,0 1);a1,b1=fi

7、nd(data(:,1)=0&data(:,2)=0); % 第一列为 0,第二列也为 0时 message(a1)=-1-j;% 输出-1-ja2,b2=find(data(:,1)=0&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=fi nd(data(:,1)=1 &data(:,2)=1); % 第一列为 1，第二列为

8、1 时 message(a4)=1+j;% 输出 1+jA=1;Tb=1;Eb=A*A*Tb;P_sig nal=Eb/Tb;N0=Eb./(10.A(SNR_DB/10);P_noise=P_sig nal*NO;sigma=sqrt(P _n oise);for Eb_N 0d=1:le ngth(sigma)n oise1=sigma(Eb_ NOd)*ra ndn (1,sum);n oise2=sigma(Eb_ NOd)*ra ndn (1,sum);receive=message+no ise1+ no ise2*j;resum=0;total=0;m1=fi nd(a ngle

9、(receive)v=pi/2&an gle(receive)>0);% 收至 U信号角度大于0 小于pi/2 时remessage(1,m1)=1+j;% 信号恢复为 1+j 即(1,1)redata(m1,1)=1;redata(m1,2)=1;m2= find( angle(receive)>pi/2&angle(receive)<=pi);% 收至U信号角度大于 pi/2小 于pi时%信号恢复为-1+jremessage(1,m2)=-1+j;即(0,1)redata(m2,1)=0;redata(m2,2)=1;m3=find( angle(rece

10、ive)>-pi&angle(receive)v=-pi/2);% 收至U信号角度大于 -pi小 于-pi/2时remessage(1,m3)=-1-j;%信号恢复为-1-j即(0,0)redata(m3,1)=0;redata(m3,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,

11、redata); pbit(Eb_N 0d)=resum/(sum*2);total,ratio2=symerr(message,remessage); pe(Eb_N Od)=total/sum;endPe=1-(1-1/2*erfc(sqrt(10.A(SNR_DB/10)/2).A2;%计算输出误码率Pbit=1/2*erfc(sqrt(10.A(SNR_DB/10)/2);%计算输出误比特率semilogy(SNR_DB,pe,':s',SNR_DB,pbit,'-o')%画出PSK仿真误码率',legend（'QPSK仿真误码率,

12、9;QPSK仿真误比特率',1）%'QPSK仿真误比特率的曲线xlabel（'信噪比 /dB'）ylabel（'概率 P'）实验一AM调制信号的MATLAB实现一、实验原理标准调幅波（AM产生原理：调制信号是来自信源的调制信号（基带信号）： 这些信号可以是模拟的，亦可以是数字的。为首调制的高频振荡信号可称为载波， 它可以是正弦波，亦可以是非正弦波（如周期性脉冲序列）。载波由高频信号源直 接产生即可，然后经过高频功率放大器进行放大， 作为调幅波的载波，调制信号 由低频信号源直接产生，二者经过乘法器后即可产生双边带的调幅波。设载波信号的表达式为cos

13、 ct，调制信号的表达式为 m t Am cos mt ， 则调幅信号的表达式为：Sam t - A mt IcQS ctCOS ctSAM (t)图1.1标准调幅波示意图所谓非相干解调（包络检波）是在接收端解调信号时不需要本地载波，而是 利用已调信号中的包络信号来恢复原基带信号。因此，非相干解调一般只适用幅度调制（AM系统。忧郁包络解调器电路简单，效率高，所以几乎所有的幅度调 制（am接收机都采用这种电路。如下为串联型包络检波器的具体电路。oXr°D 1am0）n c * 叭o11O图2.2 AM信号的非相干解调原理当RC满足条件1 Wc乞RC乞1 Wh时，包络检波器的输出基本与输

14、入信号的包络变化呈线性关系，即m°(t)二 Ao m( t)其中，Ab >|m (t)max。隔去直流后就得到原信号 m(t)二、结果分析logt趣Idsre包络检波是在A0曰m(t)max的情况下才能实现的，由上图 A=2和A=1两种 情况下可以看出当A=2时满足条件A°3|m(t)max，AM波的包络与调制信号 m(t) 的形状完全一样，可以通过隔直流的方法(包络解调)得到原信号m(t);当A=1 时，不满足条件人占向(t)max出现欠调幅，AM波的包络失真无法通过隔直流的 方法(包络解调)得到原信号 m(t)，所以在这种情况下不能用包络检波的方法， 只能通过相干

15、解调。三、源代码%2014.11.1%AM调制信号的MATLAB实现%clear all;close all;%时间采样频谱dt=O.OOO1;fc=1O;%载波中心频率T=5;%信号时长N=T/dt;t=O: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*

16、fc*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其余时间返回零 电平在单极性归零码中，t /T称为占空比单极性归零码的主要优点是可以直 接提取同步信号，因此单极性归零码常常用作其他码型提取同步信号时的过渡码 型.也就是说其他适合信道传输但不能直接提取同步信号的码型，可先

17、变换为单极性归零码，然后再提取同步信号。二、结果分析单极性归零码10111101由上图可知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=O:1/tO:le ngth(x);for i=1:le ngth(x)if x(i)=1for j=1:t0/2y(2*i-2)*t0/2+j)

18、=1;y(2*i-1)*t0/2+j)=0;endelsefor j=1:t0y(i-1)*t0+j)=0;endendendy=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