实验四-码间串扰仿真实验

实验四码间串扰仿真实验一、实验目的：1、掌握无码间串扰的传输条件；2、掌握升余弦滚降系统的MATLAB仿真；3、能够用MATLAB画出接收信号的眼图。二、实验内容：1、画出a=0.5时，升余弦滚降系统的频谱，并画出对应的时域波形。Ts=1;N_sample=17;dt=Ts/N_sample;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))/TsXf(n,k)=0;elseifabs(f(k))0.5*(1-a(n))/TsXf(n,k)=Ts;elseXf(n,k)=0.5*Ts*(1+cos(pi*Ts/(a(n)+eps)*(abs(f(k))-0.5*(1-a(n))/Ts)));end;end;xt(n,:)=sinc(t/Ts).*(cos(a(n)*pi*t/Ts))./(1-4*a(n)^2*t.^2/Ts^2+eps);endfigure(1)plot(f,Xf);axis([-1101.2]);xlabel('f/Ts');ylabel('升余弦滚降频谱');figure(2);plot(t,xt);axis([-1010-0.51.1]);xlabel('t');ylabel('升余弦滚降波形');2、例5-9，分别画出通过低通滤波器后的眼图。N=1000;N_sample=8;Ts=1;dt=Ts/N_sample;t=0:dt:(N*N_sample-1)*dt;gt=ones(1,N_sample);d=sign(randn(1,N));a=sigexpand(d,N_sample);st=conv(a,gt);ht1=5*sinc(5*t/Ts);rt1=conv(st,ht1);ht2=2*sinc(2*t/Ts);rt2=conv(st,ht2);eyediagram(rt1+j*rt2,40,5);function[out]=sigexpand(d,M)%N=length(d);out=zeros(M,N);out(1,:)=d;out=reshape(out,1,M*N);