miller.m

来自「miller码的matlab」· M 代码 · 共 50 行

clear all;
close all;
TS=1;
x=round(rand(1,1000));
L=100;
t=0:1/L:length(x);
i=1;
if(x(i)==1)
    for j=1:L/2;
        y(L/2*(2*i-2)+j)=1;
        y(L/2*(2*i-1)+j)=-1;
    end
else
    for j=1:L;
        y(L*(i-1)+j)=1;
    end
end
for i=2:length(x)
    if(x(i)==1)
        for j=1:L/2
            y(L/2*(2*i-2)+j)=y(L/2*(2*i-3)+L/4);
            y(L/2*(2*i-1)+j)=-y(L/2*(2*i-2)+j);
        end
    else
        if(x(i-1)==1)
            for j=1:L
                y(L*(i-1)+j)=y(L/2*(2*i-3)+L/4);
            end
        else
            for j=1:L
                y(L*(i-1)+j)=-y(L/2*(2*i-3)+L/4);
            end
        end
    end
end
figure(1);
y=[y x(i)];
subplot(211)
plot(t,y);
axis([0 20 -1.5 1.5]);
xlabel('Miller码');


%[f stf]=T2F(t,y);
%subplot(212)
%plot(f,10*log10(abs(stf+eps).^2/TS));grid
%axis([-5 5 -40 20]);ylabel('Miller码功率谱密度(dB/Hz)');