一维信号短时傅里叶变换域邻域降噪方法(MATLAB)

噪声在人类日常生活中无处不在,其会降低语音信号的质量和可懂度。在低信噪比的恶劣环境中,这种负面影响愈发严重。为了解决这个问题,众多研究人员在过去的几十年里提出了许多降噪算法。

根据原理的不同,降噪算法可大致分为五类:谱减法、最优滤波法、基于统计模型的方法、子空间方法、以及基于机器学习的方法。其中,最优滤波法是较为常用的一类,其将降噪问题转化为最优滤波器的设计问题。根据采用麦克风数量的不同,降噪算法又可分为单通道和多通道两类。此外,根据处理噪声信号的域,这些算法可以分为时域和变换域(如短时傅里叶变换域降噪算法。虽然这些方法已经取得了一定程度上的成功,但降噪仍然是一项具有挑战性的任务。基于相邻帧STFT系数互不相关的假设,早期的STFT域最优降噪方法绝大多数都是滤波增益,即每个频点的滤波器长度为1。然而,由于语音降噪时,语音帧需要有相互重叠,并且语音信号本身是强自相关的,所以当前帧的STFT系数与相邻帧的STFT系数高度相关。基于帧间相关性的不可忽略性,可推导出一组最优滤波器。

鉴于此,采用短时傅里叶变换域邻域降噪方法对一维信号进行降噪,运行环境为MATLAB 2021B。

Matlab 复制代码
[DD,DDT,DDTDD] = defDDt(N); % 2d Order difference operators.

    for k = 1:Nit
        u_old = u;
        
        %% u sub-problem %%
        % Solves a linear system using the conjugate gradient method.
        [u,~] = cgs(rho*DDTDD + lam*speye(N),lam*f+rho*DDT*(v - mu/rho),1e-5,100);
        
        %% v sub-problem %%
        x    = DD*u + mu/rho;
        v    = shrink(x, 1/rho);
        
       %% Update Lagrange multiplier
        mu   = mu + (DD*u - v);
        
        r1   = u - f;
        funcVal(k) = (lam/2)*norm(r1,'fro')^2 + sum(v(:));
    
    end
    
    out.sol = u;
    out.funVal = funcVal(1:k);

end

function [DD,DDT,DDTDD] = defDDt(N)
%Create a first order difference matrix D
e = ones(N,1);
B = spdiags([e -e], [1 0], N, N);
B(N,1) =  1;

D = B; 
clear B;
DD = D*D;
clear D;
% Create the transpose of D
DDT = DD'; %Remember that DT = -D, also called the backward difference.
DDTDD = DD'*DD;
end

function z = shrink(x,r)
z = sign(x).*max(abs(x)- r,0);
end


完整代码:https://mbd.pub/o/bread/ZJyTl5lq
擅长领域:现代信号处理,机器学习,深度学习,数字孪生,时间序列分析,设备缺陷检测、设备异常检测、设备智能故障诊断与健康管理PHM等。
  • 擅长领域:现代信号处理,机器学习,深度学习,数字孪生,时间序列分析,设备缺陷检测、设备异常检测、设备智能故障诊断与健康管理PHM等。
相关推荐
爱分享的飘哥6 分钟前
第七十章:告别“手写循环”噩梦!Trainer结构搭建:PyTorch Lightning让你“一键炼丹”!
人工智能·pytorch·分布式训练·lightning·accelerate·训练框架·trainer
杜子不疼.9 分钟前
《Python学习之第三方库:开启无限可能》
开发语言·python·学习
西工程小巴10 分钟前
实践笔记-VSCode与IDE同步问题解决指南;程序总是进入中断服务程序。
c语言·算法·嵌入式
阿里云大数据AI技术20 分钟前
PAIFuser:面向图像视频的训练推理加速框架
人工智能·机器学习
Tina学编程22 分钟前
48Days-Day19 | ISBN号,kotori和迷宫,矩阵最长递增路径
java·算法
盛世隐者23 分钟前
【深度学习】pytorch深度学习框架的环境配置
人工智能·pytorch·深度学习
说私域24 分钟前
基于开源链动2+1模式AI智能名片S2B2C商城小程序的流量转化策略研究
人工智能·小程序
Moonbit28 分钟前
MoonBit Perals Vol.06: MoonBit 与 LLVM 共舞 (上):编译前端实现
后端·算法·编程语言
青川入梦29 分钟前
MyBatis极速通关上篇:Spring Boot环境搭建+用户管理实战
java·开发语言·mybatis
dundunmm30 分钟前
【每天一个知识点】时间序列聚类
机器学习·支持向量机·聚类·时序·时序聚类