dmachester.m

来自「多路复用技术」· M 代码 · 共 37 行

function y=dmachester(x,grid);
t=0:1/grid:length(x);
i=1;
if(x(i)==1)
    for j=1:grid/2
        y(grid/2*(2*i-2)+j)=0;
        y(grid/2*(2*i-1)+j)=1;
    end
else
    for j=1:grid/2
        y(grid/2*(2*i-2)+j)=1;
        y(grid/2*(2*i-1)+j)=0;
    end
end
for i=2:length(x)
    if(x(i)==1)
        for j=1:grid/2
        y(grid/2*(2*i-2)+j)=1-y(grid/2*(2*i-3)+grid/4);
        y(grid/2*(2*i-1)+j)=1-y(grid/2*(2*i-2)+j);
        end
    else
        for j=1:grid/2
        y(grid/2*(2*i-2)+j)=y(grid/2*(2*i-3)+grid/4);
        y(grid/2*(2*i-1)+j)=1-y(grid/2*(2*i-2)+j);
        end
    end
end
y=[y,y(i*grid)];
M=max(y);  m=min(y);
plot(t,y); axis([0,i,m-0.1,M+0.1]);
title('1 0 0 1 1 0 0 0 0 1 0 1');