clear
clc
x=[1?0?1?0];
n=402;
for?i=1:n
???temp=xor（x（3）x（4））;
???M（i）=x（4）;
????for?j=4:-1:2
????????x（j）=x（j-1）;
????end
????x（1）=temp;
end
v=randn（1400）;%創建服從正態隨機分布矩陣
z=zeros（4021）;
z（1）=0;
z（2）=0;
for?i=3:402