资源简介

对lorenz信号的(svd)奇异值分解降噪,

资源截图

代码片段和文件信息



function jiangzao=qiyizhifenjie(datamtauEigenvaluecountK)
%相重构Eigenvaluecount
PhaseSpacedata=r(datamtau);

Length = size(PhaseSpacedata2); %相点的个数


%寻找领域 

for i=1:Length
     x=PhaseSpacedata(:i);
     dist=sum((PhaseSpacedata-repmat(x[1Length])).^21);
      [dist_sortJ]=sort(dist);
    
      Ii=J(1:K);
      P=PhaseSpacedata(:Ii)-repmat(mean(PhaseSpacedata(:Ii)2)[1K]);

%svd
   [U S V]=svd(P); 
   Temp=diag(S);
   [mn]=size(S);
   SS=zeros(mn);
   Index=1;
   for  k=1:m
       for j=1:n
          if (k==j)
             SS(kj)=SS(kj)+Temp(Index);
             Index=Index+1;
          end;
      end; 
         if(Index>Eigenvaluecount)
           break;
         end; 
       
   end;
   pp=U*SS*V‘;

%代替
  PhaseSpacedata(:Ii)

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1009  2020-09-21 02:28  ╨┬╜¿╬─╝■╝╨\qiyizhifenjie.m
     文件        1086  2020-09-21 02:28  ╨┬╜¿╬─╝■╝╨\qiyizhifenjie22.m

评论

共有 条评论