u = 0.00005;    %迭代步長
n = 64;  %WHT-LMS的抽頭數必須是2的冪
h = zeros(n ,1);    %抽頭系數，n為抽頭個數
x = zeros(1 ,n);    %濾波器抽頭輸入
%% 注意：讀者需要給出輸入信號xd, 參考信號d，運行本程序就能得到相應仿真圖
e = zeros(n ,1);       %估計誤差
y = zeros(len , 1);  %估計信號，len為數據長度
beita = 0.8;    %
p = zeros(n, 1);   %信號功率估計
for i = 1 : len
    x(1,2:end) = x(1,1:end-1);
    x(1,1) = xd(i,1);   %xd為輸入信號
    xh = WHT(x, n); %完成WHT變換
    p = p*beita + (1-beita)*abs(xh).^2;   %功率估計
    y(i) = xh.'*conj(h);       %估計信號
    e(i) = d(i) - y(i);       %d為參考信號
    h = h + 2*u.*conj(e(i))*xh./(p+0.02);   %抽頭系數更新
end
for i=1:len;
    err(i1)=abs(e(i).^2);    %統計均方誤差
end