课程设计报告模拟角度调制系统仿真.doc

青 岛 农 业 大 学信息科学与工程学院通信系统仿真课程设计报告 设 计 题 目 模拟角度调制系统仿真 学生专业班级 学生姓名（学号） 指 导 教 师 完 成 时 间 2013 实 习（设计）地点 信息楼机房 2013 年 6 月 17 日一、课程设计目的和任务1.1、设计目的(1) 、掌握调频与调相以及解调的基本原理2）、掌握模拟角度调制系统仿真，实现FM与PM调制与解调的设计方法，并比较两种调制的带宽、信噪比参数.(3)、理解模拟通信系统以及模拟调制在通信系统中的作用4）、调制与解调模拟系统，仿真实现相关功能实现相关功能设计fm、pm调制与解调模拟系统，仿真实现相关功能1.2、设计任务模拟角度调制系统仿真，实现FM与PM调制，要求信道为AWGN信道，画出调制信号、已调信号的波形图与频谱图，并比较两种调制的带宽、信噪比参数；二、分析与设计2.1.模拟通信系统简介通信系统是为了有效可靠的传输信息，信息由信源发出，以语言、图像、数据为媒体,通过电(光)信号将信息传输，由信宿接收。

通信系统又可分为数字通信与模拟通信角调制属于一类非线性调制方法，其中包括频率调制（FM）和相位调制（PM）模拟通信系统的性能指标:1:模拟信号占用的信道带宽必须等于或大于信号带宽，因此信号带宽越小，有效性越好2: 输出信噪比是按接收端输出信号的平均功率和噪声的平均功率之比在相同的条件下，系统的输出端的信噪比越大，则系统抗干扰的能力越大基于课程设计的要求，下面简要介绍模拟通信系统信源是模拟信号，信道中传输的也是模拟信号的系统为模拟通信模拟通信系统的模型如图1所示 图1 模拟通信系统模型调制器: 使信号与信道相匹配, 便于频分复用等发滤波器: 滤除调制器输出的无用信号收滤波器: 滤除信号频带以外的噪声，一般设N(t)为高斯白噪声,则Ni(t)为窄带白噪声 2.2 FM调制原理频率调制的一般表达式[1]为：（2-1）FM和PM非常相似，如果预先不知道调制信号的具体形式，则无法判断已调信号是调频信号还是调相信号 图 2-1 图 2-2图（2-1）所示的产生调频信号的方法称为直接调频法，图（2-2）所示的产生调频信号的方法称为间接调频法。

由于实际相位调制器的调节范围不可能超出 ，因而间接调频的方法仅适用于相位偏移和频率偏移不大的窄带调制情形，而直接调频则适用于宽带调制情形 为方便起见，无妨假设正弦载波的振幅A＝1，则由式（2-1）调频信号的一般表达式，得 （2-2）通过化解，利用傅立叶变化公式可得NBFM信号的频域表达式： （2-3）在NBFM中，由于下边频为负，因而合成矢量不与载波同相，而是存在相位偏移 ，当最大相位偏移满足式（2-2）时，合成矢量的幅度基本不变，这样就形成了FM信号图2-3 NBFM信号频谱2.3PM调制原理在模拟调制中，一个连续波有三个参数可以用来携带信息而构成已调信号当幅度和频率保持不变时，改变载波的相位使之随未调信号的大小而改变，这就是调相的概念角度调制信号的一般表示形式为： S (t)=Acos[ω t+φ(t)]式中，A是载波的恒定振幅；[ω t+φ(t)]是信号的瞬时相位，而φ(t)称为瞬时相位偏移；d[ω t+φ(t)]/dt为信号的瞬时频率，而dφ(t)/dt称为瞬时频率偏移，即相对于ω 的瞬时频率偏移设高频载波为u =U cosω t，调制信号为UΩ(t)，则调相信号的瞬时相位φ(t)=ω +K UΩ(t)瞬时角频率 ω(t)= =ω +K 调相信号 u =U cos［ω t+K uΩ(t)］ 将信号的信息加在载波的相位上则形成调相信号，调相的表达式为： S (t)=Acos[ω t+K f(t)+φ ]这里K 称为相移指数，这种调制方式，载波的幅度和角频率不变，而瞬时相位偏移是调制信号f(t)的线性函数，称为相位调制。

调相与调频有着相当密切的关系，我们知道相位与频率有如下关系式： ω= =ω +K f(t) φ(t)= ω t+K 所以在调相时可以先将调制信号进行微分后在进行频率调制，这样等效于调相，此方法 图2-4 PM调相信号的产生实现相位调制的基本原理是使角频率为ω 的高频载波u (t)通过一个可控相移网络, 此网络产生的相移Δφ受调制电压uΩ(t)控制, 满足Δφ=K uΩ(t)的关系, 所以网络输出就是调相信号，可控相移网络调相原理图如图2-7所示：可控相移网络调相原理图2.4解调过程分析：设输入调频信号为: 微分器的作用是把调频信号变成调幅调频波微分器输出为: 包络检波的作用是从输出信号的幅度变化中检出调制信号包络检波器输出为: 称为鉴频灵敏度（ ），是已调信号单位频偏对应的调制信号的幅度，经低通滤波器后加隔直流电容，隔除无用的直流，得 2.5高斯白噪声信道特性:设正弦波通过加性高斯白噪声信道后的信号为 其中，白噪声 的取值的概率分布服从高斯分布MATLAB本身自带了标准高斯分布的内部函数 函数产生的随机序列服从均值为 ，方差 的高斯分布正弦波通过加性高斯白噪声信道后的信号为 故其有用信号功率为: 噪声功率为: 信噪比 满足公式 则可得到公式 可以通过这个公式方便的设置高斯白噪声的方差。

设输入信号为m(t)=cos2*pi*t,载波中心频率fc=10Hz,调频器的压控震荡系数为5H／V，载波平均功率为1W三仿真结果及分析3.1仿真结果 3.2分析：调频（FM）指的是使载波的瞬时频率随着调制信号的大小而变，而幅度保持不变的调制方式由Figure1图可以得出当调制信号变小时载波的瞬时频率减小，从而得到的调频信号的频率变小，周期变大由Figure2图可以得出调频信号的频谱是对调制信号进行的频谱搬移将频谱搬移到所需较高频带fc上，从而满足信号的传输需要由Figure3图可以得出解调即调制的逆过程，从已调信号中不失真地恢复出原始基带信号，信号在信道中传输时会有噪声影响，在解调含有大信噪比高斯白噪声的已调信号比含有小信噪比高斯白噪声号的已调信号时的抗噪声性能比较好，恢复出的信号失真较小调相（PM）是用调制信号控制载波信号相位变化的一种信号变化方式当没有调制时，载波信号c(t)中的角频率w是常数4. 源代码清单Fmecho off close allclear allclcdt=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); %调制信号sn1=10; %设定信躁比(小信噪比)sn2=30; %设定信躁比(大信噪比)sn=0; %设定信躁比(无信噪比)db=am^2/(2*(10^(sn/10))); %计算对应的高斯白躁声的方差n=sqrt(db)*randn(size(t)); %生成高斯白躁声nsfm=n+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;ts=0.001; %抽样间隔fs=1/ts; %抽样频率df=0.25; 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; M=M/fs; %缩放，便于在频铺图上整体观察f=[0:df1:df1*(length(m)-1)]-fs/2; %时间向量对应的频率向量fs=1/ts;if nargin==2 n1=0;else n1=fs/df;endn2=leng。