矩阵的QR分解
#include <cstdio> #include <cstdlib> #include <algorithm> #include <cmath> #include <cassert> #include <ctime> class MclVector { public: int n; double *Mat; /** type=0: 列向量 type=1: 行向量 **/ int type; MclVector() { Mat=NULL; n=; } MclVector(int len,double initVal=0.0) { n=len; Mat=]; ;i<=n;i++) Mat[i]=initVal; type=; } double operator[](int id) const { return Mat[id]; } double& operator[](int id) { return Mat[id]; } double length() const { ; ;i<=n;i++) sum+=Mat[i]*Mat[i]; return sqrt(sum); } MclVector operator*(double val) const { MclVector ans=MclVector(n); ;i<=n;i++) ans[i]=Mat[i]*val; return ans; } MclVector operator/(double val) const { MclVector ans=MclVector(n); ;i<=n;i++) ans[i]=Mat[i]/val; return ans; } MclVector operator+(const MclVector &newVector) const { MclVector ans=MclVector(n); ;i<=n;i++) ans[i]=Mat[i]+newVector[i]; return ans; } MclVector operator-(const MclVector &newVector) const { MclVector ans=MclVector(n); ;i<=n;i++) ans[i]=Mat[i]-newVector[i]; return ans; } MclVector operator*=(double val) { ;i<=n;i++) Mat[i]=Mat[i]*val; return *this; } MclVector operator/=(double val) { ;i<=n;i++) Mat[i]=Mat[i]/val; return *this; } MclVector operator+=(const MclVector &newVector) { ;i<=n;i++) Mat[i]+=newVector[i]; return *this; } MclVector operator-=(const MclVector &newVector) { ;i<=n;i++) Mat[i]-=newVector[i]; return *this; } MclVector GetTranspose() const { MclVector ans=*this; ans.type=; return ans; } void print() const { ;i<=n;i++) printf("%8.3lf ",Mat[i]); puts(""); } }; class MclMatrix { public: int row,col; MclVector *Mat; MclMatrix() {Mat=NULL;} MclMatrix(int _row,int _col,double initVal=0.0) { row=_row; col=_col; Mat=]; ;i<=row;i++) Mat[i]=MclVector(col,initVal); } void setIdentityMatrix() { ;i<=row;i++) { ;j<=col;j++) { ; ; } } } MclMatrix GetTranspose() const { MclMatrix ans=MclMatrix(col,row); ;i<=ans.row;i++) { ;j<=ans.col;j++) { ans[i][j]=Mat[j][i]; } } return ans; } void print() const { ;i<=row;i++) Mat[i].print(); puts(""); } MclVector& operator[](int id) const { return Mat[id]; } MclVector& operator[](int id) { return Mat[id]; } MclMatrix operator*(const MclMatrix &Matrix) const { MclMatrix ans=MclMatrix(row,Matrix.col); ;i<=row;i++) { ;j<=Matrix.col;j++) { ;k<=col;k++) { ans[i][j]+=Mat[i][k]*Matrix[k][j]; } } } return ans; } MclMatrix operator+(const MclMatrix &Matrix) const { MclMatrix ans=MclMatrix(row,Matrix.col); ;i<=row;i++) { ;j<=Matrix.col;j++) { ans[i][j]=Mat[i][j]+Matrix[i][j]; } } return ans; } MclMatrix operator-(const MclMatrix &Matrix) const { MclMatrix ans=MclMatrix(row,Matrix.col); ;i<=row;i++) { ;j<=Matrix.col;j++) { ans[i][j]=Mat[i][j]-Matrix[i][j]; } } return ans; } MclVector GetCol(int colId) const { MclVector ans=MclVector(row); ;i<=row;i++) ans[i]=Mat[i][colId]; return ans; } MclVector GetRow(int rowId) const { MclVector ans=MclVector(row); ;i<=col;i++) ans[i]=Mat[rowId][i]; return ans; } MclMatrix operator*=(const MclMatrix &Matrix) { return *this=*this*Matrix; } MclMatrix operator+=(const MclMatrix &Matrix) { return *this=*this+Matrix; } MclMatrix operator-=(const MclMatrix &Matrix) { return *this=*this-Matrix; } MclMatrix operator*(double x) const { MclMatrix ans=*this; ;i<=row;i++) { ;j<=col;j++) { ans[i][j]*=x; } } return ans; } }; MclMatrix vectorMulVector(const MclVector &A,const MclVector& B) { ) { MclMatrix ans=MclMatrix(A.n,B.n); ;i<=A.n;i++) { ;j<=B.n;j++) { ans[i][j]+=A[i]*B[j]; } } return ans; } else { assert(A.n==B.n); MclMatrix ans=MclMatrix(,); ;i<=A.n;i++) { ans[][]+=A[i]*B[i]; } return ans; } } int sgn(double x) { ; ; ; } /** 将矩阵A分解为一个正交矩阵Q和一个上三角矩阵R A为任意实数矩阵 **/ std::pair<MclMatrix,MclMatrix> QRSplit(const MclMatrix &A) { assert(A.col==A.row); int n=A.row; MclMatrix Q=MclMatrix(n,n); Q.setIdentityMatrix(); MclMatrix R=A; ;i<n;i++) { MclVector s=R.GetCol(i); ;j<i;j++) s[j]=; ) continue; double c=s.length(); ) c*=-sgn(R[i][i]); MclVector u=s; u[i]-=c; MclVector uT=s.GetTranspose(); MclMatrix H=MclMatrix(n,n); H.setIdentityMatrix(); H=H-vectorMulVector(u,uT)*(2.0/(u.length()*u.length())); R=H*R; Q=Q*H; } return std::make_pair(Q,R); }
矩阵的QR分解的更多相关文章
- 矩阵的QR分解(三种方法)Python实现
1.Gram-Schmidt正交化 假设原来的矩阵为[a,b],a,b为线性无关的二维向量,下面我们通过Gram-Schmidt正交化使得矩阵A为标准正交矩阵: 假设正交化后的矩阵为Q=[A,B],我 ...
- QR 分解
将学习到什么 介绍了平面旋转矩阵,Householder 矩阵和 QR 分解以入相关性质. 预备知识 平面旋转与 Householder 矩阵是特殊的酉矩阵,它们在建立某些基本的矩阵分解过程中起着 ...
- 机器学习中的矩阵方法03:QR 分解
1. QR 分解的形式 QR 分解是把矩阵分解成一个正交矩阵与一个上三角矩阵的积.QR 分解经常用来解线性最小二乘法问题.QR 分解也是特定特征值算法即QR算法的基础.用图可以将分解形象地表示成: 其 ...
- 矩阵QR分解
1 orthonormal 向量与 Orthogonal 矩阵 orthonormal 向量定义为 ,任意向量 相互垂直,且模长为1: 如果将 orthonormal 向量按列组织成矩阵,矩阵为 ...
- QR分解
从矩阵分解的角度来看,LU和Cholesky分解目标在于将矩阵转化为三角矩阵的乘积,所以在LAPACK种对应的名称是trf(Triangular Factorization).QR分解的目的在 ...
- QR分解与最小二乘
主要内容: 1.QR分解定义 2.QR分解求法 3.QR分解与最小二乘 4.Matlab实现 一.QR分解 R分解法是三种将矩阵分解的方式之一.这种方式,把矩阵分解成一个正交矩阵与一个上三角矩阵的 ...
- QR分解与最小二乘(转载自AndyJee)
转载网址:http://www.cnblogs.com/AndyJee/p/3846455.html 主要内容: 1.QR分解定义 2.QR分解求法 3.QR分解与最小二乘 4.Matlab实现 一. ...
- QR分解迭代求特征值——原生python实现(不使用numpy)
QR分解: 有很多方法可以进行QR迭代,本文使用的是Schmidt正交化方法 具体证明请参考链接 https://wenku.baidu.com/view/c2e34678168884868762d6 ...
- MATLAB矩阵的LU分解及在解线性方程组中的应用
作者:凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ 三.实验程序 五.解答(按如下顺序提交电子版) 1.(程序) (1)LU分解源程序: function [ ...
随机推荐
- C#构造Http 破解学校教务系统学生账号密码
背景介绍 我们学校的教务系统的是以学生学号作为登陆账号,初始密码是自己的生日. 一点点想法 每次期末查成绩的时候,我都会有一个想法,要是我能跑到系统后台,把自己的成绩修改一下,那该时间多么舒坦的事情啊 ...
- Unity学习疑问记录之Apply Root Motion
Should we control the character's position from the animation itself or from script. 如果我们勾选了Animator ...
- 11高级网站构建:div和span
用<div>元素把属于一个逻辑部分的元素包围起来.可以用id属性为<div>提供一个唯一的标签. <div>的作用:1.更深一步展示页面的基本逻辑结构(相当于一个逻 ...
- ngModel
https://docs.angularjs.org/error/ngModel/numfmt?p0=sa angular.module('myApp', []) .directive('tagLis ...
- blockdev命令和blkid命令
blockdev命令和blkid命令 http://www.jb51.net/LINUXjishu/310389.html block相关的命令 这篇文章主要介绍了Linux blockdev命令设置 ...
- Python之路-python(面向对象一)
一.面向对象介绍 二.为什么要用面向对象开发 三.封装.继承.多态.类.方法 面向过程和面向对象的区别 编程范式: 一般情况下,拿到一个项目不同的两个人有不同的编写方式(相同的是固定的语法.数据结构) ...
- [Android Tips] 1. Getting StatusBar Height
public int getStatusBarHeight() { int result = 0; int resourceId = getResources().getIdentifier(&quo ...
- 数据库的Instance/Crash Recovery
crash recovery是指单实例数据库发生了failure.或者rac数据库中的所有实例都发生了failure后进行的recovery.rac数据库crash后,rac中第一个重启启动的inst ...
- ef第一次启动较慢
解决ef第一次启动较慢的问题: protected void Application_Start() { //禁用第一次ef查询对表__MigrationHistory的问题使用了ef的Code fi ...
- Android开发新手学习总结(一)——使用Android Studio搭建Android集成开发环境
[新手连载]一:使用Android Studio搭建Android集成开发环境http://bbs.itcast.cn/forum.php?mod=viewthread&tid=87055&a ...