zw_mamayu1314-9409158-Untitled.zip

zw_mamayu1314-9409158-Untitled.zip

clear all;close all;
%绘制原始心电信号的时域图和频谱图
ECG=load（‘C:\Users\Administrator\Desktop\renwushu\ecg101.txt‘）;
time=ECG（:1）;
amplitude=ECG（:2）;
figure（1）;
subplot（211）;plot（timeamplitude）;title（‘原始心电信号的时域图‘）;xlabel（‘时间/s‘）;ylabel（‘幅值/mv‘）;axis（[010-3.53.5]）;grid on;
n=length（time）;
an=fft（amplitude）;
T=ceil（max（time））;
Fs=n/T;
f=Fs/n*（0:n-1）;
subplot（212）;plot（fabs（an））;title（‘原始心电信号的频谱图‘）;xlabel（‘频率/Hz‘）;ylabel（‘|X（k）|‘）;axis（[0180-1400]）;grid on;
% 冲击响应不变法设计IIR带通滤波器设数字指标fp=[1 30]fst=[0.5 35]Rp=2As=5.
fp=[1 30];fst=[0.5 35];Rp=2;As=5;
Wp=fp*2*pi;Wst=fst*2*pi;
[N1Wc]=buttord（WpWstRpAs‘s‘）;
[zpk]=buttap（N1）;   [ba]=zp2tf（zpk）;
W1=Wp（2）-Wp（1）;   W2=sqrt（Wp（1）*Wp（2））;
[bsas]=lp2bp（baW2W1）;
[bzaz]=impinvar（bsasFs）;
[H1w1]=freqz（bzazFs）;
f1=w1/pi*Fs/2; 
figure（2）;
subplot（211）;plot（f120*log10（abs（H1）））;title（‘IIR带通滤波器幅频特性‘）;xlabel（‘频率/Hz‘）;ylabel（‘20log10|H1（K）|‘）;grid on;
subplot（212）;plot（f1angle（H1））;title（‘IIR带通滤波器相频特性‘）;xlabel（‘频率/Hz‘）;ylabel（‘arg[H1（K）]‘）;grid on;
figure（3）
h1=impz（bzaz）;
subplot（211）; plot（h1）;title（‘IIR带通滤波器的单位抽样响应‘）;xlabel（‘n‘）;ylabel（‘h1‘）;grid on;
%IIR带通滤波器滤波之后
figure（4）
y1=filter（bzazECG（:2））;
subplot（211）;plot（timey1）;title（‘IIR带通滤波器滤波后心电信号的时域图‘）;xlabel（‘时间/s‘）;ylabel（‘幅值/mv‘）;axis（[010-3.53.5]）;grid on;
subplot（212）;plot（fabs（fft（y1）））;title（‘IIR带通滤波器滤波后心电信号的频谱图‘）;xlabel（‘频率/Hz‘）;ylabel（‘|Y1（k）|‘）;axis（[0180-1400]）;grid on;
%窗函数法设计FIR高通滤波器设fp=0.3fs=0.8As=5选用矩形窗.
wp=2*pi*0.3/Fs;ws=2*pi*0.8/Fs;dw=ws-wp;
N2=ceil（1.8*pi/dw） ;
if rem（N22）==0       
    N2=N2+1;
end
wc=（ws+wp）/2;
window=boxcar（N2）;
h2=fir1（N2-1wc‘high‘window）;
figure（5）
[H2w2]=freqz（h21）; 
f2=w2/pi*Fs/2; 
subplot（211）;plot（f220*log10（abs（H2）））;title（‘FIR高通滤波器幅频特性‘）;xlabel（‘频率/Hz‘）;ylabel（‘20log10|H（K）|‘）;grid on;
subplot（212）;plot（f2angle（H2））;title（‘FIR高通滤波器相频特性‘）;xlabel（‘频率/Hz‘）;ylabel（‘arg[H（K）]‘）;grid on;
figure（6）
subplot（211）; plot（h2）;title（‘FIR高通滤波器的单位抽样响应‘）;xlabel（‘n‘）;ylabel（‘h2‘）;grid on;
%FIR高通滤波器滤波之后
figure（7）
y2=filter（h21ECG（:2））;
subplot（211）;plot（timey2）;title（‘FIR高通滤波器滤波后心电信号的时域图‘）;xlabel（‘时间/s‘）;ylabel（‘幅值/mv‘）;axis（[010-3.53.5]）;grid on;
subplot（212）;plot（fabs（fft（y2）））;title（‘FIR高通滤波器滤波后心电信号的频谱图‘）;xlabel（‘频率/Hz‘）;ylabel（‘|Y2（k）|‘）;axis（[0180-1400]）;grid on;
%对原信号2倍插值
figure（8）
x=0:n-1;
y=ECG（:2）;
xi=0:0.5:n-1;
yi_linear=interp1（xyxi）;
t=10/3600*xi;
subplot（211）;plot（tyi_linear）;title（‘2倍插值心电信号的时域图‘）;xlabel（‘时间/s‘）;ylabel（‘幅值/mv‘）;axis（[010-3.53.5]）;grid on;
n2=2*n;
y2=fft（yi_linearn2）;
Fs2=n2/T;
f3=Fs2/n2*（0:n2-1）;
subplot（212）;plot（f3abs（y2））;title（‘2倍插值心电信号的频谱图‘）;xlabel（‘频率/Hz‘）;ylabel（‘|X2（k）|‘）;axis（[0360-1800]）;grid on;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3024  2019-07-14 11:02  \mamayu1314-9409158-Untitled.m
     文件          36  2019-07-14 11:02  \no.txt