《自适应滤波器原理及Matlab仿真应用》配套matlab代码

《自适应滤波器原理及Matlab仿真应用》配套matlab代码

%clear all; clc; close all;

load 01_memory
load wo
y=filter（wo1x）;
eo=0.001*randn（size（y））;
yeo=y+eo;
r=y;
%sound（yFs）
w=zeros（Lw1）;
e=y;
zeta=zeros（size（y））;
mu=0.5;
M=5;
X=zeros（LwM）;
for n=Lw:Lw+M-2
    xtdl=x（n:-1:n-Lw+1）;
    X=[xtdl X（:1:M-1）];
end
for n=Lw+M-1:length（x）-M
    xtdl=x（n:-1:n-Lw+1）;
    X=[xtdl X（:1:M-1）];
    Yhat=X‘*w;
    yhat=Yhat（1）;
    E=yeo（n:-1:n-M+1）-Yhat;
    e（n）=E（1）;
    r（n）=y（n）-yhat;
    w=w+mu*X*inv（X‘*X+0.1*eye（M））*E;
    zeta（n）=（w-wo）‘*（w-wo）;
end
E=audioplayer（eFs）;
play（E）
ERLE=10*log10（（filter（.001[1 -.999]y.^2）+1e-12）./（filter（.001[1 -.999]r.^2）+1e-12））;
figure（2）plot（ERLE）
ERLEAPLMS=ERLE;
eAPLMS=e;
save ERLEAPLMS ERLEAPLMS eAPLMS

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        6148  2013-03-28 21:29  MATLAB\.DS_Store
     文件    30168795  2013-03-28 21:32  MATLAB\01_memory_org.mat
     文件    32676908  2013-03-28 21:32  MATLAB\03_bbc_history_org.mat
     文件         699  2013-03-28 21:32  MATLAB\APLMS_AEC_mono.m
     文件        4008  2013-03-28 21:32  MATLAB\ar_m_al1.m
     文件        2606  2013-03-28 21:32  MATLAB\ar_m_al2.m
     文件        2269  2013-03-28 21:32  MATLAB\bformer.m
     文件         658  2013-03-28 21:32  MATLAB\Blind_equalizer.m
     文件        1517  2013-03-28 21:32  MATLAB\Blind_equalizer_l1.m
     文件        1522  2013-03-28 21:32  MATLAB\Blind_equalizer_l2.m
     文件        2088  2013-03-28 21:32  MATLAB\blk_mdlg.m
     文件         407  2013-03-28 21:32  MATLAB\CORLNM2.M
     文件        3912  2013-03-28 21:32  MATLAB\CyclicEqT_eval.m
     文件         239  2013-03-28 21:32  MATLAB\CycPilot.m
     文件         308  2013-03-28 21:32  MATLAB\DCTM.M
     文件         244  2013-03-28 21:32  MATLAB\DFTM.M
     文件         281  2013-03-28 21:32  MATLAB\DHTM.M
     文件         743  2013-03-28 21:32  MATLAB\DOABartlett.m
     文件         858  2013-03-28 21:32  MATLAB\DOAMUSIC.m
     文件         742  2013-03-28 21:32  MATLAB\DOAMVDR.m
     文件         234  2013-03-28 21:32  MATLAB\DSTM.M
     文件        3782  2013-03-28 21:32  MATLAB\DTDexperiment.m
     文件        1389  2013-03-28 21:32  MATLAB\Eigenfir.m
     文件        1619  2013-03-28 21:32  MATLAB\equalizer.m
     文件        3801  2013-03-28 21:32  MATLAB\equalizer_eval.m
     文件         649  2013-03-28 21:32  MATLAB\Ex18_5.m
     文件         409  2013-03-28 21:32  MATLAB\expander.m
     文件         670  2013-03-28 21:32  MATLAB\GBTR.m
     文件        1873  2013-03-28 21:32  MATLAB\gf10_11.m
     文件        1427  2013-03-28 21:32  MATLAB\gf10_9a.m
     文件        1488  2013-03-28 21:32  MATLAB\gf10_9b.m
............此处省略43个文件信息