能量

shuizhengqi

我计算脉搏程序的谱能比，程序如下，运行结果显示 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;

shuizhengqi

[N1,N2]=period(s); 提取信号的单周期波形