clear?all;?close?all;
a=[1?0.5?0];%分子系數向量
b=[1?-1.5?0.7];%分母系數向量
d=zeros（170）;%輸出序列50個
?for?n=1:1:70
?????c=a-a（:1）./b（:1）.*b;
?????d（:n+1）=a（:1）./b（:1）;
?????a（:1）=c（:2）;
?????a（:2）=c（:3）;??
?end
plot（d‘g‘）;
grid?on;?????????%xianshi?gezi（shange）
title（‘標準單位脈沖響應與辨識結果‘）;
hold?on
?
?a=0.9;?%可變系數
?Np=63;??%M序列周期
?r=4;????%4個周期
?h=1;???????%?幅值
?L=r*Np;?????%仿真長度
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%產生e（k）
?d=[1?-1.5?0.7];?c=1?;??%?分子分母多項式系數
?nd=length（d）-1?;???%階次
?v=rand（1L+1）;
?xi=a.*v;
?ek=zeros（nd1）;
?e=zeros（1L）;%行向量
?for?k=1:L
????e（k）=-d（2:（nd+1））*ek+c*xi（k）;?
????%數據更新
????for?i=nd:-1:2
????????ek（i）=ek（i-1）;
????end
????ek（1）=e（k）;
?end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%產生Y（k）序列
M_XuLie=zeros（L1）;
%初始化M序列
M（1）=1;
M（2）=1;
M（3）=1;
M（4）=0;
M（5）=1;
M（6）=0;
for?n?=?1?:1:?L
????if?M（6）==0
??????M_XuLie（n）?=?h;
????else
??????M_XuLie（n）?=-h;
????end???????????????%輸出序列
????temp?=?xor（M（6）?M（5））;
????M（6）?=?M（5）;
????M（5）?=?M（4）;
????M（4）?=?M（3）;
????M（3）?=?M（2）;
????M（2）?=?M（1）;
????M（1）?=?temp;
end
ui=M_XuLie;???%輸入信號，M序列，列向量L行
d1=[1?-1.5?0.7];?c1=[0?1?0.5];??%分母?分子多項式系數
nd1=length（d1）-1?;nc1=length（c1）-1;???%階次
uik=zeros（nc11）;??%輸入初值
yk=zeros（nd11）;??%輸出初值
y=zeros（1L）;%行向量
for?k=1: