pres.m

北邮《通信原理基于Matlab的计算机仿真》书的源码 m文件 好不容易弄来的

%部分响应信号眼图示意
clear all; close all;
Ts=1;
N_sample=16;
eye_num = 11;

N_data=1000;

dt = Ts/N_sample;
t = -5*Ts:dt:5*Ts;

%产生双极性数字信号
d = sign(randn(1,N_data));
dd= sigexpand(d,N_sample);
%部分响应系统冲击响应
ht = sinc((t+eps)/Ts)./(1- (t+eps)./Ts);
ht( 6*N_sample+1 ) = 1;
st = conv(dd,ht);
tt = -5*Ts:dt:(N_data+5)*N_sample*dt-dt;

figure(1)
subplot(211);
plot(tt,st);
axis([0 20 -3 3]);xlabel('t/Ts');ylabel('部分响应基带信号');
subplot(212)
%画眼图
ss=zeros(1,eye_num*N_sample);
ttt = 0:dt:eye_num*N_sample*dt-dt;
for k=5:50
    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('部分响应信号眼图');