《课程设计报告：AM、FM、PM实现及性能比较.doc》由会员分享，可在线阅读，更多相关《课程设计报告：AM、FM、PM实现及性能比较.doc（21页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上课 程 设 计 说 明 书学生姓名:学 号：学 院:信息工程学院班 级:通信081题 目:AM、FM、PM实现及性能比较指导教师： 职称: 教 2011 年 7 月 10 日一、设计目的： 分别设计AM，FM，PM调制与解调通信系统，分别在理信道和非理想信道中运行，并得出仿真结果，分析比较三种不同的调制方式所存在的不同。 熟悉MATLAB文件中M文件的使用方法，包括函数、原理和方法的应用。 增强在通信原理仿真方面的动手能力与自学能力。二、设计内容和实验要求 系统建模 确定仿真算法 建立仿真模型 设计仿真程序 运行仿真程序 输出仿真结果并进行分析三、模拟通信系统简介通信
2、系统是为了有效可靠的传输信息，信息由信源发出，以语言、图像、数据为媒体,通过电(光)信号将信息传输，由信宿接收。通信系统又可分为数字通信与模拟通信。基于课程设计的要求，下面简要介绍模拟通信系统。信源是模拟信号，信道中传输的也是模拟信号的系统为模拟通信。模拟通信系统的模型如图1所示。 图1 模拟通信系统模型调制器: 使信号与信道相匹配, 便于频分复用等。发滤波器: 滤除调制器输出的无用信号。收滤波器: 滤除信号频带以外的噪声，一般设N(t)为高斯白噪声,则Ni(t)为窄带白噪声。3.1AM，FM,PM调制原理3.1.1 AM调制原理幅度调制是用调制信号去控制高频正弦载波的幅度，使其按调制信号的规
3、律变化的过程。幅度调制器的一般模型如图2.1所示。图2.1 幅度调制模型在图2-1中，若假设滤波器为全通网络（1），调制信号叠加直流后再与载波相乘，则输出的信号就是常规双边带（AM）调幅 .AM调制器模型如图2-2所示图2.2 AM调制模型AM信号波形的包络与输入基带信号成正比，故用包络检波的方法很容易恢复原始调制信号。 但为了保证包络检波时不发生失真，必须满足，否则将出现过调幅现象而带来失真。AM信号的频谱是由载频分量和上、下两个边带组成（通常称频谱中画斜线的部分为上边带，不画斜线的部分为下边带）。上边带的频谱与原调制信号的频谱结构相同，下边带是上边带的镜像。显然，无论是上边带还是下边带，都
4、含有原调制信号的完整信息。故AM信号是带有载波的双边带信号，它的带宽信号带宽的两倍。3.2.2 FM调制原理 频率调制的一般表达式1为： （2-1）FM和PM非常相似，如果预先不知道调制信号的具体形式，则无法判断已调信号是调频信号还是调相信号。 图 2-3 图 2-4 图（2-3）所示的产生调频信号的方法称为直接调频法，图（2-4）所示的产生调频信号的方法称为间接调频法4。由于实际相位调制器的调节范围不可能超出，因而间接调频的方法仅适用于相位偏移和频率偏移不大的窄带调制情形，而直接调频则适用于宽带调制情形。 根据调制后载波瞬时相位偏移的大小，可将频率调制分为宽带调频（WBFM）与窄带调频（NB
5、FM）。宽带与窄带调制的区分并无严格的界限，但通常认为由调频所引起的最大瞬时相位偏移远小于30时， （2-2）称为窄带调频。否则，称为宽带调频。 为方便起见，无妨假设正弦载波的振幅A1，则由式（2-1）调频信号的一般表达式，得= （2-3）通过化解，利用傅立叶变化公式可得NBFM信号的频域表达式：（2-4） 在NBFM中，由于下边频为负，因而合成矢量不与载波同相，而是存在相位偏移，当最大相位偏移满足式（2-2）时，合成矢量的幅度基本不变，这样就形成了FM信号。 图2-5 NBFM信号频谱3.2.3 PM调制原理在模拟调制中，一个连续波有三个参数可以用来携带信息而构成已调信号。当幅度和频率保持不
6、变时，改变载波的相位使之随未调信号的大小而改变，这就是调相的概念。角度调制信号的一般表示形式为： S (t)=Acost+(t)式中，A是载波的恒定振幅；t+(t)是信号的瞬时相位，而(t)称为瞬时相位偏移；dt+(t)/dt为信号的瞬时频率，而d(t)/dt称为瞬时频率偏移，即相对于的瞬时频率偏移。设高频载波为u=Ucost，调制信号为U(t)，则调相信号的瞬时相位(t)=+KU(t)瞬时角频率 (t)=+K调相信号 u=Ucost+Ku(t) 将信号的信息加在载波的相位上则形成调相信号，调相的表达式为： S(t)=Acost+Kf(t)+这里K称为相移指数，这种调制方式，载波的幅度和角频率
7、不变，而瞬时相位偏移是调制信号f(t)的线性函数，称为相位调制。调相与调频有着相当密切的关系，我们知道相位与频率有如下关系式： =+Kf(t) (t)=t+K所以在调相时可以先将调制信号进行微分后在进行频率调制，这样等效于调相，此方法称为间接调相，与此相对应，上述方法称为直接调相。调相信号的产生如图2-6所示：图2-6 PM调相信号的产生实现相位调制的基本原理是使角频率为的高频载波u(t)通过一个可控相移网络, 此网络产生的相移受调制电压u(t)控制, 满足=Ku(t)的关系, 所以网络输出就是调相信号，可控相移网络调相原理图如图2-7所示：图2-7 可控相移网络调相原理图四、几种调制方式间的
8、比较4.1PM与FM的比较FM(1) 瞬时频率：(2) 瞬时相位：(3) 最大频偏：(4) 最大相位：表达式：PM4.2几种不同的模拟调制方式假定所有调制系统在接收机输入端具有相等的信号功率，且加性噪声都是均值为0、双边功率谱密度为/2的高斯白噪声，基带信号带宽为，在所有系统都满足 例如，为正弦型信号。综合前面的分析，可总结各种模拟调制方式的信号带宽、制度增益、输出信噪比、设备（调制与解调）复杂程度、主要应用等如下表所示。表中还进一步假设了AM为100%调制。4.3几种模拟调制的性能比较就抗噪性能而言，WBFM最好，DSB、SSB、VSB次之，AM最差。NBFM与AM接近。图3-33示出了各种
9、模拟调制系统的性能曲线，图中的圆点表示门限点。门限点以下，曲线迅速下跌；门限点以上，DSB、SSB的信噪比比AM高4.7dB以上，而FM（6）的信噪比比AM高22dB。就频带利用率而言，SSB最好，VSB与SSB接近，DSB、AM、NBFM次之，WBFM最差。4.4几种模拟调制的特点及应用AM调制的优点是接收设备简单；缺点是功率利用率低，抗干扰能力差，信号带宽较宽，频带利用率不高。因此，AM制式用于通信质量要求不高的场合，目前主要用在中波和短波的调幅广播中。 DSB调制的优点是功率利用率高，但带宽与AM相同，频带利用率不高，接收要求同步解调，设备较复杂。只用于点对点的专用通信及低带宽信号多路复
10、用系统。 SSB调制的优点是功率利用率和频带利用率都较高，抗干扰能力和抗选择性衰落能力均优于AM，而带宽只有AM的一半；缺点是发送和接收设备都复杂。SSB制式普遍用在频带比较拥挤的场合，如短波波段的无线电广播和频分多路复用系统中。 VSB调制性能与SSB相当，原则上也需要同步解调，但在某些VSB系统中，附加一个足够大的载波，形成（VSB+C）合成信号，就可以用包络检波法进行解调。这种（VSB+C）方式综合了AM、SSB和DSB三者的优点。所以VSB在数据传输、商用电视广播等领域得到广泛使用。 FM波的幅度恒定不变，这使得它对非线性器件不甚敏感，给FM带来了抗快衰落能力。利用自动增益控制和带通限
11、幅还可以消除快衰落造成的幅度变化效应。这些特点使得NBFM对微波中继系统颇具吸引力。WBFM的抗干扰能力强，可以实现带宽与信噪比的互换，因而WBFM广泛应用于长距离高质量的通信系统中，如空间和卫星通信、调频立体声广播、短波电台等。WBFM的缺点是频带利用率低，存在门限效应，因此在接收信号弱、干扰大的情况下宜采用NBFM，这就是小型通信机常采用NBFM的原因。五、matlab生成图像AMFMPM程序清单：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);figur
12、e(1)subplot(1,1,1),plot(t,yt),title(频率为3000的调制信号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);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,y
13、c),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)subplot(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
14、(1:N/2);subplot(2,1,2),plot(q2,mx2),title(高斯白噪声频域波形)y6=y2+y4;figure(5)subplot(2,1,1),plot(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,
15、beta); %使用kaiser窗函数b=fir1(k,Wn,ftype,window,noscale); %使用标准频率响应的加窗设计函数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/r
16、o G=r/R fm：dt=0.001; %设定时间步长t=0:dt:1.5; %产生时间向量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); %调制信号%*%*添加高斯白噪
17、声*sn1=10; %设定信躁比(小信噪比)sn2=30; %设定信躁比(大信噪比)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变换，求绝对值
18、得到瞬时幅度（包络检波）zero=(max(diff_nsfmn)-min(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=f
19、ft(m,n);m=m,zeros(1,n-n2);df1=fs/n; %以上程序是对调制后的信号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; %缩放%*%*
20、disp(按任意键可以看到原调制信号、载波信号和已调信号的曲线)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(按任意键可以看到原调制信号和已调信号在频域内的图形)pause%*figure(
21、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(按任意键可以看到原调制信号、无噪声条件下已调信号和解调信号的曲线)pause%*figure(3)*figure(3)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel(时间t);title(调制信号的时域图);subp
22、lot(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(hilbert(diff_nsfm); %hilbert变换，求绝对值得到瞬时幅度（包络检波）zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn-zero;subplot(3,1,3); %
23、绘制无噪声条件下解调信号的时域图plot(1:length(diff_nsfmn1)./1000,diff_nsfmn1./400,r);xlabel(时间t); title(无噪声条件下解调信号的时域图);%*disp(按任意键可以看到原调制信号、小信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线)pause%*figure(4)*figure(4)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel(时间t);title(调制信号的时域图);db1=am2/(2*(10(sn1/10); %计算对应的小信噪比高斯白躁声的方差n1=sqrt(db1)
24、*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_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn1-zero;subplot(3,1,2);plot(1:lengt
25、h(diff_nsfm),diff_nsfm); %绘制含小信噪比高斯白噪声已调信号的时域图xlabel(时间t);title(含小信噪比高斯白噪声已调信号的时域图);subplot(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
26、,1);plot(t,mt); %绘制调制信号的时域图xlabel(时间t);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变换，求绝对值
27、得到瞬时幅度（包%络检波）zero=(max(diff_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(时间t);
28、title(含大信噪比高斯白噪声解调信号的时域图);%*结 束* pm：function 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/df;endn2=lengt
29、h(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的正交分量，t0=0.2; %信号的持续时间,用来定义时间
30、向量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_m(i+1)=int_m(i)+m(i)*ts;endM,m,df1=fftseq(m,ts,df); %对调制信号m(t)求傅里叶变换M=M/fs; %缩放，便于在频谱图上整体观察f=
31、0:df1:df1*(length(m)-1)-fs/2; %时间向量对应的频率向量u=cos(2*pi*fc*t+2*pi*kf*int_m); %调制后的信号U,u,df1=fftseq(u,ts,df); %对调制后的信号u求傅里叶变换U=U/fs; %缩放%通过调用子程序env_phas和loweq来实现解调功能v,phase=env_phas(u,ts,fc); %解调，求出u的相位phi=unwrap(phase); %校正相位角，使相位在整体上连续，便于后面对该相位角求导dem=(1/(2*pi*kf)*(diff(phi)*fs); %对校正后的相位求导 %再经一些线性变换来恢
32、复原调制信号 %乘以fs是为了恢复原信号，因为前面使用了缩放 subplot(3,2,1) %子图形式显示结果plot(t,m(1:length(t) %现在的m信号是重新构建的信号， %因为在对m求傅里叶变换时m=m,zeros(1,n-n2)axis(-0.1 0.1 -1 1) %定义两轴的刻度xlabel(时间t) title(原调制信号的时域图)subplot(3,2,2)plot(t,u(1:length(t)axis(-0.1 0.1 -1 1)xlabel(时间t)title(已调信号的时域图)subplot(3,2,3)plot(f,abs(fftshift(M) %ffts
33、hift：将FFT中的DC分量移到频谱中心axis(-600 600 0 0.04)xlabel(频率f)title(原调制信号的频谱图)subplot(3,2,4)plot(f,abs(fftshift(U)axis(-600 600 0 0.04)xlabel(频率f)title(已调信号的频谱图)subplot(3,2,5)plot(t,m(1:length(t)axis(-0.1 0.1 -1 1)xlabel(时间t)title(原调制信号的时域图)subplot(3,2,6)plot(t,dem(1:length(t)axis(-0.1 0.1 -1 1)xlabel(时间t)ti
34、tle(解调后信号的时域波形)六、心得体会通过这一次课程设计，我了解很多关于专业的知识，以前每次学这些知识时，总是不知道这些东西具体拿来有什么用，现在才知道，几个短短输入信号，在有了一个简单的电路流程后，就能仿真成我们生活中很多常见的东西。总的来说，这次课程设计过程还是比较愉快轻松的，虽然中间有过一些困难，但是在老师与同学的指点下我还是渡过了，在这里我要谢谢帮助我的老师和同学。参考文献1 樊昌信. 通信原理（第6版）.国防工业出版社，2006，092 黎洪松. 数字通信原理.西安电子系科技大学出版社，2005，073 任嘉伟. 数字频带通信系统计算机仿真J.电脑知识与技术，2008，074 吕跃广 通信系统仿真.电子工业出版社,2010.035 席在芳等 基于SIMULINK 的现代通信系统仿真分析J.系统仿真学报 2006,18(10)专心-专注-专业