模拟调制与解调算法研究

1、模拟调制与解调算法研究 科研训练论文实验总成绩： 装 订 线报告份数： 西安邮电大学 通信与信息工程学院 科研训练论文专业班级： 学生姓名： 学号(班内序号): 2013 年 9 月 22 日模拟调制与解调算法研究Research of analog modulation and demodulation algorithm摘要：对于模拟调制解调，我们选取了三种典型的方式，AM,PM,FM来进行研究。AM是通过用调制信号去控制高频正弦载波的幅度，使其按照调制信号的规律变化。其实现较为简单，但是带宽利用率不高。FM和PM非常相似，如果预先不知道调制信号的具体形式，则无法判断已调信号是调频信号还是

2、调相信号。FM是通过调制信号来控制载波频率来传递信号，FM带宽利用率高，抗干扰能力强，但实现起来比较困难，会出现门限效应，PM是通过调制信号控制载波相位来专递信号。Summary：For analog modulation and demodulation algorithm ,We chose AM, FM, PM three kinds of typical analog modulation method are discussed. AM is to control the carrier amplitude modulation signal by signal transmissi

3、on, Make it according to the law of the modulation signal to change. it is implement, but the bandwidth utilization rate is not high.FM and AM is very similiar,if you dont knoe inadvance the particular form of modulation signal is impossible to judge is modulated signal frequency modulation signal o

4、r phrase modulation signal. FM is to control the carrier frequency modulation signal to pass through the frequency of the carrier signal, FM has the high bandwidth utilization, strong anti-interference ability, but more difficult to achieve, in the demodulation occurs when the threshold effect, PM i

5、s to express signal through the modulation signal to control the carrier phase,关键词：AM FM PM 模拟 调制 解调 调制信号 载波Key word:AM FM PM analog modulation demodulation modulating signal carrier wave引言：通信系统是为了有效可靠的传输信息，信息由信源发出，以语言、图像、数据为媒体,通过电(光)信号将信息传输，由信宿接收。通信系统又可分为数字通信与模拟通信。信源是模拟信号，信道中传输的也是模拟信号的系统为模拟通信。模拟通信系

6、统的模型是模拟信源经过模拟调制器，再经过发送滤波器发送到信道，经过信道的传输和干扰到达接收端，经过接收滤波器和解调器进行解调，最后到达信宿。一、AM，FM，PM调制原理：1.AM调制原理幅度调制是用调制信号去控制高频正弦载波的幅度，使其按调制信号的规律变化的过程。幅度调制器的一般模型如图1所示。图1 幅度调制模型在图2.1中，若假设滤波器为全通网络（H（）1），调制信号叠加直流后再与载波相乘，则输出的信号就是常规双边带（AM）调幅 .AM调制器模型如图2所示：图2 AM调制模型AM信号波形的包络与输入基带信号成正比，故用包络检波的方法很容易恢复原始调制信号。 但为了保证包络检波时不发生失真，必

7、须满足，否则将出现过调幅现象而带来失真。AM信号的频谱是由载频分量和上、下两个边带组成（通常称频谱中画斜线的部分为上边带，不画斜线的部分为下边带）。上边带的频谱与原调制信号的频谱结构相同，下边带是上边带的镜像。显然，无论是上边带还是下边带，都含有原调制信号的完整信息。故AM信号是带有载波的双边带信号，它的带宽信号带宽的两倍。 2.FM调制原理 频率调制的一般表达式1为： (2-1) FM和PM非常相似，如果预先不知道调制信号的具体形式，则无法判断已调信号是调频信号还是调相信号。 图 3 直接调频法 图 4 间接调频法图3所示的产生调频信号的方法称为直接调频法，图4所示的产生调频信号的方法称为间

