课 程 设 计
课程设计名称:调频信号的解调及抗噪声性能分析
专 业 班 级 : ___ __________
学 生 姓 名 : _______________
学 号 :
_________
指 导 教 师 : ____ ___________
课程设计时间:_____
1 需求分析
应用 MATLAB 语言编写基本的通信系统的应用程序,进行模拟调制系统,数
字基带信号的传输系统的建模、设计与仿真。所有的仿真用 MATLAB 程序实现,
系统经过的信道都假设为高斯白噪声信道,并用程序画出调制信号,载波,已
调信号、解调信号的波形。
2 概要设计
通信的目的是传输信息。通信系统的作用就是将信息从信息源发送到一个
或多个目的地。对于任何一个通信系统,均可视为由发送端、信道和接收端三
大部分组成(如图 1 所示)。
发送端
信道
接收端
信息源
发送设
信 道
接受设
信息源
噪声源
图 1 通信系统一般模型
信息源(简称信源)的作用是把各种信息转换成原始信号。根据消息的种
类不同信源分为模拟信源和数字信源。发送设备的作用产生适合传输的信号,
即使发送信号的特性和信道特性相匹配,具有抗噪声的能力,并且具有足够的
功率满足原距离传输的需求。
信息源和发送设备统称为发送端。
发送端将信息直接转换得到的较低频率的原始电信号称为基带信号。通常
基带信号不宜直接在信道中传输。因此,在通信系统的发送端需将基带信号的
频谱搬移(调制)到适合信道传输的频率范围内进行传输。这就是调制的过
程。
信号通过信道传输后,具有将信号放大和反变换功能的接收端将已调制的
信号搬移(解调)到原来的频率范围,这就是解调的过程。
信号在信道中传输的过程总会受到噪声的干扰,通信系统中没有传输信号
时也有噪声,噪声永远存在于通信系统中。由于这样的噪声是叠加在信号上
的,所以有时将其称为加性噪声。噪声对于信号的传输是有害的,它能使模拟
信号失真。在本仿真的过程中我们假设信道为高斯白噪声信道。
调制在通信系统中具有十分重要的作用。一方面,通过调制可以把基带信
号的频谱搬移到所希望的位置上去,从而将调制信号转换成适合于信道传输或
便于信道多路复用的已调信号。另一方面,通过调制可以提高信号通过信道传
输时的抗干扰能力,同时,它还和传输效率有关。具体地讲,不同的调制方式
产生的已调信号的带宽不同,因此调制影响传输带宽的利用率。可见,调制方
式往往决定一个通信系统的性能。在本仿真的过程中我们选择用调频调制方法
进行调制。
调制过程是一个频谱搬移的过程,它是将低频信号的频谱搬移到载频位
置。而解调是将位于载频的信号频谱再搬回来,并且不失真地恢复出原始基带
信号。在本仿真的过程中我们选择用非相干解调方法进行解调。
3 运行环境
硬件环境:Windows XP
软件环境:MATLAB
4 详细设计
FM 调制模型的建立
图 2
FM 调制模型
其中, ( )m t 为基带调制信号,设调制信号为
)m
f t
cos(2
( )
m t
A
设正弦载波为
( )
c t
c o s ( 2
)c
f t
信号传输信道为高斯白噪声信道,其功率为 2 。
调制过程分析
在调制时,调制信号的频率去控制载波的频率的变化,载波的瞬时频偏随
d
调制信号 ( )m t 成正比例变化,即
(
d t
rad
式中, fK 为调频灵敏度( (
这时相位偏移为
s V )。
)
t
)
(
K m t
f
)
( )
t
则可得到调频信号为
( )
K m d
f
s
FM
( )
t A
cos
( )
t K m d
c
f
调制信号产生的 M 文件:
dt=0.001;
t=0:dt:1.5;
am=15;
fm=15;
%设定时间步长
%产生时间向量
%设定调制信号幅度←可更改
%设定调制信号频率←可更改
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;
end
%设定载波频率←可更改
%生成载波
%设定调频指数
%对 mt 进行积分
sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);
%调制,产生已调信号
FM 解调模型的建立
调制信号的解调分为相干解调和非相干解调两种。相干解调仅仅适用于窄
带调频信号,且需同步信号,故应用范围受限;而非相干解调不需同步信号,
且对于 NBFM 信号和 WBFM 信号均适用,因此是 FM 系统的主要解调方式。
在本仿真的过程中我们选择用非相干解调方法进行解调。
图 4 FM 解调模型
非相干解调器由限幅器、鉴频器和低通滤波器等组成,其方框图如图 5 所
示。限幅器输入为已调频信号和噪声,限幅器是为了消除接收信号在幅度上可
能出现的畸变;带通滤波器的作用是用来限制带外噪声,使调频信号顺利通
过。鉴频器中的微分器把调频信号变成调幅调频波,然后由包络检波器检出包
络,最后通过低通滤波器取出调制信号。
解调过程分析
设输入调频信号为
( )
S t
t
S
FM
( )
t
A
cos(
t K
c
f
t
( )
)
m d
微分器的作用是把调频信号变成调幅调频波。微分器输出为
( )
S t
d
c
t K
f
t
m
( )
)
d
( )
t
dS
( )
dS t
i
dt
K m t
FM
dt
( ) sin(
c
f
包络检波的作用是从输出信号的幅度变化中检出调制信号。包络检波器输
出为
o
S
( )
t
dK 称为鉴频灵敏度(V
K
d
c
( )
K m t
f
K
c
d
( )
K K m t
d
f
Hz ),是已调信号单位频偏对应的调制信号的幅
度,经低通滤波器后加隔直流电容,隔除无用的直流,得
(
m t
o
)
(
K K m t
d
f
)
微分器通过程序实现,代码如下:
%接受信号通过微分器处理
for i=1:length(t)-1
diff_nsfm(i)=(nsfm(i+1)-nsfm(i))./dt;
end
diff_nsfmn = abs(hilbert(diff_nsfm)); %hilbert 变换,求绝对值得到瞬时
幅%度(包络检波)
高斯白噪声信道特性
设正弦波通过加性高斯白噪声信道后的信号为
r
(
t
)
A
c o s (
c
t
)
(
n t
)
其中,白噪声 ( )n t 的取值的概率分布服从高斯分布。
MATLAB 本身自带了标准高斯分布的内部函数 randn 。 randn 函数产生的
随机序列服从均值为 0m ,方差 2
1 的高斯分布。
正弦波通过加性高斯白噪声信道后的信号为
)
( )
r t A
cos(
t
c
( )
n t
故其有用信号功率为
噪声功率为
信噪比 S
N 满足公式
则可得到公式
S
2
A
2
2
N
B
10log (
10
S
)
N
2
2
A
B
1 0
2
1 0
我们可以通过这个公式方便的设置高斯白噪声的方差。
添加高斯白噪声
sn1=10;
sn2=30;
sn=0;
%设定信躁比(小信噪比)
%设定信躁比(大信噪比)
%设定信躁比(无信噪比)
db=am^2/(2*(10^(sn/10)));
%计算对应的高斯白躁声的方差
n=sqrt(db)*randn(size(t));
%生成高斯白躁声
nsfm=n+sfm;
%生成含高斯白躁声的已调信号(信号通
%过信道传输)
调频系统的抗噪声性能分析
从前面的分析可知,调频信号的解调有相干解调和非相干解调两种。相干
解调仅适用于窄带调频信号,且需同步信号;而非相干解调适用于窄带和宽带
调频信号,而且不需同步信号,因而是 FM 系统的主要解调方式,所以这里仅
仅讨论非相干解调系统的抗噪声性能,其分析模型如图 9 所示。
图 9 调频系统抗噪声性能分析模型
图中带通滤波器的作用是抑制信号带宽以外的噪声。 ( )n t 是均值为零,单
边功率谱密度为 0n 的高斯白噪声,经过带通滤波器后变为窄带高斯噪声 ( )
in t 。
限幅器是为了消除接收信号在幅度上可能出现的畸变。
设调频信号为
S
F M
( )
t
A
co s(
c
t
K
f
t
m
(
)
d
)
故其输入功率为
输入噪声功率为
因此输入信噪比为
S
i
2
A
2
N n B
o FM
i
S
i
N
i
A
B
2
2
FM
在大信噪比条件下,信号和噪声的相互作用可以忽略,这时可以把信号和
噪声分开来算,这里,我们可以得到解调器的输出信噪比
3
S
N
o
o
2
2
2
f
( )
A K m t
8
3
n f
o m
2
上式中, A 为载波的振幅, fK 为调频器灵敏度, mf 为调制信号 ( )m t 的最
高频率, on 为噪声单边功率谱密度。
我们如若考虑 ( )m t 为单一频率余弦波时的情况,可得到解调器的制度增益
为
S
S
o
i
F M
G
N
N
考虑在宽带调频时,信号带宽为
FM
2(
m
f
B
3
2
m
2
f
o
i
2
A
2
n f
o m
1)
f
m
2(
f
f
m
)
则可以得到
G
FM
23
m m
f
(
f
1)
可以看出,大信噪比时宽带调频系统的信噪比增益是很高的,它与调频指
数的立方成正比。可见,加大调频指数 fm ,可使调频系统的抗噪声性能迅速改
善。
源代码
%初始化
echo
off
close
all
clear
all
clc
%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
%调制,产生已调信号