Rbo.zip_VPU_avf_心电_心电 分析

#include<iostream.h> 007vb\]Il #include<math.h> &!D}M> ^ #include<stdio.h> HIq6Gtl< int sr=256;//说明：采样率为sr，多选择II导联作为心电分析数据，该分析要用到的数据记录至少为1分钟的心电波形数据 ]L,iahaP0 void main() (E
3BGX[% { const int N=15420;//N=sr*time M./3@A const int n=(int)sr*0.04,yanshi=(int)sr*0.2,back=(int)sr*0.04;//n=sr*0.04,yanshi=sr*0.2?,back=sr*0.04 PD
l=!>u const int M=N;//负斜率找寻阵维数M=N-n "q"$bY9<3 const int R=200; UUFW=m` int Q=(int)sr*1.2; T/-%VpU;# const int RR1=1.2*sr,RR4=(int)sr*1.6;//RR1=1.2*sr,RR2=1.5*sr,RR3=0.5*sr,RR4=1.6*sr tU {1.&D float data[N][13],time,I,II,III,avr,avl,avf,v1,v2,v3,v4,v5,v6; { ft(y.#C= const int F=2; ^qH|x>Z int i=0; while(1) Hqh(-*qS { FILE *fp; if((fp=fopen("f:\\matlab\\lead12for1minute.txt","r"))==NULL) R_ ;1@{a { cout<<"can't open this file\n"; break; } } else { for(i=0;i<N;i++) { fscanf(fp,"%f %f %f %f %f %f %f %f %f %f %f %f %f\n",&time,&I,&II,&III,&avr,&avl,&avf,&v1,&v2,&v3,&v4,&v5,&v6); W6cw=<uT data[0]=time; data[1]=I; data[2]=II; data[3]=III; data[4]=avr; data[5]=avl; data[6]=avf; data[7]=v1; data[8]=v2; data[9]=v3; data[10]=v4; data[11]=v5; data[12]=v6; } } //从文本文挡中读出各个导联的心电波形数据存入data数组的各列当中 b<T>En6 fclose(fp); int note1=0,j=0,note2=0,ii; float xx,km1[10],kv1[10],kv,km=0; for(i=0;i<10;i++) { km1=0; kv1=0; } for(i=0;i<N;i++) { if(data[F]>km) { note1=1; } } if(note1==0) { cout<<"该导联心电数据不适合做心电分析！请选择其它导联数据，推荐使用II导联心电。"<<endl; \V#.lS break; } else {  double xr[M]; for(i=0,j=0;j+n<N;i++,j++) { xr=data[j+n][F]-data[j][F]; // cout<<xr<<"\t" }  for(i;i<M;i++) xr=0; for(j=0,i=0;j<10;j++) for(ii=0,kv1[j]=xr;ii<Q;ii++,i++) { if(xr<kv1[j]) { kv1[j]=xr; note2=i; } if(data[F]>km1[j]) { km1[j]=data[F]; note1=i; } } for(j=0,kv=0,km=0;j<10;j++) OO /TE { kv+=kv1[j]; km+=km1[j]; } kv/=10;//求出最大负斜率和最大幅值 Q1"m/LA* km/=10; cout<<"心电最大