8、接调频法4。由于实际相位调制器的调节范围不可能超出，因而间接调频的方法仅适用于相位偏移和频率偏移不大的窄带调制情形，而直接调频则适用于宽带调制情形。 根据调制后载波瞬时相位偏移的大小，可将频率调制分为宽带调频（WBFM）与窄带调频（NBFM）。宽带与窄带调制的区分并无严格的界限，但通常认为由调频所引起的最大瞬时相位偏移远小于30°时，   （2-2）称为窄带调频。否则，称为宽带调频。 为方便起见，无妨假设正弦载波的振幅A1，则由式调频信号的一般表达式，得= (2-3) 通过化解，利用傅立叶变化公式可得NBFM信号的频域表达式： （2-4） 在NBFM中，由于下边频为负，因而

9、合成矢量不与载波同相，而是存在相位偏移，当最大相位偏移满足式（2-2）时，合成矢量的幅度基本不变，这样就形成了FM信号。 图5 NBFM信号频谱 3.PM调制原理在模拟调制中，一个连续波有三个参数可以用来携带信息而构成已调信号。当幅度和频率保持不变时，改变载波的相位使之随未调信号的大小而改变，这就是调相的概念。角度调制信号的一般表示形式为： S (t)=Acost+(t) (2-5)式中，A是载波的恒定振幅；t+(t)是信号的瞬时相位，而(t)称为瞬时相位偏移；dt+(t)/dt为信号的瞬时频率，而d(t)/dt称为瞬时频率偏移，即相对于的瞬时频率偏移。设高频载波为u=Ucost，调制信号为U

10、(t)，则调相信号的瞬时相位(t)=+KU(t)瞬时角频率 (t)=+K调相信号 u=Ucost+Ku(t) 将信号的信息加在载波的相位上则形成调相信号，调相的表达式为： S(t)=Acost+Kf(t)+ (2-6)这里K称为相移指数，这种调制方式，载波的幅度和角频率不变，而瞬时相位偏移是调制信号f(t)的线性函数，称为相位调制。调相与调频有着相当密切的关系，我们知道相位与频率有如下关系式： =+Kf(t) (t)=t+K所以在调相时可以先将调制信号进行微分后在进行频率调制，这样等效于调相，此方法称为间接调相，与此相对应，上述方法称为直接调相。调相信号的产生如图6所示：图6 PM调相信号的产

11、生实现相位调制的基本原理是使角频率为的高频载波u(t)通过一个可控相移网络, 此网络产生的相移受调制电压u(t)控制, 满足=Ku(t)的关系, 所以网络输出就是调相信号，可控相移网络调相原理图如图7所示：图7 可控相移网络调相原理图二、几种模拟调制的性能比较：就抗噪性能而言，WBFM最好，DSB、SSB、VSB次之，AM最差。NBFM与AM接近。示出了各种模拟调制系统的性能曲线，图中的圆点表示门限点。门限点以下，曲线迅速下跌；门限点以上，DSB、SSB的信噪比比AM高4.7dB以上，而 FM（6）的信噪比比AM高22dB。就频带利用率而言，SSB最好，VSB与SSB接近，DSB、AM、NBF

12、M次之，WBFM最差。几种模拟调制的特点及应用AM调制的优点是接收设备简单；缺点是功率利用率低，抗干扰能力差，信号带宽较宽，频带利用率不高。因此，AM制式用于通信质量要求不高的场合，目前主要用在中波和短波的调幅广播中。 DSB调制的优点是功率利用率高，但带宽与AM相同，频带利用率不高，接收要求同步解调，设备较复杂。只用于点对点的专用通信及低带宽信号多路复用系统。 SSB调制的优点是功率利用率和频带利用率都较高，抗干扰能力和抗选择性衰落能力均优于AM，而带宽只有AM的一半；缺点是发送和接收设备都复杂。SSB制式普遍用在频带比较拥挤的场合，如短波波段的无线电广播和频分多路复用系统中。 VSB调制性

13、能与SSB相当，原则上也需要同步解调，但在某些VSB系统中，附加一个足够大的载波，形成（VSB+C）合成信号，就可以用包络检波法进行解调。这种（VSB+C）方式综合了AM、SSB和DSB三者的优点。所以VSB在数据传输、商用电视广播等领域得到广泛使用。 FM波的幅度恒定不变，这使得它对非线性器件不甚敏感，给FM带来了抗快衰落能力。利用自动增益控制和带通限幅还可以消除快衰落造成的幅度变化效应。这些特点使得NBFM对微波中继系统颇具吸引力。WBFM的抗干扰能力强，可以实现带宽与信噪比的互换，因而WBFM广泛应用于长距离高质量的通信系统中，如空间和卫星通信、调频立体声广播、短波电台等。WBFM的缺点

