eyespic3.m

单极性码间干扰与信道均衡

%升余弦滚降系统(2径)的眼图
clear all
close all
N=2^14;  %采样点数
L=32;  %每码元采样点数
M=N/L;  %码元数
Rb=2;   %码速率
Ts=1/Rb;  %码元间隔
dt=Ts/L;  %时域采样间隔
df=1/(N*dt);  %频域采样间隔
d=16;  %半个码元采样点数
T=N*dt;  %截短时宽
Bs=N*df/2;  %频带宽度
Na=4;   %画眼图一次四个码元
t=[-T/2+dt/2:dt:T/2];  %时域坐标
f=[-Bs+df/2:df:Bs];  %频域坐标
alpha=0.01;	 %滚降系数
g1=sin(pi*t/Ts)./(pi*t/Ts);
g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
g=g1.*g2;
Cf=1-0.5*exp(-2*i*pi*f*d); %信道传输函数
Ge=1./Cf;%信道均衡
G=t2f(g,dt).*Cf.*Ge; %系统总的频谱
%预设眼图的位置
figure(1)
set(1,'position',[150,100,350,200]);
hold on
grid
xlabel('t in us');
ylabel('s(t) in v');
for ii=1:10
    %求系统后的码元样本序列
    a=sign(randn(1,M));
    imp=zeros(1,N);
    imp(L/2:L:N)=a/dt;
    S=t2f(imp,dt).*G; %序列在频域通过系统
    s=f2t(S,dt); %转换到时域
    s=real(s)+0*randn(1,N); %可加入随机噪声
    %画眼图
    figure(1)
    tt=[0:dt:Na*L*dt];
	for jj=1:Na*L:N-Na*L
		plot(tt,s(jj:jj+Na*L))
	end
end