锐客网

  1. 首页 > 生活百科 > it技术 > >

svdjava代码 java 代码

IT技术Double

opencv中把矩阵进行奇异值分解后怎样重构矩阵？SVD相当于：
C为mxn阶矩阵，U为mxm阶酉矩阵 ， E为mxn阶实数对角矩阵，V为nxn阶酉矩阵 。
E矩阵对角线上的元素等于C的奇异值 。
在OpenCV中可以用
solve(InputArray src1, InputArray src2, OutputArray dst, intflags=DECOMP_SVD)
你是不是要用SVD求解最小二乘问题？
线性方程组Cx=b，求其最小二乘解 。
你可以先求出C的伪逆C。
C =V(E )(UT)
(E )是E的伪逆，将E主对角线上每个非零元素都求倒数之后再转置得到.
(UT)是U的转置 。
最后计算出最小二乘解为：(C )b
java 一个字符串12345678，怎样每隔2个数字 ， 用#连接起来String getStr(String str) {
if (str.length() = 2) {
return str;
}
return str.substring(0, 2)"#"getStr(str.substring(2));
}
不要太简单
求SVD算法的C实现代码/** Cfunction for SVD
函数原型:
bool svd(vectorvectordoubleA, int K, std::vectorstd::vectordoubleU, std::vectordouble S, std::vectorstd::vectordoubleV);
其中
A是输入矩阵svdjava代码，假设Asvdjava代码的维数是m*nsvdjava代码 ， 那么本函数将A分解为U diag(S) V'
其中U是m*K的列正交的矩阵. V是n*K的列正交矩阵svdjava代码，S是K维向量 。K由第二个参数指定 。
U的第i列是A的第i大奇异值对应的左歧义向量,S[i]=A的第 i大奇异值,V的第i列是A的第i大奇异值对应的右歧义响亮.
K是需要分解的rank，0K=min(m,n)
本程序采用的是最基本幂迭代算法，在linux g下编译通过
**/
#include cmath
#include iostream
#include iomanip
#include cstdlib
#include cstring
#include fstream
#include vector
using namespace std;
const int MAX_ITER=100000;
const double eps=0.0000001;
double get_norm(double *x, int n){
double r=0;
for(int i=0;in;i)
r =x[i]*x[i];
return sqrt(r);
}
double normalize(double *x, int n){
double r=get_norm(x,n);
if(reps)
return 0;
for(int i=0;in;i)
x[i]/=r;
return r;
}
inline double product(double*a, double *b,int n){
double r=0;
for(int i=0;in;i)
r =a[i]*b[i];
return r;
}
void orth(double *a, double *b, int n){//|a|=1
double r=product(a,b,n);
for(int i=0;in;i)
b[i]-=r*a[i];
}
bool svd(vectorvectordoubleA, int K, std::vectorstd::vectordoubleU, std::vectordouble S, std::vectorstd::vectordoubleV){
int M=A.size();
int N=A[0].size();
U.clear();
V.clear();
S.clear();
S.resize(K,0);
U.resize(K);
for(int i=0;iK;i)
U[i].resize(M,0);
V.resize(K);
for(int i=0;iK;i)
V[i].resize(N,0);
srand(time(0));
double *left_vector=new double[M];
double *next_left_vector=new double[M];
double *right_vector=new double[N];
double *next_right_vector=new double[N];
while(1){
for(int i=0;iM;i)
left_vector[i]= (float)rand() / RAND_MAX;
if(normalize(left_vector, M)eps)
break;
}
int col=0;
for(int col=0;colK;col){
double diff=1;
double r=-1;
for(int iter=0;diff=epsiterMAX_ITER;iter){
memset(next_left_vector,0,sizeof(double)*M);
memset(next_right_vector,0,sizeof(double)*N);
for(int i=0;iM;i)
for(int j=0;jN;j)
next_right_vector[j] =left_vector[i]*A[i][j];
r=normalize(next_right_vector,N);
if(reps) break;
for(int i=0;icol;i)
orth(V[i][0],next_right_vector,N);
normalize(next_right_vector,N);
for(int i=0;iM;i)
for(int j=0;jN;j)
next_left_vector[i] =next_right_vector[j]*A[i][j];
r=normalize(next_left_vector,M);
if(reps) break;
for(int i=0;icol;i)
orth(U[i][0],next_left_vector,M);
normalize(next_left_vector,M);
diff=0;
for(int i=0;iM;i){
double d=next_left_vector[i]-left_vector[i];
diff =d*d;
}
memcpy(left_vector,next_left_vector,sizeof(double)*M);
memcpy(right_vector,next_right_vector,sizeof(double)*N);
}
if(r=eps){
S[col]=r;
memcpy((char *)U[col][0],left_vector,sizeof(double)*M);
memcpy((char *)V[col][0],right_vector,sizeof(double)*N);
}else
break;
}
delete [] next_left_vector;
delete [] next_right_vector;
delete [] left_vector;
delete [] right_vector;
return true;
}
void print(vectorvectordoubleA){
for(int i=0;iA.size();i){
for(int j=0;jA[i].size();j){
coutsetprecision(3)A[i][j]' ';
}
coutendl;
}
}
int main(){
int m=10;
int n=5;
srand(time(0));
vectorvectordoubleA;
A.resize(m);
for(int i=0;im;i){
A[i].resize(n);
for(int j=0;jn;j)
A[i][j]=(float)rand()/RAND_MAX;
}
print(A);
coutendl;
vectorvectordoubleU;
vectordouble S;
vectorvectordoubleV;
svd(A,2,U,S,V);
cout"U="endl;
print(U);
coutendl;
cout"S="endl;
for(int i=0;iS.size();i){
coutS[i]' ';
}
coutendl;
cout"V="endl;
print(V);
return 0;
}
SVD不收敛问题，怎么解决我会
1.打开安装目录里的mods\bf2\object_server.zip(注:这是个压缩文件,先不用解压，直接浏览里面的文件就行)
2.用记事本打开Weapons\Handheld\rurif_dragunov.tweak(注:rurif_dragunov其实就是svd的代码)
然后找到ObjectTemplate.damage 45 这一行，45就是svd的威力 。改大了就可以增加威力 。
有问题再追问 。
如何将SVD算法用mapreduce实现数据挖掘比赛算法
examples/src/main/java/org/apache/mahout/cf/taste/example/kddcup/track1/svd
推荐系统中利用SVD实现降维
core/src/main/java/org/apache/mahout/cf/taste/impl/recommender/svd
【svdjava代码 java 代码】关于svdjava代码和java 代码的介绍到此就结束了，不知道你从中找到你需要的信息了吗 ？如果你还想了解更多这方面的信息，记得收藏关注本站 。

    推荐阅读


    上一篇：小米路由器1t怎么用，小米路由器1t怎么用

    下一篇：账号运营直播内容，账号运营直播内容介绍