用脉冲函数求传递函数

%第三章
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=