clear?all?;
load?foetal_ecg.dat?
abdominal?=foetal_ecg（:2:6）;?%loading??abdominal?signals??
thoraic?=foetal_ecg（:7:9）?;?%loading??thoraic?signals?
d=?mean（abdominal2）;?%Average?of?abdominal?signals
x=mean（thoraic2）;?%Average?of?thoraic?signals
M=2500;?%數(shù)據(jù)長度
N?=?8;?%order?of?the?filter
mu?=?.0000006;?%value?of?step?size
iter=length（x）;?%迭代次數(shù)默認(rèn)為輸入信號長度
w=zeros（NM）;?%每一行代表一個加權(quán)參量，每一列代表一次迭代，初始為?0
e=zeros（M1）;?%誤差序列，?e（k）?表示第?k?次迭代時預(yù)期輸出與實際輸入的誤差
%將?y?信號右移?N位，以便進(jìn)行迭代時的加權(quán)運算
x1=zeros（N1）;
for?i=N+1:（N+M-1）
x1（i）=x（i+1-N）;
end
%調(diào)整濾波器系數(shù)的?LMS算法
for?n=2:iter;
x2=x1（n:1:n+N-1）;
y（n）=w（:n-1）.‘*x2;?%每次取?N個參考信號與權(quán)系數(shù)相乘得到噪聲估計
e（n）=d（n）-y（n）;?%噪聲抵消的輸出
w（:n）=w（:n-1）+2*mu*x2*e（n）;?%權(quán)矢量迭代
end?;
%最后穩(wěn)定后的胎兒心電信號
for?n=2:iter;
x2=x1（n:1:n+N-1）;
y（n）=w（:end）.‘*x2;?%每次取?N個參考信號與最后穩(wěn)定權(quán)系數(shù)相乘得到噪聲估計
e（n）=d（n）-y（n）;?%噪聲抵消的輸出
end?
%輸出三個心電波形
figure（1）;
subplot（311）;
plot（x）;
title（?‘?母親心電信號?‘?）;
subplot（312）;
plot（d）;
title（?‘?混合心電信號?‘?）;
subplot（313）
plot（e）;
title（?‘?濾波后的胎兒心電信號?‘?）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????227500??2018-08-30?09:27??LMS?FECG\foetal_ecg.dat

?????文件???????1212??2018-08-30?09:27??LMS?FECG\yxj_lms.m

?????目錄??????????0??2018-09-21?09:45??LMS?FECG

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

???????????????228712????????????????????3