【数字信号去噪】基于matlab中值滤波+奇异值分解(SVD)数字信号降噪【含Matlab源码 1021期】
2021/6/18 22:31:45
本文主要是介绍【数字信号去噪】基于matlab中值滤波+奇异值分解(SVD)数字信号降噪【含Matlab源码 1021期】,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
一、简介
基于SVD(奇异值分解)的去噪声技术属于子空间算法的一种。简单的来说我们希望将带噪信号向量空间分解为分别由纯净信号主导和噪声信号主导的两个子空间,然后通过简单地去除落在“噪声空间”中的带噪信号向量分量来估计纯净信号。要将带噪信号向量空间分解为“信号子空间”和“噪声子空间”,可以采用线性代数中的正交矩阵分解技术,特别是奇异值分解(SVD)和特征值分解(EVD)。
二、源代码
%------------------------------中值滤波------------------------------ % 首先通过中值滤波对原数据进行处理 clear all;close all;clc; load data.mat figure(1); orig=Data(5000:9000)'; L=length(orig); plot(1:L,orig); %原始数据 N=2:12; Y1=zeros(length(N),L); Var=zeros(length(N),L); mse=zeros(1,length(N)); for i=1:length(N) %通过计算均方误差判断中值滤波的阶数最优值 Y1(i,:)=medfilt1(orig,N(i)); Var(i,:)=(Y1(i,:)-orig).^2; mse(i)=sum(Var(i,:)); end clear i; mse_best=min(mse); n_best=N(find(mse==mse_best)); x1_best=Y1(find(N==n_best),:); %----------------------奇异值分解的相关操作----------------------------- %首先是映射到相空间 % 相空间矩阵行数 %m=floor(L/2); % 最佳维数基本在P=N/2处的一个邻域内产生,并且在此邻域所取的维数值的降噪效果比较理想,而且能满足工程要求。 % 因此,重构矩阵的结构可以根据振动信号长度N来确定,工程应用中可以取P=N/2。 % 相空间间隔: 通过对时间序列计算自相关函数,取其小于某值域最下值(一般取0.1) % r_k=autocorr(orig,length(orig)-1); % m_k=find(r_k<0.1); tao=2; m=floor((L+tao)/(tao+1)); n=L+tao-m*tao; %n=L+1-m; y=zeros(1,n); X=zeros(m,n); % for i=1:m % y=x1_best(i:n+i-1); % X(i,:)=y; % end for i=1:m y=x1_best(1+(i-1)*tao:n+(i-1)*tao); X(i,:)=y; end clear i; %矩阵的奇异值分解 [U,S,V] = svd(X); sv = diag(S); newSv=sv; [a,b] = size(S); %startPosToZero=18; % 重构阶次(其确定较难) %------------重构阶次的确定方法-奇异值能力差分谱法--------------------- % for i=1:length(sv)-1 % pho(i)=(sv(i)^2-sv(i+1)^2)/(max(sv)^2-min(sv)^2); % end % figure(2) % plot(pho); startPosToZero=11; newSv(startPosToZero:end) = 0; newS = [diag(newSv),zeros(a,b-a)]; newX = U*newS*V'; % 重构相空间矩阵 %降噪后的恢复
三、运行结果
四、备注
版本:2014a
完整代码或代写加1564658423
这篇关于【数字信号去噪】基于matlab中值滤波+奇异值分解(SVD)数字信号降噪【含Matlab源码 1021期】的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-10-01基于Python+Vue开发的医院门诊预约挂号系统
- 2024-10-01基于Python+Vue开发的旅游景区管理系统
- 2024-10-01RestfulAPI入门指南:打造简单易懂的API接口
- 2024-10-01初学者指南:了解和使用Server Action
- 2024-10-01Server Component入门指南:搭建与配置详解
- 2024-10-01React 中使用 useRequest 实现数据请求
- 2024-10-01使用 golang 将ETH账户的资产平均分散到其他账户
- 2024-10-01JWT用户校验课程:从入门到实践
- 2024-10-01Server Component课程入门指南
- 2024-09-30Dnd-Kit学习:新手快速入门指南