关于数字基带传输系统的滤波器研究论文模板Word格式文档下载.docx
- 文档编号:20584833
- 上传时间:2023-01-24
- 格式:DOCX
- 页数:20
- 大小:375.36KB
关于数字基带传输系统的滤波器研究论文模板Word格式文档下载.docx
《关于数字基带传输系统的滤波器研究论文模板Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《关于数字基带传输系统的滤波器研究论文模板Word格式文档下载.docx(20页珍藏版)》请在冰豆网上搜索。
码间串扰的定义:
由于系统传输特性不良或加性噪声的影响,使信号波形发生畸变,造成收端判决
上的困难,因而造成误码,这种现象称为码间串扰。
发生码间串扰时,脉冲会被展宽,甚至重迭(串扰)到邻近时隙中去成为干扰。
假如传输的一组码元是1110、采用双极性码、经发送滤波器后变为升余弦波形所示。
经过信道后产生码间串扰,前3个“1”码的拖尾相继侵入到第4个“0”码的时隙中,如图2所示。
图2码间串扰示意
升余弦升余弦滚降特性(数字基带信号
码间串扰)
标签:
数字基带信号
码间串扰升余弦滚降系统频谱时域波形α
杂谈
%数字基带信号传输
码间串扰升余弦滚降系统的频谱及其时域波形
%文件名syx_gunjiang.m
Ts=1;
N=17;
dt=Ts/N;
df=1.0/(20.0*Ts);
t=-10*Ts:
dt:
10*Ts;
f=-2/Ts:
df:
2/Ts;
a=[0,0.5,1];
forn=1:
length(a)
fork=1:
length(f)
ifabs(f(k))>
0.5*(1+a(n))/Ts
Xf(n,k)=0;
elseifabs(f(k))<
0.5*(1-a(n))/Ts
Xf(n,k)=Ts;
else
Xf(n,k)=0.5*Ts*(1+cos(pi*Ts/(a(n)+eps)*(abs(f(k))-0.5*(1-a(n))/Ts)));
end;
xt(n,:
)=sinc(t/Ts).*(cos(a(n)*pi*t/Ts))./(1-4*a(n)^2*t.^2/Ts^2+eps);
end
subplot(211);
plot(f,Xf);
axis([-1101.2]);
xlabel('
f/Ts'
);
ylabel('
升余弦滚降系统的频谱'
legend('
α=0'
'
α=0.5'
α=1'
subplot(212);
plot(t,xt);
axis([-1010-0.51.1]);
t'
升余弦滚降系统的时域波形'
2.3无码间串扰的基带传输特性
一个好的基带传输系统,应该在传输有用信号的同时能尽量抑制码间串扰和噪声。
为便于讨论,先忽略信道噪声,同时把基带传输系统模型作一简化,如图3所示。
图3基带传输系统的简化图
图中
,为发送滤波器、信道、接收滤波器之总和,是整个系统的基带传输特性。
如果无码间串扰,系统的冲激响应满足:
即抽样时刻(
点)除当前码元有抽样值之外,其它各抽样点上的取值均应为0。
根据频谱分析,可以写出
满足此式的
)就是能实现无码间串扰的基带传输函数。
最简单的无码间串扰的基带传输函数是理想低通滤波器的传输特性
式中
为常数代表带内衰减。
波形如图4所示。
图4理想低通滤波器的传输特性
1、奈奎斯特(Nyquist)定理(奈奎斯特第一准则)
当基带传输系统具有理想低通滤波器特性时,以截至频率两倍的速率传输数字信号,使其能消除码间串扰。
2、理想低通滤波器基带传输的特征参量
(1)奈奎斯特带宽
(2)奈奎斯特速率
(3)奈奎斯特间隔
(4)无码间串扰的理想低通系统的频带利用率
3.理想低通滤波器的缺点
(1)工程不易实现,滤波器截至特性不会做得狠陡。
(2)接收时对判断要求很严。
(3)冲激响应衰减慢,拖尾长。
4.
clearall;
closeall;
aa=[1011101000011];
%选取了13个码元
ts=1;
%一个码元持续时间
nsample=17;
%一个码元由17个脉冲组成
aaa=zeros(17,13);
fori=1:
17
aaa(i,:
)=aa;
end
a=reshape(aaa,1,17*13);
dt=ts/nsample;
df=1.0/(10.0*ts);
t=-2*ts:
2*ts;
f=-2/ts:
2/ts;
alpha=1;
length(f)
0.5*(1+alpha)/ts;
f>
(1+1)/ts
xf(k)=0;
0.5*(1-alpha)/ts;
f<
xf(k)=ts;
else
xf(k)=0.5*ts*(1+cos(pi*ts/(alpha+eps)*(abs(f(k))-0.5*(1-alpha)/ts)));
xt=sinc(t/ts).*(cos(alpha*pi*t/ts))./(1-4*alpha^2*t.^2/ts^2+eps);
figure
(1)
plot(f,xf);
f/ts'
升余弦滚降系统'
grid;
figure
(2)
axis([-22-0.11.1]);
升余弦滚降波形'
figure(3)
na=length(a);
tmp=length(xt);
t1=na+tmp-1;
tt=-(na+tmp-1)/2*dt:
(na+tmp-1)/2*dt-dt;
xt3=conv(xt,a);
plot(tt,xt3)
基带波形'
figure(4),plot(a)
图形一:
设定一个随机序列的基带信号,给出的是已知基带信号序列aa=[1011101000011];
选取了13个码元,用MATLAB做出图如上。
图形二:
是升余弦滚将系统α=1的传输特性曲线图。
随着α的增大其h(t)的拖尾衰减越快,对位定时精度要求越低,但是也会使带宽增大。
图形三:
是a=1的升余弦滚将特性h(t)的波形,该h(t)满足抽样值上无串扰的传输条件,而且尾部衰减较快,这有利于减小码间串扰和位定时误差的影响。
抽样频率在1/ts的时候,满足无码间串扰。
图形四:
基带信号经过升余弦滚将系统后的输出波形,是序列和余弦滚降系统的卷积结果。
2.4滚降低通滤波器
为了解决理想低通特性存在的问题,可以使理想低通滤波器特性的边沿缓慢下降,这称为“滚降”。
因此滚降低通滤波器的设计就显得十分的有必要,这也是为此设计的原因。
1.滚降低通滤波器特性
例如图5(a),这是一个具有升余弦滚降特性传递函数的低通滤波器,
从图形上看是3个相邻段H(ω-2π/Ts)、H(ω)、H(ω+2π/Ts)分别被移到(-π/Ts,π/Ts)区间(即把图c、
d移至图b中)迭加,得到的Heq(ω)为一矩形如图e所示,此低通滤波器的带宽
当传输速率
(
)时,此基带传输系统可以实现无码间串扰。
图5升余弦滚降系统的时域响应图
在实际中传输网络不可能是理想低通,通常采用满足几对称条件的滚降低通滤波器来等效理想低通。
相应的波形如图所示。
图6升余弦滚降滤波器
2.滚降因子α
α为带宽的扩展量与奈奎斯特带宽Wc之比。
α的特性:
越大抽样函数的拖尾振荡起伏越小、衰减越快。
与理想低通相比,它付出的代价是带宽增加了一倍。
此时系统的最高传码率虽然没变,但频带宽度已被扩展
,α在0—1之间变化。
可见,图示具有升余弦滚降传输特性的滤波器满足奈氏第一准则,其带宽
传输速率
(
)
频带利用率
比理想低通滤波器的频带利用率低了一倍。
3.通过程序设计滤波器
例:
用凯塞窗设计低通滤波器,通带边界频率Ωp=0.3pi,阻带边界频率
Ωs=0.5pi,阻带衰减δs不小于50dB。
程序如下:
b=fir1(29,0.4,kaiser(30,4.55));
[h1,w1]=freqz(b,1);
plot(w1/pi,20*log10(abs(h1)));
归一化频率/p'
)
;
幅度/dB’);
波形如下:
2.5MATLAB眼图设计观察
在实际的基带传输系统中,由于难免存在滤波器的设计误差和信道特性的变化,所以无法实现理想的传输特性和滤波效果,存在码间串扰,从而导致系统性能的下降。
而这些问题的存在又是客观存在的,不可能因为试验次数的多次进行而减小,因此引入眼图的概念加以解决。
所谓眼图,是指利用实验的方法估计和改善(通过调整)传输系统性能时在示波器上观察到的一种图形。
观察眼图的方法是:
用一个示波器跨接在接收滤波器的输出端,然后调整示波器扫描周期,使示波器水平扫描周期与接收码元的周期同步,这时示波器屏幕上看到的图形像人的眼睛,故称为“眼图”。
从“眼图”上可以观察出码间串扰和噪声的影响,从而估计系统优劣程度。
另外也可以用此图形对接收滤波器的特性加以调整,以减小码间串扰和改善系统的传输性能。
下图为眼图的模型:
(1)最佳抽样时刻应在“眼睛”张开最大的时刻。
(2)对定时误差的灵敏度可由眼图斜边的斜率决定。
斜率越大,对定时误差就越灵敏。
(3)在抽样时刻上,眼图上下两分支阴影区的垂直高度,表示最大信号畸变。
(4)眼图中央的横轴位置应对应判决门限电平。
(5)在抽样时刻上,上下两分支离门限最近的一根线迹至门限的距离表示各相应电平的噪声容限,噪声瞬时值超过它就可能发生错误判决。
(6)对于利用信号过零点取平均来得到定时信息的接收系统,眼图倾斜分支与横轴相交的区域的大小,表示零点位置的变动范围,这个变动范围的大小对提取定时信息有重要的影响。
实现眼图的MATLAB编程要用到的函数有下面几个
subplotplottitlexlabelylabelaxis
下面结合具体的实例进行分析:
数字基带信号波形及其眼图:
%升余弦滚降系统眼图
%数字基带信号波形及其眼图
%文件名:
eye.m
N=15;
eye_num=6;
a=1;
N_data=1000;
t=-3*Ts:
3*Ts;
%产生双极性数字信号
d=sign(randn(1,N_data));
dd=sigexpand(d,N);
%
%基带系统冲击响应(升余弦)
ht=sinc(t/Ts).*(cos(a*pi*t/Ts))./(1-4*a^2*t.^2/Ts^2+eps);
st=conv(dd,ht);
tt=-3*Ts:
(N_data+3)*N*dt-dt;
subplot(211)
axis([020-1.21.2]);
t/Ts'
基带信号'
subplot(212)
%画眼图
ss=zeros(1,eye_num*N);
ttt=0:
eye_num*N*dt-dt;
fork=3:
50
ss=st(k*N+1:
(k+eye_num)*N);
drawnow;
plot(ttt,ss);
holdon;
基带信号眼图'
------------------------------------------------------------------------
%将输入的序列扩成间隔为N-1个0的序列
function[out]=sigexpand(d,M)
N=length(d);
out=zeros(M,N);
out(1,:
)=d;
out=reshape(out,1,M*N);
3程序总结
本程序可扩展性好。
可以多次改变参数数值,便于观察图形的变化,易于对实验原理的理解以及对实验分析的进行,本次程序编写中运用的函数较为普遍,适用性较好,且方便功能扩展;
可是也存在着一些缺点,比如仿真的效果方面,由于仿真工具低下,使一些关键部分不能很好地展现出来。
对此,是有要改进的。
4课程总结
本课程设计的主要目的是熟悉无码间串扰的时域和频域条件,掌握典型的无码间干扰基带传输系统——升余弦滚降系统。
以及用matlab显示升余弦滚降波形及频谱以及接收端的基带信号波形。
如:
设基带传输系统响应是1的升余弦滚降系统,显示升余弦滚降波形及频谱,以及接收端的数字基带信号波形。
通过实验,我们能加深对课程的理解,提高综合应用所学的知识解决实际问题的能力,得出科学研究的前期训练。
通过这次课程设计使我懂得了理论与实际相结合是很重要的,只有理论知识是远远不够的,只有把所学的理论知识与实践相结合起来,从理论中得出结论,才能真正为社会服务,从而提高自己的实际动手能力和独立思考的能力。
附录(部分程序)
clearall
globaldttfdfNT%全局变量
closeall
N=2^13;
%采样点数
L=64;
%每码元的采样点数
M=N/L;
%码元数
Rb=2;
%码速率是2Mb/s
Ts=1/Rb;
%码元间隔
dt=Ts/L;
%时域采样间隔
df=1/(N*dt);
%频域采样间隔
T=N*dt;
%截短时间
Bs=N*df/2;
%系统带宽
%alpha=0.5;
%滚降系数=0.5
Na=4;
%示波器扫描宽度为4个码元
t=linspace(-T/2,T/2,N);
%时域横坐标
f=linspace(-Bs,Bs,N);
%频域横坐标
db=input('
请选择信噪比[0-15]:
'
b=input('
采样偏差b*Ts,b=[-0.5,+0.5]'
alpha=input('
滚降系数'
Rt=input('
占空比'
ifRt==[],Rt=0.5;
hr1=sin(pi*t/Ts)./(pi*t/Ts);
hr2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
hr=hr1.*hr2;
HR=abs(t2f(hr));
GT=sqrt(HR);
GR=GT;
EP=zeros(size(f))+eps;
EPr=zeros(size(f))+eps;
['
信噪比为'
num2str(db),'
dB,采样偏差为'
num2str(b),'
*Ts'
]
forloop1=1:
16
Eb_N0(loop1)=(loop1-1);
%分贝值变为真值
eb_n0(loop1)=10^(Eb_N0(loop1)/10);
Eb=1;
n0=Eb/eb_n0(loop1);
%信道噪声谱密度
sita=n0*Bs;
%噪声功率
n_err=0;
%误码计数
forii=1:
20
code=sign(randn(1,M));
imp=zeros(1,N);
%产生冲激序列
imp(L/2:
L:
N)=code/dt;
IMP=t2f(imp);
Sa=IMP.*GT;
%升余弦信号的傅氏变换
sa=f2t(Sa);
%升余弦信号的时域波形
sa=real(sa);
P=Sa.*conj(Sa)/T;
%升余弦信号的功率谱
EP=(EP*(ii-1)+P)/ii;
n_ch=sqrt(sita)*randn(size(t));
%信道噪声
nr=real(f2t(t2f(n_ch).*GR));
%输出噪声
Sr=Sa.*GR;
%接收信号频谱
sr=real(f2t(Sr))+nr;
%接收信号
y=sign(sr(L*(.5):
N));
%抽样判别
n_err=n_err+length(find(y~=code));
%误码数
Pr=Sr.*conj(Sr)/T;
%平均功率
EPr=(EPr*(ii-1)+Pr)/ii;
tt=[0:
Na*L*dt];
ifloop1==db+1
['
画眼图'
];
forjj=1:
Na*L:
N-Na*L
figure
(2);
holdon;
subplot(2,1,1);
gridon;
plot(tt,sa(jj:
jj+Na*L));
%画发送眼图
title('
发送眼图'
xlabel('
t(us)'
)
ylabel('
s(t)(V)'
axis([0,2,-2.3,2.3]);
subplot(2,1,2);
plot(tt,sr(jj:
%画接收眼图
接收眼图'
end
画波形图'
s=zeros(1,N);
s=reshape(code(ones(1,L),:
),1,L*M);
yo=zeros(1,N);
yo=reshape(y(ones(1,L),:
figure
(1)
subplot(2,3,1);
plot(t,s,'
LineWidth'
2);
%画发送码型
axis([-10,+10,1.5*min(s),1.5*max(s)])
发送码形'
subplot(2,3,2);
plot(t,sa);
%画生成波形
axis([-10,10,1.5*min(sa),1.5*max(sa)]);
发送生成波形'
subplot(2,3,3)
plot(f,30+10*log10(EP));
%画功率谱
axis([-.05*Bs,.05*Bs,min(30+10*log10(EP))/3,1.5*max(30+10*log10(EP))]);
f(MHz)'
Ps(f)(dBm/MHz)'
发送功率谱'
subplot(2,3,4);
plot(t,yo,'
%画接收码
axis([-10,+10,1.5*min(yo),1.5*max(yo)])
接收码形'
subplot(2,3,5);
plot(t,sr);
%画接受波形
axis([-10,10,1.5*min(sr),1.5*max(sr)]);
sr(t)(V)'
接收波形'
subplot(2,3,6)
plot(f,30+10*log10(EPr+eps));
%画接收功率谱
axis([-.05*Bs,.05*Bs,min(30+10*log10(EPr+eps))/2,1.5*max(30+10*log10(EPr+eps))]);
ylabel
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 关于 数字 基带 传输 系统 滤波器 研究 论文 模板