你好,之前读到你的相关去噪文章,想自己试一下
但出了问题想咨询一下
我的数据是前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