用Android写一个ECG的小波变换

算法 //实现小波变换的ECG算法 public class ECGWavelet { // 小波变换迭代次数 private int iter; // 小波变换的基函数，默认使用Haar小波 private double[] wavelet = {1/Math.sqrt(2), 1/Math.sqrt(2)}; // 将原始信号存储到一个double数组中 private double[] signal; public ECGWavelet(int it, double[] sig) { iter = it; signal = sig; } /** * 进行小波变换，将原始信号进行分解，得到各个子带信号。 */ public void waveletTransform() { for (int i = 0; i < iter; i++) { int len = signal.length >> i + 1;// 每次迭代后信号长度减半 // 扩充基函数（wavelet）长度以适应当前信号长度 double[] extWavelet = new double[len]; for (int j = 0; j < len / 2; j++) { extWavelet[j] = wavelet[0];// 奇数位置使用基函数的左半部分 extWavelet[j + len / 2] = wavelet[1];// 偶数位置使用基函数的有半部分 } // 进行小波处理 |------------------|------------------|-------|-------|-------|-------| |-----|-----|-----|-----| |------------||------------||------------||------------||-----------||-----------| |-a-b-c-d-e-f-g-|h---------i-|j---------k-|l---------m-|n---------o-|p---------q---> a b c d e f g h i j k l m n o p q r s t u v w x y z |---------+---------+--------+--------+--------+--------+ |---w---x---y---z--+--+--+--+ 平移采样 |---------------w---------------x---------------y---------------z------> w x y z a b c d e f g h ... 通过平移采样得到斜对角形采样图 \ \ \ \ \ \ \_\_\_\_\_\_\ w x y z 通过斜对角形采样图得到子带信号 A B C D E F G H I J K L M N O P Q R S T U V W X Y Z A B C D E F G H I J K L M N O P Q R S T U V W X Y Z A B C D E F G H I J K L M N O P Q R S T U V W X Y Z A B C D E F G H I J K L M N O P Q R S T U V W X Y Z A B C D E F G H I J K L M N O P Q R S T U V W X Y Z A B C D E F G H I J K L M N O P Q R S T U V W X Y Z 获得子带A、B、C、D……Z共26个子带信号。 由于子带A、B、C……Z是相互独立的，所以它们之间不会相互干扰。 由此我们便实现了对原始ECG信号的小波处理。 最后将26个子带合成一条完整的ECG信