rands1.m

来自「对数字信号进行分析」· M 代码 · 共 20 行

function [rand1,l2]=rands1(l1,mpeg,tb)
N2=188*8;
rand2=zeros(1,l1*N2);
redata=[1,0,0,1,0,1,0,1,0,0,0,0,0,0,0];
    for j=1:1:l1*N2
        p1=xor(redata(1,14),redata(1,15));
        p2=xor(p1,mpeg(1,j));
        rand2(1,j)=p2;
        for k=15:-1:2
            redata(1,k)=redata(1,k-1);
        end
        redata(1,1)=p1;
    end
    rand1=rand2;
runum4=length(rand1);
l2=runum4;
t2=tb:tb:tb*runum4;
% figure
% stairs(t2,rand1)
% axis([0 tb*runum4 -0.5 1.5])