demo_danjinrz.m

系带传输

%单极性不归零码的时域图及功率谱密度，“0”码“1”码等概出现，只有在零频率处存在离散谱。
function demo_danjinrz
clear
close all
k=16; %取样点数＝2^k
z=8;  %每个信号取样点数＝2^z,z<k
f0=1; %f0=取1(kz)左右
N=2^k;
L=2^z;
M=N/L;
Ts=1/f0; %码元宽度是1us
dt=Ts/L; %时域采样间隔
df=1/(N*dt); %MHz
T=N*dt; %截短时间
Bs=N*df/2; %系统带宽
f=[-Bs+df/2:df:Bs]; %频域横坐标
t=[-T/2+dt/2:dt:T/2]; %时域横坐标
EP=zeros(size(f))+eps;
for ii1=1:30
    b=round(rand(1,M)); %产生随机码
    s=zeros(1,N);
    for ii=1:L;%for ii=1:L/2;码型变为占空比1/2的NRZ.
        s(ii+[0:M-1]*L)=b;
    end
    S=fftshift(fft(s))*dt; %S是s的傅氏变换
    P=S.*conj(S)/T; %功率谱
    EP=(EP*(ii1-1)+P+eps)/ii1;
end
subplot(2,1,2); %输出的功率画图
plot(f,10*log10(EP),'m');
grid
axis([-15,15,-100,20]);
xlabel('f (KHz)');ylabel('P (V/KHz)') ;title('功率谱图')
subplot(2,1,1); %输出信号波形画图
plot(t,s,'b')
grid
axis([-10,10,-0.2,1.2]);
xlabel('t (ms)');ylabel('s(t) (V)');title('单极性NRZ的时域图')