实验一-数字基带传输实验-实验总结报告

数字基带传输实验2014、3、28数字基带传输实验总结报告小组成员：所在班级：通信一班指导老师：马丕明数字基带传输实验2014、3、28目录一、实验目的................................................................................................................3二、实验原理..............................................................................................................3三、实验内容................................................................................................................4（一）因果数字升余弦滚降滤波器设计...........................................................41.窗函数法设计非匹配形式的发送滤波器..............................................42.频率抽样法设计匹配形式的发送滤波器.................................................6（二）设计无码间干扰的二进制数字基带传输系统.........................................81、子函数模块................................................................................................82、无码间干扰的数字二进制基带传输系统的模拟..................................11四、实验总结：........................................................................................................145数字基带传输实验2014、3、28一、实验目的1、提高独立学习的能力；2、培养发现问题、解决问题和分析问题的能力；3、学习Matlab的使用；4、掌握基带数字传输系统的仿真方法；5、熟悉基带传输系统的基本结构；6、掌握带限信道的仿真以及性能分析；7、通过观测眼图和星座图判断信号的传输质量。二、实验原理图1基带系统传输模型1、信源信源就是消息的源，本实验中指数字基带信号，信源序列al采用一个0、1等概率分布的二进制伪随机序列。信源序列al经在一比特周期中抽样A点，即是序列al每两点之前插A-1个零点，进行抽样，形成发送信号SigWave，即是发送滤波器模块的输入信号。2、发送滤波器匹配形式下的发送滤波器SF，通过窗函数法对模拟升余弦滚降滤波器的时域单位冲激响应hd进行时间抽样、截断、加窗、向右移位而得；非匹配形式下的发生滤波器SF，通过频率抽样法对模拟升余弦滚降滤波器的频率响应Hd进行频率抽样、离散时间傅里叶反变换、向右移位而得。发送滤波器输出SFO是由发送滤波器SF和发送信号SigWave卷积而得。3、传输信道本实验中传输信道采用理想信道，即传输信道频率响应函数为1；传输信道输出信号Co是由发送滤波器输出信号SFO和加性高斯白噪声GN叠加而成:Co=SFO+GN。4、噪声信道噪声当做加性高斯白噪声，给定标准差调用函数randn生成高斯分布随机数GN。5、接收滤波器匹配形式下，接收滤波器与发送滤波器单位冲激响应幅度相同，角度相反，均为平方根信源发送滤波器信道噪声接收滤波器抽样判决位定时提取输出数字基带传输实验2014、3、28升余弦滚降滤波器。非匹配形式下，接收滤波器为直通滤波器。6、位定时提取完成位定时信息即同步判决时间点的提取，每隔A-1个点提取一次信息，一共提取L次。注意：非匹配模式下每提取一个同步判决时间点需经过一次位延时，匹配模式下每提取一个同步判决时间点需经过两次位延时。7、抽样判决利用同步位定时信息判决接收滤波器的输出信号Rec_Sig，得到输出序列。三、实验内容（一）因果数字升余弦滚降滤波器设计1.窗函数法设计非匹配形式的发送滤波器n=-15:15;a=input('alpha=');Tc=4;hn=(sin(pi*n/Tc)./(pi*n/Tc)).*(cos(a*pi*n/Tc)./(1-4*a*a*n.*n/Tc/Tc));%升余弦滚降滤波器时域单位冲击响应hn(16)=1;%升余弦滚降滤波器中间点的校正n=0:30;wn=0.42-0.5.*cos(2.*pi.*n./(max(n)-1))+0.08.*cos(4.*pi.*n./(max(n）-1））；%Blackman窗函数表达式figure;subplot(2,2,1);stem(n,hn);xlabel('n');ylabel('hn');title('升余弦滤波器单位冲击响应时域特性');subplot(2,2,2);stem(n,wn);xlabel('n');ylabel('wn');title('布拉克曼窗单位冲击响应时域特性');h=hn.*wn;%加窗后的升余弦滚降滤波器subplot(2,2,3);stem(n,h);xlabel('n');ylabel('h');title('发送滤波器的冲击响应');nh=0:30;w=linspace(0,2*pi,512);hw=h*exp(-j*nh'*w);%频率特性subplot(2,2,4);plot(w,abs(hw));title('发送滤波器的幅频特性');axis([min(w),max(w),min(abs(hw))-0.2,max(abs(hw))+0.2])xlabel('w');ylabel('abs(hw)')g=max(hw);Figure；plot(w,10*log10(abs(hw)));xlabel('w');ylabel('Grf(dB)')title('发送滤波器的增益图形');数字基带传输实验2014、3、28Alpha=0.20102030-0.500.51nhn升余弦滤波器单位冲击响应时域特性010203000.51nwn布拉克曼窗单位冲击响应时域特性0102030-0.500.51nh发送滤波器的冲击响应024601234X:1.23Y:0.02952发送滤波器的幅频特性wabs(hw)图2窗函数法设计非匹配形式滤波器冲击响应(alpha=0.2)01234567-25-20-15-10-50510X:1.34Y:-18.96wGrf(dB)发送滤波器的增益图形图3窗函数法设计非匹配形式滤波器增益图形(alpha=0.2)Alpha=1数字基带传输实验2014、3、28-15-10-5051015-101nSF(n)非匹配滤波器时域单位冲击响应波形(1)00.10.20.30.40.50.60.70.80.9100.51/非匹配滤波器归一化幅频特性(1)00.10.20.30.40.50.60.70.80.91-200-1000/20log︳H(ej)︳(dB)非匹配滤波器增益图形（dB）(1)图3窗函数法设计非匹配滤波器冲击响应波形与增益(alpha=1)2.频率抽样法设计匹配形式的发送滤波器实验思路：匹配形式的发送滤波器的频率响应为平方根余弦滚降函数，子函数在主函数开方即可。实验程序:子函数定义的是余弦滚降函数，主函数是一个频率抽样：T=1,fs=1,Tc=4,N=31%子函数function[hn,Hf,f]=f_sampling(N,Tc,fs)alpha=0.2;K=[-(N-1)/2:(N-1)/2];n=[-(N-1)/2:(N-1)/2];f=K*fs/N;f1=(1-alpha)/(2*Tc);f2=(1+alpha)/(2*Tc);Hf=zeros(1,N);%升余弦滚降滤波器频率响应抽样函数fori=1:Nif(abs(f(i))=f1)Hf(i)=Tc;elseif(abs(f(i))=f2)Hf(i)=Tc/2*(1+cos(pi*Tc/alpha*(abs(f(i))-(1-alpha)/(2*Tc))));elseHf(i)=0;end;end;hn=1/N*Hf*exp(j*2*pi/N*K'*n);%数字升余弦滚降滤波器单位冲击响应%主函数Function[hm]=f_samp_m(N,Tc,fs)N=31;K=[-(N-1)/2:(N-1)/2];m=[-(N-1)/2:(N-1)/2];[hn,Hf,f]=f_sampling(31,4,1);HF=squrt(Hf);hm=1/N*HF*exp(j*2*pi*/N*K’*m);[Hw,w]=freqz(hm);%幅频特性；H=max(abs(Hw));plot(w,abs(Hw)/H);titel(‘归一化的幅频特性’);xlabel(‘w’);ylabei(‘Hw’);Figure；Ab=20*log(abs(Hw));plot(w,Ab);title(‘增益图形’)；xlabel(‘w’);ylabel(‘20*log(abs(Hw))(dB)’);数字基带传输实验2014、3、28Figure；subplot(2,1,1);stem(HF,’.’);title(‘匹配滤波器频域单位冲击响应波形’)；xlabel(‘f(Hz)’);ylabel(‘HF’);Figure；subplot(2,1,2);stem(real(hm),’.’);xlabel(‘m’);ylabel(‘nm’);axis([0,35,-0.15,0.7]);title(‘匹配滤波器时域单位冲击响应波形’)；Alpha=0.20102030-0.500.51nhn升余弦滤波器单位冲击响应时域特性010203000.51nwn布拉克曼窗单位冲击响应时域特性0102030-0.500.51nh发送滤波器的冲击响应0246800.51X:1.266Y:0.00168发送滤波器的归一化幅频特性wabs(hw)/max(abs(hw))图4频率抽样法设计匹配滤波器冲击响应与幅频特性(alpha=0.2)01234567-25-20-15-10-50510X:1.328Y:-18.94wGrf(dB)发送滤波器的增益图形图5频率抽样法设计匹配滤波器增益图形(alpha=0.2)Alpha=1数字基带传输实验2014、3、28-15-10-5051015-101nSF(n)非匹配滤波器时域单位冲击响应波形(1)00.10.20.30.40.50.60.70.80.9100.51/非匹配滤波器归一化幅频特性(1)00.10.20.30.40.50.60.70.80.91-100-500/20log︳H(ej)︳(dB)非匹配滤波器增益图形（dB）(1)图6、7频率抽样法设计匹配滤波器幅频特性与增益图形(alpha=1)表1窗函数法设计的非匹配滤波器表二频率抽样法设计的匹配滤波器滚降系数0.210.21第一零点带宽（Hz）0.40.450.2660.2第一旁瓣衰减（dB）18.962318.9426（二）设计无码间干扰的二进制数字基带传输系统要求传输的二进制比特个数、比特速率错误!未找到引用源。（可用与错误!未找到引用源。的关系表示）、信噪比SNR、滚降系数α是可变的．1、子函数模块1）二进制信源子函数function[al]=suiji(m)%a_n是产生的二进制随机信源序列%m是产生的序列长度%本函数功能是随机序列产生函数al=rand(1,m);fori=1:mifa(i)0.5a(i)=-1;else数字基带传输实验2014、3、2