下面是matlab程序:
clear all;
close all;
Ts = 1;%每个码元的时间长度
N_sample = 10;%每个码元内抽样点数
N = 100;%1码元数
dt = Ts/N_sample;
t = 0:dt:N_sample*N*dt-dt;
gt = ones(1,N_sample);%码型,此为不归零矩形
da = sign(randn(1,N));%产生随机序列
n = 0:N-1;
st = sigexpand(da,N_sample);%从随机序列到基带信号
st = st(1:length(st)-N_sample+1);
st = conv(st,gt);%二进制数字基带信号
h1t = gt;%匹配滤波器冲激响应的抽样
h2t = 5*sinc(5*(t-5)/Ts);%理想低通滤波器冲激响应
r1t = conv(st,h1t);%匹配滤波器输出
r1t = r1t/10;%从数字序列的卷积到对应模拟信号卷积的抽样
r2t = conv(st,h2t);%理想低通滤波器输出
r2t = r2t/10;%从数字序列的卷积到对应模拟信号卷积的抽样
figure(1)
subplot(321)
<