通信系统原理研究项目 项目2.docx

通信系统原理研究项目码间串扰与信道均衡一、 研究题目设失真信道的脉冲（冲击）响应为：其中为符号速率，对其进行速率为的采样，输入给迫零均衡器。迫零均衡器的输出为：其中2K+1为均衡器的阶数，cn为均衡器的系数。（1） 分别对于K=2,4,6，画出均衡器输入的脉冲波形；（2） 分别对于K=2,4,6，画出均衡器输出的脉冲波形，并定量比较均衡前后的效果二、 仿真结果2.1、均衡器输入的脉冲波形对于不同的K=2,4,6，输入迫零均衡器的脉冲波形是一样的，如图2.1所示，其中图中曲线各点的数据如下表2.1MATLAB中的函数代码如下：function xt=func_xt(t,T)%脉冲冲击响应函数xt=1./(1+(2*t/T).2);end%T=0.001;%周期自定义t=-6:0.1:6;k=-6:6;ys_in=func_xt(t*T,T);%连续函数波形yd_in=func_xt(k*T,T);%离散函数波形plot(t,ys_in,k,yd_in,o);图2.1迫零均衡器的输入脉冲波形k-6-5-4-3-2-10123456xkT0.00690.00990.01540.02700.05880.200010.20000.05880.02700.01540.00990.0069表2.1迫零均衡器输入波形数据2.2、均衡器参数与输出脉冲波形理论分析：根据题意：均衡器参数根据以下矩阵求出，在这里，以N=2为例其中：在MATLAB中编写程序如下：%function C=func_C(N,T) A=zeros(2*N+1,2*N+1); b=zeros(N,1);1;zeros(N,1);for i=0:2*N A(i+1,:)=func_xt(T/2)*(2*(i-N)-(-N:N),T); end C=Ab;End%function yt=func_yt(m,N,C,T)yt=0;for i=-N:N yt=yt+C(i+N+1,1)*func_xt(T/2)*(2*m-i),T);endend%function y=func_y(t,N,C,T)y = zeros(1,length(t);for i=1:length(t) y(i)=func_yt(t(i),N,C,T);endend%主程序部分%N_2=2;C_2=func_C(N_2,T);yd_2=func_y(k,N_2,C_2,T);ys_2=func_y(t,N_2,C_2,T);figure,plot(t,ys_2,k,yd_2,o);xlabel(k),ylabel(y(kT);title(N=2时迫零均衡器输出波形); N_4=4;C_4=func_C(N_4,T);yd_4=func_y(k,N_4,C_4,T);ys_4=func_y(t,N_4,C_4,T);figure,plot(t,ys_4,k,yd_4,o);xlabel(k),ylabel(y(kT);title(N=4时迫零均衡器输出波形); N_6=6;C_6=func_C(N_6,T);yd_6=func_y(k,N_6,C_6,T);ys_6=func_y(t,N_6,C_6,T);figure,plot(t,ys_6,k,yd_6,o);xlabel(k),ylabel(y(kT);title(N=6时迫零均衡器输出波形);整理仿真结果如下：N=2时迫零均衡器参数如表2.2所示C-2-1012-2.20504.8958-3.01384.8958-2.2050表2.2N=2时迫零均衡器参数N=2时迫零均衡器输出波形如图2.2图2.2N=2时迫零均衡器输出波形N=2时迫零均衡器输出的数值如表2.3k-6-5-4-3-2-10123456x(kT)0.010.010.020.030.060.210.20.060.030.020.010.01y(kT)0.020.020.030.040010000.040.030.02表2.3N=2时迫零均衡器输出N=4时迫零均衡器参数如表2.4所示C-4-3-2-1012340.0003-0.05870.5469-1.76542.5583-1.76540.5469-0.05870.0003表2.4N=4时迫零均衡器参数N=4时迫零均衡器输出波形如图2.3图2.3N=4时迫零均衡器输出波形N=4时迫零均衡器输出的数值如表2.5k-6-5-4-3-2-10123456x(kT)0.010.010.020.030.060.210.20.060.030.020.010.01y(kT)0000001000000表2.5N=4时迫零均衡器输出N=6时迫零均衡器参数如表2.6所示C-6-5-4-3-2-10-0.00020.0015-0.0121-0.03120.5519-1.80582.5926C123456-1.80580.5519-0.0312-0.01210.0015-0.0002表2.6N=6时迫零均衡器参数N=6时迫零均衡器输出波形如图2.4图2.4N=6时迫零均衡器输出波形N=6时迫零均衡器输出的数值如表2.7k-6-5-4-3-2-10123456x(kT)0.010.010.020.030.060.210.20.060.030.020.010.01y(kT)0000001000000表2.7N=6时迫零均衡器输出三、 仿真结果分析从上面的仿真结果可以看出，随着N的增大，即均衡器阶数（2N+1）的增大，迫零均衡器的调整效果越好，比如在N=2时，在k=-6处，经迫零均衡器输出的结果y=0.02，但是当