高斯消元法，是线性代数中的一个算法，可用来求解线性方程组，并可以求出矩阵的秩，以及求出可逆方阵的逆矩阵。高斯消元法的原理是：若用初等行变换将增广矩阵 化为 ，则AX = B与CX = D是同解方程组。所以我们可以用初等行变换把增广矩阵转换为行阶梯阵，然后回代求出方程的解。以上是线性代数课的回顾，下面...
分类：
其他好文 时间：
2015-02-11 18:30:17
阅读次数：
210
2.6 .1 符号矩阵分析
常用的矩阵分析指令如下：
det(A)
行列式A
diag(A)
取对角元素构成向量，或根据向量构成对角矩阵
[V,D]=eig(A)
特征值分解，使 AV=VD （没看懂）
expm(A)
矩阵指数e^A
inv(A)
A的逆矩阵
[V,J]=jordan(A)
Jordan分解，使AV=VJ
poly(A)
矩阵的特征多项式
rank...
分类：
其他好文 时间：
2015-02-08 16:56:48
阅读次数：
201
有时候我们要比较两个方法的运行效率，那么就可以利用tic、toc这两个函数去估算运行时间。使用方法如下：tic
for k=1:100
s=sum(k);
end
toc 运行结果：
Elapsed time is 0.000003 seconds.
是不是很容易呢？用tic，toc把程序围起来就可以了！下面对比一下MATLAB中求逆矩阵两种方法的效率：
tic;A=...
分类：
其他好文 时间：
2015-01-27 22:00:54
阅读次数：
323
线性规划是凸优化的一种，故线性规划只是凸优化的一个子集。是线性规划一定属于凸优化，但反之不一定成立。凸优化与线性规划定义的区别有两点：系数，等号最小二乘优化没有限制条件。但A行数不小于列数，使得转置再相乘的时候有可能得逆矩阵。最小二乘法可以通过增加权重或附加项来调整，增加灵活性。所有从n维实空间到实...
分类：
其他好文 时间：
2015-01-27 17:46:42
阅读次数：
169
本人渣学生，最近些毕设论文，看了一下德布罗意的毕设拿了炸药奖，顿时觉的自己也要得努力，不拿炸药奖，起码对得起自己。。然后我得机子又是很烂的那种，word开不了20页。然后我考虑用latex写毕设。然后写到第4章。然后第四章要写个广义逆矩阵。。然后问题来了。。。我很快就写了一个矩阵。如下：
\begin{equation}\label{equ:star_composention}
\left(...
分类：
其他好文 时间：
2015-01-17 08:53:25
阅读次数：
580
EM算法大致分为两步——E步骤和M步骤。
而在求解运算过程中，需要用到高斯分布，逆矩阵等数学知识。EM算法上篇先梳理一下基础的数学知识，具体EM算法的核心思想下篇再进行介绍。
由于公式，矩阵太多，便手写一份推导，贴在下面。...
分类：
编程语言 时间：
2015-01-09 21:00:58
阅读次数：
272
做数据结构课设时候查的资料，主要是看求逆矩阵方面的知识的。 选主元的高斯-约当（Gauss-Jordan）消元法在很多地方都会用到，例如求一个矩阵的逆矩阵、解线性方程组（插一句：LM算法求解的一个步骤），等等。它的速度不是最快的，但是它非常稳定（来自网上的定义：一个计算方法，如果在使用此方法的计算过...
分类：
其他好文 时间：
2015-01-05 20:12:12
阅读次数：
260
完成R Programming第三周这周作业有点绕，更多地是通过一个缓存逆矩阵的案例，向我们示范【词法作用域 Lexical Scopping】的功效。但是作业里给出的函数有点绕口，花费了我们蛮多心思。Lexical Scopping:The value of free variables are ...
分类：
其他好文 时间：
2014-12-20 19:39:27
阅读次数：
222
矩阵A式n阶可逆矩阵的等价条件:
1. A的行列式不等于0
2. A的秩等于n，即A为满秩矩阵
3. A的行（列）向量组线性无关
4. 齐次方程组Ax=0只有零解
5. 对于任意b属于Rn(n为上标，表示向量空间)，Ax=b总有唯一解
6. A与单位矩阵等价
7. A可表示成若干个初等矩阵的乘积
8. A的列向量可以作为n维向量空间Rn(n为上标)的一组基
9. Rn中任意一个向...
分类：
其他好文 时间：
2014-12-09 17:45:53
阅读次数：
948
高斯消元法：
数学上，高斯消元法（或译：高斯消去法）（英语：Gaussian Elimination），是线性代数中的一个算法，可用来为线性方程组求解，求出矩阵的秩，以及求出可逆方阵的逆矩阵。当用于一个矩阵时，高斯消元法会产生出一个“行梯阵式”。（来自维基百科）
构造如下方程：
a[0][0]*X0 + a[0][1] *X1 + a[0][2]*X2+...........a[0][n...
分类：
其他好文 时间：
2014-11-27 20:33:53
阅读次数：
309
