通信系统原理研究项目项目.docxVIP

通信系统原理研究项目——码间串扰与信道均衡研究题目设失真信道的脉冲（冲击）响应为：其中为符号速率，对其进行速率为的采样，输入给迫零均衡器。迫零均衡器的输出为：其中2K+1为均衡器的阶数，为均衡器的系数。分别对于，画出均衡器输入的脉冲波形；分别对于，画出均衡器输出的脉冲波形，并定量比较均衡前后的效果仿真结果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迫零均衡器的输入脉冲波形-6-5-4-3-2-101234560.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=A\b;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(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