x  = 10*randn(1,500);                     % 作為系統的輸入
    b  = fir1(63,0.5)                          %確定未知系統
    d  = filter(b,1,x)                          %期望信號
    mu =0.015;                               %固定步長值
    h  = adaptfilt.lms(64,mu);                   %模擬建模
 [y,e]  = filter(h,x,d);                           % 模擬系統輸出
     subplot(3,2,1); plot(1:500,d);axis([450,480,-25,25]);
     title(‘未知系統輸出');xlabel('采樣點'); ylabel('輸出'); 
     subplot(3,2,3); plot(1:500,y);axis([450,480,-25,25]);
     title('模擬系統輸出');xlabel('采樣點'); ylabel('輸出'); 
     subplot(3,2,5); stem(1:500,e);axis([450,480,-1,1]);
     title(‘誤差值’);xlabel(采樣點'); ylabel('誤差');