小波相关去噪读取程序

你好，之前读到你的相关去噪文章，想自己试一下
但出了问题想咨询一下
我的数据是前8192个点为I路信号，后8192个点为Q路信号
主要感觉这句有问题，我存的I路信号f应该也是个向量吧，想找您咨询一下，麻烦了

     f= data(1:8192,j);%I路 
     s_1 = xden(f,lev,wf);

我的读取程序

%相关去噪并保存去噪信号为h5文件

clc
close all
clear all

addpath('D:/test3')

savepath1 ='D:/test3/data_new.h5'; %创建h5文件(噪声)

for i=0:7

  dataset=strcat('/',num2str(i));
  sprintf('正在处理第%d个数据集',i)
  data = h5read('D:/test3/data.h5',dataset);
  [m,n]=size(data);

  data_new=zeros(m,n);%抑制主信号后的噪声信号

%开启并行加速处理
  parfor j=1:n
        disp(j)
      %分路处理
      lev = 4;
      wf='db3';

     f= data(1:8192,j);%I路 
     s_1 = xden(f,lev,wf);
   

     f=data(8193:end,j);%Q路
     s_2 = xden(f,lev,wf)
     

      %存入信号
     data_new(:,j)=[s_1';s_2'];
     

   end

%保存为h5数据
h5create(savepath1, dataset, [m,n]);
h5write(savepath1, dataset, data_new);


clear data_new
clear data

end

我的报错结果


错误使用 swt (第 62 行)
填充 Wavelet:moreMSGRF:SWT_length_MSG 的内容时出错。不允许浮点数作为填充内容。它们应转换为字符向量。

出错 xden (第 6 行)
[swa,swd]=swt(f,n,Lo_D,Hi_D);%swd是细节系数，swa是近似系数

出错 main_new (第 22 行)
  parfor j=1:n