%第三章
clear?all
close?all
clc

%周期為15的M序列
X1=1;X2=1;X3=1;X4=1;??%移位寄存器輸入Xi初T態（1111），?Yi為移位寄存器各級輸出
n=1500000;?????????????????%置M序列總長度
for?i=1:n??
???Y4=X4;??Y3=X3;??Y2=X2;???Y1=X1;
???X4=Y3;??X3=Y2;??X2=Y1;??
???X1=xor（Y3Y4）;??%異或運算
???U（i）=Y4;
???if?Y4==0
??????U（i）=1;
???else
??????U（i）=Y4*（-1）;
???end
end?
M=U‘;

%?UU（1）=1;
%?for?e=2:45
%?????UU（e）=UU（e-1）+（-1）^（e+1）;???
%?end
%?M=xor（OUU‘）;
%?輸出y
num=[1?5?7];
den=[1?8?19?12];
sys=tf（numden）;
F=c2d（sys1‘zoh‘）;
t=1:1:n;
u=M（t）;
y=