L=3000;



z（2）=0;z（1）=0;%取z的前二個初始值為零
for?k=3:L+4;?
????u（k）=randn（1）;
????z（k）=1.5*z（k-1）-0.7*z（k-2）+1.0*u（k-1）+1.5*u（k-2）+?normrnd（0111）;%給出理想的辨識輸出采樣信號
????zn（k1）=z（k）;
end

c0=[0?0?0?0?]‘;?%a1?a2?b1?b2給出被辨識參數的初始值

p0=10^6*eye（44）;

c=[c0zeros（42999）];%被辨識參數矩陣的初始值及大小


for?k=3:3000;?%開始求K?
????h1=[-z（k-1）-z（k-2）u（k-1）u（k-2）]‘;?x=h1‘*p0*h1+1;?x1=inv（x）;?%開始求K（k）
????k1=p0*h1*x1;%求出K的值
????d1=z（k）-h1‘