陈宇鹏 2023-04-21 09:34 采纳率: 0%
浏览 11

小波相关去噪读取程序

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

  • 写回答

1条回答 默认 最新

报告相同问题?

问题事件

  • 创建了问题 4月21日

悬赏问题

  • ¥15 关于smbclient 库的使用
  • ¥15 微信小程序协议怎么写
  • ¥15 c语言怎么用printf(“\b \b”)与getch()实现黑框里写入与删除?
  • ¥20 怎么用dlib库的算法识别小麦病虫害
  • ¥15 华为ensp模拟器中S5700交换机在配置过程中老是反复重启
  • ¥15 java写代码遇到问题,求帮助
  • ¥15 uniapp uview http 如何实现统一的请求异常信息提示?
  • ¥15 有了解d3和topogram.js库的吗?有偿请教
  • ¥100 任意维数的K均值聚类
  • ¥15 stamps做sbas-insar,时序沉降图怎么画