s******f
发帖数: 3984 | 1 现在在用svd的方法算一种correlation，具体方法很简单是这样的：
1）假设有2个N by M 的矩阵（暂且这样叫把！，x,y）， p1，p2分别是x，y的svd分解
的左矩阵 （[u S v] = svd(x) 中的u）。
2）再对p1'*p2 求是svd, 取其singular value。
问题： 如果在算p1,p2不做任何截取的话（因为一般情况如果N很大，我会做截取，取
前d个最大的singular value对应的column，但是N若很小，我就没有做）， 算得的
singular value全是1。为啥？
本人很土，真心求教。 |
s******f
发帖数: 3984 | |
j****s
发帖数: 156 | 3 p1是unitary，p2是unitary，所以p1*p2还是unitary的。所以singular value 都是1
【在 s******f 的大作中提到】
: 现在在用svd的方法算一种correlation，具体方法很简单是这样的：
: 1）假设有2个N by M 的矩阵（暂且这样叫把！，x,y）， p1，p2分别是x，y的svd分解
: 的左矩阵 （[u S v] = svd(x) 中的u）。
: 2）再对p1'*p2 求是svd, 取其singular value。
: 问题： 如果在算p1,p2不做任何截取的话（因为一般情况如果N很大，我会做截取，取
: 前d个最大的singular value对应的column，但是N若很小，我就没有做）， 算得的
: singular value全是1。为啥？
: 本人很土，真心求教。
|
s******f
发帖数: 3984 | 4 但是截取了以后，P1*P2不也是unitary么？ 为什么这个时候singular value就不为全1
了呢 |
s******f
发帖数: 3984 | 5 我说的截取意思是取P1的某n列， P1'*P2 不也是n by n的unitary matrix..为什么此
时singular value就不全1了
【在 j****s 的大作中提到】
: p1是unitary，p2是unitary，所以p1*p2还是unitary的。所以singular value 都是1
|
j****s
发帖数: 156 | 6 截取？怎样截取？应该不是unitary了，check in matlab.
全1
【在 s******f 的大作中提到】
: 但是截取了以后，P1*P2不也是unitary么？ 为什么这个时候singular value就不为全1
: 了呢
|
d******e
发帖数: 7844 | 7 截取之后怎么可能还是Unitary matrix啊
【在 s******f 的大作中提到】
: 我说的截取意思是取P1的某n列， P1'*P2 不也是n by n的unitary matrix..为什么此
: 时singular value就不全1了
|
