【信号去噪】基于奇异分解信号去噪matlab源码

登录后复制

%----------------------------------------------------------------------------
 
%----------------------------------------------------------------------------
clear;
 
% 调用MATLAB中含有噪声的数据文件 leleccum; 
load leleccum; 
 
index=1:3000;
x=leleccum(index);
N=8;
slength = length(x);
M=slength-100;
subplot(221);plot(x(1:M));
title('原始信号');
 
% 形成数据矩阵；
Signal=zeros(N,M);
for i=1:N
   Signal(i,:)=x(i:M+i-1);
end
 
% 对数据矩阵作特征值分解；
[U, S, V]=svd(Signal);
d=diag(S(1:N,1:N));
 
subplot(222);stem(d);
title('特征值');
for i=1:N
   if d(i)<mean(d)
      d(i)=0;
   end
end 
stemp=S;
stemp(1:8,1:8)=diag(d);
Sf=U*stemp*V';
subplot(223);plot(Sf(1,:));
title('滤波之后的信号；阈值为特征值的平均值');
d=diag(S(1:N,1:N));
for i=1:N
      if d(i)<=median(d)
      d(i)=0;
   end
end
stemp=S;
stemp(1:8,1:8)=diag(d);
Sf=U*stemp*V';
subplot(224);plot(Sf(1,:));
title('滤波之后的信号；阈值为特征值的中值');

• 1.
• 2.
• 3.
• 4.
• 5.
• 6.
• 7.
• 8.
• 9.
• 10.
• 11.
• 12.
• 13.
• 14.
• 15.
• 16.
• 17.
• 18.
• 19.
• 20.
• 21.
• 22.
• 23.
• 24.
• 25.
• 26.
• 27.
• 28.
• 29.
• 30.
• 31.
• 32.
• 33.
• 34.
• 35.
• 36.
• 37.
• 38.
• 39.
• 40.
• 41.
• 42.
• 43.
• 44.
• 45.
• 46.
• 47.
• 48.
• 49.
• 50.