|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我计算脉搏程序的谱能比,程序如下,运行结果显示 s1=[s60 s61 s51 s64 s65 s53 s54]; 有误,不知错在那里,望大家指导一二。
[N1,N2]=period(s);
N=N2-N1;
s=s(N1:N2);
t=wpdec(s',6,'coif3','shannon');
s60=wpcoef(t,[6,0]); %0-1Hz
s61=wpcoef(t,[6,1]); %1-2Hz
s51=wpcoef(t,[5,1]); %2-4HZ
s64=wpcoef(t,[6,4]); % 4-5Hz
s65=wpcoef(t,[6,5]); %5-6Hz
s53=wpcoef(t,[5,3]); % 6-8Hz
s54=wpcoef(t,[5,4]); %8-10Hz
s55=wpcoef(t,[5,5]); %10-12HZ
s43=wpcoef(t,[4,3]); % 12-16Hz
s21=wpcoef(t,[2,1]); %16-32Hz
s34=wpcoef(t,[3,4]); % 32-40Hz
s10=wpcoef(t,[1,0]); %0-32Hz
s11=wpcoef(t,[1,1]); %32-64Hz
s1=[s60 s61 s51 s64 s65 s53 s54]; %0-10Hz
s3=[s55 s43 s21 s34]; %10-40Hz
E1=s1*s1';
E3=s3*s3';
SER=E1/E3; |
|