clear all;
close all
fs=107;%采样频率
Ts=1;%基带信号码元周期
N_sample=fs;%码元采样点数
N=1000;%码元数目
% ------ 时域参数 ------ %
dt=1/N_sample;
t=0:dt:(N*N_sample-1)*dt;
% ------ 绘制升余弦滚降滤波器的冲激响应函数 ------ %
a=1;
t0=-3:dt:3;
ht=sin(pi*t0/Ts)./(pi*t0./Ts).*cos(a*pi*t0/Ts)./(1-4*a.^2*t0.^2/Ts.^2);
ht((length(t0)+1)/2)=1;
figure(1);
plot(t0,ht,'LineWidth',1.5);
axis([-3 3 -0.2 1.2]);
title('升余弦滤波器的冲击响应函数h(t)');
xlabel('时间(s)');
ylabel('电压值(V)');
grid on;
% ------ 初始化二进制数字信号码元 ------ %
bc0=[-1,zeros(1,106)];
bc1=[1,zeros(1,106)];
% ------ 生成二进制随机基带信号 ------ %
rand_code=randi([0,1],1,N);
% ------ 初始化二进制基带脉冲信号 ------ %
mt=zeros(1,N*N_sample-1);
% ------ 计算二进制基带脉冲信号 ------ %
for i=1:N
if rand_code(i)==1
mt((i-1)*N_sample+1:i*N_sample)=bc1;
else
mt((i-1)*N_sample+1:i*N_sample)=bc0;
end
end
% ------ 绘制二进制基带脉冲信号 ------ %
figure(2)
plot(t,mt,'LineWidth',1.5);
axis([0 20 -1.2 1.2]);
title('基带脉冲信号时域波形');
xlabel('时间(s)');
ylabel('电压值(V)');
grid on;
% ------ 计算并绘制二进制基带成形信号 ------ %
st=conv(ht,mt);
t1=-3:dt:(N*N_sample-1)*dt+3;
figure(3)
plot(t1,st,'LineWidth',1.5);
axis([0 20 -1.2 1.2]);
title('基带成形信号时域波形');
xlabel('时间(s)');
ylabel('电压值(V)');
grid on;
% ------ 频域参数 ------ %
df=fs/length(t1);
f=-fs/2:df:fs/2-df;
fst=fft(st);
fst=fftshift(fst);
fst=abs(fst);
fst_dB=fst.^2/Ts;
max_F=max(fst_dB);
fst_dB=fst_dB/max_F;
fst_dB=10*log10(fst_dB+eps);
figure(4)
plot(f,fst_dB);
axis([-50 50 -160 0]);
title('基带成形脉冲信号频域功率谱');
xlabel('频率(Hz)');
ylabel('功率谱幅度值(dB)');
grid on;
% ------ 绘制基带成形脉冲信号发送端眼图 ------ %
n=3;
t2=0:dt:n-dt;
figure(5);
for i=10:30
bst=st((n*(i-1)+0.5)*N_sample+1:(n*i+0.5)*N_sample);
plot(t2,bst);
title('发送端的基带信号眼图');
xlabel('时间(s)');
ylabel('电压值(V)');
grid on;
hold on;
end
% ------ 绘制基带成形脉冲信号叠加高斯白噪声后的信号 ------ %
ast=awgn(st,25);
figure(6)
plot(t1,ast,'LineWidth',1.5);
axis([0 20 -1.2 1.2]);
title('接收到的基带成形信号时域波形');
xlabel('时间(s)');
ylabel('电压值(V)');
grid on;
% ------ 绘制接收端的基带成形脉冲信号发送端眼图 ------ %
n=3;
t2=0:dt:n-dt;
figure(7);
for i=10:30
bst=ast((n*(i-1)+0.5)*N_sample+1:(n*i+0.5)*N_sample);
plot(t2,bst);
title('接受端的基带信号眼图');
xlabel('时间(s)');
ylabel('电压值(V)');
grid on;
hold on;
end
- 1
- 2
- 3
前往页