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

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/（10*N*（A^2））?%計算能夠使LMS算法收斂的delta
%delta=0.01;
u=zeros（1M）;
y=zeros（1M）;
h=zeros（1N）;?%LMS濾波器系數
y1=zeros（1N）;
a=zeros（1M）;
b=zeros（1M）;
b（N-1）=0.1;
c=zeros（1M）;
b0=0.98;??%?b=0.08??a=0.8最優
e（N-1）=d（N-1）;
for?i=N:M?%系數調整LMS算法
x1=x（i:-1:i-N+1）;
%LMS算法
y（i）=h*x1‘;
e（i）=d（i）-y（i）;
a（i）=abs（e（i）/e（i-1））;
b（i）=b0*b（i-1）+（1-b0）*abs（e（i）*e（i-1））;
c（i）=1-exp（-a（i）*abs（e（i）*e（i-1）））;
u（i）=b（i）*（1-exp（-a（i）*abs（e（i）*e（i-1））））;
h=h+u（i）*e（i）*x1;
end
error=e.^2;?%LMS算法每一步迭代的均方誤差
for?n=N:M?%利用求解得到的h，與輸入信號x做卷積，得到濾波后結果
x2=x（n:-1:n-N+1）;
y1（n）=h*x2‘;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1689??2020-11-03?05:45??╗?╙┌S║?╩2╡─╦π╖?\bijiao.asv
?????文件????????1797??2020-11-03?05:45??╗?╙┌S║?╩2╡─╦π╖?\bijiao.m
?????文件????????1125??2020-11-03?05:45??╗?╙┌S║?╩2╡─╦π╖?\EFLMSfangzhen.m
?????文件?????????922??2020-11-03?05:45??╗?╙┌S║?╩2╡─╦π╖?\GSfangzhen.m
?????文件????????1010??2020-11-03?05:45??╗?╙┌S║?╩2╡─╦π╖?\yuan.m