14、是频带利用率低，存在门限效应，因此在接收信号弱、干扰大的情况下宜采用NBFM，这就是小型通信机常采用NBFM的原因。三、AM,FM,PM的调制仿真 1.AM的调制仿真2.FM的调制仿真 3.PM的调制仿真 参考文献1 樊昌信. 通信原理（第6版）.国防工业出版社，2006，092 黎洪松. 数字通信原理.西安电子系科技大学出版社，2005，073 任嘉伟. 数字频带通信系统计算机仿真J.电脑知识与技术，2008，074 吕跃广 通信系统仿真.电子工业出版社,2010.035 席在芳等 基于SIMULINK 的现代通信系统仿真分析J.系统仿真学报 2006,18(10)6 韩利竹 王华 Matl

15、ab通信仿真与应用 国防工业出版社 20037 谢沅清 邓刚 通信电子电路 电子工业出版社 2005附 录AM ：fm=100;fc=500;fs=5000;Am=1;A=2;N=512;K=N-1;n=0:N-1;t=(0:1/fs:K/fs);yt=Am*cos(2*pi*fm*t);figure(1)subplot(1,1,1),plot(t,yt),title('调制信号f1的时时域波');y0=A+yt ;y2=y0.*cos(2*pi*fc*n/fs);y3=fft(y2,N);% fft 变换 q1=(0:N/2-1)*fs/N;mx1=abs(y3(1:N/2)

16、;figure(2)subplot(2,1,1);plot(t,y2);title('已调信号的时时域波');subplot(2,1,2);plot(q1,mx1);title('f1已调信号的频谱'); %绘图yc=cos(2*pi*fc*t);figure(3)subplot(2,1,1),plot(t,yc),title('载波fc时域波形')N=512;n=0:N-1;yc1=Am*cos(2*pi*fc*n/fs);y3=fft(yc1,N);q=(0:N/2-1)*fs/N;mx=abs(y3(1:N/2);figure(3)subp

17、lot(2,1,2),plot(q,mx),title('载波fc频谱')y4=0.01*randn(1,length(t);%用RANDN产生高斯分布序列 w=y4.2; %噪声功率figure(4)subplot(2,1,1);plot(t,y4);title('高斯白噪声时域波形')y5=fft(y4,N);q2=(0:N/2-1)*fs/N;mx2=abs(y5(1:N/2);subplot(2,1,2),plot(q2,mx2),title('高斯白噪声频域波形')y6=y2+y4;figure(5)subplot(2,1,1),plo

18、t(t,y6),title('叠加后的调制信号时域波形')q3=q1;mx3=mx1+mx2;subplot(2,1,2),plot(q3,mx3),title('叠加后的调制信号频谱波形') %调制yv=y6.*yc; %乘以载波进行解调Ws=yv.2;p1=fc-fm;k,Wn,beta,ftype=kaiserord(p1 fc,1 0,0.05 0.01,fs); %Fir数字低通滤波window=kaiser(k+1,beta); %使用kaiser窗函数b=fir1(k,Wn,ftype,window,'noscale'); %使用标

19、准频率响应的加窗设计函数yt=filter(b,1,yv);yssdb=yt.*2-2;figure(6)subplot(2,1,1),plot(t,yssdb),title('经过低通已调信号的时域波形采样')y9=fft(yssdb,N);q=(0:N/2-1)*fs/N;mx=abs(y9(1:N/2);subplot(2,1,2),plot(q,mx),title('经过低通已调信号频域波形') %解调 ro=y9-yt; W=(yt.2).*(1/2); R=W/w r=W/ro G=r/R fm：dt=0.001; %设定时间步长t=0:dt:1.5

20、; %产生时间向量am=5; %设定调制信号幅度fm=5; %设定调制信号频率mt=am*cos(2*pi*fm*t); %生成调制信号fc=50; %设定载波频率ct=cos(2*pi*fc*t); %生成载波kf=10; %设定调频指数int_mt(1)=0;for i=1:length(t)-1 int_mt(i+1)=int_mt(i)+mt(i)*dt; %求信号m(t)的积分end %调制，产生已调信号sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt); %调制信号%*%*添加高斯白噪声*sn1=10; %设定信躁比(小信噪比)sn2=30; %设定信躁比(大

21、信噪比)sn=0; %设定信躁比(无信噪比)db=am2/(2*(10(sn/10); %计算对应的高斯白躁声的方差n=sqrt(db)*randn(size(t); %生成高斯白躁声nsfm=n+sfm; %生成含高斯白躁声的已调信号（信号通%过信道传输）%*%*FM解调*for i=1:length(t)-1 %接受信号通过微分器处理 diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn = abs(hilbert(diff_nsfm); %hilbert变换，求绝对值得到瞬时幅度（包络检波）zero=(max(diff_nsfmn)-min(

22、diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn-zero;%*%*时域到频域转换*ts=0.001; %抽样间隔fs=1/ts; %抽样频率df=0.25; %所需的频率分辨率，用在求傅里叶变换%时，它表示FFT的最小频率间隔%*对调制信号m(t)求傅里叶变换*m=am*cos(2*pi*fm*t); %原调信号fs=1/ts;if nargin=2 n1=0;else n1=fs/df;endn2=length(m);n=2(max(nextpow2(n1),nextpow2(n2);M=fft(m,n);m=m,zeros(1,n-n2);df1=fs/n; %以

23、上程序是对调制后的信号u求傅里变换M=M/fs; %缩放，便于在频铺图上整体观察f=0:df1:df1*(length(m)-1)-fs/2; %时间向量对应的频率向量%*对已调信号u求傅里变换*fs=1/ts;if nargin=2 n1=0;else n1=fs/df;endn2=length(sfm);n=2(max(nextpow2(n1),nextpow2(n2);U=fft(sfm,n);u=sfm,zeros(1,n-n2);df1=fs/n; %以上是对已调信号u求傅里变换U=U/fs; %缩放%*%*disp('按任意键可以看到原调制信号、载波信号和已调信号的曲线&#

24、39;)pause%*figure(1)*figure(1)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel('时间t');title('调制信号的时域图');subplot(3,1,2);plot(t,ct); %绘制载波的时域图xlabel('时间t');title('载波的时域图');subplot(3,1,3);plot(t,sfm); %绘制已调信号的时域图xlabel('时间t');title('已调信号的时域图');%*disp('按任意

25、键可以看到原调制信号和已调信号在频域内的图形')pause%*figure(2)*figure(2)subplot(2,1,1)plot(f,abs(fftshift(M) %fftshift:将FFT中的DC分量移到频谱中心xlabel('频率f')title('原调制信号的频谱图')subplot(2,1,2)plot(f,abs(fftshift(U)xlabel('频率f')title('已调信号的频谱图')%*disp('按任意键可以看到原调制信号、无噪声条件下已调信号和解调信号的曲线')paus

26、e%*figure(3)*figure(3)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel('时间t');title('调制信号的时域图');subplot(3,1,2);plot(t,sfm); %绘制已调信号的时域图xlabel('时间t');title('无噪声条件下已调信号的时域图');nsfm=sfm; for i=1:length(t)-1 %接受信号通过微分器处理 diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn = abs(

27、hilbert(diff_nsfm); %hilbert变换，求绝对值得到瞬时幅度（包络检波）zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn-zero;subplot(3,1,3); %绘制无噪声条件下解调信号的时域图plot(1:length(diff_nsfmn1)./1000,diff_nsfmn1./400,'r');xlabel('时间t'); title('无噪声条件下解调信号的时域图');%*disp('按任意键可以看到原调制信号、小信噪比高斯白噪

28、声条件下已调信号和解调信号已调信号的曲线')pause%*figure(4)*figure(4)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel('时间t');title('调制信号的时域图');db1=am2/(2*(10(sn1/10); %计算对应的小信噪比高斯白躁声的方差n1=sqrt(db1)*randn(size(t); %生成高斯白躁声nsfm1=n1+sfm; %生成含高斯白躁声的已调信号（信号通%过信道传输）for i=1:length(t)-1 %接受信号通过微分器处理 diff_nsfm1(i

29、)=(nsfm1(i+1)-nsfm1(i)./dt;enddiff_nsfmn1 = abs(hilbert(diff_nsfm1); %hilbert变换，求绝对值得到瞬时幅度（包络检波）zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn1-zero;subplot(3,1,2);plot(1:length(diff_nsfm),diff_nsfm); %绘制含小信噪比高斯白噪声已调信号的时域图xlabel('时间t');title('含小信噪比高斯白噪声已调信号的时域图');subp

30、lot(3,1,3); %绘制含小信噪比高斯白噪声解调信号的时域图plot(1:length(diff_nsfmn1)./1000,diff_nsfmn1./400,'r');xlabel('时间t'); title('含小信噪比高斯白噪声解调信号的时域图');%*disp('按任意键可以看到原调制信号、大信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线')pause%*figure(5)*figure(5)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel('时间t')

31、;title('调制信号的时域图');db1=am2/(2*(10(sn2/10); %计算对应的大信噪比高斯白躁声的方差n1=sqrt(db1)*randn(size(t); %生成高斯白躁声nsfm1=n1+sfm; %生成含高斯白躁声的已调信号（信号通过信道传输）for i=1:length(t)-1 %接受信号通过微分器处理 diff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i)./dt;enddiff_nsfmn1 = abs(hilbert(diff_nsfm1); %hilbert变换，求绝对值得到瞬时幅度（包%络检波）zero=(max(diff_

32、nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn1-zero;subplot(3,1,2);plot(1:length(diff_nsfm1),diff_nsfm1); %绘制含大信噪比高斯白噪声已调信号%的时域图xlabel('时间t');title('含大信噪比高斯白噪声已调信号的时域图');subplot(3,1,3); %绘制含大信噪比高斯白噪声解调信号%的时域图plot(1:length(diff_nsfmn1)./1000,diff_nsfmn1./400,'r');xlabel('

33、;时间t'); title('含大信噪比高斯白噪声解调信号的时域图');PMfunction v,phi=env_phas(x,ts,f0)if nargout=2 %nargout为输出变数的个数 z=loweq(x,ts,f0); %产生调制信号的正交分量 phi=angle(z); %angle是对一个复数求相角的函数endv=abs(hilbert(x); %abs用来求复数hilbert(x)的模function M,m,df=fftseq(m,ts,df)fs=1/ts;if nargin=2 n1=0; %nargin为输入参量的个数else n1=fs/

34、df;endn2=length(m);n=2(max(nextpow2(n1),nextpow2(n2); %nextpow2(n)取n最接近的较大2次幂M=fft(m,n); %M为信号m的傅里叶变换，n为快速傅里叶变换的点数，及基n-FFT变换m=m,zeros(1,n-n2); %构建新的m信号df=fs/n; %重新定义频率分辨率function x1=loweq(x,ts,f0)t=0:ts:ts*(length(x)-1);z=hilbert(x); %希尔伯特变换对的利用-通过实部来求虚部x1=z.*exp(-j*2*pi*f0*t); %产生信号z的正交分量， %并将z信号与它的正交分量加在一起%主程序t0=0.2; %信号的持续时间,用来定义时间向量ts=0.001; %抽样间隔fs=1/ts; %抽样频率fc=300; %载波频率,fc可以任意改变t=-t0/2:ts:t0/2; %时间向量kf=100; %偏差常数df=0.25; %所需的频率分辨率，用在求傅里叶变换时，它表示FFT的最小频率间隔m=sin(100*t); %调制信号，m(t)可以任意更改int_m(1)=0; %求信号m(t)的积分for i=1:length(t)-1 int_