clc
clear
%%
fs=500;%44.1e3;%采樣率
samptime=8;%采樣時間ECG調制用4秒，其他用0.5秒，便于觀察
fm=150;%調制信號頻率
am1=1;%調制信號幅度
fc=1e3;%載波頻率
[z0t0]=ReadECGData（samptimefs）;%讀取心電信號txt文件
Amplify=800;%電路放大倍數
[~z1t]=ReadMIT_BIH（‘117‘‘.\MIT-ECG?data‘samptimefs）;
z=z1/1000*Amplify;%單位換算為V
figure（1）
plot（tz）
title（‘ECG原信號‘）;
data=z‘;
differ=zeros（40001）;
dif_max=0;
for?i=3:1:3998
????differ（i1）=data（i+11）-data（i-11）+2*data（i+21）-2*data（i-21）;%求導數
end???
dif_max=max（differ）;
thresh=dif_max*0.375;??????????%閾值為導數最大值的0.375倍
NumOfR=0;
for?i=1:4000
????if??（data（i1）>thresh）&&（differ（i1）*differ（i+1）<0）?
????????NumOfR=NumOfR+1;
????????R_position（NumOfR1）=i;
????end
end

for?i=1:NumOfR-1
????DR_R（i1）=60*500/（R_position（i+11）-R_position（i1））;
end
figure（2）
?plot（（1:4000）‘dataR_positiondata（R_position）‘rs‘）;
?title（‘117號ECG信號R波定位‘）
?%%
%Q波定位
NumOfQ=0;
for?i=1:NumOfR
????minQvalue=R_position（i）-40;
????for?m=1:40
????????if?（data（minQvalue）>data（R_position（i）-40+m））
????????minQvalue=R_position（i）-80+m;
???????end
????end
????NumOfQ=NumOfQ+1;
????Q_position（NumOfQ1）=minQvalue;
end
figure（3）
plot（（1:4000）‘dataQ_posit