matlab costas环,MATLAB写COSTAS环的问题。

X1=0;X2=0;X3=1;

m=350;               %重复50遍的7位单极性m序列

for i=1:m

Y3=X3;  Y2=X2;   Y1=X1;

X3=Y2;   X2=Y1;

X1=xor(Y3,Y1);

L(i)=Y1;

end

for i=1:m

M(i)=1-2*L(i);        %将单极性m序列变为双极性m序列

end

k=1:1:m;

% figure(1)

% subplot(3,1,1)                        %做m序列图

% stem(k-1,M);

% axis([0,7,-1,1]);

% xlabel('k');

% ylabel('M序列');

% title('移位寄存器产生的双极性7位M序列') ;

% subplot(3,1,2)

ym=fft(M,4096);

magm=abs(ym);                %求双极性m序列频谱

fm=(1:2048)*200/2048;

% plot(fm,magm(1:2048)*2/4096);

% title('双极性7位M序列的频谱')

% axis([90,140,0,0.1]);

% [a,b]=xcorr(M,'unbiased');

% subplot(3,1,3)                        %求双极性m序列自相关函数<