基带信号的眼图实验

基带信号的眼图实验

 

1. 设基带信号传输系统的响应是的升余弦滚降系统，画出在接收端的基带数字信号波形及眼图。

提示：

%需要用到sigexpand这个函数程序

% sigexpand.m

%

% Simulation program to expand the symbol sequence

%

% Programmed by Y.S. Cui and H.Harada

%

function [out]=sigexpand(d,M)  

N=length(d);             %基带信号码元长度

out=zeros(M,N);          %矩阵M为采样点  N为基带信号码元数量

out(1,:)=d;             %将零矩阵第一行换成基带信号中的8个码元

    out=reshape(out,1,M*N);  % 1行 m*n 列

 

1. 设二进制是数字基带信号，，，设加性高斯白噪声的双边功率谱密度为，画出眼图。

1. 经过理想低通滤波器后的眼图。
2. 经过理想低通后的眼图。

%基带信号眼图，

clear all;

close all;

Ts=1;

N_sample=17;

eye_num=7;

alpha=1;

N_data=1000;

dt=Ts/N_sample;

t=-3*Ts:dt:3*Ts;

%产生双极性数字信号

d=sign(randn(1,N_data));

dd=sigexpand(d,N_sample);

%基带系统冲击响应（升余弦）

ht=sinc(t/Ts).*(cos(alpha*pi*t/Ts))./(1-4*alpha^2*t.^2/Ts^2+eps);

st=conv(dd,ht);

st=conv(dd,ht);

tt=-3*Ts:dt:(N_data+3)*N_sample*dt-dt;

figure(1)

subplot(211);

plot(tt,st);

axis([0 20 -1.2 1.2]);

xlabel(t/Ts);ylabel('部分响应基带信号');

subplot(212)

%画眼图

ss=zeros(1,eye_num*N_sample);

ttt=0:dt:eye_num*N_sample*dt-dt;

for k=3:52

    ss=st(k*N_sample+1:(k+eye_num)*N_sample);

    drawnow;

    plot(ttt,ss);

    hold on;

end 

%plot(ttt,ss);

xlabel('t/Ts');ylabel('部分响应基带信号');

%需要用到sigexpand这个函数程序

% sigexpand.m

%

% Simulation program to expand the symbol sequence

%

% Programmed by Y.S. Cui and H.Harada

%

function [out]=sigexpand(d,M)  

N=length(d);             %基带信号码元长度

out=zeros(M,N);          %矩阵M为采样点  N为基带信号码元数量

out(1,:)=d;             %将零矩阵第一行换成基带信号中的8个码元

    out=reshape(out,1,M*N);  % 1行 m*n 列

 

 

 

%示意双极性NRZ基带信号经过带宽受限信号造成的码间干扰影响及其眼图

clear all;

close all;

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-5)/Ts);

rt1=conv(st,ht1);

ht2=sinc((t-5)/Ts);

rt2=conv(st,ht2);

eyediagram(rt1+j*rt2,40,5);

%调用Matlab画眼图的函数，行40点，表示5只眼

基带信号的眼图实验

 

1. 设基带信号传输系统的响应是的升余弦滚降系统，画出在接收端的基带数字信号波形及眼图。

提示：

%需要用到sigexpand这个函数程序

% sigexpand.m

%

% Simulation program to expand the symbol sequence

%

% Programmed by Y.S. Cui and H.Harada

%

function [out]=sigexpand(d,M)  

N=length(d);             %基带信号码元长度

out=zeros(M,N);          %矩阵M为采样点  N为基带信号码元数量

out(1,:)=d;             %将零矩阵第一行换成基带信号中的8个码元

    out=reshape(out,1,M*N);  % 1行 m*n 列

 

1. 设二进制是数字基带信号，，，设加性高斯白噪声的双边功率谱密度为，画出眼图。

1. 经过理想低通滤波器后的眼图。
2. 经过理想低通后的眼图。

%基带信号眼图，

clear all;

close all;

Ts=1;

N_sample=17;

eye_num=7;

alpha=1;

N_data=1000;

dt=Ts/N_sample;

t=-3*Ts:dt:3*Ts;

%产生双极性数字信号

d=sign(randn(1,N_data));

dd=sigexpand(d,N_sample);

%基带系统冲击响应（升余弦）

ht=sinc(t/Ts).*(cos(alpha*pi*t/Ts))./(1-4*alpha^2*t.^2/Ts^2+eps);

st=conv(dd,ht);

st=conv(dd,ht);

tt=-3*Ts:dt:(N_data+3)*N_sample*dt-dt;

figure(1)

subplot(211);

plot(tt,st);

axis([0 20 -1.2 1.2]);

xlabel(t/Ts);ylabel('部分响应基带信号');

subplot(212)

%画眼图

ss=zeros(1,eye_num*N_sample);

ttt=0:dt:eye_num*N_sample*dt-dt;

for k=3:52

    ss=st(k*N_sample+1:(k+eye_num)*N_sample);

    drawnow;

    plot(ttt,ss);

    hold on;

end 

%plot(ttt,ss);

xlabel('t/Ts');ylabel('部分响应基带信号');

%需要用到sigexpand这个函数程序

% sigexpand.m

%

% Simulation program to expand the symbol sequence

%

% Programmed by Y.S. Cui and H.Harada

%

function [out]=sigexpand(d,M)  

N=length(d);             %基带信号码元长度

out=zeros(M,N);          %矩阵M为采样点  N为基带信号码元数量

out(1,:)=d;             %将零矩阵第一行换成基带信号中的8个码元

    out=reshape(out,1,M*N);  % 1行 m*n 列

 

 

 

%示意双极性NRZ基带信号经过带宽受限信号造成的码间干扰影响及其眼图

clear all;

close all;

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-5)/Ts);

rt1=conv(st,ht1);

ht2=sinc((t-5)/Ts);

rt2=conv(st,ht2);

eyediagram(rt1+j*rt2,40,5);

%调用Matlab画眼图的函数，行40点，表示5只眼