m序列特性

1.移位相加后

function y=yiwei(m)

m=[1 0 0 0 1 1 1 1 0 1 0 1 1 0 0];

N=length(m);

m0=m;

x=input('请输入移位长度:');

m=[m(x+1:N) m(1:x)];

m=mod(m+m0,2);

for n=1:N

if m0==m

disp('满足移位相加后是原序列，移位：');n

break;

else

m=[m(2:N) m(1)];

end

if n==N

disp('不满足移位后是原序列');

end

end

2.游程特性

function y=youcheng(m)

m=[1 0 0 0 1 1 1 1 0 1 0 1 1 0 0];

N=length(m);

s=0;s0=1;t=0;t0=1;

m=[m m(1)];

for i=1:N

if m(i)==0&&m(i+1)==1

s=s+1;s1(s)=s0;s0=1;

else if m(i)==0&&m(i+1)==0

s0=s0+1;

end

end

if m(i)==1&&m(i+1)==0

t=t+1;t1(t)=t0;t0=1;

else if m(i)==1&&m(i+1)==1

t0=t0+1;

end

end

end

disp('零游程的个数为：');s

disp('他们的长度分别为：');s1

disp('一游程的个数为：');t

disp('他们的长度分别为：');t1

3.平衡特性

function y=pingheng(m)

m=[1 0 0 0 1 1 1 1 0 1 0 1 1 0 0];

N=length(m);

s=0;

for i=1:N

if m(i)==1

s=s+1;

end

end

s

if s==(N+1)/2

disp('满足平衡特性'); else

disp('不满足平衡特性'); end