%%%%%%%%%%??數據讀取
ECG=load（‘1record.dat‘）;
l=length（ECG）;
t=[150/l:150/l:30]
x=ECG（4:）;????%可取不同的組
y=x（1:5:end）
figure（1）
plot（t?20*log（abs（y）））;
title（‘含噪心電信號‘）;
xlabel（‘時間（s）‘）;
ylabel（‘幅度（dB）‘）;
axis（[0?15?179?183]）

%%%%%%%%%%??頻域分析
n=4000
m=abs（fft（yn））;
fs=200;
f=fs/n*（0:n-1）;
figure（2）
plot（f?m）;
title（‘心電信號的頻譜圖‘）;
xlabel（‘頻率?f/Hz‘）;ylabel（‘幅值/db‘）;
axis（[0?100?-100?40000]）


%%%%%%%%%%??加入工頻噪聲??時域及頻域分析
x=20*sin（2*50*pi*t）;
y1=y+x;
figure（3）;
subplot（211）;
plot（t?20*log（abs（y1）））;
axis（[0?15?179?183]）
title（‘加工頻噪聲時域波形‘）;
xlabel（‘時間t/s‘）;ylabel（‘幅值/db‘）;
k=abs（fft（y1n））;	
subplot（212）;
plot（fk）;
axis（[0?100?0?40000]）
title（‘加工頻噪聲頻譜圖‘）;
xlabel（‘頻率?/Hz‘）;ylabel（‘幅值‘）

%%%%%%%%%%??工頻噪聲的濾除???巴特沃思
fs=200;
wp=[47?53]*2/fs;
ws=[42?58]*2/fs;
Rp=3;Rs=45;
[Nwn]=buttord（wpwsRpRs）???????????%??N=6
[ba]=butter（Nwn‘stop‘）
figure（4）
[HW]=freqz（ba1024）;%生成濾波器的幅頻響應
k=0:1023;
plot（（fs/2）/1024?*kabs（H））;%輸出濾波器的頻率響應
axis（[0?100?0?1

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

?????文件???????7814??2013-05-30?20:14??ECG.m

?????文件????1800012??2013-04-19?11:02??1record.dat

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

??????????????1807826????????????????????2