clear?all
close?all
N=12;?????????????????????%濾波器階數
sample_N=1000;?????????????%采樣點數

A=1;?????????????????????%信號幅度
snr=20;????????????????????%信噪比
t=1:sample_N;?
length_t=100;???????%期望信號序列長度
d=A*sin（2*pi*t/length_t）;?%期望信號
M=length（d）??????????????%M為接收數據長度
x=awgn（dsnr）????????????%經過信道（加噪聲）

delta=1/（20*N*（A^2））????%計算能夠使LMS算法收斂的delta
???????????
y=zeros（1M）;
h=zeros（1N）;?????????????%LMS濾波器系數
h_normalized=zeros（1N）;??%歸一化LMS濾波器系數
y1=zeros（1N）;
u0=0.2;
r=0.8
for?n=N:M?????????????????%系數調整LMS算法
????x1=x（n:-1:n-N+1）;
????
????%LMS算法
????y（n）=h*x1‘;??????????????%濾波
????e（n）=d（n）-y（n）;??????????%誤差估計
????h=h+delta*e（n）*x1;???????%權向量更新
????
????%NLMS算法
????y_normalized（n）=h_normalized*x1‘;????????????????