| 查看: 242 | 回复: 2 | |||
| 当前主题已经存档。 | |||
[交流]
【讨论】matlab的busy问题
|
|||
|
function [Y] = lle(X,K,d) [D,N] = size(X); fprintf(1,'LLE running on %d points in %d dimensions\n',N,D); % STEP1: COMPUTE PAIRWISE DISTANCES & FIND NEIGHBORS fprintf(1,'-->Finding %d nearest neighbours.\n',K); X2 = sum(X.^2,1); distance = repmat(X2,N,1)+repmat(X2',1,N)-2*X'*X; [sorted,index] = sort(distance); neighborhood = index(2  1+K), ;% STEP2: SOLVE FOR RECONSTRUCTION WEIGHTS fprintf(1,'-->Solving for reconstruction weights.\n'); if(K>D) fprintf(1,' [note: K>D; regularization will be used]\n'); tol=1e-3; % regularlizer in case constrained fits are ill conditioned else tol=0; end W = zeros(K,N); for ii=1:N z = X(:,neighborhood(:,ii))-repmat(X(:,ii),1,K); % shift ith pt to origin C = z'*z; % local covariance C = C + eye(K,K)*tol*trace(C); % regularlization (K>D) W(:,ii) = C\ones(K,1); % solve Cw=1 W(:,ii) = W(:,ii)/sum(W(:,ii)); % enforce sum(w)=1 end; % STEP 3: COMPUTE EMBEDDING FROM EIGENVECTS OF COST MATRIX M=(I-W)'(I-W) fprintf(1,'-->Computing embedding.\n'); % M=eye(N,N); % use a sparse matrix with storage for 4KN nonzero elements M = sparse(1:N,1:N,ones(1,N),N,N,4*K*N); for ii=1:N w = W(:,ii); jj = neighborhood(:,ii); M(ii,jj) = M(ii,jj) - w'; M(jj,ii) = M(jj,ii) - w; M(jj,jj) = M(jj,jj) + w*w'; end; % CALCULATION OF EMBEDDING options.disp = 0; options.isreal = 1; options.issym = 1; [Y,eigenvals] = eigs(M,d+1,0,options); Y = Y(:,2:d+1)'*sqrt(N); % bottom evect is [1,1,1,1...] with eval 0 fprintf(1,'Done.\n'); 这是我从网上下的程序，大概意思是：输入是一个D×N的矩阵， 把每一列看作一个D维的点，然后找每个点的K个近邻点，并建立 起每个点与他的K个近邻点的关系矩阵，最后的输出是一个d×N 的矩阵（d 一个40×100的矩阵，然后运行后得到了20×100的矩阵，也就是 降了20维，但再拿另一个40×100的矩阵 来运行，matlab就busy 了，死机，有如下警告：Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.077831e-016. 请各位高手出手相救，不胜感激！ |
回复此楼» 猜你喜欢
面上本子正文33页，违规吗？会被低分嘛？
已经有10人回复
-
国自然上会要求
已经有9人回复
-
博士申请
已经有5人回复
-
评审有感
已经有6人回复
-
今年审到国自然15份，谈谈感受
已经有16人回复
-
上海大学实验技术岗位非升即走
已经有8人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
2楼2009-07-13 16:22:42
3楼2009-07-15 10:30:17
