,请自觉下载,转载请注明出处一共四个程序图片被压缩了,下载后放大会很清晰
1文件名为digit_baseband.m
%数字基带信号的功率谱密度digit_baseband.m
clear all;
close all;
Ts=1;
N_sample=8; %每个码元的抽样点数
dt=Ts/N_sample; %抽样时间间隔
N=1000; %码元数
%T=N*dt;
t=0:dt:(N*N_sample-1)*dt;
gt1=ones(1,N_sample); %NRZ非归零波形
gt2=ones(1,N_sample/2); %RZ归零波形
gt2=[gt2 zeros(1,N_sample/2)];
d=(sign(randn(1,N))+1)/2;
data=sigexpand(d,N_sample); %对序列间隔插入N_sample-1个0
st1=conv(data,gt1); %调用MATLAB的卷积函数
st2=conv(data,gt2);
d=2*d-1; %变成双极性序列
data=sigexpand(d,N_sample);
st3=conv(data,gt1);
st4=conv(data,gt2);
[f,st1f]=T2F(t,[st1(1:length(t))]);
[f,st2f]=T2F(t,[st2(1:length(t))]);
[f,st3f]=T2F(t,[st3(1:length(t))]);
[f,st4f]=T2F(t,[st4(1:length(t))]);
figure(1)
subplot(421)
plot(t,[st1(1:length(t))]);grid
axis([0 20 -1.5 1.5]);
ylabel('单极性NRZ波形');
subplot(422);
plot(f,10*log10(abs(st1f).^2/Ts));grid
axis([-5 5 -40 40]);
ylabel('单极性NRZ功率谱密度(dB/Hz)');
subplot(423)
plot(t,[st2(1:length(t))]);