资源简介
对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
评论
共有 条评论