xinlv.rar_心率_心电_心电图数据

clear clc %% 绘制图像 data=xlsread('0625.xlsx'); % x=data(:,1);%第一列的全部 y=data(:,2);%第2列的全部 % x=data(1:6778,1); %1到6778行的第一列 y=data(10:5918,2);%1到6778行的第2列 %% data=xlsread('0625.xlsx'); x=data(:,3);%第一列的全部 y1=data(:,1);%第2列的全部 x=data(1:5909,3); %1到6778行的第一列 y1=data(10:5918,1);%1到6778行的第2列 % plot(x,y,'l')%红色显示曲线 % % plot(x,y1,'b')%绘制蓝色余弦曲线 % hold on %图形保持 %% % Fs=256; % Sampling Frequency % N=32; % Order % Fc1=20; % First Cutoff Frequency % Fc2=30; % Second Cutoff Frequency % flag='b'; % Sampling Flag % win=N+1; % Create the window vector for the design algorithm. % b= fir1(N,[Fc1 Fc2]/(Fs/2),'bandpass',win,flag); % Calculate the coefficients using the FIR1 function. % filteredA = filtfilt(b,1,data); %% % time_watch =1; % RSF=data(time_watch*60*256:(time_watch+1)*60*2); % RSF=y1; % Threshold = (max(RSF) - min(RSF))*0.7 + min(RSF); % % % [R_pks,R_locs]=findpeaks(RSF,256,'MinPeakDistance',0.5,'MinPeakHeight',Threshold); % % % plot(RSF); % % title(['第',num2str(time_watch),' 分钟 EEG']) % hold on % plot(R_locs*500,R_pks,'x') %% 心率 % [maxv,maxl]=findpeaks(y1,'minpeakdistance',100); %maxv峰峰值点 ?maxl：峰峰值点对应的位置 ? 最小间隔=100 ? ? ?? [maxv1,maxl1]=findpeaks(y,'minpeakdistance',100); hold on; plot(x,y1,'b');%绘制原波形 plot(x,y,'g'); legend('疲劳驾驶','气味干预后'); % plot(maxl,maxv,'*','color','R');%绘制最大值点 % plot(maxl1,maxv1,'*','color','y'); length(maxv); length(maxv1);