m0_65124804 2023-06-06 21:15 采纳率: 0%
浏览 10

用Matlab使用wdencmp进行去燥

以下代码产生一个污染的正弦曲线。
init=1000;
[xref,x]=wnoise(5,11,7,init);
(1)对曲线使用'sym4'小波进行5层分解,使用wdencmp全局阈值对上述污染的曲线进行去噪;
(2)对曲线使用'sym4'小波进行5层分解,使用wdencmp多层阈值对上述污染的曲线进行去噪,使用软阈值。你的程序代码和运行结果截图粘贴在下面。

  • 写回答

1条回答 默认 最新

  • CSDN-Ada助手 CSDN-AI 官方账号 2023-06-07 12:06
    关注
    clc;
    clear;
    %加载信号
    load leleccum;
    %采样点
    s=leleccum(300:2000);
    %画出原始信号
    subplot(2,2,1);plot(s);title('原始信号');grid on;axis tight;
    %% 用sym4小波对原始信号进行5层小波分解并提取系数
    [c,l]=wavedec(s,5,'sym4');
    ca5=appcoef(c,l,'sym4',5);
    cd5=detcoef(c,l,5);
    cd4=detcoef(c,l,4);
    cd3=detcoef(c,l,3);
    cd2=detcoef(c,l,2);
    cd1=detcoef(c,l,1);
    %对信号进行强制性去噪处理:将每层的高频系数置零
    cdd5=zeros(1,length(cd5));
    cdd4=zeros(1,length(cd4));
    cdd3=zeros(1,length(cd3));
    cdd2=zeros(1,length(cd2));
    cdd1=zeros(1,length(cd1));
    c1=[ca5 cdd5 cdd4 cdd3 cdd2 cdd1];
    s1=waverec(c1,l,'sym4');
    subplot(2,2,2);plot(s1);title('强制去噪后信号');grid on;axis tight;
    %% 用默认阈值对信号进行去噪
    % 用ddencmp函数获得信号的默认阈值
    [THR,SORH,KEEPAPP]=ddencmp('den','wv',s);
    s2 =wdencmp('gbl' ,s,'sym4',5,THR,SORH,KEEPAPP);
    subplot(2,2,3);plot(s2);title('ddencmp/wdencmp默认阈值去噪后信号');grid on;axis tight;
    %% 用多层软阈值去噪,利用wthresh函数产生阈值
    thr=thselect(s,'rigrsure');
    softd5=wthresh(cd5,'s',thr);
    softd4=wthresh(cd4,'s',thr);
    softd3=wthresh(cd3,'s',thr);
    softd2=wthresh(cd2,'s',thr);
    softd1=wthresh(cd1,'s',thr);
    c2=[ca5 softd5 softd4 softd3 softd2 softd1];
    s3=waverec(c2,l,'sym4');
    subplot(2,2,4);plot(s3);title('多层软阈值去噪后信号');grid on;axis tight;
    

    以下为运行结果截图: 去噪结果截图

    注意这里有一个新的函数thselect,需要将softd5-softd1的阈值选择为信噪比最优阈值,可以利用该函数进行选择,具体用法可参考Matlab文档。

    评论

报告相同问题?

问题事件

  • 创建了问题 6月6日

悬赏问题

  • ¥15 App的会员连续扣费
  • ¥15 不同数据类型的特征融合应该怎么做
  • ¥15 用proteus软件设计一个基于8086微处理器的简易温度计
  • ¥15 用联想小新14Pro
  • ¥15 multisim中关于74ls192n和DSWPK开关仿真图分析(减法计数器)
  • ¥15 基于8255的交通灯设计
  • ¥15 w3wp,exe 中发生未处理的 Microsoft ,NETFramework 异常。
  • ¥20 C51单片机程序及仿真(加减器)
  • ¥15 AQWA | 水动力分析 二阶波浪力
  • ¥15 c语言题目:成绩管理系统