## 源码简介

t=[0:0.002:5.998]';      % 时间轴
d = 0:1/1.8:6;           % 尖峰信号
nNoiseECG=sin(2*pi*0.24*t)+2*pulstran(t,d,'tripuls',0.1,1)+0.1*randn(3000,1); % 叠加正弦漂移和尖峰信号以及随机白噪声
plot(t,nNoiseECG)
xlabel('t(s)')
ylabel('Voltage(v)')
title('Original Signal')

－－－－－－－

Heart Beat Period Detection

In clinic, it is frequently required to detect heart rate. The heart rate is usually computed from RR-interval, which is obtained from ECG.
However, while ECG measuring, noises are frequently superposed.

Process the noised signal to obtain ECG signal
Propose algorithm to acquire RR intervals
The Heart Rate is calculated from averaging 10 RR-intervals

The signal is given in the data file, ‘nNoiseECG.txt’, the sampling frequency is 500Hz

Finish the tasks individually and independently
Write a project report in English. Give and explain the full processes with Matlab programs
－－－－－－－－

t=[0:0.002:5.998]';      % 时间轴
d = 0:1/1.8:6;           % 尖峰信号
nNoiseECG=sin(2*pi*0.24*t)+2*pulstran(t,d,'tripuls',0.1,1)+0.1*randn(3000,1); % 叠加正弦漂移和尖峰信号以及随机白噪声
plot(t,nNoiseECG)
xlabel('t(s)')
ylabel('Voltage(v)')
title('Original Signal')

Y = fft(nNoiseECG,16384);%进行16384个点的fft变换
Pyy = Y.* conj(Y) / 16384;  %功率谱转换
f = 500*(0:8192)/16384;
plot(f,Pyy(1:8193))
title('Frequency content of y')
xlabel('frequency (Hz)')
axis([0 10 0 1500]);
grid on;

a=filter(Num,1,nNoiseECG);% 其中Num即为滤波器的极点参数，1为零点参数，a即滤波后的数据。
plot(t,a) % 滤波后画图
xlabel('t(s)')
ylabel('Voltage(v)')
title('Signal After Filtering')

for i=1:3000
if (a(i)<0.5)
a(i)=0;
end
end
b=find(diff(sign(diff(a)))<0)+1;% 求解极值点序号
plot(t, a,t(b), a(b),'r*') % 画出去噪声后的波形，由于滤波的关系，前一秒钟的信号异常，因此只画1s以后的数据
xlabel('t(s)')
ylabel('Voltage(v)')
axis([1 6 0 2]);
title('Final Signal')
grid on

for i=1:3000
if (a(i)<0.5)
a(i)=0;
end
end
b=find(diff(sign(diff(a)))<0)+1;% 求解极值点序号
plot(t, a,t(b), a(b),'r*')% 画出去噪声后的波形，由于滤波的关系，前一秒钟的信号异常，因此只画1s以后的数据
xlabel('t(s)')
ylabel('Voltage(v)')
axis([1 6 0 2]);
title('Final Signal')
grid on

(t(b(9))-t(b(1)))/10
10/(t(b(9))-t(b(1)))

## 本站长期招聘程序代写高手，欢迎加入华南地区matlab团队

 想创业却没有经验的人 无论你是否有过网上开店的经验，都可以随时联系在线客服，建立自己独立的网站 想开网店却不知道如何入手 淘宝创业成本低而且风险小，如果想开淘宝店的朋友可以联系在线客服。 想兼职创业，却不擅长交际与服务的人 在家创业月入5000元。网站程序+百套群发工具+网赚资料+域名+空间+本站终身代理资格，这样你网赚的条件全具备了。每天3小时管理、推广、收钱。 缺乏能快速赢利型产品的人 导入多种最新流行营销软件+网赚教程，让入驻者轻松加盟、复制有效成交技巧、快速赚钱。
·上一源码：信号变换