最大似然函数
我们没法直接用求导解方程的办法直接求得最大值。可以采用EM算法
换成公式3形式求最大似然函数，其中Qi(z(i))用以下求，p(xi,zi;θ)为下式分子带入3中求最大似然函数的参数值为2.
其中 ，并且 也顺理成章地可以估计为 。
重复迭代前面两步，直到似然函数的值收敛为止。
下列为GMM参考代码l连接点击...
分类：
其他好文 时间：
2016-05-12 12:37:15
阅读次数：
195
这道题是这样的： 公鸡5文钱一只，母鸡3文钱一只，小鸡三只1文钱，那怎么买，100文刚好买到100只呢？ 好！干的漂亮！所有的信息你已经提取完毕！那么接下来把它转化为python语言即可！ 你肯定知道的：电脑的计算能力比人快的多得多，人计算一次，可能电脑已经计算了成千上万次，这还算低了，对吧！这就是 ...
分类：
其他好文 时间：
2016-05-03 22:10:40
阅读次数：
250
终于又回到熟悉的Round了 数学 A - Pouring Rain 设个未知数，解方程，还好没有hack点 数学 B - Coat of Anticubism 题意：求增加最小长度的一根木棍，使得构成一个多边形。 分析：那么构成三角形，原来n条木棍分成A，B两边，A和B接近(A<=B)，那么另一条 ...
分类：
其他好文 时间：
2016-04-30 15:28:12
阅读次数：
163
高斯消元法，是线性代数中的一个算法，可用来求解线性方程组，并可以求出矩阵的秩，以及求出可逆方阵的逆矩阵。高斯消元法的原理是：若用初等行变换将增广矩阵 化为 ，则AX = B与CX = D是同解方程组。所以我们可以用初等行变换把增广矩阵转换为行阶梯阵，然后回代求出方程的解。 以上是线性代数课的回顾，下 ...
分类：
其他好文 时间：
2016-04-24 15:36:49
阅读次数：
180
A.Parallelepiped 题意：给一个六面体三面的面积，求12条边的长度和。 题解：因为都是整数，设边长分别为a,b,c，面积为xyz,那么可设x=a*b,y=b*c,z=c*a，简单解方程就可以了。 B.Array 题意：给一个长度为N的数列，求一个连续子序列，包含K个不同的数。要求这个子 ...
分类：
其他好文 时间：
2016-04-15 21:56:30
阅读次数：
250
（1）求方程ax^2+bx+0的解 方程有解有以下几种可能： 1.b^2-4ac=0,方程有两个相等实根。 2.b^2-4ac>0,方程有两个不等实根。 3.b^2-4ac<0,方程有两个共轭复根。 程序中需要判断b^2-4ac是否等于0，由于它是实数，在存储时会有一些微笑的误差，因此不能直接进行( ...
分类：
其他好文 时间：
2016-04-05 10:42:20
阅读次数：
185
方程的应用，减少了解决算术问题的思维量。有了方程，就不必砍腿去解鸡兔同笼了。列个方程，或是方程组，通过对方程进行代数变换，就可以解出未知量。 不过对方程进行代数变换是要讲逻辑的，不能随便乱变。前几天在 Quora 上看到类似这样的一道题： 已知 \( x+\frac{1}{x}=2\sqrt{5} ...
分类：
其他好文 时间：
2016-04-03 16:03:12
阅读次数：
167
题目：n个同学排成一队，一共多少种排法?设n个同学排队，共有f(n)种排法，再插入一个，对于原来已经排列的每种队型，可以放在第1位，第2位...n+1种可能，得到递归方程:f(n+1)=(n+1)f(n)，解方程如下：
分类：
其他好文 时间：
2016-03-23 15:37:40
阅读次数：
167
求方程y1=3x1+x2+x3，y2=3x1-x2-x3在x1=-2，x2=3，x3=1时的值。解：打开MATLAB输入如下指令：x1=-2;x2=3;x3=1;>>y1=3*x1+x2+x3y1=-2>>y2=3*x1-x2-x3y2=-10
分类：
其他好文 时间：
2016-03-13 18:12:38
阅读次数：
133
