yanchi2.m

延迟码的波形图。通过matlab工具仿真出了延迟码的时域和频域波形

%延迟调制码产生和功率谱密度
%时间单位 us 频率单位 MHz 码元速率 Mb/s
close all
clear all
k=input('取样点数=2^k,k=[10]');
if isempty(k),
	k=10;
end
N=2^k;	%采样点数
L=32;	%每码元采样点数
M=N/(2*L);	%码元数
Rb=2;	%码元速率
Ts=1/Rb;	%码元间隔
dt=Ts/L;	%时域取样间隔
df=1/(N*dt);	%频域取样间隔
T=N*dt;	%时域截短时间
Bs=N*df/2;	%频域系统带宽
t=[-T/2+dt/2:dt:T/2];	%时域横坐标
f=[-Bs+df/2:df:Bs];		%频域横坐标

Ep=zeros(size(f));	%定义全零距阵，将来作为求功率的累加器	
for ii=1:50
    a=zeros(1,M);   %先产生全零的数组，用于保存码元
	while a(1)==0,  %为了便于后面的差分计算，是产生的序列第一个数值为1
        a=round(rand(1,M));
    end
    %根据码元产生密勒码
    b=zeros(1,2*M); %产生全零序列，用于保存密勒码
    %产生密勒码
    for index=1:M 
        if a(index)==1,
            b(2*index-1)=0;
            b(2*index)=1;
        else
            if a(index-1)==1,
                b(2*index-1)=0;
                b(2*index)=0;
            else
                b(2*index-1)=0;
                b(2*index)=0;
            end
        end
    end
    %对密勒码进行差分编码
    b(1)=mod((b(1)+1),2);
    for tt=2:length(b)
        b(tt)=mod((b(tt)+b(tt-1)),2);
    end
    b=sign(b-0.5);%再把差分编码后的序列对应相应的波形
    %把差分编码后的序列重复L次，即相当于对每个码元采样L次
    bc=ones(1,L);
    yc=b(bc,:);
    yc=reshape(yc,1,N);    
    YC=fftshift(fft(yc)).*dt;    %进行傅立叶变换
    P=YC.*conj(YC)/T;  %求功率谱密度
    Ep=(Ep*(ii-1)+P)/ii;%累加求均值
end
%画延迟码波形
subplot(2,1,1)
plot(t,yc)
grid on
axis([-5,5,-1.2,1.2]);
title('延迟码波形');
xlabel('t(us)')
ylabel('s(t)(v)');
%画延迟码功率谱密度
aa=30+10*log10(Ep+eps);
subplot(2,1,2)
plot(f*Ts,aa)
grid on
axis([0,2,0,50])
xlabel('f*Tb')
ylabel('Ps(f)(W/Hz)')
title('延迟码功率谱密度')