%?LMS?Adaptive?Noise?cancellation

%-----Filter?Parameters-----;
M?=?20;
mu?=?0.05;
e_max?=?400;????????????????%-----maximum?of?epochs

%-----Contants-----
pi?=?3.14;
Fs?=?0.01;??????????????????%-----signal?frequency
Fn?=?0.05;??????????????????%-----noise?frequency

%-----Initialize-----
w?=?（randn（1M）?-?randn（1M））/100;
d?=?zeros（1M）;
u?=?zeros（1M）;
u_out?=?zeros（1e_max-M）;
f_out?=?zeros（1e_max-M）;

%-----Generate?desired?signal?and?input（signal+noise）-----
for?t=1:M-1
????d（t）?=?sin（2*pi*Fs*t）;??????????????????????????????????????????????????????????????????%%%???d--desired?signal;?u--mixed?signal
????u（t）?=?d（t）?+?0.5*sin（2*pi*Fn*t）?+?0.09*randn;??????????????????????????????????????????%%%???du的1~M-1節(jié)拍賦值
end
t?=?M;
epoch?=?0;

while?epoch????input?=?sin（2*pi*Fs*t）;
????for?i=2:M???????????????????????%-----shift?new?input?into?array????????????????????????%%%???移位，準(zhǔn)備下次迭代
????????d（M-i+2）?=?d（M-i+1）;
????????u（M-i+2）?=?u（M-i+1）;
????end
????d（1）?=?input;???????????????????????????????????????????????????????????????????????????%%%???d--這次的desired?signal;?u--這次的mixed?signal;（u可以得到，而d是不可得到的；）
????u（1）?=?input?+0.5*sin（2*pi*Fn*t）+0.09*randn;????????????%-----add?undesired?freq?&?random?noise
????u_out（t-M+1）?=?u（1）;???????????????????????????????????????????????????????????????????%%%???
????
????output?=?dot（wu）;??????????????%-----compute?filter?output?????????????????????????????%%%???這次的輸出僅取決于w;
????f_out（t-M+1）?=?output;
????
????%-----LMS?algorithm-----
????e?=?d（1）?-?output;??????????????%-----compute?error-----????????????????????????????????%%%???其實(shí)d（1）不可得；e也不可得；
????
????for?n=1:M???????????????????????%-----update?weights
????????w（n）?=?w（n）?+?mu*u（n）*e;????????????????????????????????????????????????????????????%%%???w的更新；需要事先知道e；
????end
????
????int（t-M+1）?=?u（1）;
????out（t-M+1）?=?output;
????err（t-M+1）?=?e;
????
????t?=?t+1;
????epoch?=?epoch?+?1;
????
????%-----plot?noise?and?filtered?signal-----
????figure（1）;
????subplot（211）plot（tu（1））axis（[0?e_max?-2.5?2.5]）title（‘LMS?Filter?Input（Siganal?+?Noise）‘）drawnowhold?on
????subplot（212）plot（toutput）axis（[0?e_max?-2.5?2.5]）title（‘LMS?Filtered?Signal‘）drawnowhold?on
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????728844??2001-01-01?01:32??LMS\LMS?DSP實(shí)現(xiàn).pdf

?????文件???????3638??2001-01-01?01:33??LMS\mls.m

?????文件???????2205??2005-03-22?16:48??LMS\LMS?源程序.txt

?????文件???????2366??2001-01-01?01:35??LMS\lms.m

?????文件???????2154??2001-01-01?01:35??LMS\lms1.m

?????文件????????830??2001-01-01?01:36??LMS\my_lms.m

?????文件???????2697??2009-06-13?21:31??LMS\Chap_5_MVDR.rar

?????文件?????766359??2009-06-13?21:32??LMS\DSP-Matlab_Programs.rar

?????文件???????7656??2009-06-13?21:32??LMS\lms?DSP.rar

?????文件?????130472??2009-06-13?21:32??LMS\matlab_algorithms.rar

?????文件?????130756??2009-06-13?21:32??LMS\MATLAB_Arithmetic.rar

?????文件???????8142??2009-06-13?21:32??LMS\mud-code.rar

?????目錄??????????0??2001-01-02?09:02??LMS

-----------?---------??----------?-----??----

??????????????1786119????????????????????13