- 博客(15)
- 收藏
- 关注
RSS订阅
原创 任意阶球谐函数可视化
先上效果图。从上到下依次为n=0:5, 从左到右依次为m=0:5, n代表阶，m代表级。上图的颜色代表的就是球谐函数的取值（后续我会把解析公式放上来），但是半径却是加工过的（具体参考下面的代码）。如果把半径设为球谐函数的绝对值，如下图所示不再过多解释，下面给出matlab制图代码% 绘制球谐函数图像% degree, Ynm的n, 阶% order, Ynm的m, 级...
2019-12-12 11:11:43
8163
10
原创 单样本KS检验，案例分析及Matlab可视化
Kolmogorov-Smirnov是比较一个经验分布与理论分布或者两个观测值分布的检验方法。其原假设H0:两个数据分布一致或者数据符合理论分布。 对于单样本KS检验，检验统计量为，当根据观测值计算出的则拒绝H0，否则接受H0假设。目录基于数据的KS检验统计量的展开式理论分布在接受域中的上下限反向验证案例分析或许可行的改进及建议附件：Kolmogorov-...
2019-11-23 14:48:46
7648
2
原创 对“陶哲轩-来自特征值的特征向量”的理解
柯西-比内(Cauchy-Binet)公式的证明 及 对《来自特征值的特征向量》的理解0.说明1. 化简的方法2. 运用行列式3. 联立可得4. 特殊情况15. 特殊情况26. 特殊情况37. 特殊情况48. 解释及备注9. 特征向量归一化的验证 据新闻报道, 3个物理学家和数学天才陶哲轩研究出一个只用特征值就可以计算矩阵特征向量的公式...
2019-11-18 09:10:02
2022
原创 Matlab 实验笔记(Live Script)保存为html格式，但Title为Untitled? 如何修改？
实验笔记是个好东西。如果你做了一堆研究，最后发现以前做过的东西忘记放在哪里了，或者说代码已经被你修改得面目全非了，你想对比一下当下的结果和当时的结果，这可能会给你带来无尽的烦恼。 Python的Jupyter确实可以解决上述的问题，它可以代码，结果，备注(说明，公式等)展示在同一个文件中，打开速度也很快，如果没有别的需求，这应该是最完美的工具。 但是我对Matlab依赖比较强，想用Matlab实现同样的功能，该怎么办呢？ Matlab的确提供了相关的功能，即实时脚本（live script）,具体使...
2020-06-30 21:53:47
3413
原创 Matlab实时函数提示
需求描述：在Matlab中自定义了一系列函数，但是这些函数的输入变量比较难以记忆，希望在输入函数时实时提示。分析：Maltab自带的有‘函数提示’功能，默认快捷键为Ctrl+F1，效果如下我们只需将Ctrl+F1修改为Shift+9，即每次输入左括号时出现函数提示。在Matlab预测选项中，具体修改办法如下如此，我们就可以畅快得敲代码了，简直不要太爽~...
2020-06-27 15:11:23
7466
2
原创 用C#编写一个图片浏览器，实现鼠标控制图片的平移缩放与图片的灰度化
1. 界面设计如图1 所示：一个名为ImView 的Form 只中包含有一个名为picturebox 的Picturebox。图2 是它的运行结果。该程序的界面设计较为简单，主要根据鼠标行为及键盘按键等事件触发行为。目前实现的功能有：图片的灰度化，鼠标滚轮控制图片缩放，鼠标控制图像平移。图1图22. 打开图像该部分的具体实现流程：双击主界面的空白部分=》弹出选择图片对话框=...
2020-05-03 12:43:11
1882
原创 一图看懂三角方程组的并行求解
三角方程组的并行求解（伪代码）算法流程的理解假设处理器个数p=3，矩阵L按列卷帘存储——即将L的第n列存放在处理器上。这样做是为了负载均衡，防止最后阶段有很多处理器空闲；一个比上述算法更好的改进方案就是将矩阵L的行和列都按列卷帘存储，当然算法会更加复杂一些。下面给出一个形象的并行计算流程图，结合上述算法和这里的流程图能够帮助理解。...
2020-05-03 11:55:11
1542
原创 基于TensorFlow.Keras实现的混合密度模型（MDN）
Note全文及代码参考自：混合密度模型Mixture Density Networks原文对模型解释的很清楚，我就不再过多解释。 如需深入，建议学习高斯混合模型(GMM)。原文的代码有些过时了，我精简了原始的代码，以MarkDown的形式展示在自己的博客里。1 The Most Common Neural Network1.1 Import Modules* TensorFlow ...
2020-04-21 12:18:54
1969
原创 LSTM反向传播算法 (基于矩阵求导的方法)
如下图所示（图片来自网络）。LSTM的核心思想是使得细胞状态Ct-1通过传送带，只经过少量的交互就输出Ct ，这使得之前的信息能够传递到当前，解决了RNN误差后向传播的梯度消失问题。它与RNN一样，参数是共享。如下图所示，这三个LSTM细胞的参数是相同的，不同的是输入值不同，因此输出值Ht和Ct也不同。对应的数学模型为如果我们使用随机梯度下降法对网络进行优化，则输入数据和输...
2020-03-31 12:20:27
1316
原创 一般线性秩统计量的方差函数 及其 极限分布
秩 及 线性秩统计量的概念见参考书1的P105或参考书2的P21-P23问题引入参考书2-P23-定理4.1下的（4.5）式给出了线性秩统计量的方差函数而对于Wilcoxon符号秩和统计量将其带入上式发现这仍然是一个随机变量，显然是与方差的性质矛盾的。我在这个问题上思考了很久，最后决定将Wilcoxon符号秩和统计量这类的统计量定义为一般线性秩统计量，然后推导...
2018-06-09 21:23:00
1429
原创 高斯-马尔科夫过程的充要条件(正态随机过程又是马尔科夫过程的充要条件)
必要性的证明充分性的证明参考参考1：《概率论与数理统计教材》（茆诗松，第二版）参考2：[公式推导]用最简洁的方法证明多元正态分布的条件分布参考3：《线性统计模型-线性回归与方差分析》（王松桂）参考3：百度文库--《随机过程-正态马尔科夫过程》。后续更新：在定理1的基础上证明的定理2：定理2就很有实用价值，由于平...
2018-06-04 22:27:00
3683
原创 用最简洁的方法证明多元正态分布的条件分布
在证明高斯-马尔科夫随机过程的性质的过程中，遇到了多元正态分布的条件分布的证明，百度发现条件分布的很多证明方法写的极其麻烦，所以自己写了一个。实际上多元随机变量的公式证明一般用矩阵的方法处理，这里就是采用这种方法，处理的结果非常好，证明过程很简洁，给大家推荐。推导的过程中，公式有可能有输错，读者有问题的地方可以留言指出。引理1引理2正式证明结论这个公式...
2018-06-03 23:00:00
3079
1
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人