13.1 极大似然估计的原理 极大似然的估计原理可以由下面的程序得到说明。我们首先生成 10 个服从 正态分布的总体,每个总体的均值都不同,依次为 0,1,2,3,4,5,6,7,8, 9。方差相同,均为 1。然后我们随机地取出一个总体,从中抽出 10 个样本,因 为事先不知道是从哪一个总体中抽出来的,所以我们分别用已知的 10 个总体参 数值代入似然函数,计算出 10 个似然函数值,取其中 大的似然值,认为该样 本是从相应的总体中取出的(从而联合概率密度也 大化)。然后我们让计算机 告诉我们它是从第几个总体中取样的,并与我们的判断进行对比。 *===========================begin================================== capt prog drop mle prog mle /*生成10个均值不同、方差均为1的正态总体,每个总体取8个样本*/ drawnorm double x0-x9,n(8) m(0,1,2,3,4,5,6,7,8,9) clear global i=int(10*uniform()) //设定一个随机数,用于随机取出一个总体 forv j=0/9 { gen lnf`j' =-0.5*ln(2*_pi)*8-sum(0.5*(x$i-`j')^2) //对取出的总体计算似然值 scalar lnf`j'=lnf`j'[_N] //最终的似然值 } scalar list // 比较10个似然值哪个最大,猜想是从第几个总体取出来的? end mle *根据10个似然值,猜想是从第几个总体取出来的? di "所抽中的样本为" as error "X"$i //显示真正的取样总体是什么 *===========================end==================================== 在现实中,我们并不知道任何一个真正的总体参数,因此,只能借助于找到 样本似然值(实际上是联合概率密度的对数值) 大的总体参数,即认为其是总 体参数。在 STATA 中实现 大似然法的估计必须自己编写程序。下面的例子说 明了如何利用 stata 编写程序来实现对模型的极大似然估计。 13.2 正态总体均值和方差的极大似然估计 *===========================begin================================== capt prog drop bb prog bb //定义程序的名称 args lnf u v //声明参数,u 为均值,v为方差 quietly replace `lnf' = -0.5*ln(2*_pi) - ln(`v') -0.5*($ML_y1-`u')^2/(`v')^2 end drawnorm x,n(100) m(10) sd(3) clear//模拟均值为10,方差为3的100个正态样本 ml model lf bb (x=) (variance:) //利用迭代法则进行极大似然估计
2022-11-05 22:27:01 2.41MB stata
1
求随机相位正弦信号的均值,比较简单的程序
2022-11-05 18:01:54 509B 随机相位
1
均值漂移聚类matlab代码全局固有对称性的组表示 代码 这些代码是本文算法的核心部分。 对于本文,您可以从我们的项目页面下载:。 该代码已在Win7 x64 + Matlab 2012a (32 bit)上进行了测试。 提供的代码仅供学术使用。 使用率 运行demo_C2Shape.m以计算和可视化马模型上的全局折返内在对称性。 运行demo_moreSymmetriesShape.m以计算和可视化表模型上的所有全局固有对称性。 引文 如果您在研究中使用我们的代码,请引用: @article {Symmetry-PG-2017, title = {全局固有对称性的组表示}, 作者= {王辉和黄辉}, journal = {计算机图形学论坛(太平洋图形学)}, 音量= {36}, 数字= {7}, 年= {2017}, } 致谢 代码的某些部分来自Internet,以下是原始URL: 测地距离: 均值漂移聚类: nn搜索: b2r.m: HKS.m: 接触 如有任何问题,请随时与(慧王)联系!
2022-11-05 00:49:15 822KB 系统开源
1
请编写MATLAB语言程序,分别计算样本函数
2022-11-03 16:57:47 3KB matlab
1
LMD_局部均值分解的源代码以及验证程序
2022-11-03 09:50:17 22KB lmd 局部均值分解
1
clear clc; %disp char "基于C 均值算法的特征点聚类" strcat datestr now ; S double imread "E: Matlab 1 jpg" ;%读取源图像 %S double imread "E: Matlab 001 jpeg" ;%读取源图像 [sr sc] find S 255 ;%sr sc存储源图像中待聚类点的坐标 len length sr ;%len存储源图像中待聚类点的个数 group zeros len 1 ;%group存储特征点与类相对应的情况">clear clc; %disp char "基于C 均值算法的特征点聚类" strcat datestr now ; S double imread "E: Matlab 1 jpg" ;%读取源图像 %S double imread "E: Matlab 001 jpeg" ;%读取源图像 [sr sc] find S 255 ;%sr sc存储源图像中待聚类点的坐标 len [更多]
2022-10-26 18:44:05 4KB C-均值聚类
1
基于模糊C均值的快速点云的聚类分析代码,可用于模糊聚类
1
局域均值分解(LMD)方法是一种新的自适应信号处理方法,Teager能量算子是求解信号瞬时能量的非线性操作算法,能有效提取信号的瞬时能量。结合这两种方法的优点,提出一种基于LMD-Teager 变换的功率谱估计方法,给出了算法原理和步骤,讨论了功率谱估计的物理意义,并在与快速傅里叶变换(FFT)方法、小波变换对比的基础上,用短数据序列和非平稳信号进行了仿真验证。结果表明:该方法突破了FFT 方法中对所分析的信号必须平稳的要求,更适合于非平稳信号的处理;且对数据长度的要求较傅里叶方法低,而其分析精度和分辨率优于传统的傅氏方法和小波变换。
2022-10-24 13:44:07 2.63MB 傅里叶光 局域均值 Teager能 功率谱估
1
基于cv2进行等间隔采样和局部均值的图像缩缩放(python版本) 下载后安装cv2包编程即可实现!操作简单!
2022-10-21 13:07:59 841B 数字图像处理 python
1