数字通信课程设计
匹配滤波器
摘要
在通信系统中，滤波器是重要的部件之一，滤波器特征的选择直接影响数字信号的恢复。
在数字信号接收中，滤波器的作用有两个方面，使滤波器输出有用信号成分尽可能强；抑制信号带外噪声，使滤波器输出噪声成分尽可能小，减少噪声对信号判决的影响。
对最佳线性滤波器的设计有一种准则是使滤波器输出信噪比在特定时刻到达最大，由此导出的最佳线性滤波器称为匹配滤波器。
在数字通信中，匹配滤波器具有广泛的应用。
因此匹配滤波器是指滤波器的性能与信号的特征取得某种一致，使滤波器输出端的信号瞬时功率与噪声平均功率的比值最大。
本文设计并仿真了一种数字基带通信系统接收端的匹配滤波器。
一、课程设计的目的
通过本次对匹配滤波器的设计，让我们对匹配滤波器的原理有更深一步的理解，掌握具体的匹配滤波器的设计方法与算法。
二、课程设计的原理
设接收滤波器的传输函数为H (f),冲击响应为h(t),滤波器输入码元s(t)的
持续时间为T s,信号和噪声之和r(t)为
r(t) s(t) n(t) 0 t T s
式中，s(t)为信号码元，n(t)为白噪声。
并设信号码元s(t)的频谱密度函数为S(f),噪声n(t)的双边功率谱密度为
P n n o/2 , n o为噪声单边功率谱密度。
假定滤波器是线性的，根据叠加定理，当滤波器输入信号和噪声两部分时，滤波器的输出也包含相应的输出信号和输出噪声两部分，即
y(t) S o(t) n o(t)
由于：&(f) H*(f)H(f )P R(f) |H(f)2P R(f)
PJf)为输出功率谱密度，P R(f)为输入功率谱密度，P R(f) n。
/2
这时的输出噪声功率N o等于
N o H⑴气f
于H
⑴2df
在抽样时刻t o上，输出信号瞬时功率与噪声平均功率之比为
为了求出「0的最大值，利用施瓦兹不等式求「0的最大值
在白噪声干扰的背景下，按上式的设计的线性滤波器，将能在给定时刻
t 0上
获得最大输出信噪比2E n 0。
是输出信噪比最大时刻。
这种滤波器就是最大信 噪比意义下的最佳线性滤波器，由于它的传输特性与信号频谱的复共轭一致， 称为匹配滤波器。
匹配滤波器的特性还可以用冲激响应函数 h(t):
h (t ) H (f )e j2 ft df kS *(f )e j 2 ft 0
e j 2 ft df
* *
k s
( ) ( t 0 t )d
ks (t
t
)
可见，匹配滤波器的冲激响应h(t)就是信号s(t)的镜像s(t 0 t),但在时 间轴上平移了 t 。
t 。
是输出信噪比最大时刻。
三、课程设计的内容
设计仿真了数字基带通信系统在接收端接收到模拟信号后通过匹配滤波 器的效果，其中匹配滤波器的输入信号我们假定为正弦信号，经过匹配滤波 器的输出信号作为抽样判决器的输入信号。
我们通过比较匹配滤波器的输入 输出信号，验证得出其性能特性。
如下图所示的仿真结果，我们假定匹配滤波器的输入信号为正弦信号与白 噪声的叠加，其频率，采样频率,延时时间，。
「0
S o (t o ) N o
H(f)S(f)e j2fto df
n o ~2
2
H(f) df 「0
2
H(f)S(f)e
j2 ft0
df
扌
H(f)2df
2 2
H(f) df S(f) df
号 |H(f)2
df
2E n 。
且当时H (f)
kS *(f)e j2ft0
，等式成立，即得到知道的信噪比为
2E n °
由实验的结果，我们可以看出匹配滤波器的单位冲激响应幅度谱与输入信号 的幅度谱一致，符合上面我们推导的公式；由滤波器的输出信号我们可以得 出，在抽样时刻处，信号有最大值，在输出信号的幅度谱中，干扰噪声白噪 声明显减小，信号的有用成分被保留下来，体现出了匹配滤波器抑制噪声， 保证在抽样时刻有最大值的特性，从而进一步保证了输出信噪比在抽样时刻 最大，为下一步抽样判决提供了有利的条件。
由此可以看出匹配滤波器对数 字基带通信系统的通信可靠性的提高有着很重要的作用。
四、附录 实验程序：
clear all ; close all ;
fc=20; %余弦信号的频率 fs=5*fc; %采样频率
ts=1/fs; k=1;
0.5
1 1.5
t/s 冲激响应
0 2
-50
0 f(Hz)
单位冲激响应幅度谱
1
|
W
0.5
50 -10
0.5
1000
f -C -
节 ---------------------- = ------------------------ ■■
O
Ik
l a
卅
I 1
-1000
0 -50
0 f(Hz)
输出信号幅度谱
100 50
-50
0 f(Hz)
50
输入信号+噪声
10 0 -10 输入信号幅度谱
1 1.5 2
10
0 t/s
匹配滤波器输出信号
12 3 4
t/s
50
to=2; %采样时刻t=[0:ts:2-ts];
t1=[0:ts:4-2*ts];
freq1=li nspace(-fs/2,fs/2,le ngth(t)); freq2=li nspace(-fs/2,fs/2,le ngth(t1));
Noise=ra ndn (1,le ngth(t)); %白噪声
sig nal=3*cos(2*pi*fc*t+pi*k*t.A2)+Noise;
h1=3*cos(2*pi*fc*(to-t)+pi*k*(to-t).A2)+Noise; So1=co nv2(sig nal,h1); %(言号卷积%信号
% 匹配滤波器
subplot(321)
plot(t,sig nal);
axis([0,2,-10,10]);
xlabel( 't/s' );
ylabel( 'S(t)');
title( '输入信号+噪声');
subplot(322) plot(freq1,fftshift(abs(fft(sig nal))/le
ngth(t))),xlabel( el( '|S(W)|' ),title( '输入信号幅度谱’);
'f(Hz)' ),ylab
subplot(323)
plot(t,h1);
axis([0,2,-10,10]);
xlabel( 't/s' );
ylabel( 'h(t)' );
title( '冲激响应');
subplot(324)
plot(freq1,fftshift(abs(fft(h1)/le ngth(t)))),xlabel(
|H1(W)|' ),title( '单位冲激响应幅度谱’);
'f(Hz)' ),ylabel(
subplot(325)
plot(t1,So1);
xlabel( 't/s' );
ylabel( 'So(t)' );
title('匹配滤波器输岀信号');
subplot(326) plot(freq2,fftshift(abs(fft(So1)/le
ngth(t1)))),xlabel( ('|So1(W)|' ),title( '输岀信号幅度
谱’);
'f(Hz)' ),ylabel。